[RegCNET] A ddf problem
Simon Krichak
shimonk at post.tau.ac.il
Mon Nov 30 15:53:30 CET 2009
Hi,
Could anybody please explain me my mistake if possible.
I created the following ddf based on the OUT_HEAD.CTL I have.
THere is a mistake somewhere however since my attemts to draw
say precipitation give me seemingly correct distribution of rain but the values
are of order of -32700.
At the same time my precipitation is very reasonable if I open
the *.nc file (see ncdump below) using sdfopen.
Thanks,
Simon
ddf file
DSET /home/kuninpav/Netta/Dead_Sea/17km_run_Sep_2009/SRFM198601_199012AVG.nc
dtype netcdf
undef -1.e34
pdef 158 148 lcc 35.00 34.00 79.00 74.00 30.00 60.00 34.00 17000. 17000.
xdef 475 linear 15.96 0.0766
ydef 318 linear 22.71 0.0766
zdef 1 levels 1008.58
tdef 59 linear 00Z01JAN1986 1mo
vars 2
RT=>rt 0 t,y,x Precipitation
TA=>ta 0 t,y,x Anemom Temp
endvars
=======================================================================================
dset ^OUT_HEAD
title RegCM domain information
options big_endian
undef -1.e34
pdef 158 148 lcc 35.00 34.00 79.00 74.00 30.00 60.00 34.00 17000. 17000.
xdef 475 linear 15.96 0.0766
ydef 318 linear 22.71 0.0766
zdef 1 levels 1000.00
tdef 1 linear 00z01Jan2001 1mo
vars 11
head 0 99 header information
ht 0 99 surface elevation
htsd 0 99 surface elevation std dev
veg2d 0 99 vegetation type in BATS
landuse 0 99 surface landuse type
xlat 0 99 latitude of cross points
xlong 0 99 longitude of cross points
xmap 0 99 map factors of cross point
dmap 0 99 map factors of dot points
coriol 0 99 coriol force
mask 0 99 land/sea mask
endvars
=========================================================================================
ncdump NOF.nc
dimensions:
lon = 158 ;
lat = 148 ;
time = UNLIMITED ; // (59 currently)
variables:
float lon(lon) ;
lon:long_name = "Longitude" ;
lon:units = "degrees_east" ;
lon:actual_range = 19.33798f, 48.84574f ;
float lat(lat) ;
lat:long_name = "Latitude" ;
lat:units = "degrees_north" ;
lat:actual_range = 23.88342f, 46.65563f ;
double time(time) ;
time:long_name = "Time" ;
time:units = "hours since 1900-1-1 00:00:0.0" ;
time:actual_range = 754200., 796560. ;
short TA(time, lat, lon) ;
TA:long_name = "Anemom Temp" ;
TA:units = "K" ;
TA:missing_value = -32767s ;
TA:add_offset = 265.f ;
TA:scale_factor = 0.002594034f ;
short RT(time, lat, lon) ;
RT:long_name = "Total Precipitation" ;
RT:units = "mm/day" ;
RT:missing_value = -32767s ;
RT:add_offset = 1249.5f ;
RT:scale_factor = 0.03816281f ;
short TAMAX(time, lat, lon) ;
TAMAX:long_name = "Max Anemom Temp" ;
TAMAX:units = "K" ;
TAMAX:missing_value = -32767s ;
TAMAX:add_offset = 275.f ;
TAMAX:scale_factor = 0.002288853f ;
short TAMIN(time, lat, lon) ;
TAMIN:long_name = "Min Anemom Temp" ;
TAMIN:units = "K" ;
TAMIN:missing_value = -32767s ;
TAMIN:add_offset = 275.f ;
TAMIN:scale_factor = 0.002288853f ;
short W10MX(time, lat, lon) ;
W10MX:long_name = "Max 10m Wind Speed" ;
W10MX:units = "m/s" ;
W10MX:missing_value = -32767s ;
W10MX:add_offset = 0.f ;
W10MX:scale_factor = 0.01525902f ;
// global attributes:
:domxmin = 19.33798f ;
:domxmax = 48.84574f ;
:domymin = 23.88342f ;
:domymax = 46.65563f ;
:domzmin = 1050.f ;
:domzmax = 1050.f ;
========================================
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20091130/3a71579b/attachment.html>
More information about the RegCNET
mailing list