Re: [ferret_users] computation and plotting of MLD from ROMS output

["William S. Kessler" <william.s.kessler@xxxxxxxx>]

One problem is your use of "z" in "let zt = (z + 0.3) * bb". You must  
tell Ferret which "z" this is (on what grid, as you do later:  
"z[gz=z100]"). You can't use x,y,z,t or i,j,k,l as simple variables.

In messages to the list, it is hard for people to help you unless you  
say specifically what kind of error was produced. Similarly, without  
knowing what your variable "h" is, it is hard to move forward.

And as a general tip, when you run into a problem, use PLOT or SHADE  
or STAT or LIST (in a limited range) on the intermediate variables to  
examine them. That will usually tell you where the problem lies. Then  
if there is still an error, you will be able to ask us a more specific  
question.

Billy K

On 07 Jun 2011, at 9:03 AM, Mahesh Shinde wrote:

> Dear ferret users,
>
> I am new to ferret. I want to compute and plot MLD from ROMS  
> output . I tried but it gives me an error.
> ROMS has sigma vertical levels. Here I am sending you a script that  
> I wrote ....
>
> cancel data/all
> cancel viewport
> !SET MODE META
> PPL SHASET RESET
>
> SET MEMORY/SIZE=2000000000
>
>
> use  "avg_rmed16_grd_e_0178.nc"
> !--------------------------------------------------------------------------------------------------------------------------
> !DEFINE AXIS/Z/UNITS=meters/DEPTH z100 = \
> ! 
> {0,2,4,6,8,10,12,15,20,25,30,35,40,45,50,60,70,80,90,100,125,150,200,250 
> ,\
> ! 
> 300,350,400,500,600,700,800,900,1000,1250,1500,2000,2500,3000,4000,5000 
> }
> !----------------------------------------------------------------------------------------------------------------------------
>   let bb = h + 0 * TEMP[l=1]
>   let zt = (z + 0.3) * bb
>   define axis/z=0:4000:10/depth z100
>   let znew = Z[gz=z100]
>   let t_nd = ZAXREPLACE(TEMP,zt,znew)
>
>    go mld_temp t_nd 0.2
>
>
> Please tell me where it is going wrong and why?
>
>
> Thanks
>
> Mahesh
>
>
>
>
>
> --
> ----------------------------------------------------------------------------------------------------------------------
> (Permanent Affiliaiton)
> Mr. Mahesh Shinde
> Climate and Global Modeling Division,
> Indian Institute of Tropical Meteorology,
> NCL POST, Dr. Homi Bhabha Road,
> Pashan, Pune-411008, INDIA
> email: samurai@xxxxxxxxxxxxxx
> cell: +91-9822050959
> skype : mshinde2k9
> ----------------------------------------------------------------------------------------------------------------------

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
William S. Kessler
NOAA / Pacific Marine Environmental Laboratory
7600 Sand Point Way NE
Seattle WA 98115 USA

william.s.kessler@xxxxxxxx
Tel: 206-526-6221
Fax: 206-526-6744
Web: http://www.pmel.noaa.gov/people/kessler