[4002] | 1 | !------------------------------------------------------------------------------- |
---|
| 2 | !-- INITIALIZATION PARAMETER NAMELIST |
---|
| 3 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar |
---|
| 4 | !------------------------------------------------------------------------------- |
---|
[4211] | 5 | ! This is a setup for generating an ocean mixed layer including the effects |
---|
| 6 | ! of Stokes drift (Langmuir circulation) and wave breaking, based on the |
---|
| 7 | ! case described in Noh, Y., H.S. Min and S. Raasch (2004): Large-eddy |
---|
| 8 | ! simulation of the ocean mixed layer: The effects of wave breaking and |
---|
| 9 | ! Langmuir circulation. J. Phys. Oceanogr., 34, 720-735. |
---|
| 10 | !------------------------------------------------------------------------------- |
---|
[4002] | 11 | &initialization_parameters |
---|
| 12 | ! |
---|
| 13 | !-- grid parameters |
---|
| 14 | !------------------------------------------------------------------------------- |
---|
[4211] | 15 | nx = 63, ! Number of gridboxes in x-direction (nx+1) |
---|
| 16 | ! For Noh et al. change to nx = 239 |
---|
| 17 | ny = 63, ! Number of gridboxes in y-direction (ny+1) |
---|
| 18 | ! For Noh et al. change to nx = 239 |
---|
| 19 | nz = 64, ! Number of gridboxes in z-direction (nz) |
---|
[4002] | 20 | |
---|
[4211] | 21 | dx = 1.25, ! Size of single gridbox in x-direction |
---|
| 22 | dy = 1.25, ! Size of single gridbox in y-direction |
---|
| 23 | dz = 1.25, ! Size of single gridbox in z-direction |
---|
[4002] | 24 | ! |
---|
| 25 | !-- initialization |
---|
| 26 | !------------------------------------------------------------------------------- |
---|
[4211] | 27 | initializing_actions = 'set_constant_profiles', ! initial conditions |
---|
[4002] | 28 | |
---|
[4211] | 29 | latitude = 55.0, ! corresponds to a Coriolis parameter |
---|
| 30 | ! f = 1.2**-4 s-1 |
---|
[4002] | 31 | |
---|
[4211] | 32 | ug_surface = 0.0, ! u-comp of geostrophic flow at surface |
---|
| 33 | vg_surface = 0.0, ! v-comp of geostrophic flow at surface |
---|
[4002] | 34 | |
---|
[4211] | 35 | pt_surface = 300.0, ! initial surface potential temp |
---|
[4002] | 36 | ! |
---|
| 37 | !-- boundary conditions |
---|
| 38 | !------------------------------------------------------------------------------- |
---|
[4211] | 39 | bc_uv_b = 'neumann', |
---|
| 40 | bc_uv_t = 'neumann', |
---|
| 41 | bc_pt_b = 'neumann', |
---|
| 42 | bc_pt_t = 'neumann', |
---|
| 43 | bc_p_b = 'neumann', |
---|
| 44 | bc_p_t = 'neumann', |
---|
[4002] | 45 | |
---|
[4211] | 46 | use_top_fluxes = .T., |
---|
| 47 | use_surface_fluxes = .F., |
---|
| 48 | constant_flux_layer = .F., |
---|
[4002] | 49 | |
---|
[4211] | 50 | top_momentumflux_u = -0.0001, ! corresponds to u* = 0.01 m/s |
---|
| 51 | top_momentumflux_v = 0.0, |
---|
[4002] | 52 | |
---|
[4211] | 53 | top_heatflux = 1.02e-4, ! gives weak cooling at ocean |
---|
| 54 | ! surface to initiate buoyancy driven |
---|
| 55 | ! turbulence, |
---|
| 56 | ! see surface_cooling_spinup_time |
---|
| 57 | ! below |
---|
[4002] | 58 | ! |
---|
| 59 | !-- numerics |
---|
| 60 | !------------------------------------------------------------------------------- |
---|
[4211] | 61 | fft_method = 'temperton-algorithm', ! build-in fft method |
---|
[4002] | 62 | |
---|
[4211] | 63 | momentum_advec = 'ws-scheme', ! attention: Noh et al. used the |
---|
| 64 | ! Piascek-Williams advection scheme |
---|
| 65 | scalar_advec = 'ws-scheme', ! attention: Noh et al. used the |
---|
| 66 | ! Piascek-Williams advection scheme |
---|
[4002] | 67 | |
---|
| 68 | / ! end of initialization parameter namelist |
---|
| 69 | |
---|
| 70 | !------------------------------------------------------------------------------- |
---|
| 71 | !-- RUNTIME PARAMETER NAMELIST |
---|
| 72 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par |
---|
| 73 | !------------------------------------------------------------------------------- |
---|
| 74 | &runtime_parameters |
---|
| 75 | ! |
---|
| 76 | !-- run steering |
---|
| 77 | !------------------------------------------------------------------------------- |
---|
[4211] | 78 | end_time = 120.0, ! simulation time of the 3D model |
---|
| 79 | ! For Noh et al. change to 28800.0 s |
---|
[4002] | 80 | |
---|
[4211] | 81 | create_disturbances = .TRUE., ! randomly perturbate horiz. velocity |
---|
| 82 | ! at beginning of the run |
---|
| 83 | disturbance_amplitude = 0.25e-3, ! maximum perturbation amplitude |
---|
| 84 | |
---|
[4002] | 85 | ! |
---|
| 86 | !-- data output |
---|
| 87 | !------------------------------------------------------------------------------- |
---|
[4211] | 88 | netcdf_data_format = 2, ! output in NetCDF3 64bit offset format |
---|
[4002] | 89 | |
---|
[4211] | 90 | dt_run_control = 0.0, ! output interval for run control |
---|
[4002] | 91 | |
---|
[4211] | 92 | ! remove the following line to create output of time series |
---|
| 93 | dt_dots = 100000.0, |
---|
| 94 | |
---|
| 95 | dt_dopr = 120.0, ! For Noh et al. change to 1800.0 |
---|
| 96 | |
---|
| 97 | averaging_interval_pr = 120.0, ! For Noh et al. change to 600.0 |
---|
| 98 | dt_averaging_input_pr = 10.0 |
---|
| 99 | data_output_pr = 'e', 'e*', '#theta', |
---|
| 100 | 'prho', 'u', 'v', 'w', 'km', 'kh', 'l', |
---|
| 101 | 'u*2', 'v*2', 'w*2', 'theta*2', 'w*3', |
---|
| 102 | 'wu', 'w*u*', 'w"u"', |
---|
| 103 | 'wv', 'w*v*', 'w"v"', |
---|
| 104 | 'w*e*', 'w*p*', |
---|
| 105 | 'w*u*u*:dz', 'w*p*:dz', 'w"e:dz', |
---|
| 106 | |
---|
| 107 | ! uncomment following lines to create some of the output that is shown in |
---|
| 108 | ! the Noh et al. paper |
---|
| 109 | ! dt_data_output = 3600.0, |
---|
| 110 | ! data_output = 'w', 'theta', |
---|
| 111 | |
---|
[4002] | 112 | / ! end of runtime parameter namelist |
---|
| 113 | |
---|
[4211] | 114 | |
---|
[4002] | 115 | !------------------------------------------------------------------------------- |
---|
| 116 | !-- OCEAN PARAMETER NAMELIST |
---|
| 117 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/oceanpar |
---|
| 118 | !------------------------------------------------------------------------------- |
---|
| 119 | &ocean_parameters |
---|
[4211] | 120 | |
---|
| 121 | stokes_waveheight = 1.0 ! results in a turbulent Langmuir |
---|
| 122 | ! number La = 0.45 |
---|
| 123 | stokes_wavelength = 40.0, ! results in a turbulent Langmuir |
---|
| 124 | ! number La = 0.45 |
---|
| 125 | wave_breaking = .F., ! The implemented parameterization is |
---|
| 126 | ! designed for a vertical grid spacing |
---|
| 127 | ! of dz = 1.25m and time steps of |
---|
| 128 | ! about 4 s. It will probably fail for |
---|
| 129 | ! other setups. |
---|
| 130 | ! Attention: |
---|
| 131 | ! For Noh et al. switch to .T. |
---|
| 132 | surface_cooling_spinup_time = 900.0, ! add a surface cooling only at start |
---|
| 133 | ! in order to initiate turbulence |
---|
| 134 | |
---|
| 135 | salinity = .TRUE., ! salinity switched on just for |
---|
| 136 | ! testing the salinity code |
---|
| 137 | ! remove this and the following two |
---|
| 138 | ! parameters for simulating the |
---|
| 139 | ! Noh et al. case |
---|
[4002] | 140 | bc_sa_t = 'neumann', |
---|
[4211] | 141 | top_salinityflux = 0.0, ! zero salinityflux at ocean surface |
---|
[4002] | 142 | |
---|
[4211] | 143 | / ! end of ocean parameters |
---|
[4002] | 144 | |
---|
| 145 | ! WARNING: do not remove this blank line from the end of this file |
---|