source: palm/trunk/TUTORIALS/cases/example_neutral/INPUT/example_neutral_p3d @ 4098

Last change on this file since 4098 was 4098, checked in by Giersch, 5 years ago

Tutorial cases updated for getting physically meaningful results (only gfortran with 4 cores is tested)

File size: 2.2 KB
Line 
1 &initialization_parameters
2
3     !-- Numerical grid   
4         nx = 49,
5         ny = 49,
6         nz = 100,
7
8         dx = 20.0,
9         dy = 20.0,
10         dz = 20.0,
11
12     !-- Pressure solver 
13         fft_method = 'temperton-algorithm',
14
15     !-- Initialization
16!         initializing_actions = 'set_constant_profiles',
17
18         initializing_actions = 'set_1d-model_profiles',
19         end_time_1d = 864000.0,
20         damp_level_1d = 1250.0,
21         mixing_length_1d = 'blackadar',
22         dissipation_1d = 'detering',
23
24         dt_run_control_1d = 3600.0,
25         dt_pr_1d = 86400.0,
26
27     !-- Rayleigh damping
28         rayleigh_damping_height = 1250.0,
29         rayleigh_damping_factor = 0.01,
30         galilei_transformation = .T.,
31
32     !-- Roughness length         
33         roughness_length = 0.1,
34
35     !-- Geostrophic wind (large-scale pressure gradient)
36         ug_surface = 5.0,
37         vg_surface = 0.0,
38
39     !-- Neutral stratification
40         pt_surface = 300.0,
41         pt_vertical_gradient = 0.0,
42         pt_vertical_gradient_level = 0.0,
43         neutral = .T.,
44
45     !-- Boundary conditions
46         bc_uv_b =  'dirichlet',
47         bc_uv_t = 'neumann',
48
49         bc_p_b = 'neumann',
50         bc_p_t = 'dirichlet',
51
52         constant_flux_layer = .T.,
53         /
54
55 &runtime_parameters 
56
57     !-- Simulation time 2 days
58         end_time = 172800.0,
59
60     !-- Disturbances
61         create_disturbances = .T.,
62         dt_disturb = 150.0,
63         disturbance_energy_limit = 0.01,
64 
65     !-- Run control
66         dt_run_control = 60.0,
67
68     !-- Time series
69         dt_dopr = 10800.0,
70         
71     !-- Vertical profiles
72         averaging_interval_pr = 10800.0,
73         dt_averaging_input_pr = 50.0,
74         data_output_pr = '#u', '#v', 'w', '#theta',
75                          'w"u"', 'w*u*', 'wu', 'w"v"', 'w*v*', 'wv','w*2', 'theta*2',
76                         
77         /
78         
79 &spectra_parameters
80 
81     !-- Output settings for spectra
82         data_output_sp =  'u', 'v', 'w',
83         dt_dosp = 10800.0,
84         comp_spectra_level =  2, 5, 10, 15, 20,
85         spectra_direction =  'x' , 'y' , 'x',
86         averaging_interval_sp = 10800.0,
87         
88         /
89 
Note: See TracBrowser for help on using the repository browser.