Changeset 1929 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
- Timestamp:
- Jun 9, 2016 4:25:25 PM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r1921 r1929 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! Bugfix: avoid segmentation fault in case one grid point is horizontally 22 ! completely surrounded by topography 22 23 ! 23 24 ! Former revisions: … … 726 727 !-- consequence would result in very large shear stresses and very 727 728 !-- small momentum fluxes (both are generally unrealistic). 728 IF ( ( z_mo / ol(j,i) ) < zeta_min ) ol(j,i) = z_mo / zeta_min 729 IF ( ( z_mo / ol(j,i) ) > zeta_max ) ol(j,i) = z_mo / zeta_max 729 IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) < zeta_min ) & 730 ol(j,i) = z_mo / zeta_min 731 IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) > zeta_max ) & 732 ol(j,i) = z_mo / zeta_max 730 733 ENDDO 731 734 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.