- Timestamp:
- May 27, 2016 11:05:02 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r1851 r1915 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! Bugfix: avoid segmentation fault in case of most_method = 'circular' at first 22 ! timestep 22 23 ! 23 24 ! Former revisions: … … 718 719 !-- consequence would result in very large shear stresses and very 719 720 !-- small momentum fluxes (both are generally unrealistic). 720 IF ( ( z_mo / ol(j,i) ) < zeta_min ) ol(j,i) = z_mo / zeta_min721 IF ( ( z_mo / ol(j,i) ) > zeta_max ) ol(j,i) = z_mo / zeta_max721 IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) < zeta_min ) ol(j,i) = z_mo / zeta_min 722 IF ( ( z_mo / ( ol(j,i) + 1E-30_wp ) ) > zeta_max ) ol(j,i) = z_mo / zeta_max 722 723 723 724 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.