[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