Changeset 3283 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM
- Timestamp:
- 2012-01-27T14:15:50+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/DIA/diafwb.F90
r3168 r3283 176 176 END SELECT 177 177 ! 178 DO ji = mi0(ii0), mi1(ii1)178 DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 179 179 DO jj = mj0(ij0), mj1(ij1) 180 180 DO jk = 1, jpk 181 181 zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 182 182 zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 183 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 183 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 184 184 185 185 IF( un(ji,jj,jk) > 0.e0 ) THEN … … 223 223 END SELECT 224 224 ! 225 DO ji = mi0(ii0), mi1(ii1)225 DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 226 226 DO jj = mj0(ij0), mj1(ij1) 227 227 DO jk = 1, jpk 228 228 zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 229 229 zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 230 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 230 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 231 231 232 232 IF( un(ji,jj,jk) > 0.e0 ) THEN … … 270 270 END SELECT 271 271 ! 272 DO ji = mi0(ii0), mi1(ii1)272 DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 273 273 DO jj = mj0(ij0), mj1(ij1) 274 274 DO jk = 1, jpk 275 275 zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 276 276 zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 277 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 277 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 278 278 279 279 IF( un(ji,jj,jk) > 0.e0 ) THEN … … 317 317 END SELECT 318 318 ! 319 DO ji = mi0(ii0), mi1(ii1)319 DO ji = mi0(ii0), MIN(mi1(ii1),jpim1) 320 320 DO jj = mj0(ij0), mj1(ij1) 321 321 DO jk = 1, jpk 322 322 zt = 0.5 * ( tsn(ji,jj,jk,jp_tem) + tsn(ji+1,jj,jk,jp_tem) ) 323 323 zs = 0.5 * ( tsn(ji,jj,jk,jp_sal) + tsn(ji+1,jj,jk,jp_sal) ) 324 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) 324 zu = un(ji,jj,jk) * fse3t(ji,jj,jk) * e2u(ji,jj) * tmask_i(ji,jj) 325 325 326 326 IF( un(ji,jj,jk) > 0.e0 ) THEN
Note: See TracChangeset
for help on using the changeset viewer.