Changeset 4185

Timestamp:

Aug 23, 2019 1:49:38 PM (5 years ago)

Author:

oliver.maas

Message:

For initializing_actions = ' cyclic_fill': Overwrite u_init, v_init, pt_init, q_init and s_init with the (temporally) and horizontally averaged vertical profiles from the end of the prerun, because these profiles shall be used as the basic state for the rayleigh damping and the pt_damping.

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! For initializing_actions = ' cyclic_fill':
+! Overwrite u_init, v_init, pt_init, q_init and s_init with the 
+! (temporally) and horizontally averaged vertical profiles from the end
+! of the prerun, because these profiles shall be used as the basic state
+! for the rayleigh damping and the pt_damping. 
+! 
+! 4182 2019-08-22 15:20:23Z scharf
 ! Corrected "Former revisions" section
 ! 
@@ -1130 +1137 @@
 !--       These profiles are the (temporally) and horizontally averaged vertical
 !--       profiles from the prerun. Alternatively, prescribed profiles
-!--       for u,v-components can be used.
+!--       for u,v-components can be used. 
+!--       Overwrite u_init, v_init, pt_init, q_init and s_init with the 
+!--       (temporally) and horizontally averaged vertical profiles from the end
+!--       of the prerun, because these profiles shall be used as the basic state
+!--       for the rayleigh damping and the pt_damping.
           ALLOCATE( mean_inflow_profiles(nzb:nzt+1,1:num_mean_inflow_profiles) )
 
@@ -1138 +1149 @@
           ELSE
              mean_inflow_profiles(:,1) = hom_sum(:,1,0)    ! u
+             u_init(:)                 = hom_sum(:,1,0)
              mean_inflow_profiles(:,2) = hom_sum(:,2,0)    ! v
+             v_init(:)                 = hom_sum(:,2,0)
           ENDIF
           mean_inflow_profiles(:,4) = hom_sum(:,4,0)       ! pt
+          pt_init(:)                = hom_sum(:,4,0)
           IF ( humidity )                                                      &
              mean_inflow_profiles(:,6) = hom_sum(:,41,0)   ! q
+             q_init(:)                 = hom_sum(:,41,0)    
           IF ( passive_scalar )                                                &
              mean_inflow_profiles(:,7) = hom_sum(:,115,0)   ! s
+             s_init(:)                 = hom_sum(:,115,0)
 !
 !--       In case of complex terrain, determine vertical displacement at inflow