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

Timestamp:

Mar 7, 2008 1:42:18 PM (17 years ago)

Author:

raasch

Message:

preliminary update for the turbulence recycling method

File:

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

palm/trunk/SOURCE/modules.f90

@@ -5 +5 @@
 ! Actual revisions:
 ! -----------------
-! +hor_index_bounds, hor_index_bounds_previous_run, numprocs_previous_run,
-! nx_on_file, ny_on_file, offset_ocean_*
+! +hor_index_bounds, hor_index_bounds_previous_run, id_inflow,
+! inflow_damping_*, mean_inflow_profiles, numprocs_previous_run, nx_on_file,
+! ny_on_file, offset_ocean_*, recycling_plane, recycling_width, turbulent_inflow
 ! -myid_char_14
 !
@@ -115 +116 @@
 
     REAL, DIMENSION(:), ALLOCATABLE ::                                         &
-          ddzu, dd2zu, dzu, ddzw, dzw, hyp, km_damp_x, km_damp_y, lad, l_grid, &
-          pt_init, q_init, rdf, sa_init, ug, u_init, u_nzb_p1_for_vfc, vg,     &
-          v_init, v_nzb_p1_for_vfc, zu, zw
+          ddzu, dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, km_damp_x,  &
+          km_damp_y, lad, l_grid, pt_init, q_init, rdf, sa_init, ug, u_init,   &
+          u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, zu, zw
 
     REAL, DIMENSION(:,:), ALLOCATABLE ::                                       &
-          c_u, c_v, c_w, dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg, pt_slope_ref,    &
-          qs, qswst_remote, ts, us, z0
+          c_u, c_v, c_w, dzu_mg, dzw_mg, f1_mg, f2_mg, f3_mg,                  &
+          mean_inflow_profiles, pt_slope_ref, qs, qswst_remote, ts, us, z0
 
     REAL, DIMENSION(:,:), ALLOCATABLE, TARGET ::                               &
@@ -309 +310 @@
                 nsor_ini = 100, n_sor, normalizing_region = 0, &
                 nz_do1d, nz_do3d = -9999, outflow_damping_width = -1, &
-                prt_time_count = 0, runnr = 0, skip_do_avs = 0, &
-                terminate_coupled = 0, terminate_coupled_remote = 0, &
-                timestep_count = 0
+                prt_time_count = 0, recycling_plane, runnr = 0, &
+                skip_do_avs = 0, terminate_coupled = 0, &
+                terminate_coupled_remote = 0, timestep_count = 0
 
     INTEGER ::  dist_nxl(0:1), dist_nxr(0:1), dist_nyn(0:1), dist_nys(0:1), &
@@ -357 +358 @@
                 random_heatflux = .FALSE., run_control_header = .FALSE., &
                 sloping_surface = .FALSE., stop_dt = .FALSE., &
-                terminate_run = .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.
+                terminate_run = .FALSE., turbulent_inflow = .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.
 
     LOGICAL ::  data_output_xy(0:1) = .FALSE., data_output_xz(0:1) = .FALSE., &
@@ -392 +394 @@
              dz_stretch_level = 100000.0, e_init = 0.0, e_min = 0.0, &
              end_time = 0.0, &
-             f = 0.0, fs = 0.0, g = 9.81, kappa = 0.4, km_constant = -1.0, &
+             f = 0.0, fs = 0.0, g = 9.81, inflow_damping_height = 9999999.9, &
+             inflow_damping_width = 9999999.9, kappa = 0.4, km_constant = -1.0,&
              km_damp_max = -1.0, lad_surface = 0.0, long_filter_factor = 0.0, &
              maximum_cpu_time_allowed = 0.0, molecular_viscosity = 1.461E-5, &
@@ -405 +408 @@
              q_surface = 0.0, q_surface_initial_change = 0.0, &
              rayleigh_damping_factor = -1.0, rayleigh_damping_height = -1.0, &
-             residual_limit = 1.0E-4, restart_time = 9999999.9, rho_reference, &
-             rho_surface, rif_max = 1.0, &
+             recycling_width = 9999999.9, residual_limit = 1.0E-4, &
+             restart_time = 9999999.9, rho_reference, rho_surface, &
+             rif_max = 1.0, &
              rif_min = -5.0, roughness_length = 0.1, sa_surface = 35.0, &
              simulated_time = 0.0, simulated_time_at_begin, sin_alpha_surface, &
@@ -961 +965 @@
 #endif
     CHARACTER(LEN=5)       ::  myid_char = ''
-    INTEGER                ::  myid=0, npex = -1, npey = -1, numprocs = 1, &
-                               numprocs_previous_run = -1,                 &
+    INTEGER                ::  id_inflow, myid=0, npex = -1, npey = -1,  &
+                               numprocs = 1, numprocs_previous_run = -1, &
                                tasks_per_node = -9999, threads_per_task = 1