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