Ignore:
Timestamp:
Aug 23, 2019 1:49:38 PM (2 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:
1 edited

Legend:

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

    r4182 r4185  
    2525! -----------------
    2626! $Id$
     27! For initializing_actions = ' cyclic_fill':
     28! Overwrite u_init, v_init, pt_init, q_init and s_init with the
     29! (temporally) and horizontally averaged vertical profiles from the end
     30! of the prerun, because these profiles shall be used as the basic state
     31! for the rayleigh damping and the pt_damping.
     32!
     33! 4182 2019-08-22 15:20:23Z scharf
    2734! Corrected "Former revisions" section
    2835!
     
    11301137!--       These profiles are the (temporally) and horizontally averaged vertical
    11311138!--       profiles from the prerun. Alternatively, prescribed profiles
    1132 !--       for u,v-components can be used.
     1139!--       for u,v-components can be used.
     1140!--       Overwrite u_init, v_init, pt_init, q_init and s_init with the
     1141!--       (temporally) and horizontally averaged vertical profiles from the end
     1142!--       of the prerun, because these profiles shall be used as the basic state
     1143!--       for the rayleigh damping and the pt_damping.
    11331144          ALLOCATE( mean_inflow_profiles(nzb:nzt+1,1:num_mean_inflow_profiles) )
    11341145
     
    11381149          ELSE
    11391150             mean_inflow_profiles(:,1) = hom_sum(:,1,0)    ! u
     1151             u_init(:)                 = hom_sum(:,1,0)
    11401152             mean_inflow_profiles(:,2) = hom_sum(:,2,0)    ! v
     1153             v_init(:)                 = hom_sum(:,2,0)
    11411154          ENDIF
    11421155          mean_inflow_profiles(:,4) = hom_sum(:,4,0)       ! pt
     1156          pt_init(:)                = hom_sum(:,4,0)
    11431157          IF ( humidity )                                                      &
    11441158             mean_inflow_profiles(:,6) = hom_sum(:,41,0)   ! q
     1159             q_init(:)                 = hom_sum(:,41,0)   
    11451160          IF ( passive_scalar )                                                &
    11461161             mean_inflow_profiles(:,7) = hom_sum(:,115,0)   ! s
     1162             s_init(:)                 = hom_sum(:,115,0)
    11471163!
    11481164!--       In case of complex terrain, determine vertical displacement at inflow
Note: See TracChangeset for help on using the changeset viewer.