[Thread Prev][Thread Next][Index]

[ferret_users] Masking of temperature to the MLD



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:
    1> ctrl_23_temp_regrid.nc  (default)
 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)

 2> ./dataatul.nc  (default)
 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.









 

[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement