Changeset 868 for trunk/NEMO/LIM_SRC_3/limdyn.F90
- Timestamp:
- 2008-03-14T19:53:00+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/limdyn.F90
r867 r868 33 33 PUBLIC lim_dyn ! routine called by ice_step 34 34 35 !! * Substitutions 36 # include "vectopt_loop_substitute.h90" 37 35 38 !! * Module variables 36 39 REAL(wp) :: rone = 1.e0 ! constant value … … 172 175 zsang = SIGN(1.e0, gphif(1,jj-1) ) * sangvg 173 176 174 DO ji = 2,jpim1177 DO ji = fs_2, fs_jpim1 175 178 ! computation of wind stress over ocean in X and Y direction 176 179 #if defined key_coupled && defined key_lim_cp1 … … 216 219 ! computation of friction velocity 217 220 DO jj = 2, jpjm1 218 DO ji = 2,jpim1221 DO ji = fs_2, fs_jpim1 219 222 220 223 zu_ice = u_ice(ji,jj) - u_io(ji,jj) … … 245 248 ! virer ca (key_lim_cp1) 246 249 DO jj = 2, jpjm1 247 DO ji = 2,jpim1250 DO ji = fs_2, fs_jpim1 248 251 #if defined key_coupled && defined key_lim_cp1 249 252 tio_u(ji,jj) = - ( gtaux(ji ,jj ) + gtaux(ji-1,jj ) &
Note: See TracChangeset
for help on using the changeset viewer.