[ferret_users] computation of MLD using Jean-Baptiste Sallée 2021, Nature

[saurabh rathore <rohitsrb2020@xxxxxxxxx>]

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

<