Hi Dr. Wittenberg and Ferreters,
Thank you for your response. However, I am still having trouble getting the deepest depth. This is most likely due to my own inexperience with Ferret.
When I tried the commands that you sent, I got an error message:
**ERROR: command syntax: "," is illegal in 'temp,r=kend'
I'm not sure how to fix the problem and again, it may be only my inexperience that is a problem. It is possible that I was supposed to insert something where you place a ' in the commands, but I'm not sure what. Also if you (or anyone else) could give me an explanation of the commands and what they are doing I would really appreciate it.
Thanks,
Jaclyn Clement
-----Original Message-----
From: Andrew Wittenberg [mailto:atw@gfdl.noaa.gov]
Sent: Tuesday, January 28, 2003 4:10 PM
To: Clement, Jaclyn USA
Cc: Ferret Users
Subject: 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 |
+--------------------------------------------------------+
<<attachment: winmail.dat>>