Changeset 4109 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Jul 22, 2019 5:00:34 PM (5 years ago)

Author:

suehring

Message:

Control discretization of advection term: separate initialization of WS advection flags for momentum and scalars. In this context, resort the bits and do some minor formatting; Make initialization of scalar-advection flags more flexible, i.e. introduce an arguemnt list to indicate non-cyclic boundaries (required for decycled scalars such as chemical species or aerosols); Introduce extended 'degradation zones', where horizontal advection of passive scalars is discretized by first-order scheme at all grid points that in the vicinity of buildings (<= 3 grid points). Even though no building is within the numerical stencil, first-order scheme is used. At fourth and fifth grid point the order of the horizontal advection scheme is successively upgraded. These extended degradation zones are used to avoid stationary numerical oscillations, which are responsible for high concentration maxima that may appear under shear-free stable conditions. Therefore, an additional 3D interger array used to store control flags is introduced; Change interface for scalar advection routine; Bugfix, avoid uninitialized value sk_num in vector version of WS scalar advection; Chemistry: Decycling boundary conditions are only set at the ghost points not on the prognostic grid points; Land-surface model: Relax checks for non-consistent initialization in case static or dynamic input is provided. For example, soil_temperature or deep_soil_temperature is not mandatory any more if dynamic input is available. Also, improper settings of x_type in namelist are only checked if no static file is available.

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (14 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Relax checks for non-consistent initialization in case static or dynamic 
+! input is provided. For example, soil_temperature or deep_soil_temperature
+! is not mandatory any more if dynamic input is available. Also, improper 
+! settings of x_type in namelist are only checked if no static file is 
+! available. 
 ! 
 ! Former revisions:
@@ -591 +595 @@
 
     USE netcdf_data_input_mod,                                                 &
-        ONLY :  building_type_f, init_3d, input_pids_static,                   &
+        ONLY :  building_type_f, init_3d,                                      &
+                input_pids_dynamic,                                            &
+                input_pids_static,                                             &
                 netcdf_data_input_interpolate, netcdf_data_input_init_lsm,     &
                 pavement_pars_f, pavement_subsurface_pars_f, pavement_type_f,  &
@@ -727 +733 @@
                                    soil_temperature = 9999999.9_wp,  & !< NAMELIST soil temperature (K) +1
                                    dz_soil  = 9999999.9_wp,          & !< (NAMELIST) soil layer depths (spacing)
-                                   zs_layer = 9999999.9_wp         !< soil layer depths (edge)
+                                   zs_layer = 9999999.9_wp             !< soil layer depths (edge)
                                  
     TYPE(surf_type_lsm), POINTER ::  t_soil_h,    & !< Soil temperature (K), horizontal surface elements
@@ -1437 +1443 @@
 !-- Check if soil types are set within a valid range.
     IF ( TRIM( surface_type ) == 'vegetation'  .OR.                            &
-         TRIM( surface_type ) == 'pavement'    .OR.                            &
-         TRIM( surface_type ) == 'netcdf' )  THEN
+         TRIM( surface_type ) == 'pavement' )  THEN
        IF ( soil_type < LBOUND( soil_pars, 2 )  .AND.                          &
             soil_type > UBOUND( soil_pars, 2 ) )  THEN
@@ -1445 +1450 @@
           CALL message( 'lsm_check_parameters', 'PA0452', 2, 2, 0, 6, 0 )
        ENDIF
+    ENDIF
+    IF ( TRIM( surface_type ) == 'netcdf' )  THEN
        IF ( soil_type_f%from_file )  THEN
           DO  i = nxl, nxr
@@ -1462 +1469 @@
 !
 !-- Check if vegetation types are set within a valid range.    
-    IF ( TRIM( surface_type ) == 'vegetation'  .OR.                            &
-         TRIM( surface_type ) == 'netcdf' )  THEN
+    IF ( TRIM( surface_type ) == 'vegetation' )  THEN
        IF ( vegetation_type < LBOUND( vegetation_pars, 2 )  .AND.              &
             vegetation_type > UBOUND( vegetation_pars, 2 ) )  THEN
@@ -1470 +1476 @@
           CALL message( 'lsm_check_parameters', 'PA0526', 2, 2, 0, 6, 0 )
        ENDIF
+    ENDIF
+    IF ( TRIM( surface_type ) == 'netcdf' )  THEN
        IF ( vegetation_type_f%from_file )  THEN
           DO  i = nxl, nxr
@@ -1488 +1496 @@
 !
 !-- Check if pavement types are set within a valid range.    
-    IF ( TRIM( surface_type ) == 'pavement'  .OR.                              &
-         TRIM( surface_type ) == 'netcdf' )  THEN
+    IF ( TRIM( surface_type ) == 'pavement' )  THEN
        IF ( pavement_type < LBOUND( pavement_pars, 2 )  .AND.                  &
             pavement_type > UBOUND( pavement_pars, 2 ) )  THEN
@@ -1496 +1503 @@
           CALL message( 'lsm_check_parameters', 'PA0527', 2, 2, 0, 6, 0 )
        ENDIF
+    ENDIF
+    IF ( TRIM( surface_type ) == 'netcdf' )  THEN
        IF ( pavement_type_f%from_file )  THEN
           DO  i = nxl, nxr
@@ -1513 +1522 @@
 !
 !-- Check if water types are set within a valid range.    
-    IF ( TRIM( surface_type ) == 'water'  .OR.                                 &
-         TRIM( surface_type ) == 'netcdf' )  THEN
+    IF ( TRIM( surface_type ) == 'water' )  THEN
        IF ( water_type < LBOUND( water_pars, 2 )  .AND.                        &
             water_type > UBOUND( water_pars, 2 ) )  THEN
@@ -1521 +1529 @@
           CALL message( 'lsm_check_parameters', 'PA0528', 2, 2, 0, 6, 0 )
        ENDIF
+    ENDIF
+    IF ( TRIM( surface_type ) == 'netcdf' )  THEN
        IF ( water_type_f%from_file )  THEN
           DO  i = nxl, nxr
@@ -1708 +1718 @@
 !
 !-- Temporary message as long as NetCDF input is not available
-    IF ( TRIM( surface_type ) == 'netcdf'  .AND.  .NOT.  input_pids_static )   &
+    IF ( TRIM( surface_type ) == 'netcdf'  .AND.  .NOT. input_pids_static )   &
     THEN
        message_string = 'surface_type = netcdf requires static input file.'
@@ -1714 +1724 @@
     ENDIF
 
-    IF ( soil_type == 0 )  THEN
+    IF ( soil_type == 0  .AND.  .NOT. input_pids_static )  THEN
 
        IF ( alpha_vangenuchten == 9999999.9_wp )  THEN
@@ -1796 +1806 @@
 !
 !-- Check whether valid soil temperatures are prescribed
-    IF ( COUNT( soil_temperature /= 9999999.9_wp ) /= nzs )  THEN
-       WRITE( message_string, * ) 'number of soil layers (', nzs, ') does not',&
+    IF ( .NOT. input_pids_dynamic )  THEN
+       IF ( COUNT( soil_temperature /= 9999999.9_wp ) /= nzs )  THEN
+          WRITE( message_string, * )                                           &
+                                  'number of soil layers (', nzs, ') does not',&
                                   ' match to the number of layers specified',  &
                                   ' in soil_temperature (', COUNT(             &
-                                   soil_temperature /= 9999999.9_wp ), ')'
-          CALL message( 'lsm_check_parameters', 'PA0471', 1, 2, 0, 6, 0 )
-    ENDIF
-
-    IF ( deep_soil_temperature == 9999999.9_wp ) THEN
-          message_string = 'deep_soil_temperature is not set but must be'//    &
-                           '/= 9999999.9'
-          CALL message( 'lsm_check_parameters', 'PA0472', 1, 2, 0, 6, 0 )
+                                  soil_temperature /= 9999999.9_wp ), ')'
+             CALL message( 'lsm_check_parameters', 'PA0471', 1, 2, 0, 6, 0 )
+       ENDIF
+
+       IF ( deep_soil_temperature == 9999999.9_wp ) THEN
+             message_string = 'deep_soil_temperature is not set but must be'// &
+                              '/= 9999999.9'
+             CALL message( 'lsm_check_parameters', 'PA0472', 1, 2, 0, 6, 0 )
+       ENDIF
     ENDIF