/[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.2 by stephd, Wed Jul 9 19:59:52 2003 UTC revision 1.7 by jmc, Sat May 14 22:51:45 2005 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
# Line 23  C  myThid               :: thread number Line 25  C  myThid               :: thread number
25  C !OUTPUT PARAMETERS: ==================================================  C !OUTPUT PARAMETERS: ==================================================
26  C  none  C  none
27    
28  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_GCHEM
29    
30  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
31    C  tIter0               :: retired parameter
32  C  iTracer              :: loop indices  C  iTracer              :: loop indices
33  C  iUnit                :: unit number for I/O  C  iUnit                :: unit number for I/O
34  C  msgBuf               :: message buffer  C  msgBuf               :: message buffer
35          INTEGER tIter0, UNSET_I
36        INTEGER iTracer        INTEGER iTracer
37        INTEGER iUnit        INTEGER iUnit
38        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
39          PARAMETER ( UNSET_I      = 123456789  )
40  CEOP  CEOP
41    
42          NAMELIST /GCHEM_PARM01/
43         &                   tIter0,
44         &                   WindFile,
45         &                   IceFile,
46         &                   AtmospFile,
47         &                   IronFile,
48         &                   SilicaFile,
49         &                   nsubtime
50    
51  C Set defaults values for parameters in GCHEM.h  C Set defaults values for parameters in GCHEM.h
52         tIter0=0         tIter0 = UNSET_I
53         nsubtime=1         nsubtime=1
54         windFile=' '         windFile=' '
55         atmospFile=' '         atmospFile=' '
# Line 44  C Set defaults values for parameters in Line 58  C Set defaults values for parameters in
58         SilicaFile=' '         SilicaFile=' '
59    
60    
61  C Open and read the data.ptracers file  C Open and read the data.gchem file
62        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
63        WRITE(msgBuf,'(A)') ' GCHEM_READPARMS: opening data.gchem'        WRITE(msgBuf,'(A)') ' GCHEM_READPARMS: opening data.gchem'
64        CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
# Line 61  C Open and read the data.ptracers file Line 75  C Open and read the data.ptracers file
75    
76  C Close the open data file  C Close the open data file
77        CLOSE(iUnit)        CLOSE(iUnit)
78    
79    C- Check for retired parameters:
80          IF ( tIter0 .NE. UNSET_I ) THEN
81    c      nRetired = nRetired+1
82           WRITE(msgBuf,'(A,A)')
83         &  'S/R GCHEM_READPARMS: Paramater "tIter0" is',
84         &  ' no longer allowed in file "data.gchem"'
85           CALL PRINT_ERROR( msgBuf , myThid)
86           WRITE(msgBuf,'(A,A)')
87         &  'S/R GCHEM_READPARMS: "tIter0" has been moved to',
88         &  ' PTRACERS_Iter0 in file "data.ptracers".'
89           CALL PRINT_ERROR( msgBuf , myThid)
90           STOP 'ABNORMAL END: S/R GCHEM_READPARMS'
91          ENDIF
92    
93        _END_MASTER(myThid)        _END_MASTER(myThid)
94    
95  C Everyone else must wait for the parameters to be loaded  C Everyone else must wait for the parameters to be loaded
96        _BARRIER        _BARRIER
97    
98  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_GCHEM */
99    
100        RETURN        RETURN
101        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22