Ignore:
Timestamp:
Mar 9, 2007 12:05:41 PM (15 years ago)
Author:
raasch
Message:

preliminary update of further changes, advec_particles is not running!

File:
1 edited

Legend:

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

    r4 r57  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Reference temperature pt_reference can be used.
    77!
    88! Former revisions:
     
    5959!
    6060!--       Normal case: horizontal surface
    61           DO  i = nxl, nxr
    62              DO  j = nys, nyn
    63                 DO  k = nzb_s_inner(j,i)+1, nzt-1
    64                     tend(k,j,i) = tend(k,j,i) + g * 0.5 * (                    &
     61          IF ( use_pt_reference )  THEN
     62             DO  i = nxl, nxr
     63                DO  j = nys, nyn
     64                   DO  k = nzb_s_inner(j,i)+1, nzt-1
     65                      tend(k,j,i) = tend(k,j,i) + g * 0.5 * (                 &
     66                        ( theta(k,j,i)   - hom(k,1,pr,0)   ) / pt_reference + &
     67                        ( theta(k+1,j,i) - hom(k+1,1,pr,0) ) / pt_reference   &
     68                                                            )
     69                   ENDDO
     70                ENDDO
     71             ENDDO
     72          ELSE
     73             DO  i = nxl, nxr
     74                DO  j = nys, nyn
     75                   DO  k = nzb_s_inner(j,i)+1, nzt-1
     76                      tend(k,j,i) = tend(k,j,i) + g * 0.5 * (                  &
    6577                        ( theta(k,j,i)   - hom(k,1,pr,0)   ) / hom(k,1,pr,0) + &
    6678                        ( theta(k+1,j,i) - hom(k+1,1,pr,0) ) / hom(k+1,1,pr,0) &
    67                                                           )
    68                 ENDDO
    69              ENDDO
    70           ENDDO
     79                                                            )
     80                   ENDDO
     81                ENDDO
     82             ENDDO
     83          ENDIF
    7184
    7285       ELSE
     
    136149!
    137150!--       Normal case: horizontal surface
    138           DO  k = nzb_s_inner(j,i)+1, nzt-1
    139               tend(k,j,i) = tend(k,j,i) + g * 0.5 * (                          &
     151          IF ( use_pt_reference )  THEN
     152             DO  k = nzb_s_inner(j,i)+1, nzt-1
     153                 tend(k,j,i) = tend(k,j,i) + g * 0.5 * (                      &
     154                        ( theta(k,j,i)   - hom(k,1,pr,0)   ) / pt_reference + &
     155                        ( theta(k+1,j,i) - hom(k+1,1,pr,0) ) / pt_reference   &
     156                                                       )
     157             ENDDO
     158          ELSE
     159             DO  k = nzb_s_inner(j,i)+1, nzt-1
     160                 tend(k,j,i) = tend(k,j,i) + g * 0.5 * (                       &
    140161                        ( theta(k,j,i)   - hom(k,1,pr,0)   ) / hom(k,1,pr,0) + &
    141162                        ( theta(k+1,j,i) - hom(k+1,1,pr,0) ) / hom(k+1,1,pr,0) &
    142                                                     )
    143           ENDDO
     163                                                       )
     164             ENDDO
     165          ENDIF
    144166
    145167       ELSE
Note: See TracChangeset for help on using the changeset viewer.