Changeset 1694 for trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
- Timestamp:
- 2009-10-30T17:06:23+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
r1693 r1694 796 796 CASE( 'U' ) 797 797 DO jj = 2, jpjm1 ! (U,V) ==> I 798 DO ji = fs_2, fs_jpim1 !vector opt.798 DO ji = 2, jpim1 ! NO vector opt. 799 799 p_taui(ji,jj) = 0.5 * ( frcv(ji-1,jj ,jpr_itx1) + frcv(ji-1,jj-1,jpr_itx1) ) 800 800 p_tauj(ji,jj) = 0.5 * ( frcv(ji ,jj-1,jpr_ity1) + frcv(ji-1,jj-1,jpr_ity1) ) … … 803 803 CASE( 'F' ) 804 804 DO jj = 2, jpjm1 ! F ==> I 805 DO ji = fs_2, fs_jpim1 !vector opt.805 DO ji = 2, jpim1 ! NO vector opt. 806 806 p_taui(ji,jj) = frcv(ji-1,jj-1,jpr_itx1) 807 807 p_tauj(ji,jj) = frcv(ji-1,jj-1,jpr_ity1) … … 810 810 CASE( 'T' ) 811 811 DO jj = 2, jpjm1 ! T ==> I 812 DO ji = fs_2, fs_jpim1 !vector opt.812 DO ji = 2, jpim1 ! NO vector opt. 813 813 p_taui(ji,jj) = 0.25 * ( frcv(ji,jj ,jpr_itx1) + frcv(ji-1,jj ,jpr_itx1) & 814 814 & + frcv(ji,jj-1,jpr_itx1) + frcv(ji-1,jj-1,jpr_itx1) ) … … 836 836 CASE( 'I' ) 837 837 DO jj = 2, jpjm1 ! I ==> F 838 DO ji = fs_2, fs_jpim1 !vector opt.838 DO ji = 2, jpim1 ! NO vector opt. 839 839 p_taui(ji,jj) = frcv(ji+1,jj+1,jpr_itx1) 840 840 p_tauj(ji,jj) = frcv(ji+1,jj+1,jpr_ity1) … … 843 843 CASE( 'T' ) 844 844 DO jj = 2, jpjm1 ! T ==> F 845 DO ji = fs_2, fs_jpim1 !vector opt.845 DO ji = 2, jpim1 ! NO vector opt. 846 846 p_taui(ji,jj) = 0.25 * ( frcv(ji,jj ,jpr_itx1) + frcv(ji+1,jj ,jpr_itx1) & 847 847 & + frcv(ji,jj+1,jpr_itx1) + frcv(ji+1,jj+1,jpr_itx1) ) … … 879 879 CASE( 'I' ) 880 880 DO jj = 2, jpjm1 ! I ==> (U,V) 881 DO ji = fs_2, fs_jpim1 !vector opt.881 DO ji = 2, jpim1 ! NO vector opt. 882 882 p_taui(ji,jj) = 0.5 * ( frcv(ji+1,jj+1,jpr_itx1) + frcv(ji+1,jj ,jpr_itx1) ) 883 883 p_tauj(ji,jj) = 0.5 * ( frcv(ji+1,jj+1,jpr_ity1) + frcv(ji ,jj+1,jpr_ity1) ) … … 1134 1134 ! 1135 1135 ! j+1 j -----V---F 1136 ! surfac zvelocity always sent from T point ! |1136 ! surface velocity always sent from T point ! | 1137 1137 ! j | T U 1138 1138 ! | | … … 1162 1162 CASE( 'I' ) ! Ocean on C grid, Ice on I-point (B-grid) ==> T 1163 1163 DO jj = 2, jpjm1 1164 DO ji = fs_2, fs_jpim1 !vector opt.1164 DO ji = 2, jpim1 ! NO vector opt. 1165 1165 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj ,1) ) * zfr_l(ji,jj) 1166 1166 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj) … … 1173 1173 CASE( 'F' ) ! Ocean on C grid, Ice on F-point (B-grid) ==> T 1174 1174 DO jj = 2, jpjm1 1175 DO ji = fs_2, fs_jpim1 !vector opt.1175 DO ji = 2, jpim1 ! NO vector opt. 1176 1176 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj ,1) ) * zfr_l(ji,jj) 1177 1177 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj) … … 1197 1197 CASE( 'I' ) ! Ocean on C grid, Ice on I-point (B-grid) ==> T 1198 1198 DO jj = 2, jpjm1 1199 DO ji = fs_2, fs_jpim1 !vector opt.1199 DO ji = 2, jpim1 ! NO vector opt. 1200 1200 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj ,1) ) * zfr_l(ji,jj) & 1201 1201 & + 0.25 * ( u_ice(ji+1,jj+1) + u_ice(ji,jj+1) & … … 1208 1208 CASE( 'F' ) ! Ocean on C grid, Ice on F-point (B-grid) ==> T 1209 1209 DO jj = 2, jpjm1 1210 DO ji = fs_2, fs_jpim1 !vector opt.1210 DO ji = 2, jpim1 ! NO vector opt. 1211 1211 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj ,1) ) * zfr_l(ji,jj) & 1212 1212 & + 0.25 * ( u_ice(ji-1,jj-1) + u_ice(ji,jj-1) &
Note: See TracChangeset
for help on using the changeset viewer.