Changeset 767 for palm/trunk/SOURCE/modules.f90

Timestamp:

Oct 14, 2011 6:39:12 AM (13 years ago)

Author:

raasch

Message:

New:
---

Flow field initialization with given (e.g. measured) profiles. Profile data
for u-,v-velocity components + respective heights are given with new
inipar-parameters u_profile, v_profile, and uv_heights. Final profiles are
calculated from these given profiles by linear interpolation.
(check_parameters, header, init_3d_model, modules, parin)

Changed:

---

ug,vg replaced by u_init,v_init as the Dirichlet top boundary condition
(boundary_conds)

dirichlet_0 conditions moved from init_3d_model to
check_parameters (check_parameters, init_3d_model)

Errors:

---

bugfix: dirichlet_0 conditions moved from init_3d_model to
check_parameters (check_parameters, init_3d_model)

File:

• palm/trunk/SOURCE/modules.f90 (modified) (3 diffs)

palm/trunk/SOURCE/modules.f90

@@ -5 +5 @@
 ! Current revisions:
 ! -----------------
-! 
+! +u_profile, v_profile, uv_heights, use_prescribed_profile_data
 !
 ! Former revisions:
@@ -536 +536 @@
                 stop_dt = .FALSE., synchronous_exchange = .FALSE., &
                 terminate_run = .FALSE., turbulent_inflow = .FALSE., &
+                use_prescribed_profile_data = .FALSE., &
                 use_prior_plot1d_parameters = .FALSE., use_reference = .FALSE.,&
                 use_surface_fluxes = .FALSE., use_top_fluxes = .FALSE., &
                 use_ug_for_galilei_tr = .TRUE., use_upstream_for_tke = .FALSE.,&
-                wall_adjustment = .TRUE., ws_scheme_sca = .FALSE.,             &
+                wall_adjustment = .TRUE., ws_scheme_sca = .FALSE., &
                 ws_scheme_mom = .FALSE.
 
@@ -636 +637 @@
              time_domask(max_masks) = 0.0, &
              tsc(10) = (/ 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /), &
+             u_profile(100) = 9999999.9, uv_heights(100) = 9999999.9, &
+             v_profile(100) = 9999999.9, &
              ug_vertical_gradient(10) = 0.0, &
              ug_vertical_gradient_level(10) = -9999999.9, &