[gradsusr] Lats4d Changing Variable and Time
Lee Byerle
lbyerle at yahoo.com
Tue Dec 4 12:51:58 EST 2012
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
More information about the gradsusr
mailing list