[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