Home

Context Navigation

warm_air_bubble_with_bcm_p3d @ 3423

Last change on this file since 3423 was 3357, checked in by knoop, 6 years ago
Added warm_air_bubble_with_bcm test case
File size: 6.6 KB

Line
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	'wpt',
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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format