Changeset 138 for palm/trunk/SOURCE/user_interface.f90

Timestamp:

Nov 28, 2007 10:03:58 AM (17 years ago)

Author:

letzel

Message:

Plant canopy model of Watanabe (2004,BLM 112,307-341) added.

File:

• palm/trunk/SOURCE/user_interface.f90 (modified) (2 diffs)

palm/trunk/SOURCE/user_interface.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! new subroutine user_init_plant_canopy
+!
 ! new subroutine user_data_output_dvrp
 ! +argument gls in user_init_grid
@@ -293 +295 @@
 
  END SUBROUTINE user_init_grid
+
+
+
+ SUBROUTINE user_init_plant_canopy
+
+!------------------------------------------------------------------------------!
+!
+! Description:
+! ------------
+! Initialisation of the leaf area density array (for scalar grid points) and 
+! the array of the canopy drag coefficient, if the user has not chosen any 
+! of the default cases 
+!------------------------------------------------------------------------------!
+
+    USE arrays_3d
+    USE control_parameters
+    USE indices
+    USE user
+
+    IMPLICIT NONE
+
+    INTEGER :: i, j 
+
+!
+!-- Here the user-defined grid initializing actions follow:
+
+!
+!-- Set the 3D-arrays lad_s and cdc for user defined canopies
+    SELECT CASE ( TRIM( canopy_mode ) )
+
+       CASE ( 'block' )
+!
+!--       Not allowed here since this is the standard case used in init_3d_model.
+
+       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
+!          DO  i = nxl-1, nxr+1
+!             IF ( i >= INT(nx+1/2) ) THEN
+!                DO  j = nys-1, nyn+1
+!                   lad_s(:,j,i) = lad(:)
+!                   cdc(:,j,i)   = drag_coefficient
+!                ENDDO
+!             ELSE
+!                lad_s(:,:,i) = 0.0
+!                cdc(:,:,i)   = 0.0
+!             ENDIF
+!          ENDDO 
+!--       After definition, please
+!--       remove the following three lines!
+          PRINT*, '+++ user_init_plant_canopy: canopy_mode "', &
+               canopy_mode, '" not available yet'
+
+       CASE DEFAULT
+!
+!--       The DEFAULT case is reached if the parameter canopy_mode contains a
+!--       wrong character string that is neither recognized in init_3d_model nor 
+!--       here in user_init_plant_canopy.
+          PRINT*, '+++ user_init_plant_canopy: unknown canopy_mode "', &
+               canopy_mode, '"'
+          CALL local_stop
+
+    END SELECT
+
+
+ END SUBROUTINE user_init_plant_canopy