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

Timestamp:

Feb 23, 2007 4:53:48 AM (18 years ago)

Author:

raasch

Message:

preliminary version of modified boundary conditions at top

File:

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

palm/trunk/SOURCE/modules.f90

@@ -5 +5 @@
 ! Actual revisions:
 ! -----------------
-! 
+! +constant_top_heatflux, top_heatflux, use_top_fluxes, +arrays for top fluxes,
+! +nzt_diff, default of bc_pt_t renamed "initial_gradient"
+! Bugfix: p is not a pointer
 !
 ! Former revisions:
@@ -75 +77 @@
 
     REAL, DIMENSION(:,:), ALLOCATABLE, TARGET ::                               &
-          qsws_1, qsws_2, rif_1, rif_2, shf_1, shf_2, usws_1, usws_2,          &
-          vsws_1, vsws_2
+          qsws_1, qsws_2, qswst_1, qswst_2, rif_1, rif_2, shf_1, shf_2,        &
+          tswst_1, tswst_2, usws_1, usws_2, vsws_1, vsws_2
 
     REAL, DIMENSION(:,:), POINTER ::                                           &
-          qsws, qsws_m, rif, rif_m, shf, shf_m, usws, usws_m, vsws, vsws_m
+          qsws, qsws_m, qswst, qswst_m, rif, rif_m, shf, shf_m, tswst,         &
+          tswst_m, usws, usws_m, vsws, vsws_m
 
     REAL, DIMENSION(:,:,:), ALLOCATABLE ::                                     &
@@ -88 +91 @@
 
     REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                             &
-          e_1, e_2, e_3, kh_1, kh_2, km_1, km_2, pt_1, pt_2, pt_3, q_1, q_2,   &
-          q_3, ql_1, ql_2, u_1, u_2, u_3, v_1, v_2, v_3, vpt_1, vpt_2, w_1,    &
-          w_2, w_3
+          e_1, e_2, e_3, kh_1, kh_2, km_1, km_2, p, pt_1, pt_2, pt_3, q_1,     &
+          q_2, q_3, ql_1, ql_2, u_1, u_2, u_3, v_1, v_2, v_3, vpt_1, vpt_2,    &
+          w_1, w_2, w_3
 
     REAL, DIMENSION(:,:,:), POINTER ::                                         &
-          e, e_m, e_p, kh, kh_m, km, km_m, p, pt, pt_m, pt_p, q, q_m, q_p, ql, &
+          e, e_m, e_p, kh, kh_m, km, km_m, pt, pt_m, pt_p, q, q_m, q_p, ql,    &
           ql_c, te_m, tpt_m, tq_m, tu_m, tv_m, tw_m, u, u_m, u_p, v, v_m, v_p, &
           vpt, vpt_m, w, w_m, w_p
@@ -201 +204 @@
                             scalar_advec = 'pw-scheme'
     CHARACTER (LEN=20)  ::  bc_e_b = 'neumann', bc_lr = 'cyclic', &
-                            bc_ns = 'cyclic', &
-                            bc_p_b = 'neumann', bc_p_t = 'dirichlet', &
-                            bc_pt_b = 'dirichlet', bc_pt_t = 'neumann', &
+                            bc_ns = 'cyclic', bc_p_b = 'neumann', &
+                            bc_p_t = 'dirichlet', bc_pt_b = 'dirichlet', &
+                            bc_pt_t = 'initial_gradient', &
                             bc_q_b = 'dirichlet', bc_q_t = 'neumann', &
                             bc_s_b = 'dirichlet', bc_s_t = 'neumann', &
@@ -263 +266 @@
                 cloud_droplets = .FALSE., cloud_physics = .FALSE., &
                 conserve_volume_flow = .FALSE., constant_diffusion = .FALSE., &
-                constant_heatflux = .TRUE., constant_waterflux = .TRUE., &
-                create_disturbances = .TRUE., cut_spline_overshoot = .TRUE., &
+                constant_heatflux = .TRUE., constant_top_heatflux = .TRUE., &
+                constant_waterflux = .TRUE., create_disturbances = .TRUE., &
+                cut_spline_overshoot = .TRUE., &
                 data_output_2d_on_each_pe = .TRUE., do2d_at_begin = .FALSE., &
                 do3d_at_begin = .FALSE., do3d_compress = .FALSE., &
@@ -270 +274 @@
                 disturbance_created = .FALSE., &
                 first_call_advec_particles = .TRUE., &
-                force_print_header = .FALSE., galilei_transformation = .FALSE., &
+                force_print_header = .FALSE., galilei_transformation = .FALSE.,&
                 inflow_l = .FALSE., inflow_n = .FALSE., inflow_r = .FALSE., &
                 inflow_s = .FALSE., iso2d_output = .FALSE., &
@@ -281 +285 @@
                 random_heatflux = .FALSE., run_control_header = .FALSE., &
                 sloping_surface = .FALSE., stop_dt = .FALSE., &
-                terminate_run = .FALSE., use_prior_plot1d_parameters = .FALSE., &
-                use_surface_fluxes = .FALSE., use_ug_for_galilei_tr = .TRUE., &
+                terminate_run = .FALSE., use_prior_plot1d_parameters = .FALSE.,&
+                use_surface_fluxes = .FALSE., use_top_fluxes = .FALSE., &
+                use_ug_for_galilei_tr = .TRUE., &
                 use_upstream_for_tke = .FALSE., wall_adjustment = .TRUE.
 
@@ -336 +341 @@
              time_do_sla = 0.0, time_dvrp = 0.0, time_prel = 0.0, &
              time_restart = 9999999.9, time_run_control = 0.0, &
-             ug_surface = 0.0, u_gtrans = 0.0, ups_limit_e = 0.0, &
-             ups_limit_pt = 0.0, ups_limit_u = 0.0, ups_limit_v = 0.0, &
-             ups_limit_w = 0.0, vg_surface = 0.0, v_gtrans = 0.0, &
-             wall_adjustment_factor = 1.8, z_max_do1d = -1.0, &
+             top_heatflux = 9999999.9, ug_surface = 0.0, u_gtrans = 0.0, &
+             ups_limit_e = 0.0, ups_limit_pt = 0.0, ups_limit_u = 0.0, &
+             ups_limit_v = 0.0, ups_limit_w = 0.0, vg_surface = 0.0, &
+             v_gtrans = 0.0, wall_adjustment_factor = 1.8, z_max_do1d = -1.0, &
              z_max_do1d_normalized = -1.0, z_max_do2d = -1.0
 
@@ -500 +505 @@
     INTEGER ::  ngp_sums, nnx, nx = 0, nxa, nxl, nxr, nxra, nny, ny = 0, nya,  &
                 nyn, nyna, nys, nnz, nz = 0, nza, nzb, nzb_diff, nzt, nzta,    &
-                uxrp = 0, vynp = 0
+                nzt_diff, uxrp = 0, vynp = 0
 
     INTEGER, DIMENSION(:), ALLOCATABLE ::                                      &