Changeset 2758 for palm/trunk


Ignore:
Timestamp:
Jan 17, 2018 12:55:21 PM (4 years ago)
Author:
suehring
Message:

Bugfix, no initial masking of wind velocity at first prognostic grid level in case of land- or urban-surface spin-up.

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r2746 r2758  
    2525! -----------------
    2626! $Id$
     27! In case of spinup of land- and urban-surface model, do not mask wind velocity
     28! at first computational grid level
     29!
     30! 2746 2018-01-15 12:06:04Z suehring
    2731! Move flag plant canopy to modules
    2832!
     
    12531257!
    12541258!--       Set velocities back to zero
    1255           DO  i = nxlg, nxrg
    1256              DO  j = nysg, nyng
    1257                 DO  k = nzb, nzt
    1258                    u(k,j,i) = MERGE( u(k,j,i), 0.0_wp,                      &
    1259                                      BTEST( wall_flags_0(k,j,i), 1 ) )
    1260                    v(k,j,i) = MERGE( v(k,j,i), 0.0_wp,                      &
    1261                                      BTEST( wall_flags_0(k,j,i), 2 ) )
    1262                 ENDDO
    1263              ENDDO
    1264           ENDDO
    1265          
     1259          u = MERGE( u, 0.0_wp, BTEST( wall_flags_0, 1 ) )
     1260          v = MERGE( v, 0.0_wp, BTEST( wall_flags_0, 2 ) )         
    12661261!
    12671262!--       WARNING: The extra boundary conditions set after running the
     
    13091304             ENDDO
    13101305          ENDDO
    1311 
     1306!
     1307!--       Mask topography
     1308          u = MERGE( u, 0.0_wp, BTEST( wall_flags_0, 1 ) )
     1309          v = MERGE( v, 0.0_wp, BTEST( wall_flags_0, 2 ) )
    13121310!
    13131311!--       Set initial horizontal velocities at the lowest computational grid
    13141312!--       levels to zero in order to avoid too small time steps caused by the
    13151313!--       diffusion limit in the initial phase of a run (at k=1, dz/2 occurs
    1316 !--       in the limiting formula!).
    1317           IF ( ibc_uv_b /= 1 )  THEN
     1314!--       in the limiting formula!).
     1315!--       Please note, in case land- or urban-surface model is used and a
     1316!--       spinup is applied, masking the lowest computational level is not
     1317!--       possible as MOST as well as energy-balance parametrizations will not
     1318!--       work with zero wind velocity.
     1319          IF ( ibc_uv_b /= 1  .AND.  .NOT.  spinup )  THEN
    13181320             DO  i = nxlg, nxrg
    13191321                DO  j = nysg, nyng
  • palm/trunk/SOURCE/time_integration_spinup.f90

    r2728 r2758  
    2525! -----------------
    2626! $Id$
     27! Comment out homogeneous setting of wind velocity as this will lead to zero
     28! friction velocity and cause problems in MOST relationships.
     29!
     30! 2728 2018-01-09 07:03:53Z maronga
    2731! Set velocity componenets to homogeneous values during spinup
    2832!
     
    170174!-- otherwise. The precise reason is still unknown. A minimum velocity of 0.1
    171175!-- m/s is used to maintain turbulent transfer at the surface.
    172     u = SIGN(1.0_wp,u) * MAX(ug_surface,0.1_wp)
    173     v = SIGN(1.0_wp,v) * MAX(vg_surface,0.1_wp)
     176!     u = SIGN(1.0_wp,u) * MAX(ug_surface,0.1_wp)
     177!     v = SIGN(1.0_wp,v) * MAX(vg_surface,0.1_wp)
    174178
    175179    dt_save = dt_3d
Note: See TracChangeset for help on using the changeset viewer.