/[MITgcm]/MITgcm/pkg/gchem/gchem_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/gchem/gchem_readparms.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by stephd, Mon Oct 6 20:04:29 2003 UTC revision 1.18 by jmc, Tue Aug 24 14:41:32 2010 UTC
# Line 1  Line 1 
1  #include "PTRACERS_OPTIONS.h"  C $Header$
2    C $Name$
3    
4  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
5    
6  CBOP  CBOP
# Line 8  C !INTERFACE: ========================== Line 10  C !INTERFACE: ==========================
10        SUBROUTINE GCHEM_READPARMS( myThid )        SUBROUTINE GCHEM_READPARMS( myThid )
11    
12  C !DESCRIPTION:  C !DESCRIPTION:
13  C     Initialize PTRACERS parameters, read in data.gchem  C     Initialize GCHEM parameters, read in data.gchem
14    
15  C !USES: ===============================================================  C !USES: ===============================================================
16        IMPLICIT NONE        IMPLICIT NONE
17  #include "SIZE.h"  #include "SIZE.h"
18  #include "EEPARAMS.h"  #include "EEPARAMS.h"
19    #include "PARAMS.h"
20  #include "GCHEM.h"  #include "GCHEM.h"
21    
22  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
# Line 23  C  myThid               :: thread number Line 26  C  myThid               :: thread number
26  C !OUTPUT PARAMETERS: ==================================================  C !OUTPUT PARAMETERS: ==================================================
27  C  none  C  none
28    
 #ifdef ALLOW_PTRACERS  
29  #ifdef ALLOW_GCHEM  #ifdef ALLOW_GCHEM
30    
31  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
32    C  tIter0               :: retired parameter
33  C  iTracer              :: loop indices  C  iTracer              :: loop indices
34  C  iUnit                :: unit number for I/O  C  iUnit                :: unit number for I/O
35  C  msgBuf               :: message buffer  C  msgBuf               :: message buffer
36          INTEGER tIter0
37        INTEGER iTracer        INTEGER iTracer
38        INTEGER iUnit        INTEGER iUnit
39        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
40    c     PARAMETER ( UNSET_I      = 123456789  )
41  CEOP  CEOP
42    
43  C Set defaults values for parameters in GCHEM.h  C- Sub-package on/off flags: not fully implemented, requires
44         tIter0=0  C  to test the flag before any corresponding pkg S/R call
45         nsubtime=1  
46         windFile=' '        NAMELIST /GCHEM_PARM01/
47         atmospFile=' '       &                   nsubtime,
48         IceFile=' '       &                   useDIC,
49         IronFile=' '       &                   useCFC,
50         SilicaFile=' '       &                   useDARWIN,
51         &                   Filename1,
52         &                   Filename2,
53         &                   Filename3,
54         &                   Filename4,
55         &                   Filename5,
56         &           gchem_int1, gchem_int2, gchem_int3,
57         &           gchem_int4, gchem_int5,
58         &           gchem_rl1, gchem_rl2, gchem_rl3,
59         &           gchem_rl4, gchem_rl5,
60         &           gchem_ForcingPeriod, gchem_ForcingCycle,
61         &           tIter0
62    
63    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
64    
 C Open and read the data.ptracers file  
65        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
66    
67    C Set defaults values for parameters in GCHEM.h
68           useDIC = .FALSE.
69           useCFC = .FALSE.
70           useDARWIN = .FALSE.
71           tIter0 = UNSET_I
72           nsubtime=1
73           Filename1=' '
74           Filename2=' '
75           Filename3=' '
76           Filename4=' '
77           Filename5=' '
78           gchem_int1=0
79           gchem_int2=0
80           gchem_int3=0
81           gchem_int4=0
82           gchem_int5=0
83           gchem_rl1=0. _d 0
84           gchem_rl2=0. _d 0
85           gchem_rl3=0. _d 0
86           gchem_rl4=0. _d 0
87           gchem_rl5=0. _d 0
88    c default periodic forcing to same as for physics
89           gchem_ForcingPeriod=externForcingPeriod
90           gchem_ForcingCycle=externForcingCycle
91    C Open and read the data.gchem file
92        WRITE(msgBuf,'(A)') ' GCHEM_READPARMS: opening data.gchem'        WRITE(msgBuf,'(A)') ' GCHEM_READPARMS: opening data.gchem'
93        CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
94       &                   SQUEEZE_RIGHT , 1)       &                   SQUEEZE_RIGHT , 1)
# Line 62  C Open and read the data.ptracers file Line 104  C Open and read the data.ptracers file
104    
105  C Close the open data file  C Close the open data file
106        CLOSE(iUnit)        CLOSE(iUnit)
107    
108    C- Check for retired parameters:
109          IF ( tIter0 .NE. UNSET_I ) THEN
110    c      nRetired = nRetired+1
111           WRITE(msgBuf,'(A,A)')
112         &  'S/R GCHEM_READPARMS: Paramater "tIter0" is',
113         &  ' no longer allowed in file "data.gchem"'
114           CALL PRINT_ERROR( msgBuf , myThid)
115           WRITE(msgBuf,'(A,A)')
116         &  'S/R GCHEM_READPARMS: "tIter0" has been moved to',
117         &  ' PTRACERS_Iter0 in file "data.ptracers".'
118           CALL PRINT_ERROR( msgBuf , myThid)
119           STOP 'ABNORMAL END: S/R GCHEM_READPARMS'
120          ENDIF
121    
122        _END_MASTER(myThid)        _END_MASTER(myThid)
123    
124  C Everyone else must wait for the parameters to be loaded  C Everyone else must wait for the parameters to be loaded
125        _BARRIER        _BARRIER
126    
127    #ifdef ALLOW_DIC
128          IF ( useDIC ) THEN
129            CALL DIC_READPARMS(myThid)
130          ENDIF
131    #endif
132    
133    #ifdef ALLOW_CFC
134          IF ( useCFC ) THEN
135            CALL CFC_READPARMS(myThid)
136          ENDIF
137  #endif  #endif
138  #endif /* ALLOW_PTRACERS */  
139    #ifdef ALLOW_DARWIN
140          IF ( useDARWIN ) THEN
141            CALL DARWIN_READPARMS(myThid)
142          ENDIF
143    #endif
144    
145    #endif /* ALLOW_GCHEM */
146    
147        RETURN        RETURN
148        END        END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22