Ignore:
Timestamp:
Nov 4, 2015 2:47:01 PM (6 years ago)
Author:
maronga
Message:

several bugfixes related to the new surface layer routine and land-surface-radiation interaction

File:
1 edited

Legend:

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

    r1701 r1709  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Updated output of Obukhov length
    2222!
    2323! Former revisions:
     
    15541554       ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)    ! v'w'    at k=0
    15551555       ts_value(21,sr) = hom(nzb,1,48,sr)           ! w"q"    at k=0
    1556  
    1557 !
    1558 !--    Calculate obukhov length
    1559        IF ( ts_value(5,sr) /= 0.0_wp )  THEN
    1560 !           IF ( most_method == 'circular' )  THEN
    1561 !              ts_value(22,sr) = ts_value(4,sr)**2 /                             &
    1562 !                           ( kappa * g * ts_value(5,sr) / ts_value(18,sr) )
    1563 !           ELSE
    1564              ts_value(22,sr) = hom(nzb,1,114,sr) 
    1565 !          ENDIF     
     1556
     1557       IF ( .NOT. neutral )  THEN
     1558          ts_value(22,sr) = hom(nzb,1,114,sr)          ! L
    15661559       ELSE
    1567           ts_value(22,sr) = 10000.0_wp
     1560          ts_value(22,sr) = 1.0E10_wp
    15681561       ENDIF
    15691562
     
    35033496       ts_value(21,sr) = hom(nzb,1,48,sr)           ! w"q"    at k=0
    35043497
    3505        IF ( ts_value(5,sr) /= 0.0_wp )  THEN
    3506           ts_value(22,sr) = hom(nzb,1,114,sr)
     3498       IF ( .NOT. neutral )  THEN
     3499          ts_value(22,sr) = hom(nzb,1,114,sr)          ! L
    35073500       ELSE
    3508           ts_value(22,sr) = 10000.0_wp
     3501          ts_value(22,sr) = 1.0E10_wp
    35093502       ENDIF
    35103503
Note: See TracChangeset for help on using the changeset viewer.