Ignore:
Timestamp:
Sep 23, 2009 9:40:33 AM (12 years ago)
Author:
raasch
Message:

in-situ AND potential density are calculated and used in the ocean version

File:
1 edited

Legend:

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

    r240 r388  
    44! Actual revisions:
    55! -----------------
     6! prho is used instead of rho in diffusion_e,
    67! external pressure gradient
    78!
     
    330331          CALL coriolis( i, j, 3 )
    331332          IF ( ocean )  THEN
    332              CALL buoyancy( i, j, rho, prho_reference, 3, 64 )
     333             CALL buoyancy( i, j, rho, rho_reference, 3, 64 )
    333334          ELSE
    334335             IF ( .NOT. humidity )  THEN
     
    725726                IF ( .NOT. humidity )  THEN
    726727                   IF ( ocean )  THEN
    727                       CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km,   &
    728                                         l_grid, rho, prho_reference, rif, tend, &
    729                                         zu, zw )
     728                      CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km,  &
     729                                        l_grid, prho, prho_reference, rif,    &
     730                                        tend, zu, zw )
    730731                   ELSE
    731                       CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, &
    732                                         l_grid, pt, pt_reference, rif, tend,  &
     732                      CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km,  &
     733                                        l_grid, pt, pt_reference, rif, tend,   &
    733734                                        zu, zw )
    734735                   ENDIF
     
    769770                      IF ( ocean )  THEN
    770771                         CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e,  &
    771                                            km, l_grid, rho, prho_reference,   &
     772                                           km, l_grid, prho, prho_reference,  &
    772773                                           rif, tend, zu, zw )
    773774                      ELSE
     
    10131014          CALL coriolis( i, j, 3 )
    10141015          IF ( ocean )  THEN
    1015              CALL buoyancy( i, j, rho, prho_reference, 3, 64 )
     1016             CALL buoyancy( i, j, rho, rho_reference, 3, 64 )
    10161017          ELSE
    10171018             IF ( .NOT. humidity )  THEN
     
    12631264                IF ( .NOT. humidity )  THEN
    12641265                   IF ( ocean )  THEN
    1265                       CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, &
    1266                                         km, l_grid, rho, prho_reference,  &
     1266                      CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e,  &
     1267                                        km, l_grid, prho, prho_reference,  &
    12671268                                        rif, tend, zu, zw )
    12681269                   ELSE
     
    15521553    CALL coriolis( 3 )
    15531554    IF ( ocean )  THEN
    1554        CALL buoyancy( rho, prho_reference, 3, 64 )
     1555       CALL buoyancy( rho, rho_reference, 3, 64 )
    15551556    ELSE
    15561557       IF ( .NOT. humidity )  THEN
     
    19671968          IF ( .NOT. humidity )  THEN
    19681969             IF ( ocean )  THEN
    1969                 CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, rho, &
    1970                                   prho_reference, rif, tend, zu, zw )
     1970                CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, &
     1971                                  prho, prho_reference, rif, tend, zu, zw )
    19711972             ELSE
    19721973                CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, pt, &
     
    20042005                IF ( ocean )  THEN
    20052006                   CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, &
    2006                                      rho, prho_reference, rif, tend, zu, zw )
     2007                                     prho, prho_reference, rif, tend, zu, zw )
    20072008                ELSE
    20082009                   CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, &
Note: See TracChangeset for help on using the changeset viewer.