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/init_3d_model.f90

    r359 r388  
    77! Current revisions:
    88! -----------------
     9! Initialization of prho added.
    910! bugfix: correction of initial volume flow for non-flat topography
    1011! bugfix: zero initialization of arrays within buildings for 'cyclic_fill'
     
    234235       ALLOCATE( saswsb_1(nys-1:nyn+1,nxl-1:nxr+1), &
    235236                 saswst_1(nys-1:nyn+1,nxl-1:nxr+1) )
    236        ALLOCATE( rho_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
    237                  sa_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
    238                  sa_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
     237       ALLOCATE( prho_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     238                 rho_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
     239                 sa_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),   &
     240                 sa_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),   &
    239241                 sa_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    240        rho => rho_1  ! routine calc_mean_profile requires density to be a
    241                      ! pointer
     242       prho => prho_1
     243       rho  => rho_1  ! routines calc_mean_profile and diffusion_e require
     244                      ! density to be apointer
    242245       IF ( humidity_remote )  THEN
    243246          ALLOCATE( qswst_remote(nys-1:nyn+1,nxl-1:nxr+1) )
     
    13281331!--    Initialize quantities needed for the ocean model
    13291332       CALL init_ocean
     1333
    13301334    ELSE
    13311335!
Note: See TracChangeset for help on using the changeset viewer.