! NOAA/PMEL TMAP
 ! FERRET v6.72  
 ! Linux 2.6.32-131.12.1.el6.x86_64 64-bit - 09/13/11
 ! 27-Apr-16 00:35     

reg/y=-30:30
use "/home/fabric/LMDZ/ETOA/un.nc"
use "/prodigfs/ipslfs/dods/fabric/lmdz/VLR/YEAR/OBS/NC/wap.nc"
let champ=wap[k=22]*86400/100/100*ts[d=1]

let masksubs=if ( champ gt 20 ) then 1.
save/clobber/file=masksubs.nc masksubs

let maskweak=if ( champ ge -20 and champ le 20  ) then 1.
save/clobber/file=maskweak.nc maskweak

let maskconv=if ( champ lt -20 ) then 1.
save/clobber/file=maskconv.nc maskconv

let maskotrop=ts[d=1]
save/clobber/file=maskotrop.nc maskotrop

quit