Changeset 14149 for NEMO

Timestamp:

2020-12-10T15:43:37+01:00 (3 years ago)

Author:

smueller

Message:

Addition of timer instructions to the subroutines of module zdfosm (ticket #2353)

File:

• NEMO/branches/2020/dev_r14122_ENHANCE-14_smueller_OSMOSIS_streamlining/src/OCE/ZDF/zdfosm.F90 (modified) (38 diffs)

NEMO/branches/2020/dev_r14122_ENHANCE-14_smueller_OSMOSIS_streamlining/src/OCE/ZDF/zdfosm.F90

@@ -67 +67 @@
    USE prtctl         ! Print control
    USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)
+   USE timing, ONLY : timing_start, timing_stop   ! Timer
 
    IMPLICIT NONE
@@ -346 +347 @@
       !!--------------------------------------------------------------------
       !
+      IF( ln_timing ) CALL timing_start('zdf_osm')
       ibld(:,:)   = 0     ; imld(:,:)  = 0
       zrad0(:,:)  = 0._wp ; zradh(:,:) = 0._wp ; zradav(:,:)    = 0._wp ; zustar(:,:)    = 0._wp
@@ -1251 +1253 @@
 
       END IF
+      IF( ln_timing ) CALL timing_stop('zdf_osm')
 
 CONTAINS
@@ -1280 +1283 @@
       REAL(wp), PARAMETER :: rn_vispyc_shr = 0.15
 
+      IF( ln_timing ) CALL timing_start('zdf_osm_dv')
       DO_2D( 0, 0, 0, 0 )
           IF ( lconv(ji,jj) ) THEN
@@ -1383 +1387 @@
           !
        END_2D
+      IF( ln_timing ) CALL timing_stop('zdf_osm_dv')
 
   END SUBROUTINE zdf_osm_diffusivity_viscosity
@@ -1421 +1426 @@
      REAL, PARAMETER :: zrot=0._wp  ! dummy rotation rate of surface stress.
 
+     IF( ln_timing ) CALL timing_start('zdf_osm_os')
 ! Determins stability and set flag lconv
      DO_2D( 0, 0, 0, 0 )
@@ -1531 +1537 @@
       ENDIF    ! lshear
      END_2D
+     IF( ln_timing ) CALL timing_stop('zdf_osm_os')
    END SUBROUTINE zdf_osm_osbl_state
 
@@ -1558 +1565 @@
 
 
+        IF( ln_timing ) CALL timing_start('zdf_osm_va')
         zt   = 0._wp
         zs   = 0._wp
@@ -1600 +1608 @@
          ENDIF
         END_2D
+     IF( ln_timing ) CALL timing_stop('zdf_osm_va')
    END SUBROUTINE zdf_osm_vertical_average
 
@@ -1619 +1628 @@
         REAL(wp) :: ztemp
 
+        IF( ln_timing ) CALL timing_start('zdf_osm_vr')
         DO_2D( 0, 0, 0, 0 )
            ztemp = zu(ji,jj)
@@ -1627 +1637 @@
            zdv(ji,jj) = zdv(ji,jj) * zsin_w(ji,jj) - ztemp * zsin_w(ji,jj)
         END_2D
+        IF( ln_timing ) CALL timing_stop('zdf_osm_vr')
     END SUBROUTINE zdf_osm_velocity_rotation
 
@@ -1651 +1662 @@
       REAL(wp)                      :: zpe_mle_ref, zwb_ent, zdbdz_mle_int
 
+      IF( ln_timing ) CALL timing_start('zdf_osm_osf')
       znd_param(:,:) = 0._wp
 
@@ -1732 +1744 @@
           ENDIF  ! lconv
         END_2D
+        IF( ln_timing ) CALL timing_stop('zdf_osm_osf')
     END SUBROUTINE zdf_osm_osbl_state_fk
 
@@ -1751 +1764 @@
 
 
+     IF( ln_timing ) CALL timing_start('zdf_osm_eg')
      DO_2D( 0, 0, 0, 0 )
         IF ( jbase(ji,jj)+1 < mbkt(ji,jj) ) THEN
@@ -1768 +1782 @@
         END IF
      END_2D
+     IF( ln_timing ) CALL timing_stop('zdf_osm_eg')
     END SUBROUTINE zdf_osm_external_gradients
 
@@ -1781 +1796 @@
      REAL(wp), PARAMETER :: zgamma_b = 2.25, zzeta_sh = 0.15
 
+     IF( ln_timing ) CALL timing_start('zdf_osm_pscp')
      DO_2D( 0, 0, 0, 0 )
         IF ( ibld(ji,jj) + jp_ext(ji,jj) < mbkt(ji,jj) ) THEN
@@ -1847 +1863 @@
         ENDIF      ! IF ( ibld(ji,jj) < mbkt(ji,jj) )
      END_2D
+     IF( ln_timing ) CALL timing_stop('zdf_osm_pscp')
 
     END SUBROUTINE zdf_osm_pycnocline_scalar_profiles
@@ -1866 +1883 @@
       REAL(wp) :: zzeta_v = 0.45
       !
+      IF( ln_timing ) CALL timing_start('zdf_osm_pshp')
       DO_2D( 0, 0, 0, 0 )
          !
@@ -1906 +1924 @@
          END IF      ! IF ( ibld(ji,jj) + ibld_ext < mbkt(ji,jj) )
       END_2D
+      IF( ln_timing ) CALL timing_stop('zdf_osm_pshp')
     END SUBROUTINE zdf_osm_pycnocline_shear_profiles
 
@@ -1929 +1948 @@
     REAL, PARAMETER :: a_ddh = 2.5, a_ddh_2 = 3.5 ! also in pycnocline_depth
 
+    IF( ln_timing ) CALL timing_start('zdf_osm_cd')
   DO_2D( 0, 0, 0, 0 )
 
@@ -2050 +2070 @@
       ENDIF ! lshear
   END_2D
+    IF( ln_timing ) CALL timing_stop('zdf_osm_cd')
     END SUBROUTINE zdf_osm_calculate_dhdt
 
@@ -2072 +2093 @@
     REAL(wp) :: zthermal, zbeta
 
+     IF( ln_timing ) CALL timing_start('zdf_osm_th')
      DO_2D( 0, 0, 0, 0 )
         IF ( ibld(ji,jj) - imld(ji,jj) > 1 ) THEN
@@ -2152 +2174 @@
         zhbl(ji,jj) = gdepw(ji,jj,ibld(ji,jj),Kmm)
      END_2D
+     IF( ln_timing ) CALL timing_stop('zdf_osm_th')
 
     END SUBROUTINE zdf_osm_timestep_hbl
@@ -2176 +2199 @@
       REAL, PARAMETER :: a_ddh_2 = 3.5 ! also in pycnocline_depth
 
+      IF( ln_timing ) CALL timing_start('zdf_osm_pt')
     DO_2D( 0, 0, 0, 0 )
 
@@ -2287 +2311 @@
       zdh(ji,jj) = zhbl(ji,jj) - zhml(ji,jj)
     END_2D
+      IF( ln_timing ) CALL timing_stop('zdf_osm_pt')
 
     END SUBROUTINE zdf_osm_pycnocline_thickness
@@ -2317 +2342 @@
 !!----------------------------------------------------------------------
       !
+      IF( ln_timing ) CALL timing_start('zdf_osm_zhg')
       !                                      !==  MLD used for MLE  ==!
 
@@ -2378 +2404 @@
              & + dbdx_mle(ji-1,jj) * dbdx_mle(ji-1,jj) + dbdy_mle(ji,jj-1) * dbdy_mle(ji,jj-1) ) )
       END_2D
+      IF( ln_timing ) CALL timing_stop('zdf_osm_zhg')
 
  END SUBROUTINE zdf_osm_zmld_horizontal_gradients
@@ -2398 +2425 @@
       REAL(wp) ::  ztmp, zdbdz, zdtdz, zdsdz, zthermal,zbeta, zbuoy, zdb_mle
 
+      IF( ln_timing ) CALL timing_start('zdf_osm_mp')
    ! Calculate vertical buoyancy, heat and salinity fluxes due to MLE.
 
@@ -2439 +2467 @@
          zhmle(ji,jj) = gdepw(ji,jj, mld_prof(ji,jj),Kmm)
       END_2D
+      IF( ln_timing ) CALL timing_stop('zdf_osm_mp')
 END SUBROUTINE zdf_osm_mle_parameters
 
@@ -2471 +2500 @@
      !!----------------------------------------------------------------------
      !
+     IF( ln_timing ) CALL timing_start('zdf_osm_init')
      READ  ( numnam_ref, namzdf_osm, IOSTAT = ios, ERR = 901)
 901  IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_osm in reference namelist' )
@@ -2672 +2702 @@
      ghamv(:,:,:) = 0.
      !
+     IF( ln_timing ) CALL timing_stop('zdf_osm_init')
    END SUBROUTINE zdf_osm_init
 
@@ -2698 +2729 @@
      !!----------------------------------------------------------------------
      !
+     IF( ln_timing ) CALL timing_start('osm_rst')
      !!-----------------------------------------------------------------------------
      ! If READ/WRITE Flag is 'READ', try to get hbl from restart file. If successful then return
@@ -2780 +2812 @@
      ww(:,:,:) = 0._wp
      WRITE(numout,*) ' ===>>>> :  wn not in restart file, set to zero initially'
+     IF( ln_timing ) CALL timing_stop('osm_rst')
    END SUBROUTINE osm_rst
 
@@ -2799 +2832 @@
       INTEGER :: ji, jj, jk
       !
+      IF( ln_timing ) CALL timing_start('tra_osm')
       IF( kt == nit000 ) THEN
          IF( ntile == 0 .OR. ntile == 1 ) THEN                    ! Do only on the first tile
@@ -2836 +2870 @@
       ENDIF
       !
+      IF( ln_timing ) CALL timing_stop('tra_osm')
    END SUBROUTINE tra_osm
 
@@ -2852 +2887 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) :: kt
+      IF( ln_timing ) CALL timing_start('trc_osm')
       WRITE(*,*) 'trc_osm: Not written yet', kt
+      IF( ln_timing ) CALL timing_stop('trc_osm')
    END SUBROUTINE trc_osm
 
@@ -2872 +2909 @@
       !!----------------------------------------------------------------------
       !
+      IF( ln_timing ) CALL timing_start('dyn_osm')
       IF( kt == nit000 ) THEN
          IF(lwp) WRITE(numout,*)
@@ -2890 +2928 @@
       ! code for saving tracer trends removed
       !
+      IF( ln_timing ) CALL timing_stop('dyn_osm')
    END SUBROUTINE dyn_osm