Changeset 4447 for branches/2011/DEV_r2739_STFC_dCSE
- Timestamp:
- 2014-02-04T13:15:59+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_iso.F90
r3432 r4447 169 169 DO jj = 2, jpjm1 ! set the slopes of iso-level 170 170 DO ji = fs_2, fs_jpim1 171 DO jk = 1, jpk171 DO jk = 1, mbkmax(ji,jj) ! jpk 172 172 #else 173 173 DO jk = 1, jpk ! set the slopes of iso-level … … 368 368 ! when jk is 1 but that doesn't matter. 369 369 !DIR$ SAFE_ADDRESS 370 DO jk = 1, jpkm1, 1370 DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1, 1 371 371 372 372 ! Vertical u- and v-shears at level jk and jk+1 … … 465 465 ! when jk is 1 but that doesn't matter. 466 466 !DIR$ SAFE_ADDRESS 467 DO jk = 1, jpkm1, 1467 DO jk = 1, mbkmax(ji,jj)-1, 1 ! jpkm1, 1 468 468 469 469 zdk1v = ( vb(ji ,jj ,jk) -vb(ji ,jj ,jk+1) ) * vmask(ji ,jj,jk+1) … … 566 566 ! when jk is 1 but that doesn't matter. 567 567 !DIR$ SAFE_ADDRESS 568 DO jk = 1, jpkm1, 1568 DO jk = 1, mbkmax(ji,jj)-1, 1 ! jpkm1, 1 569 569 570 570 ! Vertical u- and v-shears at level jk and jk+1 … … 661 661 ! when jk is 1 but that doesn't matter. 662 662 !DIR$ SAFE_ADDRESS 663 DO jk = 1, jpkm1, 1663 DO jk = 1, mbkmax(ji,jj)-1, 1 ! jpkm1, 1 664 664 665 665 zdk1v = ( vb(ji ,jj ,jk) -vb(ji ,jj ,jk+1) ) * vmask(ji ,jj,jk+1) … … 918 918 919 919 DO ji = 2, jpi 920 DO jk = 1, jpk920 DO jk = 1, mbkmax(ji,jj) ! jpk 921 921 ! i-gradient of u at jj 922 922 zdiu (ji,jk) = tmask(ji,jj ,jk) * ( ub(ji,jj ,jk) - ub(ji-1,jj ,jk) ) … … 931 931 932 932 DO ji = 1, jpim1 933 DO jk = 1, jpk933 DO jk = 1, mbkmax(ji,jj) ! jpk 934 934 ! i-gradient of v at jj 935 935 zdiv (ji,jk) = fmask(ji,jj ,jk) * ( vb(ji+1,jj,jk) - vb(ji ,jj ,jk) ) … … 951 951 ! interior (2=<jk=<jpk-1) on U field 952 952 DO ji = 2, jpim1 953 DO jk = 2, jpkm1953 DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 954 954 zcoef0= 0.5 * aht0 * umask(ji,jj,jk) 955 955 … … 977 977 ! interior (2=<jk=<jpk-1) on V field 978 978 DO ji = 2, jpim1 979 DO jk = 2, jpkm1979 DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 980 980 zcoef0= 0.5 * aht0 * vmask(ji,jj,jk) 981 981 … … 1005 1005 ! ------------------------------------------------------------------- 1006 1006 ! DO ji = 2, jpim1 1007 zfuw(ji,1) = 0.0_wp1008 zfuw(ji, jpk) = 0.0_wp1009 zfvw(ji,1) = 0.0_wp1010 zfvw(ji, jpk) = 0.0_wp1011 1012 DO jk = 1, jpkm11007 zfuw(ji,1) = 0.0_wp 1008 zfuw(ji,mbkmax(ji,jj)) = 0.0_wp 1009 zfvw(ji,1) = 0.0_wp 1010 zfvw(ji,mbkmax(ji,jj)) = 0.0_wp 1011 1012 DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1 1013 1013 ! volume elements 1014 1014 zbu = e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.