Changeset 9783 for NEMO/releases/release-3.6/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

Timestamp:

2018-06-12T17:36:37+02:00 (5 years ago)

Author:

mathiot

Message:

correcting bug introduced this morning (tmask(ji,jj) instead of tmask(ji,jj,1))

File:

• NEMO/releases/release-3.6/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90 (modified) (1 diff)

NEMO/releases/release-3.6/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

@@ -281 +281 @@
             ! ice area at u and v-points
             ! land values on u-, v- points along coastline set to adjacent t-point ocean value 
-            zat_u = ( at_i(ji,jj) * tmask(ji,jj) + at_i(ji+1,jj  ) * tmask(ji+1,jj  ) ) / MAX(1.0_wp,tmask(ji,jj)+tmask(ji+1,jj  ))
-            zat_v = ( at_i(ji,jj) * tmask(ji,jj) + at_i(ji  ,jj+1) * tmask(ji  ,jj+1) ) / MAX(1.0_wp,tmask(ji,jj)+tmask(ji  ,jj+1))
+            zat_u = ( at_i(ji,jj) * tmask(ji,jj,1) + at_i (ji+1,jj) * tmask(ji+1,jj,1) ) &
+               &    / MAX( 1.0_wp , tmask(ji,jj,1) + tmask(ji+1,jj,1))
+            zat_v = ( at_i(ji,jj) * tmask(ji,jj,1) + at_i (ji,jj+1) * tmask(ji,jj+1,1) ) &
+               &    / MAX( 1.0_wp , tmask(ji,jj,1) + tmask(ji,jj+1,1))
             !                                                   ! linearized quadratic drag formulation
             zutau_ice   = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji+1,jj) ) * ( u_ice(ji,jj) - pu_oce(ji,jj) )