Hi -
What is the error that you're seeing? You may be able to see what's
happening by using Ferret to list and give information about the
variables you've defined. Commands like
show grid zt
show grid t_nd
and various plots showing slices of these variables, may tell you if
there's an error in the definitions.
(Note by the way, that the argument to SET
MEMORY is given in mWords, so the value you've given for that
is probably larger than the memory on your machine. )
-Ansley
On 6/7/2011 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
----------------------------------------------------------------------------------------------------------------------
|