|
Hello, everyone
I want to compute Mixed Layer Depth as the
following , but I found there is no any values on the ocean.
Is there anyone know what's wrong with the
following script file?
---------------------------------------------------------------------------
!MLD=the depth of
[T(k=1)-0.5]
use
levitus_climatology
define axis/z=0:5000:2/unit=meters z1 define grid/like=temp/z=z1 gz1 let tempz1=temp[g=gz1] let mlda=tempz1[k=1]-0.5 let mldb=tempz1-mlda let mld=mldb[z=@loc:0] set var/title="Mixed Layer Depth"/units="meters" mld contour mld go fland frame/file=d:\ferret\MixedLayerDepth0.gif ---------------------------------------------------------------------------------- Gaorz
|