%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   This is an Example of an input File.               Main features:      %
%   - '+' before the Name makes the definition valid, otherwise a default  %
%     value (from your program) is used                                    %
%   - All *empty* lines and lines whithout a definition are skipped        %
%   - '%' and '#' in first column denote an explicitly commented line      %
%   - An order of definitions in the file is insignificant                 %
%   - Search goes from the BOF or Variant upto line started with 'End'     %
%   - To use this input format one should compile with -lsenq library      %
%                                                                          %
%                                                   Senya Basin, 1994.     %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  +Trace       1 
%%%% FILES: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 -Wind_file     'DATA/tau_HELm'
 -SST_file      'DATA/sst_LEVm'
 -SSS_file      'DATA/sss_LEVm'
 -Temp_file     'DATA/temp_LEVs'
 -Salt_file     'DATA/salt_LEVs'

 -Bath_file     'DATA/bath'

 -Cloud_file    'DATA/cld_ISCCPm'
 -Solar_file    'DATA/solr_ERBEm'

 -Wndspd_file   'DATA/wnsp_ECMWFm'
 -Uwnd_file     'DATA/uwnd_ECMWFm'
 -Vwnd_file     'DATA/vwnd_ECMWFm'

 -Airhum_file   'DATA/ahum_ECMWFm'
 -Airtem_file   'DATA/atem_ECMWFm'

 -EP_file       'DATA/prcp_DSilm'

%%%% 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    5 years - (1 year by default)	 

%%%% MODEL Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Salinity   4   - (1-tracer; 2-linear; 3-dens12; 4-dens17;5-UNESCO) 
 +Mixing     3   - (1-Conv.Adj; 2-Ri#; 3=1&2; 4-5 disabled)
 +Mix_step   4   - freq of mixing
 +Bathymetry 0   - (0-flat bottom, 1-read from datafile, >1-function)
 +Bath_type  0   - (0-strict z, 1-zin(k)<z<zin(k+1), 2-zin(k)<z, 3-strict
                    z except last layer)
 +Bath_min   100. 

%%%% Discretization Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 -Discretization 4

%%%% Barotropic Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Baro_solv  2    - (0: no U_BAR; 1:U = U_BAR = 0; 2:U = U_BAR(tau))
 -Baro_step  6    - (frequency of barotropic update in large time step unit)
 -Baro_rayl  .001 - (coefficient of bottom drag)
 +Baro_err   0    - (what diagnostic info to print)
 +Baro_nnonz 6    - (increase this value if ifail=5)
 -Baro_nonlin 1   - (1: turns off the nonlinear advection terms in baro mode)
                    (0: include everything/default)
 -Antarctica
 -Antarctica_lat -60 - (latitude line for line integral, default=average)
 -Channel    
 -Channel_lat 32     - (latitude line for line integral, default=average)
 -Islands
    -Iceland    1    - (0:ignore island; 1:standard island; 2:sink island)
    -Australia  1    - (0:ignore island; 1:standard island; 2:sink island)
    -Other   [-25 -12 17 24 2]      - (alon,blon,alat,blat,type(as above))

%%%% Shapiro FILTERS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Shap_vel_order 4 
 +Shap_vel_type  3
 +Shap_vel_step  8

 +Shap_scl_order 4
 +Shap_scl_type  3
 +Shap_scl_step  8

%%%% Isopycnal diffusion %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 -Diffiso_scl    1        (1-use explicit diffusion on tracers;0-don't)
 -Diffiso_vel    0        (1-use explicit diffusion on momentum;0-don't)
 -Diffiso_alpha  1.       (0.-pure horizontal diffusion; 1.-pure isopycnal)
 -Diffiso_coef   500.    (coef of isopycnal diffusion)
 -Diffiso_eps    0.00     (ratio of diapycnal/isopycnal diffusion)
 -Diffiso_cadv   500.    (coef of diffusion for Gent-McWilliams scheme)

%%%% Vertical advection of tracer %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 -Upwind_cnst    2.       (1. corresponds to original, must be >1)

%%%% FORCING - INITIALIZATION %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 +Wind_forc  5 (1-L.I.;*2-C.S.;3-AVE_cnst(TAU);5-COS(y))
 +Wind_tauxy [-0.5 0]  - taux, tauy (for COS forcing)
 +Temp_init  0 (0-tin;1-Lin.(SST:max_dep);3-Annual Levitus;4-Vary Levitus;
 +Heat_forc  0 (0-Haney;1-const(t)=temp(t=0,L=1);2-con(t)=ave_t(HTFLX)
%               3-LI_t from HTFLX;4-CU_t from HTFLX;6-no heat flux)
 +PBL_step     .5 day

 -EP_forc    6  (0-3 - like Heat_forc but with SSS, 4 - EP explicit data)
 -Ice          1   - (1-use ice;  0- no ice)
 -Ice_dynamics 1   - (1-move ice; 0- static ice)
 -Tracers      0   - (# of tracers)

 +Clim_TS       2   - (0-not used; 1-Const. in time; 2-Vary;)
 +Clim_coef     30. days
 +Clim_relax    0   - (0-no relaxation; 1-South-North; 2-Total; 3-general)
 -Pre_proc      1   - (1-convectively adjust sponge data before using)
 -Sponge_width  10

 +NZ		7
 +Nsigma	0
 +Z_profile [20 60 100 200 400 1000 3000]  - array of NZ size
 +Z_bot 5000
 +T_profile [29 28 27 16 10  5  0]  - array of NZ size
 +Temp_bot      0
 +S_profile [35 35 35 35 35 35 35]
 +Salt_bot      35

 +Bint		1.e-4
 +Cint		1.e-5   % 1e-4 equivalent to OLD CINT=0.5e-7 and 400m depth.

 +BC_type	1     - ( 1 - no slip; 2 - no normal flow/default/) 
 +BC_bot_temp	1     - ( 0 - const; 1 - no flux)
 +Grid_type	2     - ( 1 - beta plane; 2 - spherical; 3 - f plane)

%%%% GEOMETRY Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% map20
 +NX		81
 +NY		30
 +South		-29
 +North		29
 +West		124
 +East		284

 +Grid_mask [
w61,z20,
w61,z20,
w62,z19,
w63,z18,
w66,z15,
w66,z15,
w68,z13,
w72,z9,
w74,z7,
w75,z6,
w77,z4,
w79,z2,
w79,z2,
w79,z2,
z2,w77,z2,
z6,w72,z3,
z10,w68,z3,
z12,w66,z3,
z12,w66,z3,
z11,w68,z2,
z11,w69,z1,
z11,w69,z1,
z11,w70,
z12,w69,
z13,w68,
z14,w67,
z14,w67,
z15,w66,
z15,w66,
z15,w66,
]
%%%%% End Of GEOMETRY %%%%%%%%%%%%%%%%%%%%%%
End