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

Timestamp:

Oct 30, 2013 11:36:58 AM (10 years ago)

Author:

heinze

Message:

Nudging and large scale forcing from external file implemented

File:

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

palm/trunk/SOURCE/modules.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! Usage of nudging enabled 
+! +nudging, ntnudge, ptnudge, qnudge, tnudge, unudge, vnudge, wnudge
+! increase pr_palm from 80 to 90 to allow for more standard profiles
+!
+! Enable prescribed time depenend surface fluxes and geostrophic wind read in
+! from external file LSF_DATA
+! +large_scale_forcing, lsf_surf, lsf_vert, nlsf, time_surf, shf_surf, qsws_surf,
+!  pt_surf, q_surf, p_surf, time_vert, ug_vert, vg_vert, wsubs_vert
 !
 ! Former revisions:
@@ -434 +441 @@
           c_u_m, c_u_m_l, c_v_m, c_v_m_l, c_w_m, c_w_m_l, ddzu, ddzu_pres,     &
           dd2zu, dzu, ddzw, dzw, hyp, inflow_damping_factor, lad, l_grid,      &
-          nc_1d, nr_1d, ptdf_x, ptdf_y, pt_1d, pt_init, q_1d, q_init, qc_1d,   &
-          qr_1d, rdf, rdf_sc, ref_state, sa_init, ug, u_init,                  &
-          u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, w_subs, zu, zw
+          nc_1d, nr_1d, ptdf_x, ptdf_y, p_surf, pt_surf, pt_1d, pt_init,       &
+          qsws_surf, q_1d, q_init, q_surf, qc_1d, qr_1d, rdf, rdf_sc,          & 
+          ref_state, sa_init, shf_surf, timenudge, time_surf, time_vert, ug,   &
+          u_init, u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, w_subs, zu, zw
 
     REAL, DIMENSION(:,:), ALLOCATABLE ::                                       &
@@ -443 +451 @@
           flux_s_e, flux_s_nr, flux_s_pt, flux_s_q, flux_s_qr, flux_s_sa,      &
           flux_s_u, flux_s_v, flux_s_w, f1_mg, f2_mg, f3_mg,                   &
-          mean_inflow_profiles, nrs, nrsws, nrswst, pt_slope_ref, qs, qsws,    &
-          qswst, qswst_remote, qrs, qrsws, qrswst, rif, saswsb, saswst, shf,   &
-          total_2d_a, total_2d_o, ts, tswst, us, usws, uswst, vsws, vswst, z0, &
-          z0h
+          mean_inflow_profiles, nrs, nrsws, nrswst, ptnudge, pt_slope_ref,     &
+          qnudge, qs, qsws, qswst, qswst_remote, qrs, qrsws, qrswst, rif,      &
+          saswsb, saswst, shf, tnudge, total_2d_a, total_2d_o, ts, tswst,      &
+          ug_vert, unudge, us, usws, uswst, vnudge, vg_vert, vsws, vswst,      &
+          wnudge, wsubs_vert, z0, z0h
 
     REAL, DIMENSION(:,:,:), ALLOCATABLE ::                                     &
@@ -478 +487 @@
     REAL, DIMENSION(:,:,:), ALLOCATABLE :: var_x, var_y, var_z, gamma_x,       &
                                            gamma_y, gamma_z
+
 
     SAVE
@@ -715 +725 @@
                 maximum_parallel_io_streams = -1, max_pr_user = 0, &
                 mgcycles = 0, mg_cycles = -1, mg_switch_to_pe0_level = 0, mid, &
-                netcdf_data_format = 2, ngsrb = 2, nr_timesteps_this_run = 0, &
+                nlsf = 1000, ntnudge = 100, netcdf_data_format = 2, ngsrb = 2, &
+                nr_timesteps_this_run = 0, &
                 nsor = 20, nsor_ini = 100, n_sor, normalizing_region = 0, &
                 nz_do3d = -9999, pch_index = 0, prt_time_count = 0, &
@@ -772 +783 @@
                 inflow_l = .FALSE., inflow_n = .FALSE., &
                 inflow_r = .FALSE., inflow_s = .FALSE., &
-                iso2d_output = .FALSE., large_scale_subsidence = .FALSE., &
+                iso2d_output = .FALSE., large_scale_forcing = .FALSE., &
+                large_scale_subsidence = .FALSE., lsf_surf = .TRUE., &
+                lsf_vert = .TRUE., lptnudge = .FALSE., lqnudge = .FALSE., &
+                lunudge = .FALSE., lvnudge = .FALSE., lwnudge = .FALSE., &
                 masking_method = .FALSE., mg_switch_to_pe0 = .FALSE., &
-                netcdf_output = .FALSE., neutral = .FALSE., ocean = .FALSE., &
-                on_device = .FALSE., &
+                netcdf_output = .FALSE., neutral = .FALSE., nudging = .FALSE., &
+                ocean = .FALSE., on_device = .FALSE., &
                 outflow_l = .FALSE., outflow_n = .FALSE., outflow_r = .FALSE., &
                 outflow_s = .FALSE., passive_scalar = .FALSE., &
@@ -1661 +1675 @@
 
     CHARACTER (LEN=40) ::  region(0:9)
-    INTEGER ::  pr_palm = 80, statistic_regions = 0
+    INTEGER ::  pr_palm = 90, statistic_regions = 0
     INTEGER ::  u_max_ijk(3), v_max_ijk(3), w_max_ijk(3)
     LOGICAL ::  flow_statistics_called = .FALSE.