NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.

To learn about what's going on, see About the Archive Site.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20000106: decoding problem in dchrly???

>From: "David J. Knight" <address@hidden>
>Organization: .
>Keywords: 200001061655.JAA27619

>Hi,
>     There seems to be a problem in decoding many South American 
>metar observations using dchrly in gempak.
>
>Specifically it does not seem to be decoding the altimeter 
>setting when given in Hpa.
>
>For example, an observation like:
>SBSN 061230Z 25002KT 9999 FEW016 SCT100 28/24 Q1010=
>produces a missing value for alti in the surface file. As I
>understand it the Q group should be the altimeter setting in Hpa.
>
>The apparently relevant code in:
>gempak5.4/src/bridge/metar/parse_metar.c
>is:
>if( Mptr->A_altstng )
>   *alti = (float)Mptr->inches_altstng;
>else
>   {
>   if( Mptr->Q_altstng )
>      *alti = (float)Mptr->hectoPasc_altstng * 1013.25 / 29.921;
>   }
>Which appears to suggest it should be working properly. (i.e. 
>it looks like it converts the Q group to alti in inches of mercury)
>
>Suggestions - help ?
>
>Thanks
>David
>
>David Knight
>Department of Earth and Atmospheric Sciences   Tel: (518)-442-4204
>SUNYA   ES-228                                 Fax: (518)-442-4494
>Albany, NY  12222                              Email: address@hidden
>


Dave,

I fixed this in a message to Robert Mullenax at:
http://www.unidata.ucar.edu/glimpse/gempak/3160

the conversion should be:
*alti = (float)Mptr->hectoPasc_altstng * 29.921 / 1013.25;

eg, flip the conversion factors.

Recompile the bridge/metar directroy to update bridge.a, then recompile
and link $GEMPAKHOME/src/programs/dc/dchrly

Steve Chiswell