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

Timestamp:

Jan 24, 2018 3:38:29 PM (6 years ago)

Author:

kanani

Message:

Added parameter check, reduced line length, some formatting

File:

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

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added check for output quantity pcm_heatrate, some formatting
+! 
+! 2766 2018-01-22 17:17:47Z kanani
 ! Increased LEN of canopy mode to 30
 ! 
@@ -194 +197 @@
     REAL(wp) ::  leaf_scalar_exch_coeff = 0.0_wp  !< canopy scalar exchange coeff.
     REAL(wp) ::  leaf_surface_conc = 0.0_wp       !< leaf surface concentration
+    REAL(wp) ::  lsc = 0.0_wp                     !< leaf surface concentration
     REAL(wp) ::  lsec = 0.0_wp                    !< leaf scalar exchange coeff.
-    REAL(wp) ::  lsc = 0.0_wp                     !< leaf surface concentration
     REAL(wp) ::  prototype_lad                    !< prototype leaf area density for computing effective optical depth
 
@@ -204 +207 @@
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  pre_lad        !< preliminary lad
     
-    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::                                 &
-       pc_heating_rate                                    !< plant canopy heating rate
-    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  cum_lai_hf !< cumulative lai for heatflux calc.
-    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  lad_s      !< lad on scalar-grid
-
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  cum_lai_hf       !< cumulative lai for heatflux calc.
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  lad_s            !< lad on scalar-grid
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  pc_heating_rate  !< plant canopy heating rate
 
     SAVE
@@ -287 +288 @@
 
        CASE ( 'pcm_heatrate' )
+          IF ( cthf == 0.0_wp )  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 )
+          ENDIF
           unit = 'K s-1'
     
@@ -323 +329 @@
           message_string = 'plant_canopy = .TRUE. requires a non-zero drag '// &
                            'coefficient & given value is canopy_drag_coeff = 0.0'
-          CALL message( 'check_parameters', 'PA0041', 1, 2, 0, 6, 0 )
+          CALL message( 'pcm_check_parameters', 'PA0041', 1, 2, 0, 6, 0 )
        ENDIF
     
@@ -331 +337 @@
                            'of the leaf area density profile requires '    //  &
                            'both alpha_lad and beta_lad to be /= 9999999.9'
-          CALL message( 'check_parameters', 'PA0118', 1, 2, 0, 6, 0 )
+          CALL message( 'pcm_check_parameters', 'PA0118', 1, 2, 0, 6, 0 )
        ENDIF
     
@@ -339 +345 @@
                            'a non-zero lai_beta, but given value is '      //  &
                            'lai_beta = 0.0'
-          CALL message( 'check_parameters', 'PA0119', 1, 2, 0, 6, 0 )
+          CALL message( 'pcm_check_parameters', 'PA0119', 1, 2, 0, 6, 0 )
        ENDIF
 
@@ -349 +355 @@
                            'function for the construction of the leaf area '// &
                            'density profile'
-          CALL message( 'check_parameters', 'PA0120', 1, 2, 0, 6, 0 )
+          CALL message( 'pcm_check_parameters', 'PA0120', 1, 2, 0, 6, 0 )
        ENDIF 
 
@@ -355 +361 @@
           message_string = 'plant_canopy = .TRUE. requires cloud_scheme /=' // &
                           ' seifert_beheng'
-          CALL message( 'check_parameters', 'PA0360', 1, 2, 0, 6, 0 )
+          CALL message( 'pcm_check_parameters', 'PA0360', 1, 2, 0, 6, 0 )
        ENDIF
 !
@@ -365 +371 @@
                            TRIM( coupling_char ) // ' requires ' //            &
                            'canopy_mode = read_from_file_3d'
-          CALL message( 'check_parameters', 'PA0999', 1, 2, 0, 6, 0 )
+          CALL message( 'pcm_check_parameters', 'PA0999', 1, 2, 0, 6, 0 )
        ENDIF
 
@@ -844 +850 @@
 !--    When using the urban surface model (USM), canopy heating (pc_heating_rate)
 !--    by radiation is calculated in the USM.
-       IF ( cthf /= 0.0_wp  .AND. .NOT.  urban_surface)  THEN
+       IF ( cthf /= 0.0_wp  .AND. .NOT.  urban_surface )  THEN
 
           ALLOCATE( cum_lai_hf(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                 &