[RegCNET] about the post-processing

ywyang ywyang at lzb.ac.cn
Mon Jan 1 04:35:25 CET 2007


Hello,everyone:
   Happy new year!
   I am wondering the content of  OUTHEAD.ctl file.
   In my simulation  the grib number of longitude is defined 56,and that of latitude is defined 94.But in the OUTHEAD.ctl the xdef is 263,and the ydef is 125. They are not consistent,why?
   The following is the content of  OUTHEAD.ctl file:

dset ^OUT_HEAD
title RegCM domain information
options big_endian
undef -1.e34
pdef   92   54 lcc   33.00   89.00   46.00   27.00   30.00   60.00   89.00  45000.  45000.
xdef  263 linear   62.68  0.2027
ydef  125 linear   19.81  0.2027
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

After run cruPGI3.x,I got TEM2000.nc file.Then I am trying to do the post-processing.According to http://www.ictp.trieste.it/%7Epubregcm/RegCM3/postproc.html ,I converted TEM2000.nc to format grads_grib and executed a series of commands in Grads:
ga-> open tmp2000.ctl 
Scanning description file:  tmp2000.ctl
Data file tmp2000.grb is open as file 1
LON set to 67.555 110.898 
LAT set to 22.598 44.216 
LEV set to 1013 1013 
Time values set: 2000:1:1:0 2000:1:1:0 
ga-> open OUT_HEAD.CTL 
Scanning description file:  OUT_HEAD.CTL
Data file OUT_HEAD is open as file 2
ga-> set vpage 0.2 10.8 0.2 8.3
Virtual page size = 11 8.40566 
ga-> run mapoff.gs
ga-> d tmp
Contouring: -20 to 20 interval 5 
ga-> run mask2.gs 

I found the graph generated is illogical,for the lon and lat lines is located at the southwest corner.What cause the fault?How to solve it?Any help would be appreciated.Thanks in advance.

Plus:the content of tem2000.ctl file:
dset ^tmp2000.grb
title file created by grads using lats4d available from http://dao.gsfc.nasa.gov/software/grads/lats4d/
undef 1e+20
dtype grib
index ^tmp2000.gmp
xdef 92 levels
 67.555  68.010
 68.466  68.924  69.383  69.843  70.304  70.767  71.231  71.696  72.163  72.630
 73.099  73.569  74.040  74.512  74.985  75.459  75.934  76.410  76.886  77.364
 77.842  78.322  78.802  79.282  79.764  80.246  80.729  81.212  81.696  82.181
 82.666  83.151  83.637  84.123  84.610  85.097  85.584  86.072  86.559  87.047
 87.535  88.023  88.512  89.000  89.488  89.977  90.465  90.953  91.441  91.928
 92.416  92.903  93.390  93.877  94.363  94.849  95.334  95.819  96.304  96.788
 97.271  97.754  98.236  98.718  99.198  99.678  100.158  100.636  101.114  101.590
 102.066  102.541  103.015  103.488  103.960  104.431  104.901  105.370  105.837  106.304
 106.769  107.233  107.696  108.157  108.617  109.076  109.534  109.990  110.445  110.898

ydef 54 levels
 22.598  22.988
 23.378  23.770  24.162  24.555  24.950  25.345  25.741  26.137  26.535  26.933
 27.333  27.733  28.133  28.535  28.937  29.341  29.744  30.149  30.554  30.960
 31.367  31.774  32.182  32.591  33.000  33.410  33.820  34.231  34.643  35.055
 35.467  35.880  36.294  36.708  37.122  37.537  37.952  38.368  38.784  39.200
 39.617  40.034  40.451  40.868  41.286  41.704  42.122  42.541  42.959  43.378
 43.797  44.216
zdef 1 levels 1013
tdef 12 linear 0Z1jan2000  1mo
vars 1
tmp       0    1,  1,  0,  0 Mean temperature []
endvars




        ywyang
        ywyang at lzb.ac.cn
          2007-01-01


More information about the RegCNET mailing list