Changeset 1915 for palm/trunk

Timestamp:

May 27, 2016 11:05:02 AM (8 years ago)

Author:

suehring

Message:

Bugfix: avoid segmentation fault in case of most_method = 'circular' at first timstep

File:

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

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: avoid segmentation fault in case of most_method = 'circular' at first
+! timestep
 ! 
 ! Former revisions:
@@ -718 +719 @@
 !--             consequence would result in very large shear stresses and very
 !--             small momentum fluxes (both are generally unrealistic).
-                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
+                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
 
              ENDDO