Dear Ferreters
G'day
I am performing a computation of MLD using the potential density criterion i.e. potential density referenced to the surface, σ0, exceeds by a threshold of 0.03 kg m−3 the density of the water at 10 m,
σ0(z = −H) = σ0(z = −10 m) + 0.03 kg m−3, with H the mixed-layer depth.
I masked the potential density as
let dens_10m = density[z=10]
let threshold = dens_10m + 0.03
let mask = if density le threshold then density
now I want to take the deepest level of the valid density .i.e mask variable as the MLD. So how to extract this level as a variable of MLD?
let zed = z[gz=density] + 0*density
let zp = zed*mask
mld = zp[z=@max]
I like to know if this is the correct way to compute the MLD?
Cheers, Saurabh
--
REGARDS
Saurabh Rathore
Research Scholar (PhD.)
Centre For Oceans, Rivers, Atmosphere & Land Science Technology
Indian Institute Of Technology, Kharagpur
contact :- 91- 8345984434