Changeset 4365 for palm

Timestamp:

Jan 8, 2020 2:29:20 PM (5 years ago)

Author:

oliver.maas

Message:

Since revision 4185 all init-profiles (u_init, v_init, pt_init, q_init and s_init) are overwritten by the horizontally (and temporally) averaged profiles of the end of the precursor run, for initializing_actions == 'cyclic_fill' and turbulent_inflow = .TRUE.. This commit extends this functionality to all cyclic_fill-cases.

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Overwrite u_init, v_init, pt_init, q_init and s_init with hom for all
+! cyclic_fill-cases, not only for turbulent_inflow = .TRUE.
+! 
+! 4360 2020-01-07 11:25:50Z suehring
 ! Introduction of wall_flags_total_0, which currently sets bits based on static
 ! topography information used in wall_flags_static_0
@@ -1122 +1126 @@
 #endif
        ENDDO
-!
-!--    In case of cyclic fill, call init_surfaces a second time, so that 
-!--    surface properties such as heat fluxes are initialized as prescribed. 
-       IF ( TRIM( initializing_actions ) == 'cyclic_fill' )                    &
+       
+       
+       IF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
+
+!
+!--       In case of cyclic fill, call init_surfaces a second time, so that 
+!--       surface properties such as heat fluxes are initialized as prescribed. 
           CALL init_surfaces
 
+!
+!--       Overwrite u_init, v_init, pt_init, q_init and s_init with the
+!--       horizontally mean (hom) vertical profiles from the end
+!--       of the prerun, because these profiles shall be used as the reference
+!--       state for the rayleigh damping and the pt_damping. This is especially
+!--       important for the use of large_scale_subsidence, because the
+!--       reference temperature in the free atmosphere changes in time.
+          u_init(:) = hom_sum(:,1,0)
+          v_init(:) = hom_sum(:,2,0)
+          pt_init(:) = hom_sum(:,4,0)
+          IF ( humidity )                                                      &
+             q_init(:) = hom_sum(:,41,0)
+          IF ( passive_scalar )                                                &
+             s_init(:) = hom_sum(:,115,0)
+       ENDIF
 !
 !--    In case of complex terrain and cyclic fill method as initialization,
@@ -1161 +1183 @@
 !--       profiles from the prerun. Alternatively, prescribed profiles
 !--       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) )
 
@@ -1172 +1190 @@
           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