Commit aa74f04318277c689f14cd86a8299daa22237ad0

Authored by dumoda01
1 parent ade7be80
Exists in master and in 1 other branch snow

Suppression du repertoire exp

exp/airsea.nml deleted
... ... @@ -1,98 +0,0 @@
1   -!$Id: airsea.proto,v 1.1.1.1 2003/03/11 13:38:58 kbk Exp $
2   -!-------------------------------------------------------------------------------
3   -!
4   -!-------------------------------------------------------------------------------
5   -! air-sea interaction (heat, momentum and freshwater fluxes and solar radiation)
6   -!
7   -! calc_fluxes -> surface fluxes calculated by means of bulk formulae
8   -! (.true./.false.). Solar radiation is calculated from
9   -! time, latitude, longitude and clouds. If (.true.),
10   -! meteo_file must be given and wet_mode must be specified.
11   -! If (.false.), surface fluxes and solar radiation are
12   -! prescribed.
13   -!
14   -! meteo_file -> file with meteo data (for calc_fluxes=.true.) with
15   -! - date (yyyy-mm-dd hh:mm:ss)
16   -! - x-comp. of wind (10 m) in m/s
17   -! - y-comp. of wind (10 m) in m/s
18   -! - air pressure ( 2 m) in hectopascal
19   -! - dry air temp. ( 2 m) in Celsius
20   -! - relative humidity in % or wet bulb temperature in C
21   -! or dew point temperature in C (depending on wet_mode)
22   -! - cloud cover in 1/10
23   -!
24   -! wet_mode -> decides what is given in 7. column in meteo_file
25   -! 1: relative humidity
26   -! 2: wet bulb temperature
27   -! 3: dew point temperature
28   -!
29   -! heat_method -> method to provide short wave radiation (swr) and
30   -! surface heat flux (qh)
31   -! (only for calc_fluxes=.false.)
32   -! 0: heat flux not prescribed
33   -! 1: constant "const_swr" and "const_qh" given (see below)
34   -! 2: swr and qh are read from heatflux_file
35   -
36   -! const_swr -> constant value of incoming short wave radiation in W/m^2
37   -! (always positive)
38   -!
39   -! const_qh -> constant value of surface heat flux in W/m^2
40   -! (negative for heat loss)
41   -!
42   -! heatflux_file -> file with qin and qout given in W/m^2
43   -! (negative for net outgoing)
44   -!
45   -! momentum_method -> method how momentum fluxes are given
46   -! (only for calc_fluxes=.false.)
47   -! 0: momentum flux not prescribed
48   -! 1: constant surface momentum fluxes given
49   -! 2: surface momentum fluxes given from file momentumflux_file
50   -!
51   -! const_tx -> x-component of surface momentum flux in N/m^2
52   -! const_ty -> y-component of surface momentum flux in N/m^2
53   -!
54   -! momentumflux_file-> file with tx and ty given in N/m^2
55   -!
56   -! p_e_method -> method how fresh water fluxes (P-E) are given
57   -! 0: P-E not used
58   -! 1: constant value for P-E (in m/s) used
59   -! (P-E = precipitation-evaporation)
60   -! 2: values for P-E read from file
61   -!
62   -! const_p_e -> constant value for P-E in m/s (positive for P>E)
63   -!
64   -! p_e_flux_file -> file with value for P-E (positive for P>E)
65   -! used if p_e_method=2
66   -!
67   -! sst_method -> method how sea surface temperature (SST) is given
68   -! 0: no independent SST observation is read from file
69   -! 2: independent SST observation is read from file,
70   -! only for output
71   -!
72   -! sst_file -> file with independent SST observation
73   -!
74   -! sss_method -> method how sea surface salinity (SSS) is given
75   -! 0: no independent SSS observation is read from file
76   -! 2: independent SSS observation is read from file,
77   -! only for output
78   -!-------------------------------------------------------------------------------
79   - &airsea
80   - calc_fluxes= .false.
81   - meteo_file= 'meteo.dat'
82   - wet_mode= 1
83   - heat_method= 2
84   - const_swr= 100.0
85   - const_heat= -100.0
86   - heatflux_file= 'narr_daily_heatflux_ice.dat'
87   - momentum_method= 2
88   - const_tx= 0.1
89   - const_ty= 0.0
90   - momentumflux_file='narr_hourly_momentumflux.dat'
91   - p_e_method= 0
92   - const_p_e= 0.
93   - p_e_flux_file= 'p_e.dat'
94   - sst_method= 0
95   - sst_file= 'sst.dat'
96   - sss_method= 0
97   - sss_file= 'sss.dat'
98   - /
exp/bio.nml deleted
... ... @@ -1,56 +0,0 @@
1   -!$Id$
2   -!-------------------------------------------------------------------------------
3   -! Basic settings for biogeochemical model
4   -!
5   -! bio_calc -> calculation of the bio model 'bio_model' (.true./.false.)
6   -! bio_model -> choice of the bio model:
7   -! 1: NPZD (4 variables)
8   -! 2: IOW-ERGOM (9 variables)
9   -! 3: Suspended matter only (1 variable)
10   -! 4: Fasham et al. 1990 (7 variables)
11   -!
12   -! bio_eulerian -> state variables are Eulerian (.true./.false.)
13   -!
14   -! cnpar -> Cranck-Nicolson parameter for vertical diffusion
15   -!
16   -! w_adv_discr -> advection scheme for vertical motion
17   -! 1: first order upstream
18   -! 2: not coded yet
19   -! 3: third-order polynomial
20   -! 4: TVD with Superbee limiter
21   -! 5: TVD with MUSCL limiter
22   -! 6: TVD with ULTIMATE QUICKEST
23   -!
24   -! ode_method -> ODE scheme for source and sink dynamics
25   -! 1: first-order explicit (not positive)
26   -! 2: second order explicit Runge-Kutta (not positive)
27   -! 3: fourth-order explicit Runge-Kutta (not positive)
28   -! 4: Patankar (first order, not conservative)
29   -! 5: Patankar-RK (second order, not conservative)
30   -! 6: Patankar-RK (does not work, not conservative)
31   -! 7: Modified Patankar (1. order, conservat., posit.)
32   -! 8: Modified Patankar-RK (2. order, conservat., posit.)
33   -! 9: Modified Patankar-RK (does not work, conservat., posit.)
34   -! 10: Ext. Modified Patankar (1. order, conservat., posit.)
35   -! 11: Ext. Modified Patankar-RK (2. order, conservat., posit.)
36   -!
37   -! split_factor -> number of biogeochemical time steps per physical time step
38   -!
39   -! bioshade_feedback -> feedback of bio-turbidity to temp. eq. (.true./.false.)
40   -!
41   -! bio_lagrange_mean -> averaging Lagrangian conc. on output (.true./.false.)
42   -!
43   -! bio_npar -> total number of Lagrangian particles
44   -!-------------------------------------------------------------------------------
45   -&bio_nml
46   - bio_calc= .true.
47   - bio_model= 4
48   - bio_eulerian= .true.
49   - cnpar= 1.0
50   - w_adv_discr= 6
51   - ode_method= 8
52   - split_factor= 1
53   - bioshade_feedback= .true.
54   - bio_lagrange_mean= .false.
55   - bio_npar= 1000
56   - /
exp/bio_fasham.nml deleted
... ... @@ -1,91 +0,0 @@
1   -#$Id$
2   -!-------------------------------------------------------------------------------
3   -! Fasham et al. biological model with modifications by Kuehn and Radach
4   -!
5   -! numc= number of compartments for geobiochemical model
6   -!
7   -! p_initial= initial phytoplankton concentration [mmol n/m3]
8   -! z_initial= initial zooplakton concentration [mmol n/m3]
9   -! b_initial= initial bacteria concentration [mmol n/m3]
10   -! d_initial= initial detritus concentration [mmol n/m3]
11   -! n_initial= initial nitrate concentration [mmol n/m3]
12   -! a_initial= initial ammonium concentration [mmol n/m3]
13   -! l_initial= initial LDON concentration [mmol n/m3]
14   -! p0 = minimum phytoplankton concentration [mmol n/m3]
15   -! z0 = minimum zooplakton concentration [mmol n/m3]
16   -! b0 = minimum bacteria concentration [mmol n/m3]
17   -! vp = maximum phytoplankton uptake rate [1/day]
18   -! alpha = slope of the PI-curvea [m2/(W day)]
19   -! k1 = half saturation constant nitrate uptake [mmol n/m3]
20   -! k2 = half saturation constant ammonium uptake [mmol n/m3]
21   -! mu1 = phytoplankton mortality rate [1/day]
22   -! k5 = half saturation constant phytoplankton mortality [mmol n/m3]
23   -! gamma = exudation fraction [-]
24   -! w_p = phytoplankton settling velocity [m/day]
25   -! gmax = maximum ingestion rate [1/day]
26   -! k3 = half saturation constant ingestion [mmol n/m3]
27   -! beta = grazing efficiency [-]
28   -! mu2 = maximum zooplankton loss rate [1/day]
29   -! k6 = half saturation zooplankton loss [mmol n/m3]
30   -! delta = fractional zooplankton loss to LDON [-]
31   -! epsi = fractional zooplankton loss to ammonium [-]
32   -! r1 = grazing preference phytoplankton [-]
33   -! r2 = grazing preference bacteria [-]
34   -! r3 = grazing preference detritus [-]
35   -! vb = maximum bacterial uptake rate [1/day]
36   -! k4 = half saturation bacterial uptake [mmol n/m3]
37   -! mu3 = bacteria excretion rate [1/day]
38   -! eta = uptake ratio ammonium:LDON [-]
39   -! mu4 = detritus breakdown rate [1/day]
40   -! w_d = detritus settling velocity [m/day]
41   -! kc = attenuation constant for the self shading effect [m**2/mmol N]
42   -! I_min = minimum photosynthetically active radiation (PAR) [W/m**2]
43   -! I_opt = optimal photosynthetically active radiation (PAR) [W/m**2] !CHG1
44   -! inib = inhibition slope of the PI-curve (positive) [m2/(W day)] !CHG1
45   -! theta = phytoplancton buoyancy parameter [m3 day/(mmol N)] !CHG2
46   -!-------------------------------------------------------------------------------
47   - &bio_fasham_nml
48   - numc= 7
49   - p_initial= 0.012
50   - z_initial= 0.012
51   - b_initial= 0.001
52   - d_initial= 0.01
53   - n_initial= 12.0
54   - a_initial= 0.1
55   - l_initial= 0.1
56   - p0= 0.0001
57   - z0= 0.0001
58   - b0= 0.0001
59   - vp= 0.3
60   - alpha= 0.04
61   - inib= 0.06
62   - I_opt= 20.0
63   - k1= 1.0
64   - k2= 0.8
65   - mu1= 0.05
66   - k5= 0.2
67   - gamma= 0.05
68   - w_p= -0.38
69   - theta= 0.0
70   - w_pmin= -0.06
71   - w_pmax= -0.38
72   - gmax= 1.0
73   - k3= 1.0
74   - beta= 0.625
75   - mu2= 0.3
76   - k6= 0.2
77   - delta= 0.1
78   - epsi= 0.70
79   - r1= 0.55
80   - r2= 0.4
81   - r3= 0.05
82   - vb= 0.24
83   - k4= 0.5
84   - mu3= 0.03
85   - eta= 0.0
86   - mu4= 0.02
87   - w_d= -5.0
88   - kc= 0.03
89   - aa= 0.70
90   - g2= 14.0
91   - /
exp/gotmmean.nml deleted
... ... @@ -1,50 +0,0 @@
1   -!$Id: gotmmean.proto,v 1.1.1.1 2003/03/11 13:38:58 kbk Exp $
2   -!-------------------------------------------------------------------------------
3   -! The namelists 'meanflow' is read in meanflow.F90.
4   -!-------------------------------------------------------------------------------
5   -
6   -!-------------------------------------------------------------------------------
7   -! Specify variables related to the 1D meanflow model.
8   -!
9   -! h0b= bottom roughness - Note: z0b=0.03*h0b+0.1*nu/ustar [m]
10   -! z0s_min= minimum value of z0s, default value if charnock=.false. [m]
11   -! charnock= .true.: adaptation of Charnok 1955 formula used
12   -! .false.: constant surface roughness length z0s_min used
13   -! charnock_val= emp. constant in Charnok 1955 formula (default = 1400.)
14   -! ddu= grid zooming (surface), 0: no zooming; > 3 strong zooming
15   -! ddl= grid zooming (bottom), 0: no zooming; > 3 strong zooming
16   -! grid_method= 0: zooming of grid with ddl, ddu >= 0
17   -! 1: sigma grid (relative depth fractions) read from file
18   -! 2: cartesian grid (fixed layer height in m) read from file
19   -!
20   -! grid_file= file for sigma or cartesian grid. the first line gives the
21   -! number of layers, the following lines give fractions or
22   -! layer heights in m from the surface down to the bottom.
23   -! gravity= gravitational acceleration [m/s^2]
24   -! rho_0= Reference density [kg/m^3].
25   -! cp= Specific heat of sea water [J/kg/K].
26   -! avmolu= molecular viscosity for momentum [m^2/s].
27   -! avmolt= molecular diffusity for temperature [m^2/s].
28   -! avmols= molecular diffusity for salinity [m^2/s].
29   -! MaxItz0b= max # of iterations for z0b as function of u_taub.
30   -! no_shear= .true.: shear production term P is set to zero
31   -! avmoln= molecular diffusivity for nitrate [m^2/s]. !DD
32   -!-------------------------------------------------------------------------------
33   - &meanflow
34   - h0b= 0.05
35   - z0s_min= 0.02
36   - charnock= .false.
37   - charnock_val= 1400.
38   - ddu= 1.
39   - ddl= 0.
40   - grid_method= 0
41   - grid_file= 'grid.dat'
42   - gravity= 9.81
43   - rho_0= 1027.
44   - cp= 3985.
45   - avmolu= 1.3e-6
46   - avmolt= 1.4e-7
47   - avmols= 1.1e-9
48   - MaxItz0b= 1
49   - no_shear= .false.
50   - /
exp/gotmrun.nml deleted
... ... @@ -1,138 +0,0 @@
1   -!$Id: gotmrun.proto,v 1.1.1.1 2003/03/11 13:38:58 kbk Exp $
2   -!-------------------------------------------------------------------------------
3   -!
4   -!-------------------------------------------------------------------------------
5   -! general model setup
6   -!
7   -! title -> title of simulation
8   -! nlev -> number of levels
9   -! dt -> time step in seconds
10   -! cnpar -> parameter for "explicitness" of numerical scheme
11   -! (between 0.0 and 1.0)
12   -! buoy_method -> method to compute mean buoyancy
13   -! 1: from equation of state
14   -! (i.e. from potential temperature and salinity)
15   -! 2: from prognostic equation
16   -!
17   -!-------------------------------------------------------------------------------
18   - &model_setup
19   - title= "Arctic SCM"
20   - nlev= 80
21   - dt= 300.
22   - cnpar= 1.0
23   - buoy_method= 2
24   - /
25   -
26   -!-------------------------------------------------------------------------------
27   -! geographic location
28   -!
29   -! name -> name of the station
30   -! latitude -> latitude in degree (north is positive)
31   -! longitude -> longitude in degree (east is positive)
32   -! depth -> water depth in meters
33   -!
34   -!-------------------------------------------------------------------------------
35   - &station
36   - name= "Amundsen Gulf"
37   - latitude= 71.5
38   - longitude= -127.0
39   - depth= 200.0
40   - /
41   -
42   -!-------------------------------------------------------------------------------
43   -! duration of run
44   -!
45   -! timefmt -> method to specify start and duration of model run
46   -! 1: duration computed from number of time steps, MaxN
47   -! (bogus start date used)
48   -! 2: duration computed from given start and stop dates
49   -! (number of time steps MaxN computed)
50   -! 3: duration computed from number of time steps, MaxN
51   -! (start date as specified, stop date computed)
52   -!
53   -! MaxN -> nominal number of time steps (see "timefmt")
54   -! start -> nominal start date: YYYY/MM/DD HH:MM:SS (see "timefmt")
55   -! stop -> nominal stop date: YYYY/MM/DD HH:MM:SS (see "timefmt")
56   -!
57   -!-------------------------------------------------------------------------------
58   - &time
59   - timefmt= 2
60   - MaxN= 1200
61   - start= '2004-01-01 00:00:00'
62   - stop= '2004-12-31 00:00:00'
63   - /
64   -
65   -!-------------------------------------------------------------------------------
66   -! format for output and filename(s).
67   -!
68   -! out_fmt -> format for GOTM output
69   -! 1: ASCII
70   -! 2: NetCDF
71   -! 3: GrADS
72   -!
73   -! out_dir -> path to output directory (set permissions)
74   -! out_fn -> output string used to generate output file names
75   -! nsave -> save results every 'nsave' timesteps
76   -! diagnostics -> diagnostics are written to output (if .true.)
77   -!
78   -! mld_method -> how to diagnose mixed layer depth
79   -! 1: mixed layer depth computed from TKE threshold
80   -! 2: mixed layer depth from Ri threshold
81   -! diff_k -> TKE threshold [m^2/s^2] for mixed layer depth
82   -! ri_crit -> Ri threshold for mixed layer depth
83   -!
84   -! rad_corr -> correct surface buoyancy flux for solar radiation
85   -! for output (if true)
86   -!
87   -!-------------------------------------------------------------------------------
88   - &output
89   - out_fmt= 2
90   - out_dir= "."
91   - out_fn= "amdgulf"
92   - nsave= 36
93   - diagnostics= .false.
94   - mld_method= 2
95   - diff_k= 1.e-5
96   - Ri_crit= 0.5
97   - rad_corr= .true.
98   - /
99   -
100   -!-------------------------------------------------------------------------------
101   -! Specify variables related to the equation of state.
102   -!
103   -! eq_state_mode -> choice for empirical formula for equation of state
104   -! 1: UNESCO equation of state by Fofonoff and Millard (1983)
105   -! 2: equation of state according Jackett et al. (2005)
106   -!
107   -! eq_state_method -> method to compute density and buoyancy from salinity,
108   -! potential temperature and pressure
109   -! 1: full equation of state (i.e. with the LOCAL
110   -! pressure). This implies that T is NOT treated as
111   -! the potential temperature but rather as the in-situ
112   -! temperature!
113   -! 2: equation of state with pressure evaluated at the surface.
114   -! This implies that T is treated as the potential
115   -! temperature and thus rho as the potential density.
116   -! 3: linearized equation of state at T0,S0,p0
117   -! (again, use p0=p_surf to work with potential
118   -! temperature and density.)
119   -! 4: linear equation of state with T0,S0,dtr0,dsr0
120   -!
121   -! For the precise definition of the following quantities, see
122   -! GOTM documentation:
123   -!
124   -! T0 -> reference temperature (deg C) for linear equation of state
125   -! S0 -> reference salinity (psu) for linear equation of state
126   -! p0 -> reference pressure (bar) for linear equation of state
127   -! dtr0 -> thermal expansion coefficient for linear equation of state
128   -! dsr0 -> saline expansion coefficient for linear equation of state
129   -!-------------------------------------------------------------------------------
130   - &eqstate
131   - eq_state_mode = 2
132   - eq_state_method= 2
133   - T0= 10.
134   - S0= 35.
135   - p0= 0.
136   - dtr0= -0.17
137   - dsr0= 0.78
138   - /
exp/gotmturb.nml deleted
... ... @@ -1,304 +0,0 @@
1   -!$Id: gotmturb.proto,v 1.1.1.1 2003/03/11 13:38:58 kbk Exp $
2   -!-------------------------------------------------------------------------------
3   -
4   -!-------------------------------------------------------------------------------
5   -! What type of equations are solved in the turbulence model?
6   -!
7   -! turb_method -> type of turbulence closure
8   -!
9   -! 0: convective adjustment
10   -! 1: analytical eddy visc. and diff. profiles, not coded yet
11   -! 2: turbulence Model calculating TKE and length scale
12   -! (specify stability function below)
13   -! 3: second-order model (see "scnd" namelist below)
14   -! 99: KPP model (requires "kpp.inp" with specifications)
15   -!
16   -!
17   -! tke_method -> type of equation for TKE
18   -!
19   -! 1: algebraic equation
20   -! 2: dynamic equation (k-epsilon style)
21   -! 3: dynamic equation (Mellor-Yamada style)
22   -!
23   -!
24   -! len_scale_method -> type of model for dissipative length scale
25   -!
26   -! 1: parabolic shape
27   -! 2: triangle shape
28   -! 3: Xing and Davies [1995]
29   -! 4: Robert and Ouellet [1987]
30   -! 5: Blackadar (two boundaries) [1962]
31   -! 6: Bougeault and Andre [1986]
32   -! 7: Eifler and Schrimpf (ISPRAMIX) [1992]
33   -! 8: dynamic dissipation rate equation
34   -! 9: dynamic Mellor-Yamada q^2l-equation
35   -! 10: generic length scale (GLS)
36   -!
37   -!
38   -! stab_method -> type of stability function
39   -!
40   -! 1: constant stability functions
41   -! 2: Munk and Anderson [1954]
42   -! 3: Schumann and Gerz [1995]
43   -! 4: Eifler and Schrimpf [1992]
44   -!
45   -!-------------------------------------------------------------------------------
46   - &turbulence
47   - turb_method= 99
48   - tke_method= 2
49   - len_scale_method=8
50   - stab_method= 3
51   - /
52   -
53   -!-------------------------------------------------------------------------------
54   -! What boundary conditions are used?
55   -!
56   -! k_ubc, k_lbc -> upper and lower boundary conditions
57   -! for k-equation
58   -! 0: prescribed BC
59   -! 1: flux BC
60   -!
61   -! psi_ubc, psi_lbc -> upper and lower boundary conditions
62   -! for the length-scale equation (e.g.
63   -! epsilon, kl, omega, GLS)
64   -! 0: prescribed BC
65   -! 1: flux BC
66   -!
67   -!
68   -! ubc_type -> type of upper boundary layer
69   -! 0: viscous sublayer (not yet impl.)
70   -! 1: logarithmic law of the wall
71   -! 2: tke-injection (breaking waves)
72   -!
73   -! lbc_type -> type of lower boundary layer
74   -! 0: viscous sublayer (not yet impl.)
75   -! 1: logarithmic law of the wall
76   -!
77   -!-------------------------------------------------------------------------------
78   - &bc
79   - k_ubc= 1
80   - k_lbc= 1
81   - psi_ubc= 1
82   - psi_lbc= 1
83   - ubc_type= 1
84   - lbc_type= 1
85   - /
86   -
87   -!-------------------------------------------------------------------------------
88   -! What turbulence parameters have been described?
89   -!
90   -! cm0_fix -> value of cm0 for turb_method=2
91   -! Prandtl0_fix -> value of the turbulent Prandtl-number for stab_method=1-4
92   -! cw -> constant of the wave-breaking model
93   -! (Craig & Banner (1994) use cw=100)
94   -! compute_kappa -> compute von Karman constant from model parameters
95   -! kappa -> the desired von Karman constant (if compute_kappa=.true.)
96   -! compute_c3 -> compute c3 (E3 for Mellor-Yamada) for given Ri_st
97   -! Ri_st -> the desired steady-state Richardson number (if compute_c3=.true.)
98   -! length_lim -> apply length scale limitation (see Galperin et al. 1988)
99   -! galp -> coef. for length scale limitation
100   -! const_num -> minimum eddy diffusivity (only with turb_method=0)
101   -! const_nuh -> minimum heat diffusivity (only with turb_method=0)
102   -! k_min -> minimun TKE
103   -! eps_min -> minimum dissipation rate
104   -! kb_min -> minimun buoyancy variance
105   -! epsb_min -> minimum buoyancy variance destruction rate
106   -!
107   -!-------------------------------------------------------------------------------
108   - &turb_param
109   - cm0_fix= 0.5477
110   - Prandtl0_fix= 0.74
111   - cw= 100.
112   - compute_kappa= .true.
113   - kappa= 0.4
114   - compute_c3= .true.
115   - ri_st= 0.25
116   - length_lim= .false.
117   - galp= 0.53
118   - const_num= 5.e-4
119   - const_nuh= 5.e-4
120   - k_min= 1.e-10
121   - eps_min= 1.e-12
122   - kb_min= 1.e-10
123   - epsb_min= 1.e-14
124   - /
125   -
126   -!-------------------------------------------------------------------------------
127   -! The generic model (Umlauf & Burchard, J. Mar. Res., 2003)
128   -!
129   -! This part is active only, when len_scale_method=10 has been set.
130   -!
131   -! compute_param -> compute the model parameters:
132   -! if this is .false., you have to set all
133   -! model parameters (m,n,cpsi1,...) explicitly
134   -! if this is .true., all model parameters
135   -! set by you (except m) will be ignored and
136   -! re-computed from kappa, d, alpha, etc.
137   -! (see Umlauf&Burchard 2002)
138   -!
139   -! m: -> exponent for k
140   -! n: -> exponent for l
141   -! p: -> exponent for cm0
142   -!
143   -! Examples:
144   -!
145   -! k-epsilon (Rodi 1987) : m=3/2, n=-1, p=3
146   -! k-omega (Umlauf et al. 2003) : m=1/2, n=-1, p=-1
147   -!
148   -! cpsi1 -> emp. coef. in psi equation
149   -! cpsi2 -> emp. coef. in psi equation
150   -! cpsi3minus -> cpsi3 for stable stratification
151   -! cpsi3plus -> cpsi3 for unstable stratification
152   -! sig_kpsi -> Schmidt number for TKE diffusivity
153   -! sig_psi -> Schmidt number for psi diffusivity
154   -!
155   -!-------------------------------------------------------------------------------
156   - &generic
157   - compute_param= .false.
158   - gen_m= 1.0
159   - gen_n= -0.67
160   - gen_p= 3.0
161   - cpsi1= 1.
162   - cpsi2= 1.22
163   - cpsi3minus= 0.05
164   - cpsi3plus = 1.0
165   - sig_kpsi= 0.8
166   - sig_psi= 1.07
167   - gen_d= -1.2
168   - gen_alpha= -2.0
169   - gen_l= 0.2
170   - /
171   -
172   -!-------------------------------------------------------------------------------
173   -! The k-epsilon model (Rodi 1987)
174   -!
175   -! This part is active only, when len_scale_method=8 has been set.
176   -!
177   -! ce1 -> emp. coef. in diss. eq.
178   -! ce2 -> emp. coef. in diss. eq.
179   -! ce3minus -> ce3 for stable stratification, overwritten if compute_c3=.true.
180   -! ce3plus -> ce3 for unstable stratification (Rodi 1987: ce3plus=1.0)
181   -! sig_k -> Schmidt number for TKE diffusivity
182   -! sig_e -> Schmidt number for diss. diffusivity
183   -! sig_peps -> if .true. -> the wave breaking parameterisation suggested
184   -! by Burchard (JPO 31, 2001, 3133-3145) will be used.
185   -!-------------------------------------------------------------------------------
186   - &keps
187   - ce1= 1.44
188   - ce2= 1.92
189   - ce3minus= -0.4
190   - ce3plus= 1.0
191   - sig_k= 1.
192   - sig_e= 1.3
193   - sig_peps= .false.
194   - /
195   -
196   -!-------------------------------------------------------------------------------
197   -! The Mellor-Yamada model (Mellor & Yamada 1982)
198   -!
199   -! This part is active only, when len_scale_method=9 has been set!
200   -!
201   -! e1 -> coef. in MY q**2 l equation
202   -! e2 -> coef. in MY q**2 l equation
203   -! e3 -> coef. in MY q**2 l equation, overwritten if compute_c3=.true.
204   -! sq -> turbulent diffusivities of q**2 (= 2k)
205   -! sl -> turbulent diffusivities of q**2 l
206   -! my_length -> prescribed barotropic lengthscale in q**2 l equation of MY
207   -! 1: parabolic
208   -! 2: triangular
209   -! 3: lin. from surface
210   -! new_constr -> stabilisation of Mellor-Yamada stability functions
211   -! according to Burchard & Deleersnijder (2001)
212   -! (if .true.)
213   -!
214   -!-------------------------------------------------------------------------------
215   - &my
216   - e1= 1.8
217   - e2= 1.33
218   - e3= 1.8
219   - sq= 0.2
220   - sl= 0.2
221   - my_length= 3
222   - new_constr= .false.
223   - /
224   -
225   -!-------------------------------------------------------------------------------
226   -! The second-order model
227   -!
228   -! scnd_method -> type of second-order model
229   -! 1: EASM with quasi-equilibrium
230   -! 2: EASM with weak equilibrium, buoy.-variance algebraic
231   -! 3: EASM with weak equilibrium, buoy.-variance from PDE
232   -!
233   -! kb_method -> type of equation for buoyancy variance
234   -!
235   -! 1: algebraic equation for buoyancy variance
236   -! 2: PDE for buoyancy variance
237   -!
238   -!
239   -! epsb_method -> type of equation for variance destruction
240   -!
241   -! 1: algebraic equation for variance destruction
242   -! 2: PDE for variance destruction
243   -!
244   -!
245   -! scnd_coeff -> coefficients of second-order model
246   -!
247   -! 0: read the coefficients from this file
248   -! 1: coefficients of Gibson and Launder (1978)
249   -! 2: coefficients of Mellor and Yamada (1982)
250   -! 3: coefficients of Kantha and Clayson (1994)
251   -! 4: coefficients of Luyten et al. (1996)
252   -! 5: coefficients of Canuto et al. (2001) (version A)
253   -! 6: coefficients of Canuto et al. (2001) (version B)
254   -! 7: coefficients of Cheng et al. (2002)
255   -!
256   -!-------------------------------------------------------------------------------
257   - &scnd
258   - scnd_method= 1
259   - kb_method= 1
260   - epsb_method= 1
261   - scnd_coeff= 7
262   - cc1= 3.6
263   - cc2= 0.8
264   - cc3= 1.2
265   - cc4= 1.2
266   - cc5= 0.0
267   - cc6= 0.3
268   - ct1= 3.28
269   - ct2= 0.4
270   - ct3= 0.4
271   - ct4= 0.0
272   - ct5= 0.4
273   - ctt= 0.8
274   - /
275   -
276   -!-------------------------------------------------------------------------------
277   -! The internal wave model
278   -!
279   -! iw_model -> method to compute internal wave mixing
280   -! 0: no internal waves mixing parameterisation
281   -! 1: Mellor 1989 internal wave mixing
282   -! 2: Large et al. 1994 internal wave mixing
283   -!
284   -! alpha -> coeff. for Mellor IWmodel (0: no IW, 0.7 Mellor 1989)
285   -!
286   -! The following six empirical parameters are used for the
287   -! Large et al. 1994 shear instability and internal wave breaking
288   -! parameterisations (iw_model = 2, all viscosities are in m**2/s):
289   -!
290   -! klimiw -> critcal value of TKE
291   -! rich_cr -> critical Richardson number for shear instability
292   -! numshear -> background diffusivity for shear instability
293   -! numiw -> background viscosity for internal wave breaking
294   -! nuhiw -> background diffusivity for internal wave breaking
295   -!-------------------------------------------------------------------------------
296   - &iw
297   - iw_model= 0
298   - alpha= 0.7
299   - klimiw= 1e-6
300   - rich_cr= 0.7
301   - numshear= 5.e-3
302   - numiw= 1.e-4
303   - nuhiw= 1.e-5
304   - /
exp/kpp.nml deleted
... ... @@ -1,40 +0,0 @@
1   -!$Id$
2   -!-------------------------------------------------------------------------------
3   -!
4   -!-------------------------------------------------------------------------------
5   -! specifications for the KPP turbulence model
6   -!
7   -! Set "turb_method=99" in gotmturb.inp and check for the correct pre-processor
8   -! macros in cppdefs.h.
9   -! These are (see documentation at www.gotm.net):
10   -!
11   -! NONLOCAL for inclusion of non-local fluxes
12   -! KPP_SHEAR for shear instability interior mixing
13   -! KPP_INTERNAL_WAVE for internal waves interior mixing
14   -! KPP_CONVEC for convective interior mixing
15   -! KPP_DDMIX for double-diffusion interior mixing
16   -! KPP_TWOPOINT_REF for two grid points to compute reference values
17   -! KPP_IP_FC for scheme to interpolate MLD
18   -! KPP_CLIP_GS for clipping of shape function G(sigma)
19   -! KPP_SALINITY for computation of salinity diffusivity
20   -!
21   -! These pre-processor macros have been introduced for higher efficiency
22   -! of the code.
23   -!
24   -! The main flags for the KPP algorithm can be set in this file.
25   -! They are:
26   -!
27   -! kpp_sbl -> .true. for active surface boundary layer module
28   -! kpp_bbl -> .true. for active bottom boundary layer module
29   -! kpp_internal -> .true. for active interior mixing
30   -! clip_mld -> .true. for clipping of MLD at MO or Ekman scale
31   -! Ric -> critical value of bulk Richardson number
32   -!
33   -!-------------------------------------------------------------------------------
34   - &kpp
35   - kpp_sbl= .true.
36   - kpp_bbl= .true.
37   - kpp_interior= .true.
38   - clip_mld= .false.
39   - Ric= 0.3
40   - /
exp/obs.nml deleted
... ... @@ -1,564 +0,0 @@
1   -!$Id: obs.proto,v 1.1.1.1 2003/03/11 13:38:58 kbk Exp $
2   -!-------------------------------------------------------------------------------
3   -!
4   -!-------------------------------------------------------------------------------
5   -! observed or prescribed salinity profiles
6   -!
7   -! s_prof_method -> method to create initial or observed salinity profiles
8   -! 0: no initial values, S-equation is not solved
9   -! 1: use analytically prescribed initial profile
10   -! 2: read profiles at different dates from "s_prof_file"
11   -! and interpolate to GOTM timestep
12   -!
13   -! s_analyt_method -> method to create analytically precribed inital profile
14   -! 1: set profile to constant value s_1
15   -! 2: set "two layer" stratification (see user's guide)
16   -! 3: set profile with constant N^2 (see user's guide)
17   -! This option can only be used toghether with
18   -! t_analyt_method=1 (constant temperature).
19   -!
20   -! z_s1 -> upper layer thickness if s_analyt_method=2
21   -!
22   -! s_1 -> constant salinity if s_analyt_method=1
23   -! upper layer salinity if s_analyt_method=2
24   -! surface salinity if s_analyt_method=3
25   -!
26   -! z_s2 -> depth of top of lower layer if s_analyt_method=2
27   -!
28   -! s_2 -> lower layer salinity if s_analyt_method=2
29   -!
30   -! s_obs_NN -> constant value N^2 corresponding to salinity profile
31   -! if s_analyt_method=3
32   -
33   -! s_prof_file -> filename of file with salinity profiles
34   -! if s_prof_method=2
35   -!
36   -! The computed profiles can be relaxed towards observed or prescribed
37   -! profiles with a certain relaxation time. If you do not want relaxation,
38   -! set the relaxation times to 1.e15 (something large). It is possible to choose
39   -! different relaxation times in a surface and bottom layer.
40   -!
41   -! SRelaxTauM -> relaxation time for bulk of the flow
42   -! SRelaxTauB -> relaxation time for bottom layer
43   -! SRelaxTauS -> relaxation time for surface layer
44   -! SRelaxBott -> height of bottom relaxation layer
45   -! (set to 0. if not used)
46   -! SRelaxSurf -> height of surface relaxation layer
47   -! (set to 0. if not used)
48   -!
49   -!-------------------------------------------------------------------------------
50   - &sprofile
51   - s_prof_method= 2
52   - s_analyt_method= 2
53   - z_s1= 25.
54   - s_1= 31.
55   - z_s2= 35.
56   - s_2= 32.
57   - s_obs_NN= 2.56e-4
58   - s_prof_file= 'franklin_sprof_ctd.dat'
59   - SRelaxTauM= 1209600.
60   - SRelaxTauB= 1209600.
61   - SRelaxTauS= 1209600.
62   - SRelaxBott= 0.
63   - SRelaxSurf= 0.
64   - /
65   -
66   -!-------------------------------------------------------------------------------
67   -! observed or prescribed nitrate profiles
68   -!
69   -! n_prof_method -> method to create initial or observed nitrate profiles
70   -! 0: no initial values, N-equation is not solved
71   -! 1: use analytically prescribed initial profile
72   -! 2: read profiles at different dates from "n_prof_file"
73   -! and interpolate to GOTM timestep
74   -!
75   -! n_analyt_method -> method to create analytically precribed inital profile
76   -! 1: set profile to constant value n_1
77   -! 2: set "two layer" stratification (see user's guide)
78   -!
79   -! z_n1 -> upper layer thickness if n_analyt_method=2
80   -!
81   -! n_1 -> constant nitrate if n_analyt_method=1
82   -! upper layer nitrate if n_analyt_method=2
83   -!
84   -! z_n2 -> depth of top of lower layer if n_analyt_method=2
85   -!
86   -! n_2 -> lower layer salinity if n_analyt_method=2
87   -!
88   -
89   -! n_prof_file -> filename of file with nitrate profiles
90   -! if n_prof_method=2
91   -!
92   -! The computed profiles can be relaxed towards observed or prescribed
93   -! profiles with a certain relaxation time. If you do not want relaxation,
94   -! set the relaxation times to 1.e15 (something large). It is possible to choose
95   -! different relaxation times in a surface and bottom layer.
96   -!
97   -! NRelaxTauM -> relaxation time for bulk of the flow
98   -! NRelaxTauB -> relaxation time for bottom layer
99   -! NRelaxTauS -> relaxation time for surface layer
100   -! NRelaxBott -> height of bottom relaxation layer
101   -! (set to 0. if not used)
102   -! NRelaxSurf -> height of surface relaxation layer
103   -! (set to 0. if not used)
104   -!
105   -!-------------------------------------------------------------------------------
106   - &nprofile
107   - n_prof_method= 1
108   - n_analyt_method= 2
109   - z_n1= 35.
110   - n_1= 3.0
111   - z_n2= 70.
112   - n_2= 15.0
113   - n_prof_file= 'nprof_ctd.dat'
114   - NRelaxTauM= 1.e15
115   - NRelaxTauB= 1.e15
116   - NRelaxTauS= 1.e15
117   - NRelaxBott= 0.
118   - NRelaxSurf= 0.
119   - /
120   -
121   -!-------------------------------------------------------------------------------
122   -! observed or prescribed potential temperature profiles
123   -!
124   -! a_prof_method -> method to create initial or observed ammonium profiles
125   -! 0: no initial values, A-equation is not solved
126   -! 1: use analytically prescribed initial profile
127   -! 2: read profiles at different dates from "a_prof_file"
128   -! and interpolate to GOTM timestep
129   -!
130   -! a_analyt_method -> method to create analytically precribed inital profile
131   -! 1: set profile to constant value a_1
132   -! 2: set "two layer" stratification (see user's guide)
133   -! 3: set profile with constant N^2 (see user's guide)
134   -! This option can only be used toghether with
135   -! a_analyt_method=1 (constant ammonium).
136   -!
137   -! z_a1 -> upper layer thickness if a_analyt_method=2
138   -!
139   -! a_1 -> constant ammonium if a_analyt_method=1
140   -! upper layer ammonium if a_analyt_method=2
141   -! surface ammonium if a_analyt_method=3
142   -!
143   -! z_a2 -> depth of top of lower layer if a_analyt_method=2
144   -!
145   -! a_2 -> lower layer temperature if a_analyt_method=2
146   -!
147   -! a_obs_NN -> constant value N^2 corresponding to ammonium profile
148   -! if a_analyt_method=3
149   -
150   -! a_prof_file -> filename of file with ammonium profiles
151   -! if a_prof_method=2
152   -!
153   -! Computed profiles are relaxed towards observed or prescribed
154   -! profiles with a the same relaxation coefficients as for nitrate.
155   -! If you do not want relaxation, set the relaxation times to 1.e15 (something large).
156   -! It is possible to choose different relaxation times in a surface and bottom layer.
157   -!
158   -!-------------------------------------------------------------------------------
159   - &aprofile
160   - a_prof_method= 1
161   - a_analyt_method= 1
162   - z_a1= 25.
163   - a_1= 0.0
164   - z_a2= 35.
165   - a_2= 0.0
166   - a_prof_file= 'aprof.dat'
167   - /
168   -
169   -!-------------------------------------------------------------------------------
170   -! observed or prescribed potential temperature profiles
171   -!
172   -! t_prof_method -> method to create initial or observed temperature profiles
173   -! 0: no initial values, T-equation is not solved
174   -! 1: use analytically prescribed initial profile
175   -! 2: read profiles at different dates from "t_prof_file"
176   -! and interpolate to GOTM timestep
177   -!
178   -! t_analyt_method -> method to create analytically precribed inital profile
179   -! 1: set profile to constant value s_1
180   -! 2: set "two layer" stratification (see user's guide)
181   -! 3: set profile with constant N^2 (see user's guide)
182   -! This option can only be used toghether with
183   -! s_analyt_method=1 (constant salinity).
184   -!
185   -! z_t1 -> upper layer thickness if t_analyt_method=2
186   -!
187   -! t_1 -> constant temperature if t_analyt_method=1
188   -! upper layer temperature if t_analyt_method=2
189   -! surface temperature if t_analyt_method=3
190   -!
191   -! z_t2 -> depth of top of lower layer if t_analyt_method=2
192   -!
193   -! t_2 -> lower layer temperature if t_analyt_method=2
194   -!
195   -! t_obs_NN -> constant value N^2 corresponding to temperature profile
196   -! if t_analyt_method=3
197   -
198   -! t_prof_file -> filename of file with temperature profiles
199   -! if t_prof_method=2
200   -!
201   -! The computed profiles can be relaxed towards observed or prescribed
202   -! profiles with a certain relaxation time. If you do not want relaxation,
203