forked from altMITgcm/MITgcm
-
Notifications
You must be signed in to change notification settings - Fork 0
/
gchem_diagnostics_init.F
85 lines (68 loc) · 2.44 KB
/
gchem_diagnostics_init.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
C $Header: /u/gcmpack/MITgcm/pkg/gchem/gchem_diagnostics_init.F,v 1.8 2014/07/18 22:12:17 jmc Exp $
C $Name: $
#include "GCHEM_OPTIONS.h"
C !INTERFACE: ==========================================================
SUBROUTINE GCHEM_DIAGNOSTICS_INIT( myThid )
C !DESCRIPTION:
C define diagnostics for gchem package
C experiment
C !USES: ===============================================================
IMPLICIT NONE
#include "SIZE.h"
#include "PTRACERS_SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "PTRACERS_PARAMS.h"
#include "GAD.h"
C !INPUT PARAMETERS: ===================================================
C myThid :: thread number
INTEGER myThid
CEOP
#ifdef ALLOW_DIAGNOSTICS
#ifdef GCHEM_SEPARATE_FORCING
C Functions ::
INTEGER ILNBLNK
EXTERNAL ILNBLNK
c CHARACTER*4 GAD_DIAG_SUFX
c EXTERNAL GAD_DIAG_SUFX
C !LOCAL VARIABLES:
INTEGER diagNum
CHARACTER*8 diagName
CHARACTER*16 diagCode
CHARACTER*16 diagUnits
CHARACTER*(80) diagTitle
INTEGER iTrc, ilnb
CHARACTER*7 trcUnits
CHARACTER*30 locName
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
C-- Define diagnostics Names :
DO iTrc=1,PTRACERS_num
C- Set default name & tracer Units:
WRITE(locName,'(A,I4.4,A)') 'Tracer ',iTrc
trcUnits = '(kg/kg)'
C- use units & (short) name from data.ptracers :
ilnb = ILNBLNK(PTRACERS_units(iTrc))
IF ( ilnb.GE.1 ) THEN
ilnb = LEN(trcUnits)
trcUnits = PTRACERS_units(iTrc)(1:ilnb)
ENDIF
ilnb = ILNBLNK(PTRACERS_names(iTrc))
IF ( ilnb.GE.1 ) THEN
ilnb = MIN(LEN(locName),ilnb)
WRITE(locName,'(A)') PTRACERS_names(iTrc)(1:ilnb)
ENDIF
ilnb = MAX(ILNBLNK(locName),1)
WRITE(diagName,'(A5,A2,A1)') 'GC_Tr',PTRACERS_ioLabel(iTrc),' '
WRITE(diagTitle,'(2A)') locName(1:ilnb),
& ' concentration before GCHEM'
diagUnits = trcUnits//' '
diagCode = 'SMR MR '
CALL DIAGNOSTICS_ADDTOLIST( diagNum,
I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
ENDDO
C-- Other diagnostics should be initialized in each package called by gchem
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
#endif /* GCHEM_SEPARATE_FORCING */
#endif /* ALLOW_DIAGNOSTICS */
RETURN
END