[gradsusr] safe data into a NetCDF format

[Arlindo da Silva]

On Mon, Nov 28, 2011 at 1:38 PM, Muyin Wang <Muyin.Wang at noaa.gov> wrote:

> Dear Grads Users:
> I have a 40 year daily file from model outputs. Now I need to extract
> variables for limited domain from the file, and connect them into a
> single 40 year NetCDF file.
> Does any one know the trick to do this? T
> Thank you very much!
>

You can use lats4d and have this accomplished with a single command line.
The syntax is something like this:

% lats4d.sh -i yourfile.ctl -o singlefile -lat 30 40 -lon -120 -60 -v

or from inside grads

ga-> lats4d -i yourfile.ctl -o singlefile -lat 30 40 -lon -120 -60 -v


This will create a file called "singlefile.nc". There are also other
options to select variables, levels, etc. For a full list of options:

% lats4d.sh -h

Now, you can use lats4d with v1.9b4 (there is an older version version of
lats4d.gs in the script
library<http://grads.iges.org/grads/gadoc/gadocindex.html>that you can
call from inside grads 1.9), but I'd strongly recommend that
you grab the latest opengrads bundle which includes lats4d.sh:

   http://opengrads.org/wiki/index.php?title=Installing_the_OpenGrADS_Bundle

      Arlindo



> Muyin
>
>
> Here is what it looks like in the .ctl file :
>
> DSET      ^Data/2001/1dy.grd
> TITLE     20C3M+A1B
> OPTIONS   big_endian template yrev
> UNDEF     -999
> XDEF       128 LINEAR       0       2.8125
> YDEF       64 LEVELS
> -87.864 -85.097 -82.313 -79.526 -76.737 -73.948 -71.158 -68.368 -65.578
> -62.787 -59.997 -57.207 -54.416 -51.626 -48.835 -46.045 -43.254 -40.464
> -37.673 -34.883 -32.092 -29.301 -26.511 -23.720 -20.930 -18.139 -15.348
> -12.558 -9.7671 -6.9765 -4.1859 -1.3953 1.3953 4.1859 6.9765 9.7671 12.558
> 15.348 18.139 20.930 23.720 26.511 29.301 32.092 34.883 37.673 40.464
> 43.254
> 46.045 48.835 51.626 54.416 57.207 59.997 62.787 65.578 68.368 71.158
> 73.948
> 76.737 79.526 82.313 85.097 87.864
> ZDEF         1 LINEAR 1 1
> TDEF     14976 LINEAR  12:00Z01JAN2000     1dy
> VARS  8
> T2           0 99 2m temperature [K]
> q2           0 99 2m specific humidity [kg/kg]
> slp          0 99 sea level pressure [hPa]
> u10          0 99 10m zonal wind [m/s]
> v10          0 99 10m meridional wind [m/s]
> ssrd         0 99 surf.short.down [W/m**2]
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>



-- 
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111128/7e754716/attachment-0003.html