Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_3d_model.f90

    r4360 r4365  
    2525! -----------------
    2626! $Id$
     27! Overwrite u_init, v_init, pt_init, q_init and s_init with hom for all
     28! cyclic_fill-cases, not only for turbulent_inflow = .TRUE.
     29!
     30! 4360 2020-01-07 11:25:50Z suehring
    2731! Introduction of wall_flags_total_0, which currently sets bits based on static
    2832! topography information used in wall_flags_static_0
     
    11221126#endif
    11231127       ENDDO
    1124 !
    1125 !--    In case of cyclic fill, call init_surfaces a second time, so that
    1126 !--    surface properties such as heat fluxes are initialized as prescribed.
    1127        IF ( TRIM( initializing_actions ) == 'cyclic_fill' )                    &
     1128       
     1129       
     1130       IF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
     1131
     1132!
     1133!--       In case of cyclic fill, call init_surfaces a second time, so that
     1134!--       surface properties such as heat fluxes are initialized as prescribed.
    11281135          CALL init_surfaces
    11291136
     1137!
     1138!--       Overwrite u_init, v_init, pt_init, q_init and s_init with the
     1139!--       horizontally mean (hom) vertical profiles from the end
     1140!--       of the prerun, because these profiles shall be used as the reference
     1141!--       state for the rayleigh damping and the pt_damping. This is especially
     1142!--       important for the use of large_scale_subsidence, because the
     1143!--       reference temperature in the free atmosphere changes in time.
     1144          u_init(:) = hom_sum(:,1,0)
     1145          v_init(:) = hom_sum(:,2,0)
     1146          pt_init(:) = hom_sum(:,4,0)
     1147          IF ( humidity )                                                      &
     1148             q_init(:) = hom_sum(:,41,0)
     1149          IF ( passive_scalar )                                                &
     1150             s_init(:) = hom_sum(:,115,0)
     1151       ENDIF
    11301152!
    11311153!--    In case of complex terrain and cyclic fill method as initialization,
     
    11611183!--       profiles from the prerun. Alternatively, prescribed profiles
    11621184!--       for u,v-components can be used.
    1163 !--       Overwrite u_init, v_init, pt_init, q_init and s_init with the
    1164 !--       (temporally) and horizontally averaged vertical profiles from the end
    1165 !--       of the prerun, because these profiles shall be used as the basic state
    1166 !--       for the rayleigh damping and the pt_damping.
    11671185          ALLOCATE( mean_inflow_profiles(nzb:nzt+1,1:num_mean_inflow_profiles) )
    11681186
     
    11721190          ELSE
    11731191             mean_inflow_profiles(:,1) = hom_sum(:,1,0)    ! u
    1174              u_init(:)                 = hom_sum(:,1,0)
    11751192             mean_inflow_profiles(:,2) = hom_sum(:,2,0)    ! v
    1176              v_init(:)                 = hom_sum(:,2,0)
    11771193          ENDIF
    11781194          mean_inflow_profiles(:,4) = hom_sum(:,4,0)       ! pt
    1179           pt_init(:)                = hom_sum(:,4,0)
    11801195          IF ( humidity )                                                      &
    11811196             mean_inflow_profiles(:,6) = hom_sum(:,41,0)   ! q
    1182              q_init(:)                 = hom_sum(:,41,0)   
    11831197          IF ( passive_scalar )                                                &
    11841198             mean_inflow_profiles(:,7) = hom_sum(:,115,0)   ! s
    1185              s_init(:)                 = hom_sum(:,115,0)
     1199
    11861200!
    11871201!--       In case of complex terrain, determine vertical displacement at inflow
Note: See TracChangeset for help on using the changeset viewer.