This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
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
>