[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]
Andrew
+--------------------------------------------------------+
| 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