Changeset 2983 for palm/trunk
 Timestamp:
 Apr 18, 2018 10:43:40 AM (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/time_integration_spinup.f90
r2934 r2983 25 25 !  26 26 ! $Id$ 27 ! Revise limitation of walladjacent velocity. 28 ! 29 ! 2934 20180326 19:13:22Z suehring 27 30 ! Synchronize parent and child models after spinup. 28 31 ! … … 196 199 j = surf_lsm_h%j(m) 197 200 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) 200 203 ENDDO 201 204 … … 205 208 j = surf_lsm_v(l)%j(m) 206 209 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) 209 212 ENDDO 210 213 ENDDO … … 216 219 j = surf_usm_h%j(m) 217 220 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) 220 223 ENDDO 221 224 … … 225 228 j = surf_usm_v(l)%j(m) 226 229 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) 229 232 ENDDO 230 233 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.