Changeset 1920 for palm/trunk/SOURCE

Timestamp:

May 30, 2016 10:50:15 AM (9 years ago)

Author:

suehring

Message:

Avoid segmentation fault (see change in 1915) by different initialization of us instead of adding a very small number in the denominator

Location:

palm/trunk/SOURCE

Files:

• init_3d_model.f90 (modified) (2 diffs)
• surface_layer_fluxes_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Initialize us with very small number to avoid segmentation fault during 
+! calculation of Obukhov length
 !
 ! Former revisions:
@@ -951 +952 @@
           ol    = ( zu(nzb+1) - zw(nzb) ) / zeta_min
           ts    = 0.0_wp
-          us    = 0.0_wp
+!
+!--       Very small number is required for calculation of Obukhov length 
+!--       at first timestep     
+          us    = 1E-30_wp 
           usws  = 0.0_wp
           uswst = top_momentumflux_u

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Avoid segmentation fault (see change in 1915) by different initialization of 
+! us instead of adding a very small number in the denominator
 ! 
 ! Former revisions:
@@ -722 +723 @@
 !--             consequence would result in very large shear stresses and very
 !--             small momentum fluxes (both are generally unrealistic).
-                IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) < zeta_min )  ol(j,i) = z_mo / zeta_min
-                IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) > zeta_max )  ol(j,i) = z_mo / zeta_max
-
+                IF ( ( z_mo / ol(j,i) ) < zeta_min )  ol(j,i) = z_mo / zeta_min
+                IF ( ( z_mo / ol(j,i) ) > zeta_max )  ol(j,i) = z_mo / zeta_max
              ENDDO
           ENDDO