! NOAA/PMEL TMAP
 ! FERRET v6.67  
 ! Linux(gfortran) 2.6.18-238.1.1.el5 - 02/11/11
 !  9-Sep-16 10:10     

use "tavg.15500.01.01.nc"
shade f_salt[l=50,d=1]
shade f_evap[l=50,d=1]
shade f_evap[l=50,d=1]-f_precip[l=50,d=1]
shade f_precip[l=50,d=1]
shade f_evap[l=50,d=1]-f_precip[l=50,d=1]
shade f_salt[l=50,d=1]
sp grep socn code/*
sp grep "F_evap" code/*
sp grep "oevap" code/*
sp grep "F_evap" code/*
sp grep "F_precip" code/*
shade f_salt[l=50,d=1]
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1])*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^6
sp grep socn code/*
let socn = 0.03475
sp grep "F_evap" code/*
sp grep "oevap" code/*
sp grep socn code/*
shade f_salt[l=50,d=1]*10^9
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^9
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^9
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^9
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade f_salt[l=50,d=1]*10^9
shade (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^9
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade/lev=(-100,100,10) (f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
shade/lev=(-100,100,10) f_salt[l=50,d=1]*10^6/socn
shade f_salt[l=50,d=1]*10^6/socn-(f_evap[l=50,d=1]-f_precip[l=50,d=1]-f_rivdis[l=50,d=1])*10^6
list socn
exi