Mixed layer average velocity

ref

This script is useful for calculating velocity vertically averaged within MLD.

 

let v_weighted = IF zboxlo[gz=v] LT h then v*min(zbox[gz=v],h-zboxlo[gz=v])

let v_zave = v_weighted[z=@sum]/h

*v: velocity, h:mixed layer depth (MLD)

v can be u, w, or temperature.

 

So, Simple 🙂

 

 

Advertisements