Changeset 4514 for palm/trunk/SOURCE/plant_canopy_model_mod.f90

Timestamp:

Apr 30, 2020 4:29:59 PM (4 years ago)

Author:

suehring

Message:

Bugfix in plant-canopy model for output of averaged transpiration rate after a restart; Revise check for output for plant heating rate and rename error message number; Surface-data output: enable output of mixing ratio and passive scalar concentration at the surface; Surface-data input: Add possibility to prescribe surface sensible and latent heat fluxes via static input file

File:

• palm/trunk/SOURCE/plant_canopy_model_mod.f90 (modified) (10 diffs)

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! - Bugfix in output of pcm_heatrate_av in a restart run. In order to fix this, 
+!   pch_index is now output for a restart run. Therefore, define global restart
+!   routines.
+! - Error message number renamed and check for PA0505 revised in order to also
+!   consider natural surfaces with plant-canopy.
+! 
+! 4495 2020-04-13 20:11:20Z raasch
 ! restart data handling with MPI-IO added
 !
@@ -168 +175 @@
 ! @todo - precalculate constant terms in pcm_calc_transpiration_rate
 ! @todo - unify variable names (pcm_, pc_, ...)
+! @todo - get rid-off dependency on radiation model
 !------------------------------------------------------------------------------!
  MODULE plant_canopy_model_mod
@@ -187 +195 @@
               dz,                                                              &
               humidity,                                                        &
+              land_surface,                                                    &
               length,                                                          &
               message_string,                                                  &
@@ -225 +234 @@
 
     USE restart_data_mpi_io_mod,                                                                   &
-        ONLY:  wrd_mpi_io
+        ONLY:  rrd_mpi_io,                                                                         &
+               wrd_mpi_io
 
     USE surface_mod,                                                           &
@@ -294 +304 @@
            pcm_init,                                                          &
            pcm_parin,                                                         &
+           pcm_rrd_global,                                                    &
            pcm_rrd_local,                                                     &
            pcm_tendency,                                                      &
+           pcm_wrd_global,                                                    &
            pcm_wrd_local
 
@@ -349 +361 @@
     END INTERFACE pcm_rrd_local
 
+    INTERFACE pcm_rrd_global
+       MODULE PROCEDURE pcm_rrd_global_ftn
+       MODULE PROCEDURE pcm_rrd_global_mpi
+    END INTERFACE pcm_rrd_global
+
     INTERFACE pcm_tendency
        MODULE PROCEDURE pcm_tendency
@@ -357 +374 @@
        MODULE PROCEDURE pcm_wrd_local
     END INTERFACE pcm_wrd_local
+
+    INTERFACE pcm_wrd_global
+       MODULE PROCEDURE pcm_wrd_global
+    END INTERFACE pcm_wrd_global
 
 
@@ -483 +504 @@
 
        CASE ( 'pcm_heatrate' )
-          IF ( cthf == 0.0_wp  .AND. .NOT.  urban_surface )  THEN
-             message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+!
+!--       Output of heatrate can be only done if it is explicitely set by cthf,
+!--       or parametrized by absorption of radiation. The latter, however, is
+!--       only available if radiation_interactions are on. Note, these are 
+!--       enabled if land-surface or urban-surface is switched-on. Using 
+!--       radiation_interactions_on directly is not possible since it belongs
+!--       to the radition_model, which in turn depends on the plant-canopy model,
+!--       creating circular dependencies.
+          IF ( cthf == 0.0_wp  .AND. (  .NOT.  urban_surface  .AND.             &
+                                        .NOT.  land_surface ) )  THEN
+             message_string = 'output of "' // TRIM( var ) // '" requi' //      &
                               'res setting of parameter cthf /= 0.0'
-             CALL message( 'pcm_check_data_output', 'PA1000', 1, 2, 0, 6, 0 )
+             CALL message( 'pcm_check_data_output', 'PA0505', 1, 2, 0, 6, 0 )
           ENDIF
           unit = 'K s-1'
@@ -1620 +1650 @@
 ! Description:
 ! ------------
+!> Read module-specific global restart data (Fortran binary format).
+!------------------------------------------------------------------------------!
+    SUBROUTINE pcm_rrd_global_ftn( found )
+
+       LOGICAL, INTENT(OUT)  ::  found
+
+       found = .TRUE.
+
+       SELECT CASE ( restart_string(1:length) )
+
+          CASE ( 'pch_index' )
+             READ ( 13 )  pch_index
+
+          CASE DEFAULT
+
+             found = .FALSE.
+
+       END SELECT
+
+    END SUBROUTINE pcm_rrd_global_ftn
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Read module-specific global restart data (MPI-IO).
+!------------------------------------------------------------------------------!
+    SUBROUTINE pcm_rrd_global_mpi
+
+       CALL rrd_mpi_io( 'pch_index', pch_index )
+
+    END SUBROUTINE pcm_rrd_global_mpi
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
 !> Subroutine reads local (subdomain) restart data
 !------------------------------------------------------------------------------!
@@ -2437 +2502 @@
 ! Description:
 ! ------------
+!> Subroutine writes global restart data
+!------------------------------------------------------------------------------!
+    SUBROUTINE pcm_wrd_global
+
+       IF ( TRIM( restart_data_format_output ) == 'fortran_binary' )  THEN
+
+          CALL wrd_write_string( 'pch_index' )
+          WRITE ( 14 )  pch_index
+
+       ELSEIF ( TRIM( restart_data_format_output ) == 'mpi' )  THEN
+
+          CALL wrd_mpi_io( 'pch_index', pch_index )
+
+       ENDIF
+
+    END SUBROUTINE pcm_wrd_global
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
 !> Subroutine writes local (subdomain) restart data
 !------------------------------------------------------------------------------!