Ignore:
Timestamp:
Jun 29, 2020 8:49:58 AM (7 months ago)
Author:
suehring
Message:

mesoscale nesting: omit explicit pressure forcing via geostrophic wind components; chemistry: enable profile output of vertical fluxes; urban-surface: bugfix in initialization in case of cyclic_fill

File:
1 edited

Legend:

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

    r4564 r4581  
    2525! -----------------
    2626! $Id$
     27! Omit explicit pressure forcing via geostrophic wind components in case of
     28! mesoscale nesting.
     29!
     30! 4561 2020-06-12 07:05:56Z suehring
    2731! Adapt mass-flux correction also for the anelastic approximation
    2832!
     
    471475!
    472476!--    Read geostrophic wind components
    473        DO  t = nest_offl%tind, nest_offl%tind_p
    474           CALL get_variable_pr( pids_id, 'ls_forcing_ug', t+1,                 &
    475                                 nest_offl%ug(t-nest_offl%tind,nzb+1:nzt) )
    476           CALL get_variable_pr( pids_id, 'ls_forcing_vg', t+1,                 &
    477                                 nest_offl%vg(t-nest_offl%tind,nzb+1:nzt) )
    478        ENDDO
     477!        DO  t = nest_offl%tind, nest_offl%tind_p
     478!           CALL get_variable_pr( pids_id, 'ls_forcing_ug', t+1,                 &
     479!                                 nest_offl%ug(t-nest_offl%tind,nzb+1:nzt) )
     480!           CALL get_variable_pr( pids_id, 'ls_forcing_vg', t+1,                 &
     481!                                 nest_offl%vg(t-nest_offl%tind,nzb+1:nzt) )
     482!        ENDDO
    479483!
    480484!--    Read data at lateral and top boundaries. Please note, at left and
     
    15671571! Description:
    15681572!------------------------------------------------------------------------------!
    1569 !>  Update of the geostrophic wind components.
    1570 !>  @todo: update geostrophic wind also in the child domains (should be done
    1571 !>         in the nesting.
     1573!>  Update of the geostrophic wind components. Note, currently this routine is
     1574!>  not invoked.
    15721575!------------------------------------------------------------------------------!
    15731576    SUBROUTINE nesting_offl_geostrophic_wind
     
    21352138          ENDIF
    21362139!
    2137 !--       Initialize geostrophic wind components. Actually this is already done in
    2138 !--       init_3d_model when initializing_action = 'inifor', however, in speical
    2139 !--       case of user-defined initialization this will be done here again, in
    2140 !--       order to have a consistent initialization.
    2141           ug(nzb+1:nzt) = nest_offl%ug(0,nzb+1:nzt)
    2142           vg(nzb+1:nzt) = nest_offl%vg(0,nzb+1:nzt)
    2143 !
    2144 !--       Set bottom and top boundary condition for geostrophic wind components
    2145           ug(nzt+1) = ug(nzt)
    2146           vg(nzt+1) = vg(nzt)
    2147           ug(nzb)   = ug(nzb+1)
    2148           vg(nzb)   = vg(nzb+1)
     2140!--       In case of offline nesting the pressure forms itself based on the prescribed lateral
     2141!--       boundary conditions. Hence, explicit forcing by pressure gradients via geostrophic
     2142!--       wind components is not necessary and would be canceled out by the perturbation
     2143!--       pressure otherwise. For this reason, set geostrophic wind components to zero.
     2144          ug(nzb+1:nzt) = 0.0_wp
     2145          vg(nzb+1:nzt) = 0.0_wp
     2146
    21492147       ENDIF
    21502148!
Note: See TracChangeset for help on using the changeset viewer.