Hi -
Some things to look into: Are the units of the the variables in the
datset correct for the calculation? Does the original Z axis have
units of meters? Did you look at the intermediate variables salt_z
and temp_z to make sure they look correct?
-Ansley
On 1/5/2012 5:20 AM, Mahesh Shinde wrote:
Dear friends,
I am trying to plot MLD in ROMS output. I have plot it for single
day using density criteria ..please have a look at fig attached
herewith.
a script is also attached here....there is something wrong in the
plot......
!the script
can data/all
SET MEMORY/SIZE=100
use "../his_rmed16_grd_gls_2506.nc"
DEFINE AXIS/Z/UNITS=meters/DEPTH znew = \
{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 depth =zeta+z[gz=s_rho]*h
let zaxnew = z[gz=znew]
let Zdepth = depth * Z[g=GPV16] ! GPV16 comes
from the sh grid var
let mask_temp = IF mask_rho EQ 1 THEN temp
let mask_salt = IF mask_rho EQ 1 THEN salt
let temp_z = ZAXREPLACE(temp, Zdepth, zaxnew)
let salt_z = ZAXREPLACE(salt, Zdepth, zaxnew)
!-----------------------------------------------------------------
!based on density criteria......
go mld_dens salt_z,temp_z, 0.5
fill mld_d, lon_rho, lat_rho
go land
frame/file=r_mld_dens.gif
!end
Please tell why i am getting such a plot ...
thanking you
regards
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 : mashinde2k9
----------------------------------------------------------------------------------------------------------------------
|