- Timestamp:
- 2016-12-01T18:46:13+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_closea.F90
r7200 r7426 214 214 ! 215 215 ! ! surface of closed seas 216 IF( lk_mpp_rep ) THEN ! MPP reproductible calculation 217 DO jc = 1, jpncs 218 ctmp = CMPLX( 0.e0, 0.e0, wp ) 219 DO jj = ncsj1(jc), ncsj2(jc) 220 DO ji = ncsi1(jc), ncsi2(jc) 221 ztmp = e1e2t(ji,jj) * tmask_i(ji,jj) 222 CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 223 END DO 216 DO jc = 1, jpncs 217 ctmp = CMPLX( 0.e0, 0.e0, wp ) 218 DO jj = ncsj1(jc), ncsj2(jc) 219 DO ji = ncsi1(jc), ncsi2(jc) 220 ztmp = e1e2t(ji,jj) * tmask_i(ji,jj) 221 CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 224 222 END DO 225 IF( lk_mpp ) CALL mpp_sum( ctmp ) 226 surf(jc) = REAL(ctmp,wp) 227 END DO 228 ELSE ! Standard calculation 229 DO jc = 1, jpncs 230 DO jj = ncsj1(jc), ncsj2(jc) 231 DO ji = ncsi1(jc), ncsi2(jc) 232 surf(jc) = surf(jc) + e1e2t(ji,jj) * tmask_i(ji,jj) ! surface of closed seas 233 END DO 234 END DO 235 END DO 236 IF( lk_mpp ) CALL mpp_sum ( surf, jpncs ) ! mpp: sum over all the global domain 237 ENDIF 223 END DO 224 IF( lk_mpp ) CALL mpp_sum( ctmp ) 225 surf(jc) = REAL(ctmp,wp) 226 END DO 238 227 239 228 IF(lwp) WRITE(numout,*)' Closed sea surfaces'
Note: See TracChangeset
for help on using the changeset viewer.