On 16/12/2015 13:08, Delclaux Francois
wrote:
Hi Ferret users,
Sorry for my last mail, but I did not change the object of the
mail. So I re-post my question with a new subject
I have a netcdf file with a 3D variable according to the
following:
1> ./DRAIN.nc
(default)
name title I
J K L M N
DR2 Averaged_Soil_Drainage_Flux 1:9 1:13
1:26280 ... ... ...
kg/m2s on grid GJI1 with 1.E+20 for missing data
X=0.5E:9.5E Y=0.5N:13.5N Z=0.5:26281
LAT ... 1:13
... ... ... ...
degrees_north on grid GJI2 with 1.E+20 for
missing data
Y=0.5N:13.5N
LON 1:9 ...
... ... ... ...
degrees_east on grid GJI3 with 1.E+20 for missing
data
X=0.5E:9.5E
TIME ... ...
1:26280 ... ... ...
hours since 2013-01-01 00:00:00 on grid GJI4 with
1.E+20 for missing data
Z=0.5:26281
Now I would like to transform this file in a new one DR2 =
DR2(x,y,t) where the axes LAT, LON and TIME axis are directly
included in the file
What's the best way for doing that ?
Hi François,
Try these lines:
def axis/X/units="degrees_east" lonaxis=lon
def axis/Y/units="degrees_north" lataxis=lat
let dr2_georeferenced=dr2[gx=lonaxis, gy=lataxis]
then you should be able to do a:
shade dr2_georeferenced
Regards
Patrick
--
LSCE/IPSL, Laboratoire CEA-CNRS-UVSQ
Data Analysis and Visualization Engineer
ICMC - IPSL Climate Modelling Centre
--
|