Changeset 2696 for palm/trunk/SOURCE/parin.f90

Timestamp:

Dec 14, 2017 5:12:51 PM (7 years ago)

Author:

kanani

Message:

Merge of branch palm4u into trunk

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/parin.f90 (modified) (13 diffs)

palm/trunk/SOURCE/parin.f90

@@ -1 +1 @@
 !> @file parin.f90
 !------------------------------------------------------------------------------!
-! This file is part of PALM.
+! This file is part of the PALM model system.
 !
 ! PALM is free software: you can redistribute it and/or modify it under the
@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of uv exposure model (FK)
+! Added rans_mode and turbulence_closure to inipar (TG)
+! Implementation of chemistry module
+! Sorting of USE list (FK)
+! Forcing implemented, and initialization with inifor (MS)
+! 
+! 2600 2017-11-01 14:11:20Z raasch
 ! some comments added and variables renamed concerning r2599
 ! 
@@ -316 +323 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  pt_init, q_init, ref_state, s_init, sa_init, ug, u_init, v_init,&
-               vg
+        ONLY:  pt_init, q_init, ref_state, s_init, sa_init,                    &     
+               ug, u_init, v_init, vg
+
+#if defined( __chem )
+    USE chemistry_model_mod,                                                   &
+        ONLY:  chem_parin
+        
+    USE chem_modules
+#endif
+
+    USE control_parameters
+
+    USE cpulog,                                                                &
+        ONLY:  cpu_log_barrierwait
 
     USE date_and_time_mod,                                                     &
         ONLY:  day_of_year_init, time_utc_init
-               
-    USE plant_canopy_model_mod,                                                &
-         ONLY: pcm_parin  
-
-    USE control_parameters
-
-    USE cpulog,                                                                &
-        ONLY:  cpu_log_barrierwait
 
     USE dvrp_variables,                                                        &
@@ -353 +364 @@
                sigma_bulk, ventilation_effect
 
-
     USE model_1d_mod,                                                          &
         ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
 
-    USE pegrid
-
     USE netcdf_interface,                                                      &
         ONLY:  netcdf_data_format, netcdf_deflate, netcdf_precision
 
+    USE pegrid
+               
+    USE plant_canopy_model_mod,                                                &
+         ONLY: pcm_parin
+
     USE pmc_interface,                                                         &
         ONLY:  nested_run, nesting_mode
@@ -386 +399 @@
         ONLY: usm_parin
 
+    USE uv_exposure_model_mod,                                                 &
+        ONLY:  uvem_parin
+
+    USE vertical_nesting_mod,                                                  &
+        ONLY:  vnest_start_time
+
     USE wind_turbine_model_mod,                                                &
         ONLY:  wtm_parin
 
-    USE vertical_nesting_mod,                                                  &
-        ONLY:  vnest_start_time
 
     IMPLICIT NONE
@@ -422 +439 @@
              dz_stretch_factor, dz_stretch_level, end_time_1d,                 &
              ensemble_member_nr, e_init, e_min, fft_method,                    &
-             flux_input_mode, flux_output_mode,                                &
+             flux_input_mode, flux_output_mode, forcing,                       &
              galilei_transformation, humidity,                                 &
              inflow_damping_height, inflow_damping_width,                      &
@@ -439 +456 @@
              pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
              q_vertical_gradient, q_vertical_gradient_level,                   &
-             random_generator, random_heatflux,                                &
+             random_generator, random_heatflux, rans_mode,                     &
              rayleigh_damping_factor, rayleigh_damping_height,                 &
              recycling_width, recycling_yshift,                                &
@@ -456 +473 @@
              top_scalarflux, transpose_compute_overlap,                        &
              tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
-             tunnel_wall_depth,                                                &
+             tunnel_wall_depth, turbulence_closure,                            &
              turbulent_inflow, turbulent_outflow,                              &
              use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
@@ -617 +634 @@
           ENDIF
 
+          IF ( forcing )  THEN
+             bc_lr   = 'forcing'
+             bc_ns   = 'forcing'
+             bc_uv_t = 'forcing'
+             bc_pt_t = 'forcing'
+             bc_q_t  = 'forcing'
+             bc_s_t  = 'forcing'  ! scalar boundary condition is not clear
+             bc_p_t  = 'neumann'
+
+          ENDIF
+
 !         
 !--       In case of nested runs, make sure that initializing_actions =
@@ -625 +653 @@
 !--       are set properly. An exception is made in case of restart runs and
 !--       if user decides to do everything by its own.
-!--       MS: is this setting really necessary?
-          IF ( nest_domain  .AND.  (                                           &
-               TRIM( initializing_actions ) /= 'read_restart_data'  .OR.       &
-               TRIM( initializing_actions ) /= 'by_user' ) )  THEN
+!--       MS: is this really necessary?
+          IF ( nest_domain  .AND.  .NOT. (                                     &
+               TRIM( initializing_actions ) == 'read_restart_data'  .OR.       &
+               TRIM( initializing_actions ) == 'by_user'            .OR.       &
+               TRIM( initializing_actions ) == 'inifor' ) )  THEN
              initializing_actions = 'set_constant_profiles'
           ENDIF
@@ -637 +666 @@
 !--       therefore cannot be check in check_parameters
           IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
-               bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested' )  THEN
+               bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested'  .AND. &
+               bc_lr /= 'forcing' )  THEN
              message_string = 'unknown boundary condition: bc_lr = "' // &
                               TRIM( bc_lr ) // '"'
@@ -643 +673 @@
           ENDIF
           IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
-               bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested' )  THEN
+               bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested'  .AND. &
+               bc_ns /= 'forcing' )  THEN
              message_string = 'unknown boundary condition: bc_ns = "' // &
                               TRIM( bc_ns ) // '"'
              CALL message( 'check_parameters', 'PA0050', 1, 2, 0, 6, 0 )
           ENDIF
-
 !
 !--       Set internal variables used for speed optimization in if clauses
@@ -707 +737 @@
 !--       required
           CALL stg_parin
+!
+!--       Read chemistry variables
+#if defined( __chem )
+          CALL chem_parin
+#endif
+!
+!--       Check if uv exposure model is used and read &uvexposure_par
+          CALL uvem_parin
 !
 !--       Read user-defined variables