Changeset 1484 for palm/trunk/SOURCE/package_parin.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/package_parin.f90 (modified) (9 diffs)

palm/trunk/SOURCE/package_parin.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Changes due to new module structure of the plant canopy model:
+!   module plant_canopy_model_mod added,
+!   new package/namelist canopy_par added, i.e. the canopy model is no longer 
+!   steered over the inipar-namelist,
+!   drag_coefficient, leaf_surface_concentration and scalar_exchange_coefficient
+!   renamed to canopy_drag_coeff, leaf_surface_conc and leaf_scalar_exch_coeff.
+! Changed statement tags in CONTINUE-statement
 ! 
 ! Former revisions:
@@ -103 +109 @@
                write_particle_statistics
 
+    USE plant_canopy_model_mod,                                                &
+        ONLY:  alpha_lad, beta_lad, calc_beta_lad_profile, canopy_drag_coeff,  &
+               canopy_mode, cthf, lad_surface,                                 &
+               lad_vertical_gradient, lad_vertical_gradient_level, lai_beta,   &
+               leaf_scalar_exch_coeff, leaf_surface_conc, pch_index,           &
+               plant_canopy
+
     USE spectrum,                                                              &
         ONLY:  comp_spectra_level, data_output_sp, plot_spectra_level,         &
@@ -110 +123 @@
 
     CHARACTER (LEN=80) ::  line  !:
+
+    NAMELIST /canopy_par/         alpha_lad, beta_lad, canopy_drag_coeff,      &
+                                  canopy_mode, cthf,                           &
+                                  lad_surface,                                 & 
+                                  lad_vertical_gradient,                       &
+                                  lad_vertical_gradient_level,                 &
+                                  lai_beta,                                    &
+                                  leaf_scalar_exch_coeff,                      &
+                                  leaf_surface_conc, pch_index
 
     NAMELIST /dvrp_graphics_par/  clip_dvrp_l, clip_dvrp_n, clip_dvrp_r,       &
@@ -159 +181 @@
     line = ' '
 
+!
+!-- Try to find canopy package
+    REWIND ( 11 )
+    line = ' '
+    DO   WHILE ( INDEX( line, '&canopy_par' ) == 0 )
+       READ ( 11, '(A)', END=10 )  line
+    ENDDO
+    BACKSPACE ( 11 )
+
+!
+!-- Read user-defined namelist
+    READ ( 11, canopy_par )
+
+!
+!-- Set flag that indicates that canopy model is switched on
+    plant_canopy = .TRUE.
+
+!
+!-- Set flag that indicates that the lad-profile shall be calculated by using
+!-- a beta probability density function
+    IF ( alpha_lad /= 9999999.9_wp  .AND.  beta_lad /= 9999999.9_wp )          &
+       calc_beta_lad_profile = .TRUE.
+
+ 10 CONTINUE
+
+
 #if defined( __dvrp_graphics )
     REWIND ( 11 )
     line = ' '
     DO   WHILE ( INDEX( line, '&dvrp_graphics_par' ) == 0 )
-       READ ( 11, '(A)', END=10 )  line
+       READ ( 11, '(A)', END=20 )  line
     ENDDO
     BACKSPACE ( 11 )
@@ -171 +219 @@
     READ ( 11, dvrp_graphics_par )
 
- 10 CONTINUE
+ 20 CONTINUE
 #endif
 
@@ -179 +227 @@
     line = ' '
     DO   WHILE ( INDEX( line, '&particles_par' ) == 0 )
-       READ ( 11, '(A)', END=20 )  line
+       READ ( 11, '(A)', END=30 )  line
     ENDDO
     BACKSPACE ( 11 )
@@ -191 +239 @@
     particle_advection = .TRUE.
 
- 20 CONTINUE
+ 30 CONTINUE
 
 
@@ -198 +246 @@
     line = ' '
     DO   WHILE ( INDEX( line, '&spectra_par' ) == 0 )
-       READ ( 11, '(A)', END=30 )  line
+       READ ( 11, '(A)', END=40 )  line
     ENDDO
     BACKSPACE ( 11 )
@@ -211 +259 @@
     IF ( dt_dosp == 9999999.9_wp )  dt_dosp = dt_data_output
 
- 30 CONTINUE
+ 40 CONTINUE
 #endif