Changeset 9772 for NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/src/OCE/nemogcm.F90

Timestamp:

2018-06-08T14:52:48+02:00 (6 years ago)

Author:

smasson

Message:

dev_r9759_HPC09_ESIWACE: add benchmark features

File:

• NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/src/OCE/nemogcm.F90 (modified) (4 diffs)

NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/src/OCE/nemogcm.F90

@@ -99 +99 @@
    CHARACTER(lc) ::   cform_aaa="( /, 'AAAAAAAA', / ) "     ! flag for output listing
 
+#if defined key_mpp_mpi
+   INCLUDE 'mpif.h'
+#endif
+
    !!----------------------------------------------------------------------
    !! NEMO/OCE 4.0 , NEMO Consortium (2018)
@@ -121 +125 @@
       !!----------------------------------------------------------------------
       INTEGER ::   istp   ! time step index
+      INTEGER ::   i, ik, id, iend, iseq   ! time step index
+      REAL(wp) ::   tic  
+
       !!----------------------------------------------------------------------
       !
@@ -178 +185 @@
          !
          DO WHILE( istp <= nitend .AND. nstop == 0 )
-            CALL stp        ( istp ) 
+#if   defined key_mpp_mpi
+            ncom_stp = istp
+#endif
+            IF ( istp == ( nit000 + 1 ) ) tic = MPI_Wtime()
+            IF ( ln_comm_only .AND. istp >= ( nit000  + 1 )) THEN
+               DO i = 1, n_sequence
+                 IF ( icomm_sequence(i,1) == 1 ) THEN
+                    SELECT CASE ( icomm_sequence(i,2) )
+                      CASE (1)
+                      CALL lbc_lnk_multi( un(:,:,1)  , 'U', -1. )
+                      CASE (2)
+                      CALL lbc_lnk_multi( un(:,:,1)  , 'U', -1., vn(:,:,1)  , 'V', -1.)
+                      CASE (3)
+                      CALL lbc_lnk_multi( un(:,:,1)  , 'U', -1., vn(:,:,1)  , 'V', -1., ua(:,:,1), 'U', -1.)
+                      CASE (4)
+                      CALL lbc_lnk_multi( un(:,:,1)  , 'U', -1., vn(:,:,1)  , 'V', -1., ua(:,:,1), 'U', -1., va(:,:,1)  , 'V', -1.)
+                      CASE (5)
+                      CALL lbc_lnk_multi( un(:,:,1)  , 'U', -1., vn(:,:,1)  , 'V', -1., ua(:,:,1), 'U', -1., va(:,:,1)  , 'V', -1., ub(:,:,1)  , 'U', -1.)
+                      CASE (6)
+                      CALL lbc_lnk_multi( un(:,:,1)  , 'U', -1., vn(:,:,1)  , 'V', -1., ua(:,:,1), 'U', -1., va(:,:,1)  , 'V', -1., ub(:,:,1)  , 'U', -1., vb(:,:,1)  , 'V', -1. )
+                    END SELECT
+                 ELSE
+                    SELECT CASE ( icomm_sequence(i,2) )
+                      CASE (1)
+                      CALL lbc_lnk_multi( un  , 'U', -1. )
+                      CASE (2)
+                      CALL lbc_lnk_multi( un  , 'U', -1., vn  , 'V', -1.)
+                      CASE (3)
+                      CALL lbc_lnk_multi( un  , 'U', -1., vn  , 'V', -1., ua, 'U', -1.)
+                      CASE (4)
+                      CALL lbc_lnk_multi( un  , 'U', -1., vn  , 'V', -1., ua, 'U', -1., va  , 'V', -1.)
+                      CASE (5)
+                      CALL lbc_lnk_multi( un  , 'U', -1., vn  , 'V', -1., ua, 'U', -1., va  , 'V', -1., ub  , 'U', -1.)
+                      CASE (6)
+                      CALL lbc_lnk_multi( un  , 'U', -1., vn  , 'V', -1., ua, 'U', -1., va  , 'V', -1., ub  , 'U', -1., vb  , 'V', -1. )
+                    END SELECT
+                 ENDIF
+               ENDDO
+            ELSE
+              CALL stp        ( istp ) 
+            ENDIF
             istp = istp + 1
+            IF ( istp == ( nitend - 1) ) tic = MPI_Wtime() - tic
          END DO
+         WRITE(6,'(A20,F11.6,A15,I8)') 'timeloop duration : ',tic,' on MPI rank : ', narea
          !
       ELSE                                            !==  diurnal SST time-steeping only  ==!
@@ -596 +645 @@
       !
       ierr =        oce_alloc    ()    ! ocean 
-      ierr = ierr + dia_wri_alloc()
+      !EM useless alloc for BENCH
+      !EM ierr = ierr + dia_wri_alloc()
       ierr = ierr + dom_oce_alloc()    ! ocean domain
       ierr = ierr + zdf_oce_alloc()    ! ocean vertical physics