Changeset 2207 for branches/DEV_r2191_3partymerge2010/NEMO/OPA_SRC/SOL
- Timestamp:
- 2010-10-11T18:51:54+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2191_3partymerge2010/NEMO/OPA_SRC/SOL/solmat.F90
r2031 r2207 15 15 !! 3.2 ! 2009-06 (S. Masson) distributed restart using iom 16 16 !! - ! 2009-07 (R. Benshila) suppression of rigid-lid option 17 !! 3.3 ! 2010-09 (D. Storkey) update for BDY module. 17 18 !!---------------------------------------------------------------------- 18 19 … … 26 27 USE phycst ! physical constants 27 28 USE obc_oce ! ocean open boundary conditions 29 USE bdy_oce ! unstructured open boundary conditions 28 30 USE lbclnk ! lateral boudary conditions 29 31 USE lib_mpp ! distributed memory computing … … 80 82 ENDIF 81 83 82 #if defined key_dynspg_flt 84 #if defined key_dynspg_flt && ! defined key_bdy 83 85 # if ! defined key_obc 84 86 … … 158 160 ENDIF 159 161 # endif 162 163 # elif defined key_dynspg_flt && defined key_bdy 164 165 ! defined gcdmat in the case of unstructured open boundaries 166 DO jj = 2, jpjm1 167 DO ji = 2, jpim1 168 zcoef = z2dt * z2dt * grav * bmask(ji,jj) 169 170 ! south coefficient 171 zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1) 172 zcoefs = zcoefs * bdyvmask(ji,jj-1) 173 gcp(ji,jj,1) = zcoefs 174 175 ! west coefficient 176 zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj) 177 zcoefw = zcoefw * bdyumask(ji-1,jj) 178 gcp(ji,jj,2) = zcoefw 179 180 ! east coefficient 181 zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj) 182 zcoefe = zcoefe * bdyumask(ji,jj) 183 gcp(ji,jj,3) = zcoefe 184 185 ! north coefficient 186 zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj) 187 zcoefn = zcoefn * bdyvmask(ji,jj) 188 gcp(ji,jj,4) = zcoefn 189 190 ! diagonal coefficient 191 gcdmat(ji,jj) = e1t(ji,jj)*e2t(ji,jj)*bmask(ji,jj) & 192 - zcoefs -zcoefw -zcoefe -zcoefn 193 END DO 194 END DO 195 160 196 #endif 161 197
Note: See TracChangeset
for help on using the changeset viewer.