\cancel mode verify
use gasx_ocmip2.nc
!
! >> Plot the sea/land mask (1 if by sea; 2 if by land)
say "Sea/land mask (press any key to visualize)"
message/quiet
shade tmask
!
! >> Plot the term Kw (excluding Fice and Sc terms):
say "Monthly piston velocity (press any key to visualize)"
message/quiet
let/quiet b=if tmask eq 1 then (1-fice)*xkw
set variable/title="Xconv * a * (u2 + v)"/units="m/s" b 
repeat/l=1:12 (shade b;go land;say "Kw: press any key for next month";message/quiet)
!
! >> Plot the Atmospheric Pressure:
say "Monthly atmospheric pressure [atm] (press any key to visualize)"
message/quiet
let/quiet c=if tmask eq 1 then P
set variable/title="Dry Atmospheric Pressure (sea level)"/units="atm" c 
repeat/l=1:12 (shade c;go land;say "P [atm]: press any key for next month";message/quiet)
!
! >> Plot the Fractional seaice
say "Monthly fractional sea ice (press any key to visualize)"
message/quiet
let/quiet d=if tmask eq 1 then fice
set variable/title="Fractional sea ice (fice)" d 
repeat/l=1:12 (shade d;go land;say "Ice: press any key for next month";message/quiet)