%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  +Trace       1 

%%%% TRACER PARAMETERS  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Tracers      1        - (# of tracers, if != 0 - additional data needed)
 +Tr_name      ['IAGE'] - name here must be consistent with
.ti
os file)
 +Trac_forc   [52]
 +Trac_init    [2]
 +Trac_init_value [1]

%%%% FILES: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
 +Wind_file     '/c4/naomi/DATA/tau_HEL'        

 +SST_file      '/c4/naomi/DATA/sst_LEV'
 +SSS_file      '/c4/naomi/DATA/sss_LEV'
 +Temp_file     '/c4/naomi/DATA/potemp_LEVs'
 +Salt_file     '/c4/naomi/DATA/salt_LEVs'

 +Cloud_file    '/c4/naomi/DATA/cld_ISCCP'
 +Solar_file    '/c4/naomi/DATA/solr_ISCCP'

 +Wndspd_file   '/c4/naomi/DATA/wnsp_ECMWF'
 +Uwnd_file     '/c4/naomi/DATA/uwnd_ECMWF'
 +Vwnd_file     '/c4/naomi/DATA/vwnd_ECMWF'


 +Precip_file   '/c4/naomi/DATA/prcp_Oberh'

 +Airhum_file   '/c4/naomi/DATA/ahum_ECMWFm'
 +Airtem_file   '/c4/naomi/DATA/atem_ECMWFm'

 +Bath_file     '/c4/naomi/DATA/bath_max'
 +River_file     '/c4/naomi/DATA/river.data'


 -Save_file    '/home/naomi/atl/output/atl.save'
 +Restart_file 'output/relax_spin.80y.restart'

%%%% CONTROL parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Restart 3    0:new run;1:continue;2:run "+Run_time" more;3:restart w/ new TIOS
%%%% TIME Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Time_step   1. hour  
 +Run_time    20 years (1 year by default)	 
%%%% MODEL Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Baro_solv  1    - (0: no U_BAR; 1:U = U_BAR = 0; 2:U = U_BAR(tau))
 +Baro_err   0

%%%% Island info for barotropic solve%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +New_island  [1 2 344 65 ]
 +Island_poly1  [340 63 340 67 348 67 348 63]
 +Sink_region [286 290 12 20 2000]

    -Iceland    1
    -Other_island  [340 348 63 67 2]


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

 +Bathymetry 1   - (0-flat bottom, 1-read from datafile, 2-function)
 +Bath_type  1   - (0-strict z, 1-zin(k)<z<zin(k+1), 2-zin(k)<z, 3-strict
                    z except last layer)
 +Bath_min   100. 

%%%% 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  0  (0-don't use;3-usual choice)
 +Shap_scl_step  8
 +Shap_scl_damp  1.
 
%%%% Isopycnal diffusion %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Diffiso_vel       0        (1-use explicit diffusion on momentum;0-don't)
 +Diffiso_coef_mo 1000.0

 +Diffiso_scl       1        (1-use explicit diffusion on tracers;0-don't)
 +Diffiso_coef   1000

 +Diffiso_alpha  1.      (0.-pure horizontal diffusion; 1.-pure isopycnal)
 +Diffiso_eps    0.00     (ratio of diapycnal/isopycnal diffusion)

%%%% partial upwind advection of scalars %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Upwind_cnst_z       1000.   (this is only for momentum advection)
 +Upwind_cnst_xy      1000.   (this is only for momentum advection)
 +TVD_vert   
 +TVD_horz   
 
%%%% Vertical diffusion %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Bint 1.e-4  vert diffusion coef. for momentum
 +Cint 1.e-5  vert diffusion coef. for heat & salt

 -Vert_kbot 1e-3   vert diffusion coef in bottom cells, (implicit)

%%%% Bottom drag %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Bottom_draglaw   2
 +Bottom_friction  0.002

%%%% FORCING - INITIALIZATION %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Temp_init    3  (0:T_profile;3-from data)
 +Salt_init    3  ( -"- )
 +Wind_forc    1  (1-L.I.;*2-C.S.;3-AVE_cnst(TAU);5-COS(y))
 -Wind_tauxy [-0. 0]
 +Heat_forc    8  (0-3:Haney;0-30deg,1-T(k=1);2-annual SST;3-climatological SST;
 +EP_forc      3  (0-3 - like Heat_forc but with SSS, 4 - EP explicit data)
 -EP_coef_precip  1.1
 -EP_budget    1
 -EP_bud_time   5 year
 -Rivers       1
 -River_coef   1.5
 -River_dmin   10.
 +Ice          1
 +Ice_dynamics 1
     +ICE_dyice_diff 500
 +PBL_step     4 hours
 +PBL_eddy     180000

 +Rlx_time_sst  0 year
 +Rlx_time_sss  1 month

 +Solar_penrad 1 (0:don't use; 1:use for interior)

%%%% SPONGE INFO %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Clim_H        0   - (0-initialize with hin; 1-Const. in time; 2-Vary;)  
 +Clim_TS       2   - (0-initialize with hin; 1-Const. in time; 2-Vary;)  
 -Clim_nz       10
 +Sponge_time    1 month
 +Clim_relax     3  - (0-no relaxation; 1-South-North; 2-Total; 3-Generalized) 
 +Clim_psi      0    - (0-no sponge for psi; 1-Const. in time; 2-Vary;)
 +Sponge_width  1   - (number of grid points)
 +Relax_width   1   - (number of grid points)

 +NZ 30
 +Z_profile [12 36 68 109 153.6 221.4 300 389.4 489.6 600.6 722.4 855.0 
998.4 1250 1500 1750 2000 2250 2500 2750 3000 3250 3500 3750 4000 4250
4500 4750 5000 5250]
 +Z_bot 5375

 +T_profile [24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1
1 1 1 1 1 1]
 +Temp_bot 1.
 +S_profile [35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35] 
 +Salt_bot 35

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

%% GEOMETRY Parameters %%%%%%%%%%%%
 +NX    60
 +NY    56
 +South    -35.0
 +North     75.0
 +West     262.0
 +East     380.0

 +Grid_mask [
z16,s5,z18,s19,z2,
z16,w6,z17,w18,w1,z2,
z16,w6,z17,w18,w1,z2,
z16,w7,z14,w20,w1,z2,
z17,w6,z10,w24,z3,
z17,w7,z5,w11,z3,w12,z5,
z17,w8,z4,w25,z6,
z17,w9,z2,w24,z8,
z18,w34,z8,
z19,w27,z2,w4,z8,
z20,w25,z15,
z21,w24,z15,
z22,w25,z13,
z23,w24,z13,
z19,w28,z13,
z19,w28,z13,
z14,w31,z15,
z13,w32,z15,
z12,w33,z15,
z11,w34,z15,
z11,w34,z15,
z10,w35,z15,
z9,w36,z15,
z2,w6,z1,w35,z16,
z2,w6,z1,w34,z17,
z2,w40,z18,
z2,w39,z19,
z6,w35,z19,
z6,w35,z19,
z6,w35,z19,
z8,w33,z19,
z8,w33,z19,
z8,w4,z3,w27,z18,
z19,w24,z17,
z20,w24,z16,
z23,w22,z15,
z24,w30,z6,
z24,w30,z6,
z26,w28,z6,
z30,w25,z5,
z31,w25,z4,
z32,w24,z4,
z32,w24,z4,
z31,w25,z4,
z30,w26,z4,
z30,w26,z4,
z30,w26,z4,
z30,w26,z4,
z29,w27,z4,
z27,w30,z3,
z26,w31,z3,
z25,w32,z3,
z25,w33,z2,
z25,w33,z2,
z25,w33,z2,
z25,s33,z2,
]
%%%%
 Of GEOMETRY %%%%%%%%%%%%%%%%%%%%%%
End