Changeset 4186 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
- Timestamp:
- Aug 23, 2019 4:06:14 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r4182 r4186 26 26 ! ----------------- 27 27 ! $Id$ 28 ! - To enable limitation of Obukhov length, move it before exit-cycle construct. 29 ! Further, change the limit to 10E-5 in order to get rid-off unrealistic 30 ! peaks in the heat fluxes during nighttime 31 ! - Unused variable removed 32 ! 33 ! 4182 2019-08-22 15:20:23Z scharf 28 34 ! Corrected "Former revisions" section 29 35 ! … … 129 135 INTEGER(iwp) :: l !< loop index for surf type 130 136 131 INTEGER(iwp), PARAMETER :: num_steps = 15000 !< number of steps in the lookup table132 133 137 LOGICAL :: coupled_run !< Flag for coupled atmosphere-ocean runs 134 138 LOGICAL :: downward = .FALSE.!< Flag indicating downward-facing horizontal surface … … 1008 1012 ENDIF 1009 1013 ! 1014 !-- Assure that Obukhov length does not become zero. If the limit is 1015 !-- reached, exit the loop. 1016 IF ( ABS( surf%ol(m) ) < 1E-5_wp ) THEN 1017 surf%ol(m) = SIGN( 1E-5_wp, surf%ol(m) ) 1018 EXIT 1019 ENDIF 1020 ! 1010 1021 !-- Check for convergence 1011 1022 IF ( ABS( ( surf%ol(m) - ol_m ) / surf%ol(m) ) < 1.0E-4_wp ) THEN … … 1013 1024 ELSE 1014 1025 CYCLE 1015 ENDIF 1016 ! 1017 !-- Assure that Obukhov length does not become zero 1018 IF ( ABS( surf%ol(m) ) < 1E-6_wp ) & 1019 surf%ol(m) = SIGN( 10E-6_wp, surf%ol(m) ) 1026 ENDIF 1020 1027 1021 1028 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.