Changeset 402 for palm/trunk


Ignore:
Timestamp:
Oct 21, 2009 11:59:41 AM (15 years ago)
Author:
maronga
Message:

bugfix: calculation of time-averaged surface-heatfluxes

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r398 r402  
    1414Bugfix: exchange of ghost points for prho included (time_integration)
    1515
    16 time_integration
     16Bugfix: calculation of time-averaged surface heatfluxes (sum_up_3d_data)
     17
     18sum_up_3d_data, time_integration
    1719
    1820
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r392 r402  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix in calculation of shf*_av, qsws*_av
    77!
    88! Former revisions:
     
    321321                ENDDO
    322322             ENDDO
    323              
     323
    324324          CASE ( 'ql' )
    325325             DO  i = nxl-1, nxr+1
     
    358358             ENDDO
    359359
     360          CASE ( 'qsws*' )
     361             DO  i = nxl-1, nxr+1
     362                DO  j = nys-1, nyn+1
     363                   qsws_av(j,i) = qsws_av(j,i) + qsws(j,i)
     364                ENDDO
     365             ENDDO
     366
    360367          CASE ( 'qv' )
    361368             DO  i = nxl-1, nxr+1
     
    375382                ENDDO
    376383             ENDDO
    377              
     384
    378385          CASE ( 's' )
    379386             DO  i = nxl-1, nxr+1
     
    384391                ENDDO
    385392             ENDDO
    386              
     393
    387394          CASE ( 'sa' )
    388395             DO  i = nxl-1, nxr+1
     
    393400                ENDDO
    394401             ENDDO
    395              
     402
     403          CASE ( 'shf*' )
     404             DO  i = nxl-1, nxr+1
     405                DO  j = nys-1, nyn+1
     406                   shf_av(j,i) = shf_av(j,i) + shf(j,i)
     407                ENDDO
     408             ENDDO
     409
    396410          CASE ( 't*' )
    397411             DO  i = nxl-1, nxr+1
Note: See TracChangeset for help on using the changeset viewer.