[3357] | 1 | !------------------------------------------------------------------------------- |
---|
| 2 | !-- INITIALIZATION PARAMETER NAMELIST |
---|
| 3 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar |
---|
| 4 | !------------------------------------------------------------------------------- |
---|
| 5 | &initialization_parameters |
---|
| 6 | ! |
---|
| 7 | !-- grid parameters |
---|
| 8 | !------------------------------------------------------------------------------- |
---|
| 9 | nx = 23, ! Number of gridboxes in x-direction (nx+1) |
---|
| 10 | ny = 23, ! Number of gridboxes in y-direction (ny+1) |
---|
| 11 | nz = 40, ! Number of gridboxes in z-direction (nz) |
---|
| 12 | |
---|
| 13 | dx = 50.0, ! Size of single gridbox in x-direction |
---|
| 14 | dy = 50.0, ! Size of single gridbox in y-direction |
---|
| 15 | dz = 50.0, ! Size of single gridbox in z-direction |
---|
| 16 | |
---|
| 17 | dz_stretch_level = 1000.0, ! Height (in m) where stretching starts |
---|
| 18 | dz_stretch_factor = 1.08, ! dz(k+1) = dz(k) * dz_stretch_factor |
---|
| 19 | |
---|
| 20 | ! |
---|
| 21 | !-- initialization and vertical profiles |
---|
| 22 | !------------------------------------------------------------------------------- |
---|
| 23 | initializing_actions = 'set_constant_profiles initialize_bubble', |
---|
| 24 | ! initial conditions (with warm air bubble) |
---|
| 25 | |
---|
| 26 | ug_surface = 0.0, ! u-comp of geostrophic wind at surface |
---|
| 27 | surface_pressure = 1015.4, ! surface pressure |
---|
| 28 | |
---|
| 29 | pt_surface = 297.9, ! temperature at surface |
---|
| 30 | pt_vertical_gradient = 0.0, 0.58588957, ! vertical gradient of temperature |
---|
| 31 | pt_vertical_gradient_level = 0.0, 740.0, ! height level of temp gradients |
---|
| 32 | |
---|
| 33 | q_surface = 0.016, ! mixing ratio at surface |
---|
| 34 | q_vertical_gradient = -2.97297E-4, -4.5238095E-4, -8.108108E-5, ! gradient for mix. ratio |
---|
| 35 | q_vertical_gradient_level = 0.0, 740.0, 3260.0, ! height lev. for gradients |
---|
| 36 | |
---|
| 37 | ! |
---|
| 38 | !-- passive scalar |
---|
| 39 | !------------------------------------------------------------------------------- |
---|
| 40 | passive_scalar = .TRUE., ! enable passive scalar |
---|
| 41 | s_surface = 0.0008, ! equivalent to 420ppm background |
---|
| 42 | surface_scalarflux = 0.002, ! kg emission per square meter and second |
---|
| 43 | bc_s_b = 'neumann', ! Neumann boundary conditions |
---|
| 44 | scalar_advec = 'bc-scheme', ! Bott chlond scheme for scalar advection |
---|
| 45 | loop_optimization = 'vector', ! loops over k,j,i |
---|
| 46 | |
---|
| 47 | ! |
---|
| 48 | !-- bulk microphysics |
---|
| 49 | !------------------------------------------------------------------------------- |
---|
| 50 | humidity = .TRUE., ! prog. equation for total water mixing ratio |
---|
| 51 | |
---|
| 52 | ! |
---|
| 53 | !-- numerics |
---|
| 54 | !------------------------------------------------------------------------------- |
---|
| 55 | fft_method = 'temperton-algorithm', ! build-in fft method |
---|
| 56 | |
---|
| 57 | / ! end of initialization parameter namelist |
---|
| 58 | |
---|
| 59 | |
---|
| 60 | !------------------------------------------------------------------------------- |
---|
| 61 | !-- RUNTIME PARAMETER NAMELIST |
---|
| 62 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par |
---|
| 63 | !------------------------------------------------------------------------------- |
---|
| 64 | &runtime_parameters |
---|
| 65 | ! |
---|
| 66 | !-- run steering |
---|
| 67 | !------------------------------------------------------------------------------- |
---|
| 68 | end_time = 1200.0, ! simulation time of the 3D model |
---|
| 69 | |
---|
| 70 | create_disturbances = .TRUE., ! randomly perturbate horiz. velocity |
---|
| 71 | dt_disturb = 150.0, ! interval for random perturbations |
---|
| 72 | disturbance_energy_limit = 0.01, ! upper limit for perturbation energy |
---|
| 73 | |
---|
| 74 | data_output_2d_on_each_pe = .FALSE., ! don't do 2D output on each MPI rank |
---|
| 75 | |
---|
| 76 | ! |
---|
| 77 | !-- data output |
---|
| 78 | !------------------------------------------------------------------------------- |
---|
| 79 | netcdf_data_format = 2, ! use NetCDF3 |
---|
| 80 | dt_run_control = 0.0, ! output interval for run control |
---|
| 81 | dt_data_output = 400.0, ! output interval for general data |
---|
| 82 | dt_dopr = 400.0, ! output interval for profile data |
---|
| 83 | |
---|
| 84 | data_output = 'w_xy', |
---|
| 85 | 'w_xz', |
---|
| 86 | 'pt_xy', |
---|
| 87 | 'pt_xz', |
---|
| 88 | 'ql_yz', |
---|
| 89 | 'w_yz', |
---|
| 90 | 'u_yz', |
---|
| 91 | 'v_yz', |
---|
| 92 | 's_yz', |
---|
| 93 | 'nc_yz', |
---|
| 94 | |
---|
| 95 | data_output_pr = 'u', |
---|
| 96 | 'v', |
---|
| 97 | 'w', |
---|
| 98 | 'ql', |
---|
| 99 | '#pt', |
---|
| 100 | 'w"pt"', |
---|
| 101 | 'w*pt*', |
---|
| 102 | 'wpt', |
---|
| 103 | 'w*2', |
---|
| 104 | 'pt*2', |
---|
| 105 | |
---|
| 106 | section_xy = 2, 10, ! grid index for 2D XY cross sections |
---|
| 107 | section_xz = 20, ! grid index for 2D XZ cross sections |
---|
| 108 | section_yz = -1, 10, ! grid index for 2D YZ cross sections |
---|
| 109 | |
---|
| 110 | averaging_interval_pr = 100.0, ! averaging interval profile data |
---|
| 111 | dt_averaging_input_pr = 10.0, ! averaging profile data sampling rate |
---|
| 112 | |
---|
| 113 | / ! end of runtime parameter namelist |
---|
| 114 | |
---|
| 115 | |
---|
| 116 | !------------------------------------------------------------------------------- |
---|
| 117 | !-- BULK CLOUD PARAMETER NAMELIST |
---|
| 118 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/bcmpar |
---|
| 119 | !------------------------------------------------------------------------------- |
---|
| 120 | &bulk_cloud_parameters |
---|
| 121 | |
---|
| 122 | bulk_cloud_model = .TRUE., ! enable bulk cloud model |
---|
| 123 | cloud_scheme = 'morrison', ! two-moment liquid water microphysics |
---|
| 124 | ! including activation, diffusion, autoconversion, |
---|
| 125 | ! accretion, sedimentation, precipitation |
---|
| 126 | |
---|
| 127 | aerosol_bulk = 'nacl', ! chemical composition of aerosol |
---|
| 128 | |
---|
| 129 | na_init = 1.5E8, ! background aerosol concetration |
---|
| 130 | dry_aerosol_radius = 0.15E-6, ! mean dry aerosol radius |
---|
| 131 | sigma_bulk = 1.61, ! width of dry aerosol spectrum |
---|
| 132 | |
---|
| 133 | curvature_solution_effects_bulk = .TRUE., ! considering aerosol properties for activation |
---|
| 134 | collision_turbulence = .TRUE., ! parametize turbulence effects for collision |
---|
| 135 | cloud_water_sedimentation = .TRUE., ! enable cloud water sedimentation |
---|
| 136 | |
---|
| 137 | call_microphysics_at_all_substeps = .TRUE., ! microphysics every timestep |
---|
| 138 | |
---|
| 139 | / ! end of bulk cloud parameter namelist |
---|