Changeset 5974 for branches/2015/dev_r5072_UKMO2_OBS_simplification/NEMOGCM/NEMO/NST_SRC/agrif_opa_sponge.F90
- Timestamp:
- 2015-12-02T11:52:05+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5072_UKMO2_OBS_simplification/NEMOGCM/NEMO/NST_SRC/agrif_opa_sponge.F90
r5682 r5974 210 210 DO jj = j1,j2-1 211 211 DO ji = i1,i2-1 212 zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk)213 zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk)212 zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * e2_e1u(ji,jj) * fse3u_n(ji,jj,jk) 213 zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * e1_e2v(ji,jj) * fse3v_n(ji,jj,jk) 214 214 ztu(ji,jj,jk) = zabe1 * ( tsbdiff(ji+1,jj ,jk,jn) - tsbdiff(ji,jj,jk,jn) ) 215 215 ztv(ji,jj,jk) = zabe2 * ( tsbdiff(ji ,jj+1,jk,jn) - tsbdiff(ji,jj,jk,jn) ) … … 239 239 240 240 IF (.NOT. tabspongedone_tsn(ji,jj)) THEN 241 zbtr = r1_e1 2t(ji,jj) / fse3t_n(ji,jj,jk)241 zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) 242 242 ! horizontal diffusive trends 243 243 ztsa = zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji ,jj-1,jk) ) … … 290 290 DO jj = j1,j2 291 291 DO ji = i1+1,i2 ! vector opt. 292 zbtr = r1_e1 2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj)292 zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 293 293 hdivdiff(ji,jj,jk) = ( e2u(ji ,jj)*fse3u_n(ji ,jj,jk) * ubdiff(ji ,jj,jk) & 294 294 & -e2u(ji-1,jj)*fse3u_n(ji-1,jj,jk) * ubdiff(ji-1,jj,jk) ) * zbtr … … 298 298 DO jj = j1,j2-1 299 299 DO ji = i1,i2 ! vector opt. 300 zbtr = r1_e1 2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj)300 zbtr = r1_e1e2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 301 301 rotdiff(ji,jj,jk) = (-e1u(ji,jj+1) * ubdiff(ji,jj+1,jk) & 302 302 +e1u(ji,jj ) * ubdiff(ji,jj ,jk) & … … 396 396 DO jj = j1+1,j2 397 397 DO ji = i1,i2 ! vector opt. 398 zbtr = r1_e1 2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj)398 zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 399 399 hdivdiff(ji,jj,jk) = ( e1v(ji,jj ) * fse3v(ji,jj ,jk) * vbdiff(ji,jj ,jk) & 400 400 & -e1v(ji,jj-1) * fse3v(ji,jj-1,jk) * vbdiff(ji,jj-1,jk) ) * zbtr … … 403 403 DO jj = j1,j2 404 404 DO ji = i1,i2-1 ! vector opt. 405 zbtr = r1_e1 2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj)405 zbtr = r1_e1e2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 406 406 rotdiff(ji,jj,jk) = ( e2v(ji+1,jj) * vbdiff(ji+1,jj,jk) & 407 407 & -e2v(ji ,jj) * vbdiff(ji ,jj,jk) &
Note: See TracChangeset
for help on using the changeset viewer.