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

palm/trunk/SOURCE/init_3d_model.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Changes due to new module structure of the plant canopy model:
+!   canopy-related initialization (e.g. lad and canopy_heat_flux) moved to new 
+!   subroutine init_plant_canopy within the module plant_canopy_model_mod,
+!   call of subroutine init_plant_canopy added.
 ! 
 ! Former revisions:
@@ -217 +220 @@
     USE pegrid
     
+    USE plant_canopy_model_mod,                                                &
+        ONLY:  init_plant_canopy, plant_canopy
+    
     USE random_function_mod 
     
@@ -517 +523 @@
     
 !
-!-- 3D-arrays for the leaf area density and the canopy drag coefficient
-    IF ( plant_canopy ) THEN
-       ALLOCATE ( lad_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
-                  lad_u(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
-                  lad_v(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
-                  lad_w(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
-                  cdc(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-
-       IF ( passive_scalar ) THEN
-          ALLOCATE ( sls(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
-                     sec(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 
-       ENDIF
-
-       IF ( cthf /= 0.0_wp ) THEN
-          ALLOCATE ( lai(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
-                     canopy_heat_flux(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-       ENDIF
-
-    ENDIF
-
-!
 !-- 4D-array for storing the Rif-values at vertical walls
     IF ( topography /= 'flat' )  THEN
@@ -1571 +1556 @@
 
 !
-!-- Initialization of the leaf area density
-    IF ( plant_canopy )  THEN
- 
-       SELECT CASE ( TRIM( canopy_mode ) )
-
-          CASE( 'block' )
-
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   lad_s(:,j,i) = lad(:)
-                   cdc(:,j,i)   = drag_coefficient
-                   IF ( passive_scalar )  THEN
-                      sls(:,j,i) = leaf_surface_concentration
-                      sec(:,j,i) = scalar_exchange_coefficient
-                   ENDIF
-                ENDDO
-             ENDDO
-
-          CASE DEFAULT
-
-!
-!--          The DEFAULT case is reached either if the parameter 
-!--          canopy mode contains a wrong character string or if the 
-!--          user has coded a special case in the user interface. 
-!--          There, the subroutine user_init_plant_canopy checks 
-!--          which of these two conditions applies.
-             CALL user_init_plant_canopy
- 
-          END SELECT
-
-       CALL exchange_horiz( lad_s, nbgp )
-       CALL exchange_horiz( cdc, nbgp )
-
-       IF ( passive_scalar )  THEN
-          CALL exchange_horiz( sls, nbgp )
-          CALL exchange_horiz( sec, nbgp )
-       ENDIF
-
-!
-!--    Sharp boundaries of the plant canopy in horizontal directions
-!--    In vertical direction the interpolation is retained, as the leaf
-!--    area density is initialised by prescribing a vertical profile 
-!--    consisting of piecewise linear segments. The upper boundary 
-!--    of the plant canopy is now defined by lad_w(pch_index,:,:) = 0.0.
-
-       DO  i = nxl, nxr
-          DO  j = nys, nyn
-             DO  k = nzb, nzt+1 
-                IF ( lad_s(k,j,i) > 0.0_wp )  THEN
-                   lad_u(k,j,i)   = lad_s(k,j,i) 
-                   lad_u(k,j,i+1) = lad_s(k,j,i)
-                   lad_v(k,j,i)   = lad_s(k,j,i)
-                   lad_v(k,j+1,i) = lad_s(k,j,i)
-                ENDIF
-             ENDDO
-             DO  k = nzb, nzt
-                lad_w(k,j,i) = 0.5_wp * ( lad_s(k+1,j,i) + lad_s(k,j,i) )
-             ENDDO
-          ENDDO
-       ENDDO
-
-       lad_w(pch_index,:,:) = 0.0_wp
-       lad_w(nzt+1,:,:)     = lad_w(nzt,:,:)
-
-       CALL exchange_horiz( lad_u, nbgp )
-       CALL exchange_horiz( lad_v, nbgp )
-       CALL exchange_horiz( lad_w, nbgp )
-
-!
-!--    Initialisation of the canopy heat source distribution
-       IF ( cthf /= 0.0_wp )  THEN
-!
-!--       Piecewise evaluation of the leaf area index by 
-!--       integration of the leaf area density
-          lai(:,:,:) = 0.0_wp
-          DO  i = nxlg, nxrg
-             DO  j = nysg, nyng
-                DO  k = pch_index-1, 0, -1
-                   lai(k,j,i) = lai(k+1,j,i) +                   &
-                                ( 0.5_wp * ( lad_w(k+1,j,i) +    &
-                                          lad_s(k+1,j,i) ) *     &
-                                  ( zw(k+1) - zu(k+1) ) )  +     &
-                                ( 0.5_wp * ( lad_w(k,j,i)   +    &
-                                          lad_s(k+1,j,i) ) *     &
-                                  ( zu(k+1) - zw(k) ) )
-                ENDDO
-             ENDDO
-          ENDDO
-
-!
-!--       Evaluation of the upward kinematic vertical heat flux within the 
-!--       canopy
-          DO  i = nxlg, nxrg
-             DO  j = nysg, nyng
-                DO  k = 0, pch_index
-                   canopy_heat_flux(k,j,i) = cthf *                    &
-                                             exp( -0.6_wp * lai(k,j,i) )
-                ENDDO
-             ENDDO
-          ENDDO
-
-!
-!--       The near surface heat flux is derived from the heat flux 
-!--       distribution within the canopy
-          shf(:,:) = canopy_heat_flux(0,:,:)
-
-       ENDIF
-
-    ENDIF
+!-- If required, initialize quantities needed for the plant canopy model
+    IF ( plant_canopy )  CALL init_plant_canopy
 
 !