[gradsusr] Regridding MERRA data with lats4d

[Bill Reilly]

Maura:

There's another recipe that may help you:

http://cookbooks.opengrads.org/index.php?title=Recipe-018:_Automatic_regridding_of_MERRA_data_with_Lats4D

Bill

--
Bill Reilly
(+34) 686-110-748

bill at passageweather.com
www.passageweather.com

bill at brwebworks.com
www.brwebworks.com


On 12/6/2011 3:21 PM, Tokay, Maura M. (GSFC-618.0)[SCIENCE SYSTEMS 
APPLICATIONS] wrote:
> Hi  GrADS users,
>
> I need to regrid MERRA data from .5x.6666 to .1x.1. My file is within the following boundaries: 27N-30N and 84E-87.333E. I read the documentation on how to use lats4d but it is not working and I am not sure what I might be missing. I have tried two different approaches with and without the de option:
>
> lats4d.sh -i MERRA300.prod.assim.inst6_3d_ana_Nv.20011012.SUB.nc -o MERRA300.inst6_3d_ana_Nv.20011012 -format netcdf -func "re(@,33,linear,84,0.1,30,linear,27,0.1,bl)"
>
> lats4d.sh -i MERRA300.prod.assim.inst6_3d_ana_Nv.20011012.SUB.nc -o MERRA300.inst6_3d_ana_Nv.20011012 -format netcdf -de MERRA300.prod.assim.inst6_3d_ana_Nv -func "re(@,33,linear,84,0.1,30,linear,27,0.1,bl)"
>
> My control file MERRA300.prod.assim.inst6_3d_ana_Nv.ctl is described below:
>
> dset nofile
> title MERRA data
> options template
> undef 1e+15
> xdef 33 linear  84 0.1
> ydef 30 linear 27 0.1
> zdef 72 levels 976.62 961.913 947.258 932.602 917.946 903.289 888.632 873.975 859.316 844.658 829.999 815.34 798.236 776.243 751.808 727.369 702.925 678.481 647.923 611.244 574.557 537.86 501.148 464.421 427.673 390.897 354.091 311.373 265.745 226.234 192
> 485 163.662 139.115 118.25 100.514 85.439 72.5579 61.4957 52.0159 43.9097 36.9927 31.0889 26.0491 21.761 18.1244 15.0503 12.4602 10.2849 8.45639 6.91832 5.6318 4.56169 3.6765 2.94832 2.35259 1.86788 1.47565 1.15998 0.907287 0.705957 0.546293 0.420424 0.32
> 784 0.244938 0.185422 0.139599 0.104524 0.0776725 0.0567925 0.0401425 0.02635 0.015
> tdef 4 linear 00Z12OCT2001 6hr
> vars 5
> ps  0  t,y,x  Surface pressure
> t  72  t,z,y,x  Air Temperature
> u  72  t,z,y,x  Eastward wind component
> v  72  t,z,y,x  Northward wind component
> qv  72  t,z,y,x  Specific humidity
> endvars
>
> Does anybody have any idea of what I might be missing?
>
> Thanks,
> Maura Tokay
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>