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

Timestamp:

Jul 27, 2018 1:36:03 PM (6 years ago)

Author:

suehring

Message:

New Inifor features: grid stretching, improved command-interface, support start dates in different formats in both YYYYMMDD and YYYYMMDDHH, Ability to manually control input file prefixes (--radiation-prefix, --soil-preifx, --flow-prefix, --soilmoisture-prefix) for compatiblity with DWD forcast naming scheme; GNU-style short and long option; Prepared output of large-scale forcing profiles (no computation yet); Added preprocessor flag netcdf4 to switch output format between netCDF 3 and 4; Updated netCDF variable names and attributes to comply with PIDS v1.9; Inifor bugfixes: Improved compatibility with older Intel Intel compilers by avoiding implicit array allocation; Added origin_lon/_lat values and correct reference time in dynamic driver global attributes; corresponding PALM changes: adjustments to revised Inifor; variables names in dynamic driver adjusted; enable geostrophic forcing also in offline nested mode; variable names in LES-LES and COSMO offline nesting changed; lateral boundary flags for nesting, in- and outflow conditions renamed

File:

• palm/trunk/SOURCE/parin.f90 (modified) (12 diffs)

palm/trunk/SOURCE/parin.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables and boundary conditions in mesoscale-offline nesting mode
 ! 
 ! Former revisions:
@@ -529 +529 @@
              dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
              dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
-             e_min, fft_method, flux_input_mode, flux_output_mode, forcing,    &
+             e_min, fft_method, flux_input_mode, flux_output_mode,             &
              galilei_transformation, humidity,                                 &
              inflow_damping_height, inflow_damping_width,                      &
@@ -538 +538 @@
              loop_optimization, lsf_exception, masking_method, mg_cycles,      &
              mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
-             most_method, na_init, nc_const, netcdf_precision, neutral, ngsrb, &
+             most_method, na_init, nc_const, nesting_offline,                  &
+             netcdf_precision, neutral, ngsrb,                                 &
              nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
              outflow_source_plane, passive_scalar,                             &
@@ -601 +602 @@
              dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
              dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
-             e_min, fft_method, flux_input_mode, flux_output_mode, forcing,    &
+             e_min, fft_method, flux_input_mode, flux_output_mode,             &
              galilei_transformation, humidity,                                 &
              inflow_damping_height, inflow_damping_width,                      &
@@ -610 +611 @@
              loop_optimization, lsf_exception, masking_method, mg_cycles,      &
              mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
-             most_method, na_init, nc_const, netcdf_precision, neutral, ngsrb, &
+             most_method, na_init, nc_const, nesting_offline,                  &
+             netcdf_precision, neutral, ngsrb,                                 &
              nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
              outflow_source_plane, passive_scalar,                             &
@@ -924 +926 @@
                    bc_ns   = 'cyclic'
                 ENDIF
-                IF ( nest_domain )  THEN
+                IF ( child_domain )  THEN
                    bc_uv_t  = 'nested'
                    bc_pt_t  = 'nested'
@@ -936 +938 @@
 !--          nested domains.
              ELSE 
-                IF ( nest_domain )  THEN
+                IF ( child_domain )  THEN
                    bc_lr    = 'nested'
                    bc_ns    = 'nested'
@@ -948 +950 @@
              ENDIF
           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_cs_t  = 'forcing'  ! same for chemical species
+!
+!--       Set boundary conditions also in case the model is offline-nested in
+!--       larger-scale models.
+          IF ( nesting_offline )  THEN
+             bc_lr    = 'nesting_offline'
+             bc_ns    = 'nesting_offline'
+             bc_uv_t  = 'nesting_offline'
+             bc_pt_t  = 'nesting_offline'
+             bc_q_t   = 'nesting_offline'
+             bc_s_t   = 'nesting_offline'  ! scalar boundary condition is not clear
+             bc_cs_t  = 'nesting_offline'  ! same for chemical species
              bc_p_t   = 'neumann'
           ENDIF
@@ -968 +972 @@
 !--       are set properly. An exception is made in case of restart runs and
 !--       if user decides to do everything by its own.
-          IF ( nest_domain  .AND.  .NOT. (                                     &
+          IF ( child_domain  .AND.  .NOT. (                                    &
                TRIM( initializing_actions ) == 'read_restart_data'      .OR.   &
                TRIM( initializing_actions ) == 'set_constant_profiles'  .OR.   &
@@ -979 +983 @@
 
              initializing_actions = 'set_constant_profiles'
-          ENDIF
-            
+          ENDIF            
 !
 !--       Check validity of lateral boundary conditions. This has to be done
@@ -987 +990 @@
           IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
                bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested'  .AND. &
-               bc_lr /= 'forcing' )  THEN
+               bc_lr /= 'nesting_offline' )  THEN
              message_string = 'unknown boundary condition: bc_lr = "' // &
                               TRIM( bc_lr ) // '"'
@@ -994 +997 @@
           IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
                bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested'  .AND. &
-               bc_ns /= 'forcing' )  THEN
+               bc_ns /= 'nesting_offline' )  THEN
              message_string = 'unknown boundary condition: bc_ns = "' // &
                               TRIM( bc_ns ) // '"'