source: NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfutils.F90 @ 11852

Last change on this file since 11852 was 11852, checked in by mathiot, 12 months ago

ENHANCE-02_ISF_nemo: fix WED025 restartability, finish removing useless USE, remove useless lbc_lnk

File size: 4.3 KB
Line 
1MODULE isfutils
2   !!======================================================================
3   !!                       ***  MODULE  isfutils  ***
4   !! istutils module : miscelenious useful routines
5   !!======================================================================
6   !! History :  4.1  !  2019-09  (P. Mathiot) original code
7   !!----------------------------------------------------------------------
8
9   !!----------------------------------------------------------------------
10   !!   isfutils       : - read_2dcstdta to read a constant input file with iom_get
11   !!                    - debug to print array sum, min, max in ocean.output
12   !!----------------------------------------------------------------------
13
14   USE iom           , ONLY: iom_open, iom_get, iom_close, jpdom_data ! read input file
15   USE lib_fortran   , ONLY: glob_sum, glob_min, glob_max             ! compute global value
16   USE dom_oce       , ONLY: jpi,jpj,jpk                              ! domain size
17   USE in_out_manager, ONLY: wp, lwp, numout                          ! miscelenious
18
19   IMPLICIT NONE
20
21   PRIVATE
22
23   INTERFACE debug
24      MODULE PROCEDURE debug2d, debug3d
25   END INTERFACE debug
26
27   PUBLIC read_2dcstdta, debug
28
29CONTAINS
30
31   SUBROUTINE read_2dcstdta(cdfile, cdvar, pvar)
32      !!--------------------------------------------------------------------
33      !!                  ***  ROUTINE read_2dcstdta  ***
34      !!
35      !! ** Purpose : read input file
36      !!
37      !!--------------------------------------------------------------------
38      !!-------------------------- OUT -------------------------------------
39      REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) :: pvar          ! output variable
40      !!-------------------------- IN  -------------------------------------
41      CHARACTER(len=256)          , INTENT(in   ) :: cdfile   ! input file name
42      CHARACTER(len=34)           , INTENT(in   ) :: cdvar    ! variable name
43      !!--------------------------------------------------------------------
44      INTEGER :: inum
45      !!--------------------------------------------------------------------
46
47      CALL iom_open( TRIM(cdfile), inum )
48      CALL iom_get( inum, jpdom_data, TRIM(cdvar), pvar)
49      CALL iom_close(inum)
50
51   END SUBROUTINE read_2dcstdta
52
53   SUBROUTINE debug2d(cdtxt,pvar)
54      !!--------------------------------------------------------------------
55      !!                  ***  ROUTINE isf_debug2d  ***
56      !!
57      !! ** Purpose : add debug print
58      !!
59      !!--------------------------------------------------------------------
60      !!-------------------------- OUT -------------------------------------
61      !!-------------------------- IN  -------------------------------------
62      CHARACTER(LEN=256)          , INTENT(in   ) :: cdtxt
63      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pvar
64      !!--------------------------------------------------------------------
65      REAL(wp) :: zmin, zmax, zsum
66      !!--------------------------------------------------------------------
67      !
68      zsum = glob_sum( 'debug', pvar(:,:) )
69      zmin = glob_min( 'debug', pvar(:,:) )
70      zmax = glob_max( 'debug', pvar(:,:) )
71      !
72      IF (lwp) WRITE(numout,*) TRIM(cdtxt),' (min, max, sum) : ',zmin, zmax, zsum
73      !
74      CALL FLUSH(numout)
75      !
76   END SUBROUTINE debug2d
77
78   SUBROUTINE debug3d(cdtxt,pvar)
79      !!--------------------------------------------------------------------
80      !!                  ***  ROUTINE isf_debug3d  ***
81      !!
82      !! ** Purpose : add debug print
83      !!
84      !!--------------------------------------------------------------------
85      !!-------------------------- OUT -------------------------------------
86      !!-------------------------- IN  -------------------------------------
87      CHARACTER(LEN=256)              , INTENT(in   ) :: cdtxt
88      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in   ) :: pvar
89      !!--------------------------------------------------------------------
90      REAL(wp) :: zmin, zmax, zsum
91      !!--------------------------------------------------------------------
92      !
93      zsum = glob_sum( 'debug', pvar(:,:,:) )
94      zmin = glob_min( 'debug', pvar(:,:,:) )
95      zmax = glob_max( 'debug', pvar(:,:,:) )
96      !
97      IF (lwp) WRITE(numout,*) TRIM(cdtxt),' (min, max, sum) : ',zmin, zmax, zsum
98      !
99      CALL FLUSH(numout)
100      !
101   END SUBROUTINE debug3d
102
103END MODULE isfutils
104
Note: See TracBrowser for help on using the repository browser.