gdvint and dewpoint temperatures question

• To: gembud@xxxxxxxxxxxxxxxx
• Subject: gdvint and dewpoint temperatures question
• From: David Ovens <ovens@xxxxxxxxxxxxxxxxxxxx>
• Date: Tue, 15 Nov 2005 16:54:42 -0800

Gembuds,

I have converted some sigma level data (from MM5 output) into pressure
level data using gdvint as follows:

GDFILE = outfile2
GDOUTF = outfile2
GDATTIM = LAST
GLEVEL = 1000-900-10
GVCORD = sgma:100/pres
MAXGRD = 600
GAREA = grid
VCOORD = pres;none
$respond=yes
run
GLEVEL = 875-750-25
run
GLEVEL = 700-100-50
run
GLEVEL = 1050-1010-10
run

This mostly works great with SPFH (specific humidity) appearing as the
moisture value -- the original sigma levels had TMPK and DWPK at all

In gdprof, I can plot TMPK (UWND, VWND, PRES, etc.) all the way down

YAXIS    = 1050/400/50/;1
The text output snippet for TMPK is

                 PRES                  TMPK
               1020.28              278.679
               1000.00              277.595
                990.00              277.818

but for DWPK, it stops at 1000 mb

                PRES                  DWPK
               1000.00              270.792
                990.00              269.367

I think there is some problem in calculating DWPK from TMPK and SPFH
below 1000 mb and that is why gdprof does not want to display this.  I
think this is probably a failing in gdvint since I do not see
DWPK@0%PRES -- I only see TMPK, HGHT, PRES, UWND, UREL, SPFH, WWND
@0 %PRES -- even though I had TMPK and DWPK at 0%NONE.

Has anyone figured out how to work around this either in gdvint or

David
--
David Ovens              e-mail: ovens@xxxxxxxxxxxxxxxxxxxx
Research Meteorologist    phone: (206) 685-8108
Dept of Atm. Sciences      plan: Real-time MM5 forecasting for the
Box 351640                        Pacific Northwest
University of Washington          http://www.atmos.washington.edu/mm5rt
Seattle, WA  98195               Weather Graphics and Loops
                                 http://www.atmos.washington.edu/~ovens/loops