<div dir="ltr">Dear GrADS Users,<br><br>I have downloaded ECMWF (ERA40) data for surface temperature at 2.5 by 2.5 resolution for years 1958-2001. My intention was to see the trend for the month of January only. So I extracted all the Januaries between 1958-2001 in a binary/dat file by using the following commands in GrADS. When I tried to open this file by making a CTL (choosing the same x & y dimentions) and display January of any one year, the spatial picture is not the same as if I display the same January in original data file. Could you please help me?<br>
<br><b>Here is the script to create a binary/dat file</b><br clear="all">'reinit'<br>'open EraT2M.ctl'<br>'set gxout fwrite'<br>'set fwrite January.dat'<br>t=5<br>*tmax=6<br>tmax=532<br>while(t<=tmax)<br>
'set t 't<br>'d no2tsfc-273.15'<br>t=t+12<br>endwhile<br>'disable fwrite'<br>'reinit'<br><br><b>Here is the control file to open the above January.dat file</b><br>DSET january.dat<br>TITLE Sample Data Set<br>
UNDEF -9.999E33<br>*OPTIONS big_endian<br>OPTIONS little_endian<br>*OPTIONS template<br>xdef 144 linear 0 2.5<br>ydef 73 linear -90 2.5<br>ZDEF 1 LINEAR 1 1<br>TDEF 44 LINEAR 1JAN1958 1mon<br>VARS 1<br>
no2tsfc 0 99 Surface Temperature_2m<br>ENDVARS<br><br><br>Regards,<br><br>Mubashar<br>
</div>