Changeset 3706 for LMDZ6/branches/Optimisation_LMDZ/libf/phylmd

Timestamp:

Jun 11, 2020, 11:09:38 AM (4 years ago)

Author:

adurocher

Message:

Added timers for physiq and display physic profiling

File:

• LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/physiq_mod.F90 (modified) (34 diffs)

LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/physiq_mod.F90

@@ -15 +15 @@
        flxmass_w, &
        d_u, d_v, d_t, d_qx, d_ps)
-
+    USE profiling_physic_mod, only : enter_profile, exit_profile, print_profile
     USE assert_m, only: assert
     USE ioipsl, only: histbeg, histvert, histdef, histend, histsync, &
@@ -1182 +1182 @@
     REAL pi
 
+    integer, save :: itau_profiling_physiq ! Print frequency for physiq profiling
+
     pi = 4. * ATAN(1.)
 
@@ -1226 +1228 @@
     CALL assert(paprs(:, nbp_lev + 1) < paprs(:, nbp_lev), &
             "physiq_mod paprs bad order")
-
-    IF (first) THEN 
+    IF (first) THEN      
+       call enter_profile("phy_init")
+
+       itau_profiling_physiq=-1 ! Default is -1 : never
+       CALL getin_p('itau_profiling_physiq', itau_profiling_physiq )
+
        CALL init_etat0_limit_unstruct
        IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
@@ -1284 +1290 @@
        first=.FALSE.
 
+       call exit_profile("phy_init")
     ENDIF  ! first
+    
+    call enter_profile("physiq")
 
     !ym => necessaire pour iflag_con != 2    
@@ -1304 +1313 @@
 
     IF (debut) THEN
+       call enter_profile("phy_init")
+
        CALL suphel ! initialiser constantes et parametres phys.
 ! tau_gl : constante de rappel de la temperature a la surface de la glace - en
@@ -1368 +1379 @@
        WRITE(lunout,*) 'Call to infocfields from physiq'
        CALL infocfields_init
-
     ENDIF
 
@@ -2022 +2032 @@
       sollwdown(:)= sollwdown(:) + betalwoff *(-1.*ZFLDN0(:,1) - &
                     sollwdown(:))
-
-
-    ENDIF
+      call exit_profile("phy_init")
+   ENDIF
     !
     !   ****************     Fin  de   IF ( debut  )   ***************
@@ -2445 +2454 @@
 
     IF (iflag_pbl/=0) THEN
+      
+       call enter_profile("phy_pbl")
 
        !jyg+nrlmd<
@@ -2623 +2634 @@
             t_seri(:,1),q_seri(:,1),u_seri(:,1),v_seri(:,1),evap_pot)
 
+      call exit_profile("phy_pbl")
     ENDIF
     ! =================================================================== c
@@ -2697 +2709 @@
 !!                       itapcv, cvpas, itap-1, cvpas_0
     IF (MOD(itapcv,cvpas).EQ.0 .OR. MOD(itap-1,cvpas_0).EQ.0) THEN
+      
+      call enter_profile("phy_convection")
 
     !
@@ -3018 +3032 @@
     proba_notrig(:) = 1.
     itapcv = 0
+    
+    call exit_profile("phy_convection")
     ENDIF !  (MOD(itapcv,cvpas).EQ.0 .OR. MOD(itapcv,cvpas_0).EQ.0)
 !
@@ -3108 +3124 @@
        !
        IF (MOD(itapwk,wkpas).EQ.0) THEN
-          !
+          call enter_profile("phy_wake")
+          
           DO k=1,klev
              DO i=1,klon
@@ -3166 +3183 @@
           !jyg    Reinitialize itapwk when wakes have been called
           itapwk = 0
+          call exit_profile("phy_wake")
        ENDIF !  (MOD(itapwk,wkpas).EQ.0)
        !
@@ -3233 +3251 @@
 
     ELSE
-
+       call enter_profile("phy_thermique")
        !  Thermiques
        !  ==========
@@ -3402 +3420 @@
        ENDIF
 
+       call exit_profile("phy_thermique")
     ENDIF
     !
@@ -3423 +3442 @@
        print *,'itap, ->fisrtilp ',itap
     ENDIF
-    !
+    call enter_profile("phy_ls_condens")
+    
     CALL fisrtilp(phys_tstep,paprs,pplay, &
          t_seri, q_seri,ptconv,ratqs, &
@@ -3436 +3456 @@
     WHERE (rain_lsc < 0) rain_lsc = 0.
     WHERE (snow_lsc < 0) snow_lsc = 0.
+    
+    call exit_profile("phy_ls_condens")
 
 !+JLD
@@ -3699 +3721 @@
 
     IF (type_trac == 'inca') THEN
+       call enter_profile("phy_inca")
 #ifdef INCA
        CALL VTe(VTphysiq)
@@ -3753 +3776 @@
        CALL VTb(VTphysiq)
 #endif 
+       call exit_profile("phy_inca")
     ENDIF !type_trac = inca
 
@@ -3760 +3784 @@
     !
     IF (MOD(itaprad,radpas).EQ.0) THEN
-
+       call enter_profile("phy_rayonnement")
        !
        !jq - introduce the aerosol direct and first indirect radiative forcings
        !jq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
        IF (flag_aerosol .GT. 0) THEN
+          call enter_profile("read_aerosol")
           IF (iflag_rrtm .EQ. 0) THEN !--old radiation
              IF (.NOT. aerosol_couple) THEN
@@ -3842 +3867 @@
              ENDIF
           ENDIF
+          call exit_profile("read_aerosol")
        ELSE   !--flag_aerosol = 0 
           tausum_aero(:,:,:) = 0.
@@ -3865 +3891 @@
        !--updates tausum_aero,tau_aero,piz_aero,cg_aero
        IF (flag_aerosol_strat.GT.0) THEN
+          call enter_profile("read_aerosol")
           IF (prt_level .GE.10) THEN
              PRINT *,'appel a readaerosolstrat', mth_cur
@@ -3895 +3922 @@
 #endif
           ENDIF
+          call exit_profile("read_aerosol")
        ELSE
           tausum_aero(:,:,id_STRAT_phy) = 0. 
@@ -4220 +4248 @@
           zxtsol(:) = zsav_tsol (:)
        ENDIF
+       call exit_profile("phy_rayonnement")
     ENDIF ! MOD(itaprad,radpas)
     itaprad = itaprad + 1
@@ -4502 +4531 @@
 ! Inititialization
 !------------------
-
+       call enter_profile("phy_init")
        addtkeoro=0   
        CALL getin_p('addtkeoro',addtkeoro) 
@@ -4514 +4543 @@
 
        smallscales_tkeoro=.FALSE.   
-       CALL getin_p('smallscales_tkeoro',smallscales_tkeoro) 
+       CALL getin_p('smallscales_tkeoro',smallscales_tkeoro)
+       call exit_profile("phy_init")
 
 
@@ -4605 +4635 @@
 
     IF (ok_cosp) THEN
+       call enter_profile("phy_cosp")
        ! adeclarer 
 #ifdef CPP_COSP
@@ -4684 +4715 @@
        ENDIF
 #endif
-
+       call exit_profile("phy_cosp")
     ENDIF  !ok_cosp
 
@@ -4726 +4757 @@
 
     IF (iflag_phytrac == 1 ) THEN
-
+      call enter_profile("phy_phytrac")
 #ifdef CPP_Dust
       CALL       phytracr_spl ( debut,lafin , jD_cur,jH_cur,iflag_con,       &  ! I
@@ -4770 +4801 @@
          tr_seri, init_source)
 #endif
+      call exit_profile("phy_phytrac")
     ENDIF    ! (iflag_phytrac=1)
 
@@ -5027 +5059 @@
     !   Ecriture des sorties
     !=============================================================
+    call enter_profile("phy_output")
+
+    if( itau_profiling_physiq>0 .and. 0 == mod( itap, itau_profiling_physiq ) ) call print_profile()
+
 #ifdef CPP_IOIPSL
 
@@ -5091 +5127 @@
 ! On remet des variables a .false. apres un premier appel
     IF (debut) THEN
+      call enter_profile("phy_init")
 #ifdef CPP_XIOS
       swaero_diag=.FALSE.
@@ -5137 +5174 @@
 !      write (lunout,*)'ok_4xCO2atm= ',swaero_diag, swaerofree_diag, dryaod_diag, ok_4xCO2atm
 #endif
-    ENDIF
+      call exit_profile("phy_init")
+   ENDIF
 
     !====================================================================
@@ -5177 +5215 @@
        WRITE(lunout,*) ' physiq fin, nombre de steps ou cvpas = 1 : ', Ncvpaseq1
     ENDIF
-
+    
+    call exit_profile("phy_output")
+    call exit_profile("physiq")
     !      first=.false.