[gradsusr] merge data

[Ricardo Hallak]

Hi Ugo,

you extracted only one level (1000 hPa), but your .ctl to plot data is expecting 17 levels.

In the .ctl, try:
zdef 1 levels 1000

It must work.
Ricardo 

On Mon, 18 Jun 2012 10:26:08 -0700 (PDT), ugo okoro wrote
> hello all,
> 
> please i am working with ncep data in .nc format and had merged it with the help of this forum. however i am faced with some challenges as when i subtract the merged data with the original data for some variables, i dont get zero.
>  
> 
> 
> for instance for u-wind, i merged it using the following .gs script;
> 
> *script to merge data from multiple netcdf files, e.g. uwnd.1979.nc, uwnd.1980.nc, etc
> *edit start and end dates as needed
> *view output with Mergedncep_uwnd1000.ctl (edit tdef for longer times)
> 'reinit'
> iyear=1979
> 'set gxout fwrite'
> 'set fwrite F:\merged_uwnd1000.dat'
> while (iyear <= 2011)
> 'sdfopenF:\ncep_uwind\uwnd.'iyear'.nc'
> 'q file'
> res=sublin(result,5)
> tsize=subwrd(res,12)
> say tsize
> tim = 1
> 'set x 1 144'
> 'set y 1 73'
> while (tim <= tsize)
> 'set t 'tim 
> 'd uwnd(lev=1000)'
> tim=tim+1
> endwhile
> 'close 1'
> iyear=iyear+1
> endwhile 
> 'disable fwrite'
> 'reinit'
> 
> and ran it with the following .ctl script;
> 
> dset F:\merged_uwnd1000.dat
> title daily NCEP/DOE Reanalysis 2
> undef -9.99e+33
> xdef 144 linear 0 2.5
> ydef 73 linear -90 2.5
> zdef 17 levels 1000 925 850 700 600 500 400 300
> 250 200 150 100 70 50 30 20 10
> tdef 12053 linear 00Z01JAN1979 1dy
> vars 1
> uwnd 0 t,z,y,x Daily U-wind on Pressure Levels
> endvars
> 
> please were did i get it wrong?
> 
> ugo okoro
> cmsr, iap-cas,
> beijing.
> 
>   
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120618/f6a4d5f3/attachment-0003.html