Changeset 389 for trunk/NEMO/OPA_SRC/SOL/solmat.F90
- Timestamp:
- 2006-03-09T18:22:04+01:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SOL/solmat.F90
r359 r389 182 182 183 183 !!cr ENDIF 184 #endif 185 #if defined key_AGRIF 186 IF (.NOT.AGRIF_ROOT()) THEN 187 188 IF ( (nbondi == -1) .OR. (nbondi == 2) ) bmask(2,:)=0. 189 IF ( (nbondi == 1) .OR. (nbondi == 2) ) bmask(nlci-1,:)=0. 190 IF ( (nbondj == -1) .OR. (nbondj == 2) ) bmask(:,2)=0. 191 IF ( (nbondj == 1) .OR. (nbondj == 2) ) bmask(:,nlcj-1)=0. 192 193 DO jj = 2, jpjm1 194 DO ji = 2, jpim1 195 zcoef = z2dt * z2dt * grav * rnu * bmask(ji,jj) 196 ! south coefficient 197 IF( ((nbondj == -1) .OR. (nbondj == 2)) .AND. ( jj == 3 ) ) THEN 198 zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1)*(1.-vmask(ji,jj-1,1)) 199 ELSE 200 zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1) 201 END IF 202 gcp(ji,jj,1) = zcoefs 203 204 ! west coefficient 205 IF( ( (nbondi == -1) .OR. (nbondi == 2) ) .AND. ( ji == 3 ) ) THEN 206 zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj)*(1.-umask(ji-1,jj,1)) 207 ELSE 208 zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj) 209 END IF 210 gcp(ji,jj,2) = zcoefw 211 212 ! east coefficient 213 IF( ((nbondi == 1) .OR. (nbondi == 2)) .AND. ( ji == nlci-2 ) ) THEN 214 zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj)*(1.-umask(ji,jj,1)) 215 ELSE 216 zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj) 217 END IF 218 gcp(ji,jj,3) = zcoefe 219 220 ! north coefficient 221 IF( ((nbondj == 1) .OR. (nbondj == 2)) .AND. ( jj == nlcj-2 ) ) THEN 222 zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj)*(1.-vmask(ji,jj,1)) 223 ELSE 224 zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj) 225 END IF 226 gcp(ji,jj,4) = zcoefn 227 228 ! diagonal coefficient 229 gcdmat(ji,jj) = e1t(ji,jj)*e2t(ji,jj)*bmask(ji,jj) & 230 - zcoefs -zcoefw -zcoefe -zcoefn 231 END DO 232 END DO 233 234 ENDIF 184 235 #endif 185 236
Note: See TracChangeset
for help on using the changeset viewer.