Changeset 1709 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Nov 4, 2015 2:47:01 PM (8 years ago)

Author:

maronga

Message:

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

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (3 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Updated output of Obukhov length
 ! 
 ! Former revisions:
@@ -1554 +1554 @@
        ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)    ! v'w'    at k=0
        ts_value(21,sr) = hom(nzb,1,48,sr)           ! w"q"    at k=0
- 
-!
-!--    Calculate obukhov length
-       IF ( ts_value(5,sr) /= 0.0_wp )  THEN
-!           IF ( most_method == 'circular' )  THEN
-!              ts_value(22,sr) = ts_value(4,sr)**2 /                             &
-!                           ( kappa * g * ts_value(5,sr) / ts_value(18,sr) ) 
-!           ELSE
-             ts_value(22,sr) = hom(nzb,1,114,sr)  
-!          ENDIF     
+
+       IF ( .NOT. neutral )  THEN
+          ts_value(22,sr) = hom(nzb,1,114,sr)          ! L
        ELSE
-          ts_value(22,sr) = 10000.0_wp
+          ts_value(22,sr) = 1.0E10_wp
        ENDIF
 
@@ -3503 +3496 @@
        ts_value(21,sr) = hom(nzb,1,48,sr)           ! w"q"    at k=0
 
-       IF ( ts_value(5,sr) /= 0.0_wp )  THEN
-          ts_value(22,sr) = hom(nzb,1,114,sr)
+       IF ( .NOT. neutral )  THEN
+          ts_value(22,sr) = hom(nzb,1,114,sr)          ! L
        ELSE
-          ts_value(22,sr) = 10000.0_wp
+          ts_value(22,sr) = 1.0E10_wp
        ENDIF