[ferret_users] computation and plotting of MLD from ROMS output

[Mahesh Shinde <mashinde05@xxxxxxxxx>]

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
----------------------------------------------------------------------------------------------------------------------