Changeset 1484 for palm/trunk/SOURCE/user_init_plant_canopy.f90

Timestamp:

Oct 21, 2014 10:53:05 AM (10 years ago)

Author:

kanani

Message:

New:
---
Subroutine init_plant_canopy added to module plant_canopy_model_mod. (plant_canopy_model)
Alternative method for lad-profile construction added, also, new parameters added.
(header, package_parin, plant_canopy_model, read_var_list, write_var_list)
plant_canopy_model-dependency added to several subroutines. (Makefile)
New package/namelist canopy_par for canopy-related parameters added. (package_parin)

Changed:
---
Code structure of the plant canopy model changed, all canopy-model related code
combined to module plant_canopy_model_mod. (check_parameters, init_3d_model,
modules, timestep)
Module plant_canopy_model_mod added in USE-lists of some subroutines. (check_parameters,
header, init_3d_model, package_parin, read_var_list, user_init_plant_canopy, write_var_list)
Canopy initialization moved to new subroutine init_plant_canopy. (check_parameters,
init_3d_model, plant_canopy_model)
Calculation of canopy timestep-criterion removed, instead, the canopy
drag is now directly limited in the calculation of the canopy tendency terms.
(plant_canopy_model, timestep)
Some parameters renamed. (check_parameters, header, init_plant_canopy,
plant_canopy_model, read_var_list, write_var_list)
Unnecessary 3d-arrays removed. (init_plant_canopy, plant_canopy_model, user_init_plant_canopy)
Parameter checks regarding canopy initialization added. (check_parameters)
All canopy steering parameters moved from namelist inipar to canopy_par. (package_parin, parin)
Some redundant MPI communication removed. (init_plant_canopy)

Bugfix:
---
Missing KIND-attribute for REAL constant added. (check_parameters)
DO-WHILE-loop for lad-profile output restricted. (header)
Removed double-listing of use_upstream_for_tke in ONLY-list of module
control_parameters. (prognostic_equations)

File:

• palm/trunk/SOURCE/user_init_plant_canopy.f90 (modified) (5 diffs)

palm/trunk/SOURCE/user_init_plant_canopy.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Changes in the course of the canopy-model modularization:
+!   module plant_canopy_model_mod added,
+!   definition of array cdc (canopy drag coefficient) removed, since it is now
+!   defined purely as a single constant value (see module plant_canopy_model)
 ! 
 ! Former revisions:
@@ -57 +60 @@
     
     USE kinds
+
+    USE plant_canopy_model_mod
     
     USE user
@@ -62 +67 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) :: i   !: 
-    INTEGER(iwp) :: j   !: 
+    INTEGER(iwp) :: i   !: running index
+    INTEGER(iwp) :: j   !: running index
 
 !
@@ -69 +74 @@
 
 !
-!-- Set the 3D-arrays lad_s and cdc for user defined canopies
+!-- Set the 3D-array lad_s for user defined canopies
     SELECT CASE ( TRIM( canopy_mode ) )
 
@@ -78 +83 @@
        CASE ( 'user_defined_canopy_1' )
 !
-!--       Here the user can define his own topography. 
-!--       The following lines contain an example in that the 
-!--       plant canopy extends only over the second half of the 
-!--       model domain along x
+!--       Here the user can define his own forest topography. 
+!--       The following lines contain an example, where the plant canopy extends
+!--       only over the second half of the model domain along x.
+!--       Attention: DO-loops have to include the ghost points (nxlg-nxrg, 
+!--       nysg-nyng), because no exchange of ghost point information is intended,
+!--       in order to minimize communication between CPUs
 !          DO  i = nxlg, nxrg
 !             IF ( i >= INT(nx+1/2) ) THEN
 !                DO  j = nysg, nyng
 !                   lad_s(:,j,i) = lad(:)
-!                   cdc(:,j,i)   = drag_coefficient
 !                ENDDO
 !             ELSE
 !                lad_s(:,:,i) = 0.0_wp
-!                cdc(:,:,i)   = 0.0_wp
 !             ENDIF
 !          ENDDO 
+!
 !--       After definition, please
 !--       remove the following three lines!