obs.nml 17.9 KB
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453
!$Id: obs.proto,v 1.1.1.1 2003/03/11 13:38:58 kbk Exp $
!-------------------------------------------------------------------------------
! 
!-------------------------------------------------------------------------------
! observed or prescribed salinity profiles
!
! s_prof_method    -> method to create initial or observed salinity profiles
!                  0: no initial values, S-equation is not solved
!                  1: use analytically prescribed initial profile
!                  2: read profiles at different dates from "s_prof_file"
!                     and interpolate to GOTM timestep
!
! s_analyt_method  -> method to create analytically precribed inital profile
!                  1: set profile to constant value s_1
!                  2: set "two layer" stratification (see user's guide)
!                  3: set profile with constant N^2 (see user's guide)
!                     This option can only be used toghether with
!                     t_analyt_method=1 (constant temperature).
!
! z_s1             -> upper layer thickness if s_analyt_method=2
!
! s_1              -> constant salinity     if s_analyt_method=1
!                     upper layer salinity  if s_analyt_method=2
!                     surface salinity      if s_analyt_method=3
!
! z_s2             -> depth of top of lower layer if s_analyt_method=2
!
! s_2              -> lower layer salinity if s_analyt_method=2
!
! s_obs_NN         -> constant value N^2 corresponding to salinity profile
!                     if s_analyt_method=3

! s_prof_file      -> filename of file with salinity profiles
!                     if s_prof_method=2
!
! The computed profiles can be relaxed towards observed or prescribed
! profiles with a certain relaxation time. If you do not want relaxation,
! set the relaxation times to 1.e15 (something large). It is possible to choose
! different relaxation times in a surface and bottom layer.
!
! SRelaxTauM       -> relaxation time for bulk of the flow 
! SRelaxTauB       -> relaxation time for bottom layer 
! SRelaxTauS       -> relaxation time for surface layer
! SRelaxBott       -> height of bottom relaxation layer
!                     (set to 0. if not used)
! SRelaxSurf       -> height of surface relaxation layer
!                     (set to 0. if not used)
!
!-------------------------------------------------------------------------------
 &sprofile
  s_prof_method=   1
  s_analyt_method= 1
  z_s1=            5.
  s_1=             15.
  z_s2=            5.
  s_2=             15.
  s_obs_NN=        2.56e-4
  s_prof_file=     'sprof.dat'
  SRelaxTauM=      44714.
  SRelaxTauB=      44714.
  SRelaxTauS=      44714.
  SRelaxBott=      0.
  SRelaxSurf=      0.
 /

!-------------------------------------------------------------------------------
! observed or prescribed potential temperature profiles
!
! t_prof_method    -> method to create initial or observed temperature profiles
!                  0: no initial values, T-equation is not solved
!                  1: use analytically prescribed initial profile
!                  2: read profiles at different dates from "t_prof_file"
!                     and interpolate to GOTM timestep
!
! t_analyt_method  -> method to create analytically precribed inital profile
!                  1: set profile to constant value s_1
!                  2: set "two layer" stratification (see user's guide)
!                  3: set profile with constant N^2 (see user's guide)
!                     This option can only be used toghether with
!                     s_analyt_method=1 (constant salinity).
!
! z_t1             -> upper layer thickness if t_analyt_method=2
!
! t_1              -> constant temperature     if t_analyt_method=1
!                     upper layer temperature  if t_analyt_method=2
!                     surface temperature      if t_analyt_method=3
!
! z_t2             -> depth of top of lower layer if t_analyt_method=2
!
! t_2              -> lower layer temperature if t_analyt_method=2
!
! t_obs_NN         -> constant value N^2 corresponding to temperature profile
!                     if t_analyt_method=3

! t_prof_file      -> filename of file with temperature profiles
!                     if t_prof_method=2
!
! The computed profiles can be relaxed towards observed or prescribed
! profiles with a certain relaxation time. If you do not want relaxation,
! set the relaxation times to 1.e15 (something large). It is possible to choose
! different relaxation times in a surface and bottom layer.
!
! TRelaxTauM       -> relaxation time for bulk of the flow 
! TRelaxTauB       -> relaxation time for bottom layer 
! TRelaxTauS       -> relaxation time for surface layer
! TRelaxBott       -> height of bottom relaxation layer
!                     (set to 0. if not used)
! TRelaxSurf       -> height of surface relaxation layer
!                     (set to 0. if not used)
!
!-------------------------------------------------------------------------------
 &tprofile
  t_prof_method=   0
  t_analyt_method= 3
  z_t1=            30.
  t_1=             20.
  z_t2=            40.
  t_2=             15.
  t_obs_NN=        1.0e-4
  t_prof_file=     'tprof.dat'
  TRelaxTauM=      1.e15
  TRelaxTauB=      1.e15
  TRelaxTauS=      1.e15
  TRelaxBott=      0.
  TRelaxSurf=      0.
 /

!-------------------------------------------------------------------------------
! external pressure gradients
!
! ext_press_method -> method to compute external pressure gradient from data
!                  0: constant external pressure gradient
!                  1: external pressure gradient from tidal constituents
!                  2: external pressure gradient from data
!                     found in "ext_press_file". 
!                  The data supplied from 0-2 are interpreted differently
!                  depending on the value of "ext_press_mode" !
!
!
! ext_press_mode   -> how to interprete ANY prescribed data for 
!                     the computation of the external pressure gradient
!                  0: interprete all data as surface elevation gradients
!                  1: interprete all data as current meter measurements
!                     at a given height "PressHeight" 
!                     (see documentation)
!                  2: interprete all data as vertically averaged current speeds
!                     (see documentation)
!
! ext_press_file   -> filename of input file for ext_press_method=2. 
!                     The data in the file are interpreted according to the 
!                     value of "ext_press_mode".
!
! PressConstU      -> constant pressure gradient data for x-direction for 
!                     ext_press_method=0. The value is interpreted according to
!                     the value of "ext_press_mode".
!
! PressConstV      -> constant pressure gradient data for y-direction for 
!                     ext_press_method=0. The value is interpreted according to
!                     the value of "ext_press_mode".
!
! PressHeight      -> height above bottom for current meter observations
!                     (only used for ext_press_mode=1)
!
! The following data specify the tidal constituents for ext_press_method=1. The
! tidal amplitudes are interpreted according to the value of "ext_press_mode".
!
! PeriodM          -> period of 1. harmonic (eg. M2-tide)
! AmpMu            -> u amplitude of 1. harmonic
! AmpMv            -> v amplitude of 1. harmonic
! PhaseMu          -> u phase of 1. harmonic
! PhaseMv          -> v phase of 1. harmonic
! PeriodS          -> period of 2. harmonic (eg. S2-tide)
! AmpSu            -> v amplitude of 2. harmonic
! AmpSv            -> v amplitude of 2. harmonic
! PhaseSu          -> v phase of 2. harmonic
! PhaseSv          -> v phase of 2. harmonic
!
!-------------------------------------------------------------------------------
 &ext_pressure
  ext_press_method=1
  ext_press_mode=  2
  ext_press_file=  'pressure.dat'
  PressConstU=     -0.02
  PressConstV=     0.0
  PressHeight=     0.0
  PeriodM=         44714.0
  AmpMu=           0.5
  AmpMv=           0.0
  PhaseMu=         0.0
  PhaseMv=         0.0
  PeriodS=         43200.0
  AmpSu=           0.0
  AmpSv=           0.0
  PhaseSu=         0.0
  PhaseSv=         0.0
 /

!-------------------------------------------------------------------------------
! internal pressure gradients
!
!
! int_press_method -> method to compute internal pressure gradient from data
!                  0: no internal pressure gradient
!                  1: vertically and temporally constant internal pressure 
!                     gradient computed from fixed horizontal salinity and 
!                     temperature gradients
!                  2: vertically and temporally variable int. pressure gradient
!                     computed from data in file "int_press_file"
!
! int_press_file   -> file with profiles of dsdx,dsdy,dtdx,dtdy 
!                     (only read for int_press_method=2)
!
! const_dsdx       -> x-gradient of S [psu/m] for int_press_method=1
! const_dsdy       -> y-gradient of S [psu/m] for int_press_method=1
! const_dtdx       -> x-gradient of T [K/m]   for int_press_method=1
! const_dtdy       -> y-gradient of T [K/m]   for int_press_method=1
!
! For all values of "int_press_method", it can be specified if 
! the given temperature and salinity gradients should also be taken into 
! account in the lateral advection terms of temperature and salinity.
!
! s_adv             -> advection of salinity    (.true./.false.)
! t_adv             -> advection of temperature (.true./.false.)
!
!-------------------------------------------------------------------------------
 &int_pressure
  int_press_method=1
  int_press_file=  'intern_press.dat'
  const_dsdx=      -0.0002
  const_dsdy=      0.0
  const_dtdx=      0.0
  const_dtdy=      0.0
  s_adv=           .true.
  t_adv=           .false.
 /

!-------------------------------------------------------------------------------
! Light extinction - Jerlov type or from file
!
! extinct_method   -> which method used to compute the extinction coefficient
!                  0: from file
!                  1: Jerlov type I
!                  2: Jerlov type 1 (upper 50 m)
!                  3: Jerlov type IA
!                  4: Jerlov type IB
!                  5: Jerlov type II
!                  6: Jerlov type III
!                  7: Adolf Stips, Lago Maggiore 
! 
! extinct_file     -> name of file used if extinct_method=0
!
!-------------------------------------------------------------------------------
 &extinct
  extinct_method=  1
  extinct_file=    'extinction.dat'
 /

!-------------------------------------------------------------------------------
! vertical advection
!
! It is possible to specify a vertical advection velocity. All fields (except 
! the turbulence quantities) are then vertically advected with this velocity,
! for example to mimic the effect of horizontal divergence leading to vertical
! velocities. The profiles of vertical velocity are determined by two values,
! the height of maximum absolute value of vertical velocity, w_height, and the
! vertical velocity at this height, w_adv. From w_height, the vertical velocity
! is linearly decreasing towards the surface and the bottom, where is value
! is zero.
!  
! w_adv_method     -> method to specify the vertical advection velocity
!                  0: zero vertical advection velocity
!                  1: temporally constant vertical advection velocity profile
!                  2: temporally varying vertical advection velocity profile
!                     with "w_adv" and "w_height" both read in from file 
!                     "w_adv_file")
!
! w_adv_height0    -> constant height of maximum absolute value of vertical 
!                     advection velocity
!
! w_adv0           -> constant vertical advection velocity for 
!                     w_adv_method=1

! w_adv_file       -> filename for file with "w_height" and "w_adv"
!                     (for w_adv_method=2)
!
! w_adv_discr      -> method to discretize the vertical advection term    
!                  1: first order upstream
!                  2: second-order polynomial (P2)
!                  3: third-order polynomial  (P3)
!                  4: TVD with Superbee limiter
!                  5: TVD with MUSCL limiter
!                  6: TVD with ULTIMATE QUICKEST
!
!-------------------------------------------------------------------------------
 &w_advspec
  w_adv_method=    0
  w_adv_height0=   0.
  w_adv0=          0.
  w_adv_file=      'w_adv.dat'
  w_adv_discr=     4
 /

!-------------------------------------------------------------------------------
! sea surface elevations
!
! Varying sea surface elevations can be specified, where the depth of the water
! column in GOTM will vary, for example to account for water level 
! fluctuations in tidally influenced estuaries.
!
! zeta_method      -> method to prescribe the sea surface elevation
!                  0: constant sea surface elevation 
!                  1: sea surface elevation from tidal constituents
!                  2: sea surface elevation from file
!
! zeta_file        -> filename for file with sea surface elevations
!                     (if zeta_method=2)
!
! zeta0            -> constant value for sea surface elevation
!                     (if zeta_method=0)
!
! The following 6 variables are used only if zeta_method=1
! period1          -> period of 1. harmonic (eg. M2-tide) - [s]
! amp1             -> amplitude of 1. harmonic  - [m]
! phase1           -> phase of 1. harmonic - [s]
! period2          -> period lsof 2. harmonic (eg. S2-tide) - [s]
! amp2             -> amplitude of 2. harmonic  - [m]
! phase2           -> phase of 2. harmonic - [s]
!
!-------------------------------------------------------------------------------
 &zetaspec
  zeta_method=     1
  zeta_file=       'zeta.dat'
  zeta_0=          0.00000
  period_1=        44714.0
  amp_1=           1.00000
  phase_1=         0.00000
  period_2=        43200.0
  amp_2=           0.50000
  phase_2=         0.00000
 /

!-------------------------------------------------------------------------------
! wind waves
!
! Wind induced waves can be specified. Used in e.g. the spm module to account 
! for wave-current interaction.
!
! wave_method      -> method to prescribe the wind waves
!                  0: nothing is done - all varibles equal 0
!                  1: constant values
!                  2: time varying variables read from file
!
! wave_file        -> filename for file wind waves
!                     (if wave_method=2)
!
! The following variables are used if wave_method=1
! Hs               -> constant value for significant wave-height
! Tz               -> constant value for mean zero-crossing period
! phiw             -> constant value for mean zero-crossing period
!-------------------------------------------------------------------------------
 &wave_nml
  wave_method=     0
  wave_file=       'wave.dat'
  Hs=              0.
  Tz=              0.
  phiw=            0.
 /

!-------------------------------------------------------------------------------
! observed or prescribed velocity profiles
!
! vel_prof_method  -> method to create initial or observed velocity profiles
!                  0: start from zero inital velocities
!                  1: (not implemented)
!                  2: read profiles at different dates from "vel_prof_file"
!                     and interpolate to GOTM timestep
!
! vel_prof_file    -> filename of file with u- and v profiles
!                     if vel_prof_method=2
!
! vel_relax_tau    -> relaxation time for velocity in
!                      (set to 1.e15 s for no relaxation)
!
!  vel_relax_ramp  -> duration of initial relaxation in [s]
!                     (set to 1.e15 for constant relaxation)
!
!-------------------------------------------------------------------------------
 &velprofile
  vel_prof_method= 0
  vel_prof_file=   'velprof.dat'
  vel_relax_tau=   1.e15
  vel_relax_ramp=  1.e15
 /

!-------------------------------------------------------------------------------
! turbulent dissipation rate profiles
!
! e_prof_method    -> method to create initial or observed dissipation profiles
!                  0: no observed dissipation rate profiles
!                  1: (not implemented)
!                  2: read profiles at different dates from "e_prof_file"
!                     and interpolate to GOTM timestep
!
! e_obs_const      -> a constant filling value [W/kg]
!
! e_prof_file      -> filename of file with dissipation rate profiles
!                     if e_prof_method=2
!
!-------------------------------------------------------------------------------
 &eprofile
  e_prof_method=   0
  e_obs_const=     1.e-12
  e_prof_file=     'eprof.dat'
 /

!-------------------------------------------------------------------------------
! prescibred inital profile and inital surface boundary conditions for the 
! dynamic buoyancy equation (only used if buoy_method=2 in gotmrun.inp) 
!
! b_obs_surf       -> prescribed initial buoyancy at the surface
! b_obs_NN         -> prescribed initial value of N^2 (=buoyancy gradient)
! b_obs_sbf        -> prescribed constant surface buoyancy flux
!
!-------------------------------------------------------------------------------
 &bprofile
  b_obs_surf=      0.0
  b_obs_NN=        1.e-4
  b_obs_sbf=       0.0
 /

!-------------------------------------------------------------------------------
! observed oxygen profiles
!
! o2_prof_method   -> method to create initial or observed oxygen profiles
!                  0: no observed oxygen rate profiles
!                  1: (not used)
!                  2: read profiles at different dates from "o2_prof_file"
!                     and interpolate to GOTM timestep
!
! o2_unit          -> which unit are the observations reported in
!                  0: saturation (%)
!                  1: mmol/m3 ???
!                  2: ml/m3 ???
!
! o2_prof_file     -> filename of file with oxygen profiles
!                     if o2_prof_method=2
!
!-------------------------------------------------------------------------------
 &o2_profile
  o2_prof_method=0
  o2_units=1
  o2_prof_file='o2.dat'
 /