doc/tut/examples/2013_heinze_DYCOMS: dycoms_p3d

File dycoms_p3d, 2.5 KB (added by heinze, 12 years ago)

Input parameter file

Line 
1 &inipar nx = 159, ny = 159, nz = 320,
2         dx = 35.0, dy = 35.0, dz = 5.0,
3
4         fft_method = 'temperton-algorithm',
5
6         humidity = .T., cloud_physics = .T.,
7
8         momentum_advec = 'ws-scheme',
9         scalar_advec   = 'ws-scheme',
10
11         galilei_transformation = .F.,
12
13         initializing_actions = 'by_user',
14
15         bc_p_t = 'dirichlet',
16
17         ug_surface = 7.0, vg_surface = -5.5,
18
19         pt_reference = 300.0,
20       
21         pt_surface = 289.0,
22         !pt-gradient is set in User Interface (UI)
23
24         q_surface = 0.009, 
25         !q-gradient is set in UI
26
27         !ql is set in UI
28
29         subs_vertical_gradient       = -3.75E-04, ! in (m/s)/100m
30         subs_vertical_gradient_level = 0.0,
31
32         surface_heatflux = 0.0121, bc_pt_b = 'neumann',   ! = 15W/m2
33         surface_waterflux = 3.816E-05, bc_q_b = 'neumann',! = 115W/m2
34
35         surface_pressure = 1017.8,
36
37         phi = 31.5,
38
39         rayleigh_damping_factor = 0.05,
40         rayleigh_damping_height = 1300.0,
41         scalar_rayleigh_damping = .F.,
42
43         netcdf_precision = 'pr_NF90_REAL8', /
44
45
46 &d3par  end_time = 14400.0,
47         termination_time_needed = 900.0,
48
49         npex = 16,
50         npey = 8,
51
52         create_disturbances = .T.,
53         dt_disturb = 150.0, disturbance_energy_limit = 0.01,
54
55         dt_run_control = 0.0,
56
57         netcdf_data_format = 2,
58
59         nz_do3d = 200,
60
61         section_xy = 1,
62
63         dt_data_output = 1800.0,
64         data_output    = 'u','v','w','lpt','q','ql','p','e',
65                          'lwp*_xy',
66                   
67         dt_dopr = 1800.0,
68         averaging_interval_pr = 1800.0,
69         dt_averaging_input_pr = 20.0,
70
71         data_output_pr = '#u','#v','w','#pt','#vpt','#lpt','#q','#qv','ql',
72                          'e','e*','p','#km','#kh','#l',
73                          'w"u"','w*u*','wu', 'w"v"','w*v*','wv',
74                          'w"pt"', 'w*pt*', 'wpt',
75                          'w"vpt"', 'w*vpt*', 'wvpt',
76                          'w"q"', 'w*q*', 'wq',
77                          'w"qv"', 'w*qv*', 'wqv',
78                          'w*e*', 'u*2', 'v*2', 'w*2','pt*2',
79                          'w*3', 'Sw', 'w*2pt*', 'w*pt*2',
80                          'w*u*u*:dz','w*p*:dz', 'w"e:dz','hyp','q*2',
81         
82         dt_dots = 0.0,   
83
84         data_output_2d_on_each_pe = .T., /
85
86&userpar radiation_dycoms = .T.,
87 
88         data_output_pr_user = 'TCC_pr', 'rho_ref', 'w"ql"', 'w*ql*',
89                                'wql', 'F_rad', 'Q_r', /