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

script question

Steve,

I have been tinkering around with the script you gave me
a while back. I have attached it below. It seems to run
fine with no errors however, I notice when I display the
new calculated grdnam (stbc) in gdprof it outputs a
level of 975mb for den. Using garp and looking at
the pres@0%none (station pressure), I can see that the
station pressure for den is less than 900mb. I would think
that the first level displayed in the output file from gdprof
would less than the station pressure. I am thinking about
this correctly? I have attached my gdprof.csh and its output
as well.

Thank you again,
Matt



SCRIPT

#!/bin/csh

set GDFILE=/data1/gempak/grids/99061412_eta.gem
#set FHRS=("f00" "f06" "f12" "f18" "f24" "f30" "f36" "f42" "f48")

  $GEMDATA/gddiag << EOF1
  gdfile  = $GDFILE
  gdoutf  = $GDFILE
  gdattim = f00
  gfunc   = pres
  glevel  = 0
  gvcord  = none
  grdnam  = pres@0%pres
  r

EOF1

set GLEV1 = 1000
while($GLEV1 > 100)
  @ GLEV2 = $GLEV1 - 50
  @ GLEVA = $GLEV1 + $GLEV2
  @ GLEVA = $GLEVA / 2

  $GEMEXE/gddiag << EOF2
  gdfile  = $GDFILE
  gdoutf  = $GDFILE
  gfunc   = stab(tmpc)
  gdattim = f00
  glevel  = ${GLEV1}:${GLEV2}
  gvcord  = pres
  grdnam  = stbc@${GLEVA}
  r

  e
EOF2
  @ GLEV1 = $GLEV1 - 50

end   #end while GLEV
#end  #foreach FHR

$GEMEXE/gpend
rm *.nts


GDPROF
#!/bin/csh -v


   $GEMEXE/gdprof <<EOF
   gpoint  = den
   gdattim = f00
   gvcord  = pres
   gfunc   = stbc
   gvect   = wnd
   gdfile  = /data1/gempak/grids/99072012_eta.gem
   line    = 5
   marker  = 0
   border  = 1
   ptype   = log
   scale   = 0
   xaxis   = -15/10/5
   yaxis   = 1000/100/100
   wind    = bm6//2
   refvec  =
   winpos  = 1
   filter  = 0
   title   = 1/3/~ @ Stability Profile     ! 0
   panel   = 0
   clear   = yes
   text    = 1
   device  = xw
   output  = f/dumb
   thtaln  =
   thteln  =
   mixrln  =
   r

   exit
EOF

$GEMEXE/gpend
rm *.nts


OUTPUT FROM GDPROF

                  PRES                  STBC
                 975.00               -6.470
                 925.00               -5.993
                 875.00               -7.307
                 825.00               -3.994
                 775.00               -4.346
                 725.00               -5.697
                 675.00               -5.752
                 625.00               -6.031
                 575.00               -6.694
                 525.00               -6.726
                 475.00               -6.902
                 425.00               -6.533
                 375.00               -6.664
                 325.00               -7.174
                 275.00               -7.095
                 225.00               -7.937
                 175.00               -6.247
                 125.00               -0.153



From address@hidden  Mon Jul 26 14:33:30 1999
Received: from screen.dri.edu (screen.dri.edu [134.197.100.251])
        by unidata.ucar.edu (8.8.8/8.8.8) with ESMTP id OAA01567
        for <address@hidden>; Mon, 26 Jul 1999 14:33:25 -0600 (MDT)
Keywords: 199907262033.OAA01567
Received: (from fw_mx@localhost)
        by screen.dri.edu (8.9.1b+Sun/8.9.1) id NAA00933
        for <address@hidden>; Mon, 26 Jul 1999 13:33:24 -0700 (PDT)
Received: from rotor.dri.edu(10.10.9.96) by screen.dri.edu via smap (V2.1.dri)
        id xma000930; Mon, 26 Jul 99 13:32:28 -0700
Received: from dri.edu (fropa.dri.edu [10.10.9.93])
        by rotor.dri.edu (8.9.1b+Sun/8.9.1) with ESMTP id NAA24585
        for <address@hidden>; Mon, 26 Jul 1999 13:32:54 -0700 (PDT)
Message-ID: <address@hidden>
Date: Mon, 26 Jul 1999 13:32:41 -0800
From: Matt Fearon <address@hidden>
Organization: Desert Research Institute
X-Mailer: Mozilla 4.5 (Macintosh; U; PPC)
X-Accept-Language: en
MIME-Version: 1.0
To: address@hidden
Subject: script assistance
Content-Type: text/plain; charset=us-ascii; x-mac-type="54455854"; 
x-mac-creator="4D4F5353"
Content-Transfer-Encoding: 7bit

Steve,

You mentioned in a prior email if I explained to you how
I wanted to calculate mixing depth, you may be able to
give me some assistance with scripting. An idea I have
been considering is comparing two forecast hours.
I am using strictly the eta model. My thoughts were
if I could somehow take the 2m hagl tmpk say at the
f06 and raise it adiabatically and as that tmpk is being
raised, I would like to compare those temp values to the
tmpk values from f00 at various altitudes above ground.
If the lifted adiabatic air from f06 intersects or becomes
cooler than the tmpk from the f00, I could say that
intersection (altitude) is a tentative mixing height.
Would you know of way do this calculation for a chosen
grid area? Let me know when you have a chance.

Thank you,
Matt Fearon
address@hidden