- Timestamp:
- 2014-07-23T18:27:21+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO2_ice_shelves/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90
r4666 r4726 186 186 187 187 ! (ISF) define barotropic mask and mask the ice shelf point 188 lmask(:,:)=tmask(:,:,1) ! at this stage ice shelf is not masked188 ssmask(:,:)=tmask(:,:,1) ! at this stage ice shelf is not masked 189 189 190 190 DO jk = 1, jpk 191 191 DO jj = 1, jpj 192 192 DO ji = 1, jpi 193 IF( REAL( mi cedep(ji,jj) - jk, wp ) - 0.1_wp >= 0._wp ) THEN193 IF( REAL( misfdep(ji,jj) - jk, wp ) - 0.1_wp >= 0._wp ) THEN 194 194 tmask(ji,jj,jk) = 0._wp 195 195 END IF … … 220 220 ! Interior domain mask (used for global sum) 221 221 ! -------------------- 222 tmask_i(:,:) = lmask(:,:) ! (ISH) tmask_i = 1 even on the ice shelf222 tmask_i(:,:) = ssmask(:,:) ! (ISH) tmask_i = 1 even on the ice shelf 223 223 iif = jpreci ! ??? 224 224 iil = nlci - jpreci + 1 … … 266 266 DO jj = 1, jpjm1 267 267 DO ji = 1, fs_jpim1 ! vector loop 268 umask_i(ji,jj) = lmask(ji,jj) * lmask(ji+1,jj ) * MIN(1._wp,SUM(umask(ji,jj,:)))269 vmask_i(ji,jj) = lmask(ji,jj) * lmask(ji ,jj+1) * MIN(1._wp,SUM(vmask(ji,jj,:)))268 umask_i(ji,jj) = ssmask(ji,jj) * ssmask(ji+1,jj ) * MIN(1._wp,SUM(umask(ji,jj,:))) 269 vmask_i(ji,jj) = ssmask(ji,jj) * ssmask(ji ,jj+1) * MIN(1._wp,SUM(vmask(ji,jj,:))) 270 270 END DO 271 271 DO ji = 1, jpim1 ! NO vector opt. 272 fmask_i(ji,jj) = lmask(ji,jj ) * lmask(ji+1,jj ) &273 & * lmask(ji,jj+1) * lmask(ji+1,jj+1) * MIN(1._wp,SUM(fmask(ji,jj,:)))272 fmask_i(ji,jj) = ssmask(ji,jj ) * ssmask(ji+1,jj ) & 273 & * ssmask(ji,jj+1) * ssmask(ji+1,jj+1) * MIN(1._wp,SUM(fmask(ji,jj,:))) 274 274 END DO 275 275 END DO … … 284 284 ! 4. ocean/land mask for the elliptic equation 285 285 ! -------------------------------------------- 286 bmask(:,:) = lmask(:,:) ! elliptic equation is written at t-point286 bmask(:,:) = ssmask(:,:) ! elliptic equation is written at t-point 287 287 ! 288 288 ! ! Boundary conditions
Note: See TracChangeset
for help on using the changeset viewer.