%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  +Trace       1 
%%%% FILES: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%% CONTROL parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Restart 0    0:new run;1:continue;2:run "+Run_time" more;3:restart w/ new TIOS
%%%% TIME Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Time_step   1.  hour  
 +Run_time    1 month - (1 year by default)	 

%%%% Barotropic Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Baro_solv  2    - (0: no U_BAR; 1:U = U_BAR = 0; 2:U = U_BAR(tau))
 +Baro_step  1
 +Baro_rayl  .0001
 +Baro_err   0
 +Baro_nnonz 6

%%%% MODEL Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Salinity   1   - (1-tracer; 2-linear; 3-dens12; 4-dens17;5-UNESCO) 
 +Mixing     1   - (1-Conv.Adj; 2-Ri#; 3=1&2; 4-5 disabled)
 +Mix_step   4   - freq of mixing

%%%% Topography%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Bathymetry 4   - (0-flat bottom, 1-read from datafile, >1-function)
 +Bath_type  2
 +Bath_fac   .1 
 +Ridge_min  5
 +Ridge_max  6
 -Bath_min   700. 

%%%% Shapiro FILTERS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% NOTE THAT it only matters whether Shap_vel_type is zero or non-zero
 +Shap_vel_order 4 
 +Shap_vel_type  3
 +Shap_vel_step  8
 +Shap_vel_damp  1.

 +Shap_scl_order 4 
 +Shap_scl_type  3
 +Shap_scl_step  8
 +Shap_scl_damp  1

%%%% Vertical diffusion %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Bint 1.e-4  vert diffusion coef. for momentum
 +Cint 1.e-5  vert diffusion coef. for heat & salt

%%%% Bottom friction %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Bottom_friction .000

%%%% Isopycnal diffusion %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Diffiso_vel       0     (1-use explicit diffusion on momentum;0-don't)
 +Diffiso_coef_mo 100.

 +Diffiso_scl       0     (1-use explicit diffusion on tracers;0-don't)
 +Diffiso_coef    100.    (coef of isopycnal diffusion)

 +Diffiso_alpha  0.      (0.-pure horizontal diffusion; 1.-pure isopycnal)
 +Diffiso_eps    0.00    (ratio of diapycnal/isopycnal diffusion)
 +Diffiso_cadv   100.    (coef of diffusion for Gent-McWilliams scheme)

%%%% Vertical advection of tracer %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Upwind_cnst    1000.

%%%% FORCING - INITIALIZATION %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Temp_init    0  (0:T_profile;3-from data)
 -Temp_coef    0.01  (percent random noise added to temp)
 +Salt_init    0  ( -"- )
 +Wind_forc    5  (1-L.I.;*2-C.S.;3-AVE_cnst(TAU);5-COS(y))
 +Wind_tauxy [1. 0.]  (in dynes/cm^2)
 +Wind_cos     1
 +Heat_forc    6  (0-3:Haney;0-30deg,1-T(k=1);2-annual SST;3-climatological SST
                    4- Q explicit data; 6- none; 8-AML)
 +EP_forc      6  ( as above)

 +NZ 4
 +H_profile [1000 1000 1000 1000]
 +T_profile [30 25 20 15 ]
 +S_profile [35 35 35 35 ]

 +Periodic	0
 +BC_type	2     - ( 1 - no slip; 2 - no normal flow/default/) 
 +BC_coef       10000 days
 +BC_bot_temp	1     - ( 0 - Const; 1 - No Flux)
 +Grid_type	1     - ( 1 - beta plane; 2 - spherical; 3 - f plane)

%% GEOMETRY Parameters %%%%%%%%%%%%
 +NX		10
 +NY	        10
 +South		20
 +North		60
 +West		 0
 +East		40

 +Grid_mask [
w10,r10,
]
%%%%% End Of GEOMETRY %%%%%%%%%%%%%%%%%%%%%%
End
