[gradsusr] Query re working with lats4d...

[Filipkowski, John D]

Hmmn, is there a way of telling lats4d which grib table to use?

-----Original Message-----
From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Wesley Ebisuzaki
Sent: Wednesday, June 23, 2010 4:27 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] Query re working with lats4d...

John,

   I am not a user of lat4d but my understanding is that lats4d creates a custom grib table for every file.  wgrib doesn't know about the grib table, so it uses the operational NCEP table.

     Wesley Ebisuzaki


Filipkowski, John D wrote:
>
> Hi GrADS users!
>
>  - I'm looking for a link  or  reference do the following:
>
> For a GRIB-1 file, for a selected variable, I want to:
>   Smooth the data, output smoothed field to GRIB-1 format   ..
>
> Test data:
> Input files:       *TSTM.ctl*
> *                      TSTM.grib                      - GRIB-1 data
> file *
> *                      TSTM.grib.idx*
>
>
> *wgrib  TSTM.grib*
> 1:0:d=10060912:TSTM:kpds5=60:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sf
> c:0-0hr
> acc:NAve=0
>
> *Contents ot TSTM.ctl*
> dset  ^TSTM.grib
> index ^TSTM.grib.idx
> undef 9.999E+20
> title  TSTM.grib
> * produced by grib2ctl v0.9.12.5p32
> dtype grib 212
> pdef 185 129 lcc 35.0 -95.0 105 49 25 25 -95 40635 40635 xdef 181
> linear -140. 0.5 ydef  91 linear  15  0.5 tdef 1 linear 12Z09jun2010
> 1mo zdef 1 linear 1 1 vars 1
> TSTMsfc  60,1, 0      ** surface Thunderstorm probability [%]
> ENDVARS
>
> I've used the "lats4d" utility of gradsc/gradsnc as follows:
>> *gradsnc  -blc 'lats4d -v -q -i TSTM.ctl  -ftype ctl -func
>> smth9(tstmsfc) -o test  -format grads_grib'*
>>
>>
>> *Output:  (excerpted)...*
>> lats4d: vertical levels 1
>> lats4d: using PRESSURE for vertical coordinate ...
>> lats4d: Function expression smth9(tstmsfc)
>> lats4d: surface  variables  tstmsfc
>> lats4d: writing to grads_grib file test.grb on 12Z09JUN2010  ...
>>
>> *wgrib  test.grb*
>> Undefined parameter table (center 100-2 table 128), using NCEP-opn
>> 1:0:d=10060912:*PRES*:kpds5=1:kpgs6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:
>> sfc:anl:Nave=0
>>
>>
>> I don't understand why the data written to the new GRIB file is now
>> called "PRES"...
>>
>> Can anyone shed some light on this?
>>
>> Cheers,
>>
>>
>>
>
>
>
> The information contained in this e-mail is intended only for the
> individual or entity to whom it is addressed.
> Its contents (including any attachments) may contain confidential
> and/or privileged information.
> If you are not an intended recipient you must not use, disclose,
> disseminate, copy or print its contents.
> If you receive this e-mail in error, please notify the sender by reply
> e-mail and delete and destroy the message.
> ----------------------------------------------------------------------
> --
>
> _______________________________________________
> 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

The information contained in this e-mail is intended only for the individual or entity to whom it is addressed.
Its contents (including any attachments) may contain confidential and/or privileged information.
If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents.
If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message.