Changeset 1386 for palm

Timestamp:

May 5, 2014 1:55:30 PM (11 years ago)

Author:

boeske

Message:

bugfix: simulated time before last timestep must be used for the profile data output of large scale forcing tendencies

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -21 +21 @@
 ! Current revisions:
 ! -----------------
-! 
+! bugfix: simulated time before the last timestep is needed for the correct 
+! calculation of the profiles of large scale forcing tendencies 
 ! 
 ! Former revisions:
@@ -1009 +1010 @@
 !--       Interpolation in time of LSF_DATA 
           nt = 1
-          DO WHILE ( simulated_time > time_vert(nt) )
+          DO WHILE ( simulated_time - dt_3d > time_vert(nt) )
              nt = nt + 1
           ENDDO
-          IF ( simulated_time /= time_vert(nt) )  THEN
+          IF ( simulated_time - dt_3d /= time_vert(nt) )  THEN
             nt = nt - 1
           ENDIF
 
-          fac = ( simulated_time-time_vert(nt) )                               &
+          fac = ( simulated_time - dt_3d - time_vert(nt) )                     &
                 / ( time_vert(nt+1)-time_vert(nt) )
 
@@ -2764 +2765 @@
 !--       Interpolation in time of LSF_DATA 
           nt = 1
-          DO WHILE ( simulated_time > time_vert(nt) )
+          DO WHILE ( simulated_time - dt_3d > time_vert(nt) )
              nt = nt + 1
           ENDDO
-          IF ( simulated_time /= time_vert(nt) )  THEN
+          IF ( simulated_time - dt_3d /= time_vert(nt) )  THEN
             nt = nt - 1
           ENDIF
 
-          fac = ( simulated_time-time_vert(nt) )                               &
+          fac = ( simulated_time - dt_3d - time_vert(nt) )                     &
                 / ( time_vert(nt+1)-time_vert(nt) )