Changeset 4127 for palm/trunk/SOURCE/biometeorology_mod.f90

Timestamp:

Jul 30, 2019 2:47:10 PM (5 years ago)

Author:

suehring

Message:

Merge with branch resler: biomet- output of bio_mrt added; plant_canopy - separate vertical dimension for 3D output (to save disk space); radiation - remove unused plant canopy variables; urban-surface model - do not add anthropogenic heat during wall spin-up

File:

• palm/trunk/SOURCE/biometeorology_mod.f90 (modified) (16 diffs)

palm/trunk/SOURCE/biometeorology_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Output for bio_mrt added (merge from branch resler)
+! 
+! 4126 2019-07-30 11:09:11Z gronemeier
 ! renamed vitd3_exposure_av into vitd3_dose,
 ! renamed uvem_calc_exposure into bio_calculate_uv_exposure
@@ -226 +229 @@
     LOGICAL ::  average_trigger_utci   = .FALSE.  !< update averaged input on call to bio_utci?
     LOGICAL ::  average_trigger_pet    = .FALSE.  !< update averaged input on call to bio_pet?
+    LOGICAL ::  average_trigger_mrt    = .FALSE.  !< update averaged input on call to bio_pet?
     LOGICAL ::  do_calculate_perct     = .FALSE.  !< Turn index PT (instant. input) on or off
     LOGICAL ::  do_calculate_perct_av  = .FALSE.  !< Turn index PT (averaged input) on or off
@@ -232 +236 @@
     LOGICAL ::  do_calculate_utci      = .FALSE.  !< Turn index UTCI (instant. input) on or off
     LOGICAL ::  do_calculate_utci_av   = .FALSE.  !< Turn index UTCI (averaged input) on or off
+    LOGICAL ::  do_calculate_mrt2d     = .FALSE.  !< Turn index MRT 2D (averaged or inst) on or off
 
 !
@@ -432 +437 @@
 
 
-          CASE ( 'bio_perct*', 'bio_utci*', 'bio_pet*' )
+          CASE ( 'bio_perct*', 'bio_utci*', 'bio_pet*', 'bio_mrt*' )
 
 !
@@ -447 +452 @@
              IF ( .NOT. average_trigger_perct  .AND.                           &
                   .NOT. average_trigger_utci   .AND.                           &
-                  .NOT. average_trigger_pet )  THEN
+                  .NOT. average_trigger_pet    .AND.                           &
+                  .NOT. average_trigger_mrt )  THEN
 !
 !--             Memorize the first index called to control averaging
@@ -458 +464 @@
                 IF ( TRIM( variable ) == 'bio_pet*' )  THEN
                     average_trigger_pet = .TRUE.
+                ENDIF
+                IF ( TRIM( variable ) == 'bio_mrt*' )  THEN
+                    average_trigger_mrt = .TRUE.
                 ENDIF
              ENDIF
@@ -499 +508 @@
              ENDIF
 
-          CASE ( 'bio_perct*', 'bio_utci*', 'bio_pet*' )
+          CASE ( 'bio_perct*', 'bio_utci*', 'bio_pet*', 'bio_mrt*' )
 !
 !--          Only continue if the current index is the one to trigger the input 
@@ -509 +518 @@
              IF ( average_trigger_pet    .AND.  TRIM( variable ) /=            &
                 'bio_pet*')    RETURN
+             IF ( average_trigger_mrt    .AND.  TRIM( variable ) /=            &
+                'bio_mrt*')    RETURN
 !
 !--          Now memorize which of the input grids are not averaged by other 
@@ -637 +648 @@
              ENDIF
 
-          CASE ( 'bio_perct*', 'bio_utci*', 'bio_pet*' )
+          CASE ( 'bio_perct*', 'bio_utci*', 'bio_pet*', 'bio_mrt*' )
 !
 !--          Only continue if update index, see above
@@ -646 +657 @@
              IF ( average_trigger_pet   .AND.                                  &
                 TRIM( variable ) /= 'bio_pet*' )  RETURN
+             IF ( average_trigger_mrt   .AND.                                  &
+                TRIM( variable ) /= 'bio_mrt*' )  RETURN
 
              IF ( ALLOCATED( pt_av )  .AND.  do_average_theta )  THEN
@@ -748 +761 @@
 !--    derived in a single step based on priorly averaged arrays (see
 !--    bio_calculate_thermal_index_maps).
-       CASE ( 'bio_mrt' )
+       CASE ( 'bio_mrt', 'bio_mrt*' )
           unit = 'degree_C'
           thermal_comfort = .TRUE.  !< enable thermal_comfort if user forgot to do so
@@ -755 +768 @@
              tmrt_grid = REAL( bio_fill_value, KIND = wp )
           ENDIF
+          IF ( TRIM( var ) == 'bio_mrt*' )  THEN
+             do_calculate_mrt2d = .TRUE.
+          END IF
 
        CASE ( 'bio_perct*' )
@@ -961 +977 @@
            ENDDO
 
+        CASE ( 'bio_mrt*_xy' )        ! 2d-array
+           grid = 'zu1'
+           two_d = .TRUE.
+           IF ( av == 0 )  THEN
+              DO  i = nxl, nxr
+                 DO  j = nys, nyn
+                    local_pf(i,j,nzb+1) = tmrt_grid(j,i)
+                 ENDDO
+              ENDDO
+           ELSE
+              DO  i = nxl, nxr
+                 DO  j = nys, nyn
+                    local_pf(i,j,nzb+1) = tmrt_av_grid(j,i)
+                 ENDDO
+              ENDDO
+           ENDIF
+
 
         CASE ( 'bio_perct*_xy' )        ! 2d-array
@@ -1372 +1405 @@
           READ ( 13 )  average_trigger_pet
 
+       CASE ( 'average_trigger_mrt' )
+          READ ( 13 )  average_trigger_mrt
+
 
        CASE DEFAULT
@@ -1451 +1487 @@
     CALL wrd_write_string( 'average_trigger_pet' )
     WRITE ( 14 )  average_trigger_pet
+    CALL wrd_write_string( 'average_trigger_mrt' )
+    WRITE ( 14 )  average_trigger_mrt
 
  END SUBROUTINE bio_wrd_global
@@ -1694 +1732 @@
     IF ( do_calculate_perct  .OR.  do_calculate_perct_av  .OR.                 &
        do_calculate_utci  .OR.  do_calculate_utci_av  .OR.                     &
-       do_calculate_pet  .OR.  do_calculate_pet_av )  THEN
+       do_calculate_pet  .OR.  do_calculate_pet_av  .OR.                       &
+       do_calculate_mrt2d )  THEN
 
 !