Changeset 95 for palm/trunk/SOURCE/boundary_conds.f90
- Timestamp:
- Jun 2, 2007 4:48:38 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/boundary_conds.f90
r77 r95 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! Boundary conditions for salinity added 7 7 ! 8 8 ! Former revisions: … … 128 128 129 129 ! 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 ! 130 152 !-- Boundary conditions for total water content or scalar, 131 !-- bottom and surfaceboundary (see also temperature)153 !-- bottom and top boundary (see also temperature) 132 154 IF ( humidity .OR. passive_scalar ) THEN 133 155 !
Note: See TracChangeset
for help on using the changeset viewer.