source: palm/trunk/TESTS/cases/oceanml/INPUT/oceanml_p3d

Last change on this file was 4211, checked in by raasch, 5 years ago

ocean mixed layer test case modified; can now be used as a basis to reproduce the Noh et.al. (2004) results

File size: 7.0 KB
RevLine 
[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
Note: See TracBrowser for help on using the repository browser.