[gradsusr] Problem using lats4d.sh to clip a file

[Lee Byerle]

Hi Maura,
i have note tried to reproduce your result, but have clipped regions using lats4d and a re-grid function, and output to grib.  Maybe if you specify the grid using "-de" and "-func re().." it will return the same result.  From the opengrads lats4d documentation, I include a dummy ctl file that gives the dimensions of the clipped file.  Here is how I list the lats4d command using an example for a clipped region:

lats4d... .... -de clipped_ncep1x1 -func 're(@,20,linear,0,1,20,linear,35,1,bl)' -format grads_grib...etc


clipped_ncep1x1.ctl is:

DSET nofile
title Template GrADS for regridding
options template
undef 1e+20
xdef 20 linear   0 1.
ydef 20 linear 35  1.
zdef   1 levels  1000
tdef 1 linear 0Z1jan1900 1dy
vars 1
var       0 0        generic sfc variable
endvars

good luck,
Lee

On Nov 7, 2012, at 12:42 PM, Tokay, Maura M. (GSFC-610.4)[SCIENCE SYSTEMS AND APPLICATIONS INC] wrote:

> Hello -,
> 
> I am using the lats4d.sh  to clip the MERRA file that have the physical constants. The file is global but I need a small region and I want to convert from HDF to netcdf format. My command line is
> 
> lats4d.sh -i MERRA300.prod.assim.const_2d_slv_Nx.ctl -lon 84.0 87.0 -lat 27.0 30.0 -o merra.prod.assim.const.2d.slv.Nx.nc
> 
> Because I am having some problems with temperature when I downscale the MERRA temperature data using the micromet technique one of suggestions I got to troubleshoot the temperature problem was to do the downscale to a larger area and see if the results would be the same. I clipped the same original MERRA file using the same command line for a different area
> 
> lats4d.sh -i MERRA300.prod.assim.const_2d_slv_Nx.ctl -lon 80.0 90.0 -lat 26.0 33.0 -o merra.prod.assim.const.2d.slv.Nx.nc
> 
> To make sure the data was correct since I saw a significant difference after I applied the downscale procedure, I loaded my 2 clipped MERRA files on ArcMap and when I overlaped them and checked the surface geopotential value for the same cell I get two distinct results. I don't understand how this can be possible and I would like to know if anybody has any idea of what I might be doing worng. Any help will be greatly appreciated.
> 
> Thanks,
> Maura
> 
> Maura Tokay, PMP - SSAI, Inc
> Lead Software Programmer
> NASA Goddard Space Flight Center
> Office of Applied Sciences
> Phone: 301-614-6334
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr