Changeset 11738 for branches/UKMO/dev_r5518_GO6_under_ice_relax_dr_hook/NEMOGCM/NEMO/OPA_SRC/STO/stopar.F90
- Timestamp:
- 2019-10-21T09:34:11+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_under_ice_relax_dr_hook/NEMOGCM/NEMO/OPA_SRC/STO/stopar.F90
r9321 r11738 24 24 USE timing 25 25 26 USE yomhook, ONLY: lhook, dr_hook 27 USE parkind1, ONLY: jprb, jpim 28 26 29 IMPLICIT NONE 27 30 PRIVATE … … 159 162 INTEGER :: ji, jj, jk, jsto, jflt 160 163 REAL(wp) :: stomax 164 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 165 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 166 REAL(KIND=jprb) :: zhook_handle 167 168 CHARACTER(LEN=*), PARAMETER :: RoutineName='STO_PAR' 169 170 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 171 161 172 162 173 ! … … 236 247 END DO 237 248 249 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 238 250 END SUBROUTINE sto_par 239 251 … … 259 271 REAL(wp) :: rinflate 260 272 INTEGER :: ios ! Local integer output status for namelist read 273 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 274 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 275 REAL(KIND=jprb) :: zhook_handle 276 277 CHARACTER(LEN=*), PARAMETER :: RoutineName='STO_PAR_INIT' 278 279 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 280 261 281 262 282 ! Read namsto namelist : stochastic parameterization … … 666 686 IF( ln_rststo ) CALL sto_rst_read 667 687 688 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 668 689 END SUBROUTINE sto_par_init 669 690 … … 683 704 CHARACTER(LEN=9) :: clsto3d='sto3d_000' ! stochastic parameter variable name 684 705 CHARACTER(LEN=10) :: clseed='seed0_0000' ! seed variable name 706 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 707 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 708 REAL(KIND=jprb) :: zhook_handle 709 710 CHARACTER(LEN=*), PARAMETER :: RoutineName='STO_RST_READ' 711 712 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 713 685 714 686 715 IF ( jpsto2d > 0 .OR. jpsto3d > 0 ) THEN … … 725 754 ENDIF 726 755 756 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 727 757 END SUBROUTINE sto_rst_read 728 758 … … 744 774 CHARACTER(LEN=9) :: clsto3d='sto3d_000' ! stochastic parameter variable name 745 775 CHARACTER(LEN=10) :: clseed='seed0_0000' ! seed variable name 776 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 777 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 778 REAL(KIND=jprb) :: zhook_handle 779 780 CHARACTER(LEN=*), PARAMETER :: RoutineName='STO_RST_WRITE' 781 782 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 783 746 784 747 785 IF ( jpsto2d > 0 .OR. jpsto3d > 0 ) THEN … … 808 846 ENDIF 809 847 848 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 810 849 END SUBROUTINE sto_rst_write 811 850 … … 821 860 INTEGER :: ji, jj 822 861 REAL(KIND=8) :: gran ! Gaussian random number (forced KIND=8 as in kiss_gaussian) 862 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 863 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 864 REAL(KIND=jprb) :: zhook_handle 865 866 CHARACTER(LEN=*), PARAMETER :: RoutineName='STO_PAR_WHITE' 867 868 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 869 823 870 824 871 DO jj = 1, jpj … … 829 876 END DO 830 877 878 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 831 879 END SUBROUTINE sto_par_white 832 880 … … 841 889 !! 842 890 INTEGER :: ji, jj 891 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 892 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 893 REAL(KIND=jprb) :: zhook_handle 894 895 CHARACTER(LEN=*), PARAMETER :: RoutineName='STO_PAR_FLT' 896 897 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 898 843 899 844 900 DO jj = 2, jpj-1 … … 850 906 END DO 851 907 908 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 852 909 END SUBROUTINE sto_par_flt 853 910 … … 869 926 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pflta 870 927 REAL(wp) :: ratio 928 INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 929 INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 930 REAL(KIND=jprb) :: zhook_handle 931 932 CHARACTER(LEN=*), PARAMETER :: RoutineName='STO_PAR_FLT_FAC' 933 934 IF (lhook) CALL dr_hook(RoutineName,zhook_in,zhook_handle) 935 871 936 872 937 pflt0(-1,-1) = 0 ; pflt0(-1,0) = 1 ; pflt0(-1,1) = 0 … … 901 966 sto_par_flt_fac = ratio 902 967 968 IF (lhook) CALL dr_hook(RoutineName,zhook_out,zhook_handle) 903 969 END FUNCTION sto_par_flt_fac 904 970
Note: See TracChangeset
for help on using the changeset viewer.