Hi, all of you,
From the Ferret Email Achives I found the
following scripts used for calculating thermocline depth, but I have some
questions about it.
Please see below, thanks very
much.
! calculation of thermocline depth
(by computing the depth of maximum vertical gradient)
use levitus_climatology let dtdz = temp[z=@ddc] let dtdz_min = dtdz[z=@min] ! Why is "@min" here, not "@max" ???????? let zero_at_min = dtdz - dtdz_min let thermocline_rough = zero_at_min[z=@loc:0] let thermocline = thermocline_rough[x=@sbx:11,y=@sbx:11] set var/title="Thermocline Depth"/units="meters" thermocline fill/lev=30 thermocline go fland ---------------
Gaorz
|