Ignore:
Timestamp:
Jun 4, 2007 8:07:41 AM (14 years ago)
Author:
raasch
Message:

more preliminary uncomplete changes for ocean version

File:
1 edited

Legend:

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

    r77 r96  
    44! Actual revisions:
    55! -----------------
    6 !
     6! +sum-up of density and salinity
    77!
    88! Former revisions:
     
    132132                qv_av = 0.0
    133133
     134             CASE ( 'rho' )
     135                IF ( .NOT. ALLOCATED( rho_av ) )  THEN
     136                   ALLOCATE( rho_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
     137                ENDIF
     138                rho_av = 0.0
     139
    134140             CASE ( 's' )
    135141                IF ( .NOT. ALLOCATED( s_av ) )  THEN
     
    137143                ENDIF
    138144                s_av = 0.0
     145
     146             CASE ( 'sa' )
     147                IF ( .NOT. ALLOCATED( sa_av ) )  THEN
     148                   ALLOCATE( sa_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
     149                ENDIF
     150                sa_av = 0.0
    139151
    140152             CASE ( 't*' )
     
    337349             ENDDO
    338350
     351          CASE ( 'rho' )
     352             DO  i = nxl-1, nxr+1
     353                DO  j = nys-1, nyn+1
     354                   DO  k = nzb, nzt+1
     355                      rho_av(k,j,i) = rho_av(k,j,i) + rho(k,j,i)
     356                   ENDDO
     357                ENDDO
     358             ENDDO
     359             
    339360          CASE ( 's' )
    340361             DO  i = nxl-1, nxr+1
     
    342363                   DO  k = nzb, nzt+1
    343364                      s_av(k,j,i) = s_av(k,j,i) + q(k,j,i)
     365                   ENDDO
     366                ENDDO
     367             ENDDO
     368             
     369          CASE ( 'sa' )
     370             DO  i = nxl-1, nxr+1
     371                DO  j = nys-1, nyn+1
     372                   DO  k = nzb, nzt+1
     373                      sa_av(k,j,i) = sa_av(k,j,i) + sa(k,j,i)
    344374                   ENDDO
    345375                ENDDO
Note: See TracChangeset for help on using the changeset viewer.