Changeset 1920 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
- Timestamp:
- May 30, 2016 10:50:15 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r1916 r1920 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! Avoid segmentation fault (see change in 1915) by different initialization of 22 ! us instead of adding a very small number in the denominator 22 23 ! 23 24 ! Former revisions: … … 722 723 !-- consequence would result in very large shear stresses and very 723 724 !-- small momentum fluxes (both are generally unrealistic). 724 IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) < zeta_min ) ol(j,i) = z_mo / zeta_min 725 IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) > zeta_max ) ol(j,i) = z_mo / zeta_max 726 725 IF ( ( z_mo / ol(j,i) ) < zeta_min ) ol(j,i) = z_mo / zeta_min 726 IF ( ( z_mo / ol(j,i) ) > zeta_max ) ol(j,i) = z_mo / zeta_max 727 727 ENDDO 728 728 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.