[Thread Prev][Thread Next][Index]
[las_users] wrong plots of data with curvilinear grids
LAS 7.0.2 tomcat 6 ferret 6.1
problem: A netcdf file with curvilinear coordinates is properly
plotted by Ferret, but wrongly by LAS.
LAS plots only the lower half of the (lat/lon) data, leaving the upper
half blank, as if there were missing values.
Probably I wrote a wrong xml file, but I did my best. If the xml is
right, then may be there are problems in LAS ferret/scripts dealing
with curvilinear transformations (may be LAS_transforms or
transform_curvi_xy ?).
I'm not able to solve the problem on my own, thus I ask for your help.
You can access LAS at
https://clima.casaccia.enea.it:8443/las
and then select the file "ernesto"
I prefere not to attach any file, thus please find at
ftp://clima.casaccia.enea.it/emanuele.lombardi
the following files:
ferret.gif = the output of ferret shade/k=1/l=1
votemper,nav_lon,nav_lat
LAS.gif = the outut of LAS
ernesto.nc = the netcdf file
ernesto.xml = the xml I use in LAS
ernesto.cdl = the "ncdump -h" of the netcdf file
Thanks for any help,
Emanuele
Here it is the "ncdump -h" of the netcdf file
netcdf ernesto {
dimensions:
deptht = 2 ;
y = 376 ;
x = 360 ;
time_counter = UNLIMITED ; // (2 currently)
variables:
float deptht(deptht) ;
deptht:long_name = "Vertical levels" ;
deptht:units = "m" ;
deptht:title = "deptht" ;
float nav_lat(y, x) ;
nav_lat:units = "degrees_north" ;
nav_lat:long_name = "Latitude" ;
nav_lat:nav_model = "Default grid" ;
float nav_lon(y, x) ;
nav_lon:units = "degrees_east" ;
nav_lon:long_name = "Longitude" ;
nav_lon:nav_model = "Default grid" ;
float time_counter(time_counter) ;
time_counter:units = "days since 2009-01-20 00:00:00"
;
time_counter:calendar = "gregorian" ;
time_counter:title = "Time" ;
time_counter:long_name = "Time axis" ;
time_counter:time_origin = "2009-01-20 00:00:00" ;
float votemper(time_counter, deptht, y, x) ;
votemper:units = "C" ;
votemper:missing_value = 0.f ;
votemper:valid_min = 10.f ;
votemper:valid_max = 35.f ;
votemper:long_name = "Temperat" ;
votemper:short_name = "votemper" ;
votemper:online_operation = "ave(only(x))" ;
votemper:axis = "TZYX" ;
votemper:interval_operation = 600.f ;
votemper:interval_write = 3600.f ;
votemper:associate = "time_counter deptht nav_lat
nav_lon" ;
// global attributes:
:Conventions = "GDT 1.3" ;
:file_name = " " ;
:production = "ENEA CLIMMOD" ;
:TimeStamp = " " ;
:associate_file = " " ;
:history = "Fri Jan 23 22:24:52 2009: ncks -v
votemper,nav_lon,nav_lat -d time_counter,0,1 -d deptht,0,1
TIRRfpE090120_01_grid_T.nc ernesto.nc" ;
}
And this is the xml I wrote for the above netcdf file:
<datasets>
<ciccio name="ernesto1" url="/home/bancadat/ernesto.nc">
<properties>
<ferret>
<curvi_coord_lon>nav_lon</curvi_coord_lon>
<curvi_coord_lat>nav_lat</curvi_coord_lat>
</ferret>
</properties>
<variables>
<votemper name="Potential Temperature" units="C">
<link match="/lasdata/grids/grid_xyzt-ciccio"/>
</votemper>
</variables>
</ciccio>
</datasets>
<grids>
<grid_xyzt-ciccio>
<link match="/lasdata/axes/x-ciccio"/>
<link match="/lasdata/axes/y-ciccio"/>
<link match="/lasdata/axes/z-ciccio"/>
<link match="/lasdata/axes/t-ciccio"/>
</grid_xyzt-ciccio>
</grids>
<axes>
<x-ciccio type="x" units="degrees_east">
<arange start="8.8125" size="360" step="0.0217" />
</x-ciccio>
<y-ciccio type="y" units="degrees_north">
<arange start="36.688" size="376" step="0.0208" />
</y-ciccio>
<z-ciccio type="z" units="m">
<v>1.47</v>
<v>4.55</v>
</z-ciccio>
<t-ciccio type="t" units="t">
<v>27-jan-2009:00:00:00</v>
<v>27-jan-2009:00:59:59</v>
</t-ciccio>
</axes>
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement