Changeset 2782 for palm

Timestamp:

Feb 2, 2018 11:51:10 AM (6 years ago)

Author:

maronga

Message:

bugfix/adjusment in spinup mechanism

File:

• palm/trunk/SOURCE/time_integration_spinup.f90 (modified) (2 diffs)

palm/trunk/SOURCE/time_integration_spinup.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix and re-activation of homogeneous setting of velocity components 
+! during spinup
+! 
+! 2758 2018-01-17 12:55:21Z suehring
 ! Comment out homogeneous setting of wind velocity as this will lead to zero
 ! friction velocity and cause problems in MOST relationships.
@@ -174 +178 @@
 !-- otherwise. The precise reason is still unknown. A minimum velocity of 0.1
 !-- m/s is used to maintain turbulent transfer at the surface.
-!     u = SIGN(1.0_wp,u) * MAX(ug_surface,0.1_wp)
-!     v = SIGN(1.0_wp,v) * MAX(vg_surface,0.1_wp)
+    IF ( land_surface )  THEN
+       DO  m = 1, surf_lsm_h%ns
+          i   = surf_lsm_h%i(m)            
+          j   = surf_lsm_h%j(m)
+          k   = surf_lsm_h%k(m)
+          u(k,j,i) =  SIGN(1.0_wp,u(k,j,i)) * MAX(ug_surface,0.1_wp)
+          v(k,j,i) =  SIGN(1.0_wp,v(k,j,i)) * MAX(vg_surface,0.1_wp)
+       ENDDO
+
+       DO  l = 0, 3
+          DO  m = 1, surf_lsm_v(l)%ns
+             i   = surf_lsm_v(l)%i(m)            
+             j   = surf_lsm_v(l)%j(m)
+             k   = surf_lsm_v(l)%k(m)
+             u(k,j,i) =  SIGN(1.0_wp,u(k,j,i)) * MAX(ug_surface,0.1_wp)
+             v(k,j,i) =  SIGN(1.0_wp,v(k,j,i)) * MAX(vg_surface,0.1_wp)
+          ENDDO
+       ENDDO
+    ENDIF
+
+    IF ( urban_surface )  THEN
+       DO  m = 1, surf_usm_h%ns
+          i   = surf_usm_h%i(m)            
+          j   = surf_usm_h%j(m)
+          k   = surf_usm_h%k(m)
+          u(k,j,i) =  SIGN(1.0_wp,u(k,j,i)) * MAX(ug_surface,0.1_wp)
+          v(k,j,i) =  SIGN(1.0_wp,v(k,j,i)) * MAX(vg_surface,0.1_wp)
+       ENDDO
+
+       DO  l = 0, 3
+          DO  m = 1, surf_usm_v(l)%ns
+             i   = surf_usm_v(l)%i(m)            
+             j   = surf_usm_v(l)%j(m)
+             k   = surf_usm_v(l)%k(m)
+             u(k,j,i) =  SIGN(1.0_wp,u(k,j,i)) * MAX(ug_surface,0.1_wp)
+             v(k,j,i) =  SIGN(1.0_wp,v(k,j,i)) * MAX(vg_surface,0.1_wp)
+          ENDDO
+       ENDDO
+    ENDIF
 
     dt_save = dt_3d