Changeset 1386


Ignore:
Timestamp:
May 5, 2014 1:55:30 PM (7 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:
1 edited

Legend:

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

    r1383 r1386  
    2121! Current revisions:
    2222! -----------------
    23 !
     23! bugfix: simulated time before the last timestep is needed for the correct
     24! calculation of the profiles of large scale forcing tendencies
    2425!
    2526! Former revisions:
     
    10091010!--       Interpolation in time of LSF_DATA
    10101011          nt = 1
    1011           DO WHILE ( simulated_time > time_vert(nt) )
     1012          DO WHILE ( simulated_time - dt_3d > time_vert(nt) )
    10121013             nt = nt + 1
    10131014          ENDDO
    1014           IF ( simulated_time /= time_vert(nt) )  THEN
     1015          IF ( simulated_time - dt_3d /= time_vert(nt) )  THEN
    10151016            nt = nt - 1
    10161017          ENDIF
    10171018
    1018           fac = ( simulated_time-time_vert(nt) )                               &
     1019          fac = ( simulated_time - dt_3d - time_vert(nt) )                     &
    10191020                / ( time_vert(nt+1)-time_vert(nt) )
    10201021
     
    27642765!--       Interpolation in time of LSF_DATA
    27652766          nt = 1
    2766           DO WHILE ( simulated_time > time_vert(nt) )
     2767          DO WHILE ( simulated_time - dt_3d > time_vert(nt) )
    27672768             nt = nt + 1
    27682769          ENDDO
    2769           IF ( simulated_time /= time_vert(nt) )  THEN
     2770          IF ( simulated_time - dt_3d /= time_vert(nt) )  THEN
    27702771            nt = nt - 1
    27712772          ENDIF
    27722773
    2773           fac = ( simulated_time-time_vert(nt) )                               &
     2774          fac = ( simulated_time - dt_3d - time_vert(nt) )                     &
    27742775                / ( time_vert(nt+1)-time_vert(nt) )
    27752776
Note: See TracChangeset for help on using the changeset viewer.