Ignore:
Timestamp:
May 30, 2016 10:50:15 AM (5 years ago)
Author:
suehring
Message:

Avoid segmentation fault (see change in 1915) by different initialization of us instead of adding a very small number in the denominator

File:
1 edited

Legend:

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

    r1916 r1920  
    1919! Current revisions:
    2020! ------------------
    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
    2223!
    2324! Former revisions:
     
    722723!--             consequence would result in very large shear stresses and very
    723724!--             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
    727727             ENDDO
    728728          ENDDO
Note: See TracChangeset for help on using the changeset viewer.