Changeset 2983


Ignore:
Timestamp:
Apr 18, 2018 10:43:40 AM (3 years ago)
Author:
suehring
Message:

Revise limitation of wall-adjacent velocity in spinup

File:
1 edited

Legend:

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

    r2934 r2983  
    2525! -----------------
    2626! $Id$
     27! Revise limitation of wall-adjacent velocity.
     28!
     29! 2934 2018-03-26 19:13:22Z suehring
    2730! Synchronize parent and child models after spinup.
    2831!
     
    196199          j   = surf_lsm_h%j(m)
    197200          k   = surf_lsm_h%k(m)
    198           u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX(u_init(k),0.1_wp)
    199           v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX(v_init(k),0.1_wp)
     201          u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX( ABS( u_init(k) ),0.1_wp)
     202          v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX( ABS( v_init(k) ),0.1_wp)
    200203       ENDDO
    201204
     
    205208             j   = surf_lsm_v(l)%j(m)
    206209             k   = surf_lsm_v(l)%k(m)
    207              u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX(u_init(k),0.1_wp)
    208              v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX(v_init(k),0.1_wp)
     210             u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX( ABS( u_init(k) ),0.1_wp)
     211             v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX( ABS( v_init(k) ),0.1_wp)
    209212          ENDDO
    210213       ENDDO
     
    216219          j   = surf_usm_h%j(m)
    217220          k   = surf_usm_h%k(m)
    218           u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX(u_init(k),0.1_wp)
    219           v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX(v_init(k),0.1_wp)
     221          u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX( ABS( u_init(k) ),0.1_wp)
     222          v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX( ABS( v_init(k) ),0.1_wp)
    220223       ENDDO
    221224
     
    225228             j   = surf_usm_v(l)%j(m)
    226229             k   = surf_usm_v(l)%k(m)
    227              u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX(u_init(k),0.1_wp)
    228              v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX(v_init(k),0.1_wp)
     230             u(k,j,i) = SIGN(1.0_wp,u_init(k)) * MAX( ABS( u_init(k) ),0.1_wp)
     231             v(k,j,i) = SIGN(1.0_wp,v_init(k)) * MAX( ABS( v_init(k) ),0.1_wp)
    229232          ENDDO
    230233       ENDDO
Note: See TracChangeset for help on using the changeset viewer.