20010518: 20010517: cape, bulk richardson

Diana,

Since you sent me the STNDEX values for:
STID = TOP           STNM =    72456   TIME = 930723/0000

Just run snlist with:
SNPARM=dset
LEVELS=all
VCOORD=pres
AREA=@72456
DATTIM=930723/0000
SNFILE=1993_72456.gem

The SNPARM=dset will list out the all the data values for
each pressure level.
With that data, I can see why CAPE is arriving at a value of 0
when you expect something much larger.

Some possibilities are that the sounding prematurely terminated
so that there aren't enough data levels to calculate integrated
parameters, or that there are problems with duplicate levels.

Steve Chiswell
Unidata User Support



>From: address@hidden
>Organization: UCAR/Unidata
>Keywords: 200105182156.f4ILuZp08228

>hmm, i haven't been able to get a list of all levels.  I remember the first
>script I ran gave me all of that info, but I changed some part of the script
>to condense the output.  Can you tell what is missing in the script, so 
>I can get the info you need?  
>
>#!/bin/csh -f
>
>set STNS=  ("72456")
>
>@ YEAR = 97
>
>while ($YEAR <= 97)
>  
> foreach STN ($STNS)
>       snlist << EOF
>               snfile = 19${YEAR}_${STN}.gem
>               dattim = 970621/0000
>               output = f/stat-run/out_19${YEAR}_${STN}.dat
>               snparm = 
>               stndex = lift; show; totl; kinx; swet; cape; brch
>               levels = all
>               vcoord = pres
>               mrgdat = y
>               area   = @${STN}
>               r
>               
>               e
>EOF
>       end
>       @ YEAR = $YEAR + 1
>end
>echo "Finished"
>
>
>
>
>thanks, Diana
>