[gradsusr] Lats4d Changing Variable and Time
Fishel, Greg
gfishel at wral.com
Tue Dec 4 13:43:05 EST 2012
Lee,
Sorry I did overlook your email yesterday. However, I saved your file and
added it to my lats4d command. Unfortunately, I am getting this:
LATS error (lats_parmtab): LATS (latsint.c) --> Cannot open parameter file
greg.lats.table
No such file or directory
When I put a @ in front of my table name, I don't get the error but the
parameter name still changes to PRES. Boy am I confused.
Greg
On 12/4/12 12:51 PM, "Lee Byerle" <lbyerle at yahoo.com> wrote:
>Greg
>Are you able to view the output correctly in GrADs?
>
>Did you receive the parameter table file example? I mailed yesterday. Try
>adding to your lats4d command the option that points to an edited version
>of that parameter table that includes snodsfc. (see my previous post).
>
>-table table-name
>
>This should label the data properly as snow depth. Still to be addressed
>is the issue of specifying your forecast times correctly...
>
>If you can view the output in grads, you can use g2grb in the session to
>control grib header info (
>http://www.cpc.ncep.noaa.gov/products/wesley/g2grb.html ).
>
>Lee
>
>
>Sent from my iPhone
>
>On Dec 4, 2012, at 11:05 AM, "Fishel, Greg" <gfishel at wral.com> wrote:
>
>> OK, I have a cfs grib2 file, as evidenced by the following inventory
>> produced by wgrib2:
>>
>> 1:0:d=2012120312:SNOD:surface:6 hour fcst:
>> 2:33208:d=2012120312:SNOD:surface:12 hour fcst:
>> 3:66380:d=2012120312:SNOD:surface:18 hour fcst:
>> 4:99643:d=2012120312:SNOD:surface:24 hour fcst:
>> 5:132801:d=2012120312:SNOD:surface:30 hour fcst:
>> 6:165539:d=2012120312:SNOD:surface:36 hour fcst:
>> 7:198305:d=2012120312:SNOD:surface:42 hour fcst:
>> 8:231208:d=2012120312:SNOD:surface:48 hour fcst:
>>
>>
>>-------------------------------------------------------------------------
>>--
>>
>>-------------------------------------------------------------------------
>>--
>> --------
>>
>> The data is on a gaussian grid, again as shown by wgrib2:
>>
>> 1:0:vt=2012120318:surface:6 hour fcst:SNOD Snow Depth [m]:
>> ndata=72960:undef=48334:mean=0.442751:min=0:max=1.99616
>> grid_template=40:winds(N/S):
>> Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN
>> number of latitudes between pole-equator=95 #points=72960
>> lat 89.277000 to -89.277000
>> lon 0.000000 to 359.062000 by 0.938000
>>
>>
>>-------------------------------------------------------------------------
>>--
>>
>>-------------------------------------------------------------------------
>>--
>> --------
>>
>> I then create a control file for the gfs file:
>>
>>
>> dset ^gfs.snow.grib2
>> index ^gfs.snow.grib2.idx
>> undef 9.999E+20
>> title gfs.snow.grib2
>> * produced by g2ctl v0.0.8.8
>> * command line options: -verf gfs.snow.grib2
>> * griddef=1:0:(720 x 361):grid_template=0:winds(N/S): lat-lon grid:(720
>>x
>> 361) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to
>> -90.000000 by 0.500000 lon 0.000000 to 359.500000 by 0.500000
>> #points=259920:winds(N/S)
>> dtype grib2
>> ydef 361 linear -90.000000 0.5
>> xdef 720 linear 0.000000 0.500000
>> tdef 8 linear 18Z03dec2012 6hr
>> zdef 1 linear 1 1
>> vars 1
>> SNODsfc 0,1,0 0,1,11 ** surface Snow Depth [m]
>> ENDVARS
>>
>>
>>-------------------------------------------------------------------------
>>--
>>
>>-------------------------------------------------------------------------
>>--
>> --------
>>
>> And then a control file for the cfs file
>>
>> dset ^cfs.snow.grib2
>>
>>
>> index ^cfs.snow.grib2.idx
>> undef 9.999E+20
>> title cfs.snow.grib2
>> * produced by g2ctl v0.0.8.8
>> * command line options: -verf cfs.snow.grib2
>> * griddef=1:0:(384 x 190):grid_template=40:winds(N/S): Gaussian grid:
>>(384
>> x 190) units 1e-06 input WE:NS output WE:SN number of latitudes between
>> pole-equator=95 #points=72960 lat 89.277000 to -89.277000 lon 0.000000
>>to
>> 359.062000 by 0.938000:winds(
>> dtype grib2
>> * ny=190 nlat=190 lat0=-89.277000 lat1=89.277000
>> xdef 384 linear 0.000000 0.937498694516971
>> ydef 190 levels
>> -89.277 -88.340 -87.397 -86.454 -85.509 -84.565 -83.620 -82.676 -81.731
>> -80.786 -79.841 -78.897 -77.952 -77.007
>> -76.062 -75.117 -74.173 -73.228 -72.283 -71.338 -70.393 -69.448 -68.503
>> -67.559 -66.614 -65.669 -64.724 -63.779
>> -62.834 -61.889 -60.945 -60.000 -59.055 -58.110 -57.165 -56.220 -55.275
>> -54.330 -53.386 -52.441 -51.496 -50.551
>> -49.606 -48.661 -47.716 -46.771 -45.827 -44.882 -43.937 -42.992 -42.047
>> -41.102 -40.157 -39.212 -38.268 -37.323
>> -36.378 -35.433 -34.488 -33.543 -32.598 -31.653 -30.709 -29.764 -28.819
>> -27.874 -26.929 -25.984 -25.039 -24.094
>> -23.150 -22.205 -21.260 -20.315 -19.370 -18.425 -17.480 -16.535 -15.590
>> -14.646 -13.701 -12.756 -11.811 -10.866
>> -9.921 -8.976 -8.031 -7.087 -6.142 -5.197 -4.252 -3.307 -2.362
>> -1.417 -0.472 0.472 1.417 2.362
>> 3.307 4.252 5.197 6.142 7.087 8.031 8.976 9.921 10.866
>> 11.811 12.756 13.701 14.646 15.590
>> 16.535 17.480 18.425 19.370 20.315 21.260 22.205 23.150 24.094
>> 25.039 25.984 26.929 27.874 28.819
>> 29.764 30.709 31.653 32.598 33.543 34.488 35.433 36.378 37.323
>> 38.268 39.212 40.157 41.102 42.047
>> 42.992 43.937 44.882 45.827 46.771 47.716 48.661 49.606 50.551
>> 51.496 52.441 53.386 54.330 55.275
>> 56.220 57.165 58.110 59.055 60.000 60.945 61.889 62.834 63.779
>> 64.724 65.669 66.614 67.559 68.503
>> 69.448 70.393 71.338 72.283 73.228 74.173 75.117 76.062 77.007
>> 77.952 78.897 79.841 80.786 81.731
>> 82.676 83.620 84.565 85.509 86.454 87.397 88.340 89.277
>> tdef 8 linear 18Z03dec2012 6hr
>> zdef 1 linear 1 1
>> vars 1
>> SNODsfc 0,1,0 0,1,11 ** surface Snow Depth [m]
>> ENDVARS
>>
>>
>>-------------------------------------------------------------------------
>>--
>>
>>-------------------------------------------------------------------------
>>--
>> --------
>>
>> I want to convert this to a standard GFS .5 degree grid. I use the
>> following lats4d command:
>>
>>
>>
>> lats4d.sh -i cfs.snow.grib2.ctl -precision 16 -vars snodsfc -center NCEP
>> -time 18z3dec12 12z5dec12 1 -func
>> 're(@,720,linear,0,.5,361,linear,-90,.5,bl)' -de gfs.snow.grib2.ctl
>> -format grads_grib -o cfs.snow.newgrid.grib2
>>
>> The grib file is created but the inventory indicates the parameter has
>> changed from snodsfc to pres and all times are analyses, as shown by
>>wgrib:
>>
>> 1:0:d=12120318:PRES:sfc:anl:NAve=0
>> 2:189878:d=12120400:PRES:sfc:anl:NAve=0
>> 3:379756:d=12120406:PRES:sfc:anl:NAve=0
>> 4:569634:d=12120412:PRES:sfc:anl:NAve=0
>> 5:759512:d=12120418:PRES:sfc:anl:NAve=0
>> 6:949390:d=12120500:PRES:sfc:anl:NAve=0
>> 7:1139268:d=12120506:PRES:sfc:anl:NAve=0
>> 8:1329146:d=12120512:PRES:sfc:anl:NAve=0
>>
>> I have googled until I'm blue in the face. Any suggestions would be
>> extremely welcome!
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
>_______________________________________________
>gradsusr mailing list
>gradsusr at gradsusr.org
>http://gradsusr.org/mailman/listinfo/gradsusr
More information about the gradsusr
mailing list