Ignore:
Timestamp:
Jun 2, 2007 4:48:38 PM (17 years ago)
Author:
raasch
Message:

further preliminary uncomplete changes for ocean version

File:
1 edited

Legend:

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

    r77 r95  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Boundary conditions for salinity added
    77!
    88! Former revisions:
     
    128128
    129129!
     130!--    Boundary conditions for salinity
     131       IF ( ocean )  THEN
     132!
     133!--       Bottom boundary: Neumann condition because salinity flux is always
     134!--       given
     135          DO  i = nxl-1, nxr+1
     136             DO  j = nys-1, nyn+1
     137                sa_p(nzb_s_inner(j,i),j,i) = sa_p(nzb_s_inner(j,i)+1,j,i)
     138             ENDDO
     139          ENDDO
     140
     141!
     142!--       Top boundary: Dirichlet or Neumann
     143          IF ( ibc_sa_t == 0 )  THEN
     144             sa_p(nzt+1,:,:) = sa(nzt+1,:,:)
     145          ELSEIF ( ibc_sa_t == 1 )  THEN
     146             sa_p(nzt+1,:,:) = sa_p(nzt,:,:)
     147          ENDIF
     148
     149       ENDIF
     150
     151!
    130152!--    Boundary conditions for total water content or scalar,
    131 !--    bottom and surface boundary (see also temperature)
     153!--    bottom and top boundary (see also temperature)
    132154       IF ( humidity  .OR.  passive_scalar )  THEN
    133155!
Note: See TracChangeset for help on using the changeset viewer.