Changeset 11738 for branches/UKMO/dev_r5518_GO6_under_ice_relax_dr_hook/NEMOGCM/NEMO/TOP_SRC/MEDUSA/mocsy_phsolvers.F90

Timestamp:

2019-10-21T09:34:11+02:00 (5 years ago)

Author:

marc

Message:

The Dr Hook changes from my perl code.

File:

• branches/UKMO/dev_r5518_GO6_under_ice_relax_dr_hook/NEMOGCM/NEMO/TOP_SRC/MEDUSA/mocsy_phsolvers.F90 (modified) (20 diffs)

branches/UKMO/dev_r5518_GO6_under_ice_relax_dr_hook/NEMOGCM/NEMO/TOP_SRC/MEDUSA/mocsy_phsolvers.F90

@@ -17 +17 @@
 
 USE mocsy_singledouble
+USE yomhook, ONLY: lhook, dr_hook
+USE parkind1, ONLY: jprb, jpim
+
 IMPLICIT NONE
 
@@ -54 +57 @@
 
 USE mocsy_singledouble
+USE yomhook, ONLY: lhook, dr_hook
+USE parkind1, ONLY: jprb, jpim
+
 IMPLICIT NONE
 
@@ -67 +73 @@
 REAL(KIND=wp), INTENT(OUT) :: p_alknw_inf
 REAL(KIND=wp), INTENT(OUT) :: p_alknw_sup
+INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0
+INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1
+REAL(KIND=jprb)               :: zhook_handle
+
+CHARACTER(LEN=*), PARAMETER :: RoutineName='ANW_INFSUP'
+
+IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle)
+
 
 p_alknw_inf =  -p_po4tot - p_so4tot - p_flutot
@@ -73 +87 @@
 !             + p_nh4tot + p_h2stot
 
+IF (lhook) CALL dr_hook(RoutineName,zhook_out ,zhook_handle)
 RETURN
+IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle)
 END SUBROUTINE anw_infsup
 
@@ -85 +101 @@
 
 USE mocsy_singledouble
+USE yomhook, ONLY: lhook, dr_hook
+USE parkind1, ONLY: jprb, jpim
+
 IMPLICIT NONE
 REAL(KIND=wp) :: equation_at
@@ -115 +134 @@
 REAL(KIND=wp) ::                                      zalk_wat, zdalk_wat
 REAL(KIND=wp) :: aphscale
+INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0
+INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1
+REAL(KIND=jprb)               :: zhook_handle
+
+CHARACTER(LEN=*), PARAMETER :: RoutineName='EQUATION_AT'
+
+IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle)
+
 
 ! TOTAL H+ scale: conversion factor for Htot = aphscale * Hfree
@@ -212 +239 @@
                 - Kw/p_h**2 - 1._wp/aphscale
 ENDIF
+IF (lhook) CALL dr_hook(RoutineName,zhook_out ,zhook_handle)
 RETURN
+IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle)
 END FUNCTION equation_at
 
@@ -231 +260 @@
 
 USE mocsy_singledouble
+USE yomhook, ONLY: lhook, dr_hook
+USE parkind1, ONLY: jprb, jpim
+
 IMPLICIT NONE
 
@@ -244 +276 @@
 REAL(KIND=wp)  ::  zd, zsqrtd, zhmin
 REAL(KIND=wp)  ::  za2, za1, za0
+INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0
+INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1
+REAL(KIND=jprb)               :: zhook_handle
+
+CHARACTER(LEN=*), PARAMETER :: RoutineName='AHINI_FOR_AT'
+
+IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle)
+
 
 IF (p_alkcb <= 0._wp) THEN
@@ -274 +314 @@
 
 ENDIF
+IF (lhook) CALL dr_hook(RoutineName,zhook_out ,zhook_handle)
 RETURN
+IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle)
 END SUBROUTINE ahini_for_at
 
@@ -289 +331 @@
 
 USE mocsy_singledouble
+USE yomhook, ONLY: lhook, dr_hook
+USE parkind1, ONLY: jprb, jpim
+
 IMPLICIT NONE
 REAL(KIND=wp) :: SOLVE_AT_GENERAL
@@ -318 +363 @@
 LOGICAL        :: l_exitnow
 REAL(KIND=wp), PARAMETER :: pz_exp_threshold = 1.0_wp
+INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0
+INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1
+REAL(KIND=jprb)               :: zhook_handle
+
+CHARACTER(LEN=*), PARAMETER :: RoutineName='SOLVE_AT_GENERAL'
+
+IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle)
+
 
 ! TOTAL H+ scale: conversion factor for Htot = aphscale * Hfree
@@ -477 +530 @@
    ENDIF
 ENDIF
+IF (lhook) CALL dr_hook(RoutineName,zhook_out ,zhook_handle)
 RETURN
+IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle)
 END FUNCTION solve_at_general
 
@@ -493 +548 @@
 !USE MOD_CHEMCONST, ONLY: api1_wat, aphscale
 USE mocsy_singledouble
+USE yomhook, ONLY: lhook, dr_hook
+USE parkind1, ONLY: jprb, jpim
+
 IMPLICIT NONE
 REAL(KIND=wp) :: SOLVE_AT_GENERAL_SEC
@@ -519 +577 @@
 REAL(KIND=wp)  ::  aphscale
 LOGICAL        ::  l_exitnow
+INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0
+INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1
+REAL(KIND=jprb)               :: zhook_handle
+
+CHARACTER(LEN=*), PARAMETER :: RoutineName='SOLVE_AT_GENERAL_SEC'
+
+IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle)
+
 
 ! TOTAL H+ scale: conversion factor for Htot = aphscale * Hfree
@@ -590 +656 @@
    solve_at_general_sec = zh_2
    IF(PRESENT(p_val)) p_val = zeqn_2
+IF (lhook) CALL dr_hook(RoutineName,zhook_out ,zhook_handle)
    RETURN
 ENDIF
@@ -747 +814 @@
 ENDIF
 
+IF (lhook) CALL dr_hook(RoutineName,zhook_out ,zhook_handle)
 RETURN
+IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle)
 END FUNCTION SOLVE_AT_GENERAL_SEC
 
@@ -761 +830 @@
 
 USE mocsy_singledouble
+USE yomhook, ONLY: lhook, dr_hook
+USE parkind1, ONLY: jprb, jpim
+
 IMPLICIT NONE
 REAL(KIND=wp) :: SOLVE_AT_FAST
@@ -786 +858 @@
 LOGICAL        :: l_exitnow
 REAL(KIND=wp), PARAMETER :: pz_exp_threshold = 1.0_wp
+INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0
+INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1
+REAL(KIND=jprb)               :: zhook_handle
+
+CHARACTER(LEN=*), PARAMETER :: RoutineName='SOLVE_AT_FAST'
+
+IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle)
+
 
 ! TOTAL H+ scale: conversion factor for Htot = aphscale * Hfree
@@ -854 +934 @@
 ENDIF
 
+IF (lhook) CALL dr_hook(RoutineName,zhook_out ,zhook_handle)
 RETURN
+IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle)
 END FUNCTION solve_at_fast
 !===============================================================================