[Thread Prev][Thread Next][Index]
Misuse (?) of "MAX" command?
I'm trying to generate some nice timeseries PLOTs of several fields plus a set
of fixed values, and I want to set the VLIMITS properly so all the quantities
are plotted, ie, I want all of the fixed value lines on the plot as well.
I used the following commands:
[...]
LET fix_lo = 17.82 ; LET fix_hi = 18.39
LET var6 = var[d=data6.nc] ; LET mx6 = var6[l=@max] ; LET min6 = var6[l=@min]
LET var7 = var[d=data7.nc] ; LET mx7 = var7[l=@max] ; LET min7 = var7[l=@min]
LET var8 = var[d=data8.nc] ; LET mx8 = var8[l=@max] ; LET min8 = var8[l=@min]
LET var9 = var[d=data9.nc] ; LET mx9 = var9[l=@max] ; LET min9 = var9[l=@min]
!
LET varmx = max(max(max(mx6,mx7),mx8),mx9)
LET varmn = min(min(min(mn6,mn7),mn8),mn9)
LET vltop = max(varmx,fix_hi)*1.01
LET vlbot = min(varmn,fix_lo)*0.99
!
PLOT/VL=`vlbot`:`vltop` var6,var7,var8,var9
[...]
Should work, right?
Well, doing a LIST gives the following:
LIST mx6,mx7,mx8,mx9,varmx,vltop
MX6 MX7 MX8 MX9 VARMX VLTOP
I / *: 19.09 22.09 20.09 20.16 19.16 19.35
LIST mn6,mn7,mn8,mn9,varmn,vlbot
MN6 MN7 MN8 MN9 VARMN VLBOT
I / *: 16.85 17.28 17.55 17.24 16.85 16.68
How come VARMX is wrong?
--
/\ Gary Strand (303) 497-1336 NCAR CG2260
\_][ www.cgd.ucar.edu/ccr/strandwg 3080 Center Green Dr
\___strandwg@ucar.edu Boulder, Colorado, USA 80301
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement