[ferret_users] about mld plots ....

[Mahesh Shinde <mashinde05@xxxxxxxxx>]

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

Attachment: r_mld_dens.gif
Description: GIF image