[Thread Prev][Thread Next][Index]

Re: deepest depth

On Tue, 28 Jan 2003, Clement, Jaclyn USA wrote:
> I'm trying to shade a variable (e.g. temp) at the deepest depth (highest
> k level) not including missing values.  I've been able to define and
> shade a variable (kmax) by using the following: let kmax3D = k[g=temp] +
> 0*temp ! an XYZ field of K masked by bottom let kmax =
> kmax3D[k=1:33@max] ! if 33 is number of vertical levels This variable
> (kmax) correctly represents the deepest depth at each x,y data point,
> but how do I tell Ferret to use kmax as the k when shading temp?

Hi Jaclyn,

The @WEQ ("weighted equal") transformation does the trick (see the User's
Guide, Ch3 Sec2.4.27, Ex. 5):

let kvals = k[g=temp] + 0*temp
let kzero = kvals - kvals[k=1:`temp,r=kend`@max]
let integrand = temp * kzero[k=@weq:0]
let temp_bottom = integrand[k=1:`temp,r=kend`@sum]


|   Dr. Andrew T. Wittenberg   |        GFDL/NOAA        |
|      atw@gfdl.noaa.gov       |      Princeton, NJ      |

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement