Changeset 2618 for palm/trunk/SOURCE/user_init_3d_model.f90

Timestamp:

Nov 16, 2017 3:37:30 PM (6 years ago)

Author:

suehring

Message:

Reorder calls for initializing surface quantities and give example of user-defined initialization of surface quantities

File:

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

palm/trunk/SOURCE/user_init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Provide example for user-defined initialization of surface-related quantities
+! 
+! 2233 2017-05-30 18:08:54Z suehring
 !
 ! 2232 2017-05-30 17:47:52Z suehring
@@ -68 +71 @@
     IMPLICIT NONE
 
+    INTEGER(iwp) ::  l !< running index surface orientation
+    INTEGER(iwp) ::  m !< running index surface elements
+
+!
+!-- Initialization of surface-related quantities.
+!-- The following example shows required initialization of surface quantitites
+!-- at default-type upward-facing surfaces.  
+!   DO  m = 1, surf_def_h(0)%ns
+!      surf_def_h(0)%ol(m)   = ...    ! Obukhov length
+!      surf_def_h(0)%us(m  ) = ...    ! friction velocity
+!      surf_def_h(0)%usws(m) = ...    ! vertical momentum flux, u-component
+!      surf_def_h(0)%vsws(m) = ...    ! vertical momentum flux, v-component
+!      surf_def_h(0)%z0(m)   = ...    ! roughness length for momentum
+!      IF ( .NOT. neutral )  THEN
+!         surf_def_h(0)%ts(m)   = ... ! scaling parameter
+!         surf_def_h(0)%shf(m)  = ... ! surface sensible heat flux
+!         surf_def_h(0)%z0h(m)  = ... ! roughness length for heat
+!      ENDIF
+!      IF ( humditiy )  THEN
+!         surf_def_h(0)%qs(m)   = ... ! scaling parameter
+!         surf_def_h(0)%qsws(m) = ... ! surface latent heat flux
+!         surf_def_h(0)%z0q(m)  = ... ! roughness length for moisture
+!      ENDIF
+!      IF ( passive_scalar )  THEN
+!         surf_def_h(0)%ss(m)   = ... ! scaling parameter
+!         surf_def_h(0)%ssws(m) = ... ! surface latent heat flux
+!      ENDIF
+!   ENDDO 
+!
+!-- Same for natural and urban type surfaces
+!   DO  m = 1, surf_lsm_h%ns
+!      ...
+!   ENDDO 
+!   DO  m = 1, surf_usm_h%ns
+!      ...
+!   ENDDO
+!
+!-- Also care for vertically aligned surfaces (default-, natural-, and 
+!-- urban-type).
+!   DO  l = 0, 3
+!      DO  m = 1, surf_def_v(l)%ns
+!         ...
+!      ENDDO
+!      DO  m = 1, surf_lsm_v(l)%ns
+!         ...
+!      ENDDO
+!      DO  m = 1, surf_usm_v(l)%ns
+!         ...
+!      ENDDO
+!   ENDDO
+!
+!
+!-- In the following, initialize 3D quantities, e.g. u, v, w, pt, etc..
 
  END SUBROUTINE user_init_3d_model