Changeset 1707

Timestamp:

Nov 2, 2015 3:24:52 PM (9 years ago)

Author:

maronga

Message:

bugfix: initialization of Obukhov length from 1d precursor in case of neutral stratification

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: transfer of Richardson number from 1D model to Obukhov length caused
+! devision by zero in neutral stratification
 ! 
 ! Former revisions:
@@ -721 +722 @@
 
              IF ( constant_flux_layer )  THEN
-                ol   = ( zu(nzb+1) - zw(nzb) ) / rif1d(nzb+1)
+                ol   = ( zu(nzb+1) - zw(nzb) ) / ( rif1d(nzb+1) + 1.0E-20_wp )
                 ts   = 0.0_wp  ! could actually be computed more accurately in the
                                ! 1D model. Update when opportunity arises.
@@ -760 +761 @@
 
 !
-!--       inside buildings set velocities back to zero
           IF ( topography /= 'flat' )  THEN
              DO  i = nxl-1, nxr+1