Ignore:
Timestamp:
Aug 23, 2019 4:06:14 PM (22 months ago)
Author:
suehring
Message:

Enable limitation of Obukhov length so that it does not become zero; to read input data from netcdf in init_3d_model use provided module variables instead of defining local ones; tests updated because changes in Obukhov lengths causes small differences during the initial phase of the run

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r4182 r4186  
    2626! -----------------
    2727! $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
    2834! Corrected "Former revisions" section
    2935!
     
    129135    INTEGER(iwp) ::  l              !< loop index for surf type
    130136
    131     INTEGER(iwp), PARAMETER ::  num_steps = 15000  !< number of steps in the lookup table
    132 
    133137    LOGICAL      ::  coupled_run       !< Flag for coupled atmosphere-ocean runs
    134138    LOGICAL      ::  downward = .FALSE.!< Flag indicating downward-facing horizontal surface
     
    10081012             ENDIF
    10091013!
     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!
    10101021!--          Check for convergence
    10111022             IF ( ABS( ( surf%ol(m) - ol_m ) /  surf%ol(m) ) < 1.0E-4_wp )  THEN
     
    10131024             ELSE
    10141025                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             
    10201027
    10211028          ENDDO
Note: See TracChangeset for help on using the changeset viewer.