Dear ferret users please help me how to mask the temperature upto MLD. I have two data files one is temperature file and other one is MLD file
currently SET data sets:
name title I J K L M N
DEPTHT1_1_BNDS
1:2 ... 1:75 ... ... ...
(invalid axis bounds)
TEMPE TEMP[X=40E:110E,Y=20S:30N] 1:71 1:51 1:75 1:73 ... ...
(T=01-JAN-2001 00:00:01-JAN-2002 00:00)
name title I J K L M N
LATITUDE Station latitude ... 1:51 ... ... ... ...
LONGITUDE
Station longitude 1:71 ... ... ... ... ...
MLD mld at 1 degree grid 1:71 1:51 ... 1:146 ... ...
I am using following command to mask my temperature file but I got this error
define grid/like=tempe[d=1] gcommon
let temp1=tempe[d=1]
let ZMIX=mld[d=2]
let mask = IF Z[g=temp1] LT zmix[g=gcommon] THEN 1
LET masked_temp = mask*temp1
when I am plotting the data
fill/k=1/l=1 masked_temp
it show no valid data are found
Please tell me how to remove this problem.