[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