Changeset 2126


Ignore:
Timestamp:
Jan 20, 2017 3:54:21 PM (5 years ago)
Author:
raasch
Message:

for ocean simulations, density in SGS-TKE-production is replaced by potential density

File:
1 edited

Legend:

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

    r2119 r2126  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! density in ocean case replaced by potential density
    2323!
    2424! Former revisions:
     
    138138
    139139       USE arrays_3d,                                                          &
    140            ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho_ocean, shf,       &
     140           ONLY:  ddzw, dd2zu, kh, km, prho, pt, q, ql, qsws, qswst, shf,      &
    141141                  tend, tswst, u, v, vpt, w
    142142
     
    537537                      DO  j = nys, nyn
    538538                         DO  k = nzb_s_inner(j,i)+1, nzt
    539                             tend(k,j,i) = tend(k,j,i) +                     &
    540                                           kh(k,j,i) * g / rho_reference *   &
    541                                           ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
     539                            tend(k,j,i) = tend(k,j,i) +                        &
     540                                          kh(k,j,i) * g / rho_reference *      &
     541                                          ( prho(k+1,j,i) - prho(k-1,j,i) ) * &
    542542                                          dd2zu(k)
    543543                         ENDDO
     
    577577                      DO  j = nys, nyn
    578578                         DO  k = nzb_s_inner(j,i)+1, nzt
    579                             tend(k,j,i) = tend(k,j,i) +                     &
    580                                           kh(k,j,i) * g / rho_ocean(k,j,i) *      &
    581                                           ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
     579                            tend(k,j,i) = tend(k,j,i) +                        &
     580                                          kh(k,j,i) * g / prho(k,j,i) *        &
     581                                          ( prho(k+1,j,i) - prho(k-1,j,i) ) * &
    582582                                          dd2zu(k)
    583583                         ENDDO
     
    744744
    745745       USE arrays_3d,                                                          &
    746            ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho_ocean, shf,       &
     746           ONLY:  ddzw, dd2zu, kh, km, prho, pt, q, ql, qsws, qswst, shf,      &
    747747                  tend, tswst, u, v, vpt, w
    748748
     
    11001100                      tend(k,j,i) = tend(k,j,i) +                   &
    11011101                                    kh(k,j,i) * g / rho_reference * &
    1102                                     ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * dd2zu(k)
     1102                                    ( prho(k+1,j,i) - prho(k-1,j,i) ) * dd2zu(k)
    11031103                   ENDDO
    11041104
     
    11301130!--                the bottom and top surface layer
    11311131                   DO  k = nzb_s_inner(j,i)+1, nzt
    1132                       tend(k,j,i) = tend(k,j,i) +                &
    1133                                     kh(k,j,i) * g / rho_ocean(k,j,i) * &
    1134                                     ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * dd2zu(k)
     1132                      tend(k,j,i) = tend(k,j,i) +                              &
     1133                                    kh(k,j,i) * g / prho(k,j,i) *              &
     1134                                    ( prho(k+1,j,i) - prho(k-1,j,i) ) * dd2zu(k)
    11351135                   ENDDO
    11361136
Note: See TracChangeset for help on using the changeset viewer.