Commit 72e7304e6365838df8f17a78c8170fbbb5f095ee

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

Ajout du repertoire nml contenant les namelists

nml/airsea.nml 0 → 100644
... ... @@ -0,0 +1,98 @@
  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 + /
... ...
nml/bio.nml 0 → 100644
... ... @@ -0,0 +1,56 @@
  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 + /
... ...
nml/bio_fasham.nml 0 → 100644
... ... @@ -0,0 +1,91 @@
  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 + /
... ...
nml/gotmmean.nml 0 → 100644
... ... @@ -0,0 +1,50 @@
  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 + /
... ...
nml/gotmrun.nml 0 → 100644
... ... @@ -0,0 +1,138 @@
  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 + /
... ...
nml/gotmturb.nml 0 → 100644
... ... @@ -0,0 +1,304 @@
  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 + /
... ...
nml/kpp.nml 0 → 100644
... ... @@ -0,0 +1,40 @@
  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 + /
... ...
nml/obs.nml 0 → 100644
... ... @@ -0,0 +1,564 @@
  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