Changeset 2618 for palm/trunk

Timestamp:

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

Author:

suehring

Message:

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

Location:

palm/trunk/SOURCE

Files:

• init_3d_model.f90 (modified) (5 diffs)
• user_init_3d_model.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Reorder calls of init_surfaces.
+! 
+! 2564 2017-10-19 15:56:56Z Giersch
 ! Variable wind_turbine was added to control_parameters.
 ! 
@@ -1135 +1138 @@
 
           ENDIF
+!
+!--       Initialize surface variables, e.g. friction velocity, momentum 
+!--       fluxes, etc. 
+          CALL init_surfaces
 
           CALL location_message( 'finished', .TRUE. )
@@ -1247 +1254 @@
 !--       of a sloping surface
           IF ( sloping_surface )  CALL init_slope
+!
+!--       Initialize surface variables, e.g. friction velocity, momentum 
+!--       fluxes, etc. 
+          CALL init_surfaces
 
           CALL location_message( 'finished', .TRUE. )
@@ -1254 +1265 @@
 
           CALL location_message( 'initializing by user', .FALSE. )
+!
+!--       Pre-initialize surface variables, i.e. setting start- and end-indices
+!--       at each (j,i)-location. Please note, this does not supersede 
+!--       user-defined initialization of surface quantities. 
+          CALL init_surfaces
 !
 !--       Initialization will completely be done by the user
@@ -1924 +1940 @@
     ENDIF
 !
-!-- Initialize surface elements and its attributes, e.g. heat- and 
-!-- momentumfluxes, roughness, scaling parameters. 
-!-- This is already done in case of restart data.  
+!-- Finally, if random_heatflux is set, disturb shf at horizontal
+!-- surfaces. Actually, this should be done in surface_mod, where all other 
+!-- initializations of surface quantities are done. However, this
+!-- would create a ring dependency, hence, it is done here. Maybe delete
+!-- disturb_heatflux and tranfer the respective code directly into the 
+!-- initialization in surface_mod.         
     IF ( TRIM( initializing_actions ) /= 'read_restart_data'  .AND.            &
          TRIM( initializing_actions ) /= 'cyclic_fill' )  THEN
-       CALL init_surfaces
-!
-!--    Finally, if random_heatflux is set, disturb shf at horizontal
-!--    surfaces. Actually, this should be done in surface_mod, where all other 
-!--    initializations of surface quantities are done. However, this
-!--    would create a ring dependency, hence, it is done here. Maybe delete
-!--    disturb_heatflux and tranfer the respective code directly into the 
-!--    initialization in surface_mod.         
+ 
        IF ( use_surface_fluxes  .AND.  constant_heatflux  .AND.                &
             random_heatflux )  THEN

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