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 >