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

Timestamp:

Apr 17, 2018 10:27:57 AM (6 years ago)

Author:

kanani

Message:

Fixes for radiative transfer model

File:

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

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -19 +19 @@
 !
 ! Current revisions:
-! -----------------
+! ------------------
 ! 
 ! 
@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implement changes from branch radiation (r2948-2971) with minor modifications,
+! plus some formatting.
+! (moh.hefny):
+! Add plant canopy type to account for changes in LAD (based on the changes
+! done by Resler & Pavel) and correct the error message to PALM Standard.
+! 
+! 2932 2018-03-26 09:39:22Z maronga
 ! renamed canopy_par to plant_canopy_parameters
 ! 
@@ -213 +220 @@
     REAL(wp) ::  lad_vertical_gradient(10) = 0.0_wp              !< lad gradient
     REAL(wp) ::  lad_vertical_gradient_level(10) = -9999999.9_wp !< lad-prof. levels (in m)
+
+    REAL(wp) ::  lad_type_coef(0:10) = 1.0_wp   !< multiplicative coeficients for particular types
+                                                !< of plant canopy (e.g. deciduous tree during winter)
 
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  lad            !< leaf area density
@@ -962 +972 @@
                                   alpha_lad, beta_lad, canopy_drag_coeff,      &
                                   canopy_mode, cthf,                           &
-                                  lad_surface,                                 & 
+                                  lad_surface, lad_type_coef,                  & 
                                   lad_vertical_gradient,                       &
                                   lad_vertical_gradient_level,                 &
@@ -971 +981 @@
        NAMELIST /canopy_par/      alpha_lad, beta_lad, canopy_drag_coeff,      &
                                   canopy_mode, cthf,                           &
-                                  lad_surface,                                 & 
+                                  lad_surface, lad_type_coef,                  & 
                                   lad_vertical_gradient,                       &
                                   lad_vertical_gradient_level,                 &
@@ -1034 +1044 @@
 !>
 !> num_levels
-!> dtype,x,y,value(nzb),value(nzb+1), ... ,value(nzb+num_levels-1)
-!> dtype,x,y,value(nzb),value(nzb+1), ... ,value(nzb+num_levels-1)
-!> dtype,x,y,value(nzb),value(nzb+1), ... ,value(nzb+num_levels-1)
+!> dtype,x,y,pctype,value(nzb),value(nzb+1), ... ,value(nzb+num_levels-1)
+!> dtype,x,y,pctype,value(nzb),value(nzb+1), ... ,value(nzb+num_levels-1)
+!> dtype,x,y,pctype,value(nzb),value(nzb+1), ... ,value(nzb+num_levels-1)
 !> ...
 !>
@@ -1060 +1070 @@
 
        INTEGER(iwp)                        ::  dtype     !< type of input data (1=lad)
+       INTEGER(iwp)                        ::  pctype    !< type of plant canopy (deciduous,non-deciduous,...)
        INTEGER(iwp)                        ::  i, j      !< running index
        INTEGER(iwp)                        ::  nzp       !< number of vertical layers of plant canopy
@@ -1073 +1084 @@
 !
 !--    Open and read plant canopy input data
-       OPEN(152, file='PLANT_CANOPY_DATA_3D' // TRIM( coupling_char ),         &
-                 access='SEQUENTIAL', action='READ', status='OLD',             &
-                 form='FORMATTED', err=515)
-       READ(152, *, err=516, end=517) nzp   !< read first line = number of vertical layers
+       OPEN(152, FILE='PLANT_CANOPY_DATA_3D' // TRIM( coupling_char ),         &
+                 ACCESS='SEQUENTIAL', ACTION='READ', STATUS='OLD',             &
+                 FORM='FORMATTED', ERR=515)
+       READ(152, *, ERR=516, END=517)  nzp   !< read first line = number of vertical layers
        
-       ALLOCATE(col(0:nzp-1))
+       ALLOCATE( col(0:nzp-1) )
 
        DO
-          READ(152, *, err=516, end=517) dtype, i, j, col(:)
-          IF ( i < nxlg .or. i > nxrg .or. j < nysg .or. j > nyng ) CYCLE
-
-             SELECT CASE (dtype)
-                CASE( 1 )   !< leaf area density
-!
-!--                 This is just the pure canopy layer assumed to be grounded to
-!--                 a flat domain surface. At locations where plant canopy sits 
-!--                 on top of any kind of topography, the vertical plant column 
-!--                 must be "lifted", which is done in SUBROUTINE pcm_tendency.
-                    lad_s(0:nzp-1, j, i) = col(0:nzp-1)
-                    
-                CASE DEFAULT
-                   WRITE(message_string, '(a,i2,a)')   &
-                      'Unknown record type in file PLANT_CANOPY_DATA_3D: "', dtype, '"'
-                   CALL message( 'pcm_read_plant_canopy_3d', 'PA0530', 1, 2, 0, 6, 0 )
-             END SELECT
+          READ(152, *, ERR=516, END=517) dtype, i, j, pctype, col(:)
+          IF ( i < nxlg  .OR.  i > nxrg  .OR.  j < nysg  .OR.  j > nyng )  CYCLE
+          
+          SELECT CASE (dtype)
+             CASE( 1 )   !< leaf area density
+!
+!--             This is just the pure canopy layer assumed to be grounded to
+!--             a flat domain surface. At locations where plant canopy sits 
+!--             on top of any kind of topography, the vertical plant column 
+!--             must be "lifted", which is done in SUBROUTINE pcm_tendency.           
+                IF ( pctype < 0  .OR.  pctype > 10 )  THEN   !< incorrect plant canopy type
+                   WRITE( message_string, * ) 'Incorrect type of plant canopy. '   //  &
+                                              'Allowed values 0 <= pctype <= 10, ' //  &
+                                              'but pctype is ', pctype
+                   CALL message( 'pcm_read_plant_canopy_3d', 'PA0349', 1, 2, 0, 6, 0 )
+                ENDIF
+                lad_s(0:nzp-1,j,i) = col(0:nzp-1) * lad_type_coef(pctype)
+                
+             CASE DEFAULT
+                WRITE(message_string, '(a,i2,a)')   &
+                     'Unknown record type in file PLANT_CANOPY_DATA_3D: "', dtype, '"'
+                CALL message( 'pcm_read_plant_canopy_3d', 'PA0530', 1, 2, 0, 6, 0 )
+          END SELECT
        ENDDO
 
@@ -1107 +1124 @@
 
 517    CLOSE(152)
-       DEALLOCATE(col)
+       DEALLOCATE( col )
        
        CALL exchange_horiz( lad_s, nbgp )