Changeset 10425 for NEMO/trunk/src/OCE/DOM/closea.F90
- Timestamp:
- 2018-12-19T22:54:16+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DOM/closea.F90
r10069 r10425 118 118 ! number of closed seas = global maximum value in closea_mask field 119 119 jncs = maxval(closea_mask(:,:)) 120 IF( lk_mpp ) CALL mpp_max(jncs)120 CALL mpp_max('closea', jncs) 121 121 IF( jncs > 0 ) THEN 122 122 IF( lwp ) WRITE(numout,*) 'Number of closed seas : ',jncs … … 146 146 ! number of closed seas rnf mappings = global maximum in closea_mask_rnf field 147 147 jncsr = maxval(closea_mask_rnf(:,:)) 148 IF( lk_mpp ) CALL mpp_max(jncsr)148 CALL mpp_max('closea', jncsr) 149 149 IF( jncsr > 0 ) THEN 150 150 IF( lwp ) WRITE(numout,*) 'Number of closed seas rnf mappings : ',jncsr … … 166 166 ! number of closed seas empmr mappings = global maximum value in closea_mask_empmr field 167 167 jncse = maxval(closea_mask_empmr(:,:)) 168 IF( lk_mpp ) CALL mpp_max(jncse)168 CALL mpp_max('closea', jncse) 169 169 IF( jncse > 0 ) THEN 170 170 IF( lwp ) WRITE(numout,*) 'Number of closed seas empmr mappings : ',jncse … … 237 237 surfe(:) = 0.e0_wp 238 238 ! 239 surf(jncs+1) = glob_sum( e1e2t(:,:) ) ! surface of the global ocean239 surf(jncs+1) = glob_sum( 'closea', e1e2t(:,:) ) ! surface of the global ocean 240 240 ! 241 241 ! ! surface areas of closed seas … … 243 243 ztmp2d(:,:) = 0.e0_wp 244 244 WHERE( closea_mask(:,:) == jc ) ztmp2d(:,:) = e1e2t(:,:) * tmask_i(:,:) 245 surf(jc) = glob_sum( ztmp2d(:,:) )245 surf(jc) = glob_sum( 'closea', ztmp2d(:,:) ) 246 246 END DO 247 247 ! … … 254 254 ztmp2d(:,:) = 0.e0_wp 255 255 WHERE( closea_mask_rnf(:,:) == jcr .and. closea_mask(:,:) == 0 ) ztmp2d(:,:) = e1e2t(:,:) * tmask_i(:,:) 256 surfr(jcr) = glob_sum( ztmp2d(:,:) )256 surfr(jcr) = glob_sum( 'closea', ztmp2d(:,:) ) 257 257 END DO 258 258 ENDIF … … 263 263 ztmp2d(:,:) = 0.e0_wp 264 264 WHERE( closea_mask_empmr(:,:) == jce .and. closea_mask(:,:) == 0 ) ztmp2d(:,:) = e1e2t(:,:) * tmask_i(:,:) 265 surfe(jce) = glob_sum( ztmp2d(:,:) )265 surfe(jce) = glob_sum( 'closea', ztmp2d(:,:) ) 266 266 END DO 267 267 ENDIF … … 301 301 ztmp2d(:,:) = 0.e0_wp 302 302 WHERE( closea_mask(:,:) == jc ) ztmp2d(:,:) = e1e2t(:,:) * ( emp(:,:)-rnf(:,:) ) * tmask_i(:,:) 303 zfwf(jc) = glob_sum( ztmp2d(:,:) )303 zfwf(jc) = glob_sum( 'closea', ztmp2d(:,:) ) 304 304 END DO 305 305 zfwf_total = SUM(zfwf) … … 316 316 ztmp2d(:,:) = 0.e0_wp 317 317 WHERE( closea_mask_rnf(:,:) == jcr .and. closea_mask(:,:) > 0 ) ztmp2d(:,:) = e1e2t(:,:) * ( emp(:,:)-rnf(:,:) ) * tmask_i(:,:) 318 zfwfr(jcr) = glob_sum( ztmp2d(:,:) )318 zfwfr(jcr) = glob_sum( 'closea', ztmp2d(:,:) ) 319 319 ! 320 320 ! The following if avoids the redistribution of the round off … … 345 345 ztmp2d(:,:) = 0.e0_wp 346 346 WHERE( closea_mask_empmr(:,:) == jce .and. closea_mask(:,:) > 0 ) ztmp2d(:,:) = e1e2t(:,:) * ( emp(:,:)-rnf(:,:) ) * tmask_i(:,:) 347 zfwfe(jce) = glob_sum( ztmp2d(:,:) )347 zfwfe(jce) = glob_sum( 'closea', ztmp2d(:,:) ) 348 348 ! 349 349 ! The following if avoids the redistribution of the round off … … 396 396 emp (:,:) = emp (:,:) * tmask(:,:,1) 397 397 ! 398 CALL lbc_lnk( emp , 'T', 1._wp )398 CALL lbc_lnk( 'closea', emp , 'T', 1._wp ) 399 399 ! 400 400 END SUBROUTINE sbc_clo
Note: See TracChangeset
for help on using the changeset viewer.