Changeset 11738 for branches/UKMO/dev_r5518_GO6_under_ice_relax_dr_hook/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
- Timestamp:
- 2019-10-21T09:34:11+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_under_ice_relax_dr_hook/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r6487 r11738 45 45 USE timing ! Timing 46 46 47 USE yomhook, ONLY: lhook, dr_hook 48 USE parkind1, ONLY: jprb, jpim 49 47 50 IMPLICIT NONE 48 51 PRIVATE … … 103 106 ! 104 107 NAMELIST/namzgr/ ln_zco, ln_zps, ln_sco, ln_isfcav 108 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 109 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 110 REAL(KIND=jprb) :: zhook_handle 111 112 CHARACTER(LEN=*), PARAMETER :: RoutineName='DOM_ZGR' 113 114 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 115 105 116 !!---------------------------------------------------------------------- 106 117 ! … … 173 184 IF( nn_timing == 1 ) CALL timing_stop('dom_zgr') 174 185 ! 186 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 175 187 END SUBROUTINE dom_zgr 176 188 … … 203 215 REAL(wp) :: zrefdep ! depth of the reference level (~10m) 204 216 REAL(wp) :: za2, zkth2, zacr2 ! Values for optional double tanh function set from parameters 217 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 218 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 219 REAL(KIND=jprb) :: zhook_handle 220 221 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_Z' 222 223 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 224 205 225 !!---------------------------------------------------------------------- 206 226 ! … … 346 366 IF( nn_timing == 1 ) CALL timing_stop('zgr_z') 347 367 ! 368 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 348 369 END SUBROUTINE zgr_z 349 370 … … 388 409 INTEGER , ALLOCATABLE, DIMENSION(:,:) :: idta ! global domain integer data 389 410 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zdta ! global domain scalar data 411 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 412 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 413 REAL(KIND=jprb) :: zhook_handle 414 415 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_BAT' 416 417 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 418 390 419 !!---------------------------------------------------------------------- 391 420 ! … … 600 629 IF( nn_timing == 1 ) CALL timing_stop('zgr_bat') 601 630 ! 631 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 602 632 END SUBROUTINE zgr_bat 603 633 … … 615 645 !!---------------------------------------------------------------------- 616 646 INTEGER :: ii0, ii1, ij0, ij1 ! temporary integers 647 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 648 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 649 REAL(KIND=jprb) :: zhook_handle 650 651 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_BAT_ZOOM' 652 653 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 654 617 655 !!---------------------------------------------------------------------- 618 656 ! … … 652 690 ENDIF 653 691 ! 692 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 654 693 END SUBROUTINE zgr_bat_zoom 655 694 … … 681 720 INTEGER :: icompt, ibtest, ikmax ! temporary integers 682 721 REAL(wp), POINTER, DIMENSION(:,:) :: zbathy 722 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 723 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 724 REAL(KIND=jprb) :: zhook_handle 725 726 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_BAT_CTL' 727 728 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 729 683 730 684 731 !!---------------------------------------------------------------------- … … 791 838 IF( nn_timing == 1 ) CALL timing_stop('zgr_bat_ctl') 792 839 ! 840 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 793 841 END SUBROUTINE zgr_bat_ctl 794 842 … … 809 857 INTEGER :: ji, jj ! dummy loop indices 810 858 REAL(wp), POINTER, DIMENSION(:,:) :: zmbk 859 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 860 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 861 REAL(KIND=jprb) :: zhook_handle 862 863 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_BOT_LEVEL' 864 865 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 866 811 867 !!---------------------------------------------------------------------- 812 868 ! … … 836 892 IF( nn_timing == 1 ) CALL timing_stop('zgr_bot_level') 837 893 ! 894 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 838 895 END SUBROUTINE zgr_bot_level 839 896 … … 853 910 INTEGER :: ji, jj ! dummy loop indices 854 911 REAL(wp), POINTER, DIMENSION(:,:) :: zmik 912 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 913 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 914 REAL(KIND=jprb) :: zhook_handle 915 916 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_TOP_LEVEL' 917 918 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 919 855 920 !!---------------------------------------------------------------------- 856 921 ! … … 882 947 IF( nn_timing == 1 ) CALL timing_stop('zgr_top_level') 883 948 ! 949 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 884 950 END SUBROUTINE zgr_top_level 885 951 … … 893 959 !!---------------------------------------------------------------------- 894 960 INTEGER :: jk 961 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 962 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 963 REAL(KIND=jprb) :: zhook_handle 964 965 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_ZCO' 966 967 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 968 895 969 !!---------------------------------------------------------------------- 896 970 ! … … 912 986 IF( nn_timing == 1 ) CALL timing_stop('zgr_zco') 913 987 ! 988 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 914 989 END SUBROUTINE zgr_zco 915 990 … … 970 1045 REAL(wp) :: zrefdep ! temporary scalar 971 1046 REAL(wp), POINTER, DIMENSION(:,:,:) :: zprt 1047 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 1048 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 1049 REAL(KIND=jprb) :: zhook_handle 1050 1051 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_ZPS' 1052 1053 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 1054 972 1055 !!--------------------------------------------------------------------- 973 1056 ! … … 1252 1335 IF( nn_timing == 1 ) CALL timing_stop('zgr_zps') 1253 1336 ! 1337 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 1254 1338 END SUBROUTINE zgr_zps 1255 1339 … … 1282 1366 REAL(wp), POINTER, DIMENSION(:,:) :: zrisfdep, zbathy, zmask ! 2D workspace (ISH) 1283 1367 INTEGER , POINTER, DIMENSION(:,:) :: zmbathy, zmisfdep ! 2D workspace (ISH) 1368 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 1369 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 1370 REAL(KIND=jprb) :: zhook_handle 1371 1372 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_ISF' 1373 1374 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 1375 1284 1376 !!--------------------------------------------------------------------- 1285 1377 ! … … 1758 1850 IF( nn_timing == 1 ) CALL timing_stop('zgr_isf') 1759 1851 1852 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 1760 1853 END SUBROUTINE 1761 1854 … … 1816 1909 NAMELIST/namzgr_sco/ln_s_sh94, ln_s_sf12, ln_sigcrit, rn_sbot_min, rn_sbot_max, rn_hc, rn_rmax,rn_theta, & 1817 1910 rn_thetb, rn_bb, rn_alpha, rn_efold, rn_zs, rn_zb_a, rn_zb_b 1911 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 1912 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 1913 REAL(KIND=jprb) :: zhook_handle 1914 1915 CHARACTER(LEN=*), PARAMETER :: RoutineName='ZGR_SCO' 1916 1917 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 1918 1818 1919 !!---------------------------------------------------------------------- 1819 1920 ! … … 2240 2341 IF( nn_timing == 1 ) CALL timing_stop('zgr_sco') 2241 2342 ! 2343 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 2242 2344 END SUBROUTINE zgr_sco 2243 2345 … … 2261 2363 REAL(wp), POINTER, DIMENSION(:,:,:) :: z_gsigw3, z_gsigt3, z_gsi3w3 2262 2364 REAL(wp), POINTER, DIMENSION(:,:,:) :: z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 2365 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 2366 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 2367 REAL(KIND=jprb) :: zhook_handle 2368 2369 CHARACTER(LEN=*), PARAMETER :: RoutineName='S_SH94' 2370 2371 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 2372 2263 2373 2264 2374 CALL wrk_alloc( jpi, jpj, jpk, z_gsigw3, z_gsigt3, z_gsi3w3 ) … … 2339 2449 CALL wrk_dealloc( jpi, jpj, jpk, z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 2340 2450 2451 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 2341 2452 END SUBROUTINE s_sh94 2342 2453 … … 2365 2476 REAL(wp), POINTER, DIMENSION(:,:,:) :: z_gsigw3, z_gsigt3, z_gsi3w3 2366 2477 REAL(wp), POINTER, DIMENSION(:,:,:) :: z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 2478 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 2479 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 2480 REAL(KIND=jprb) :: zhook_handle 2481 2482 CHARACTER(LEN=*), PARAMETER :: RoutineName='S_SF12' 2483 2484 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 2485 2367 2486 2368 2487 ! … … 2476 2595 CALL wrk_dealloc( jpi, jpj, jpk, z_esigt3, z_esigw3, z_esigtu3, z_esigtv3, z_esigtf3, z_esigwu3, z_esigwv3 ) 2477 2596 2597 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 2478 2598 END SUBROUTINE s_sf12 2479 2599 … … 2495 2615 REAL(wp), POINTER, DIMENSION(:) :: z_gsigw, z_gsigt, z_gsi3w 2496 2616 REAL(wp), POINTER, DIMENSION(:) :: z_esigt, z_esigw 2617 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 2618 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 2619 REAL(KIND=jprb) :: zhook_handle 2620 2621 CHARACTER(LEN=*), PARAMETER :: RoutineName='S_TANH' 2622 2623 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 2624 2497 2625 2498 2626 CALL wrk_alloc( jpk, z_gsigw, z_gsigt, z_gsi3w ) … … 2550 2678 CALL wrk_dealloc( jpk, z_esigt, z_esigw ) 2551 2679 2680 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 2552 2681 END SUBROUTINE s_tanh 2553 2682 … … 2565 2694 REAL(wp), INTENT(in) :: pk ! continuous "k" coordinate 2566 2695 REAL(wp) :: pf ! sigma value 2696 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 2697 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 2698 REAL(KIND=jprb) :: zhook_handle 2699 2700 CHARACTER(LEN=*), PARAMETER :: RoutineName='FSSIG' 2701 2702 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 2703 2567 2704 !!---------------------------------------------------------------------- 2568 2705 ! … … 2573 2710 & / ( 2._wp * SINH( rn_theta ) ) 2574 2711 ! 2712 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 2575 2713 END FUNCTION fssig 2576 2714 … … 2590 2728 REAL(wp), INTENT(in) :: pbb ! Stretching coefficient 2591 2729 REAL(wp) :: pf1 ! sigma value 2730 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 2731 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 2732 REAL(KIND=jprb) :: zhook_handle 2733 2734 CHARACTER(LEN=*), PARAMETER :: RoutineName='FSSIG1' 2735 2736 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 2737 2592 2738 !!---------------------------------------------------------------------- 2593 2739 ! … … 2600 2746 ENDIF 2601 2747 ! 2748 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 2602 2749 END FUNCTION fssig1 2603 2750 … … 2629 2776 REAL(wp) :: za,zb,zx ! local variables 2630 2777 integer :: jk 2778 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 2779 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 2780 REAL(KIND=jprb) :: zhook_handle 2781 2782 CHARACTER(LEN=*), PARAMETER :: RoutineName='FGAMMA' 2783 2784 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 2785 2631 2786 !!---------------------------------------------------------------------- 2632 2787 ! … … 2652 2807 2653 2808 ! 2809 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 2654 2810 END FUNCTION fgamma 2655 2811
Note: See TracChangeset
for help on using the changeset viewer.