/[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.8 by stephd, Thu Aug 18 18:23:02 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         &                   Filename1,
50         &                   Filename2,
51         &                   Filename3,
52         &                   Filename4,
53         &                   Filename5,
54         &                   nsubtime
55    
56  C Set defaults values for parameters in GCHEM.h  C Set defaults values for parameters in GCHEM.h
57         tIter0=0         tIter0 = UNSET_I
58         nsubtime=1         nsubtime=1
59         windFile=' '         windFile=' '
60         atmospFile=' '         atmospFile=' '
61         IceFile=' '         IceFile=' '
62         IronFile=' '         IronFile=' '
63         SilicaFile=' '         SilicaFile=' '
64           Filename1=' '
65           Filename2=' '
66  C Open and read the data.ptracers file         Filename3=' '
67           Filename4=' '
68           Filename5=' '
69    C Open and read the data.gchem file
70        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
71        WRITE(msgBuf,'(A)') ' GCHEM_READPARMS: opening data.gchem'        WRITE(msgBuf,'(A)') ' GCHEM_READPARMS: opening data.gchem'
72        CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
# Line 61  C Open and read the data.ptracers file Line 83  C Open and read the data.ptracers file
83    
84  C Close the open data file  C Close the open data file
85        CLOSE(iUnit)        CLOSE(iUnit)
86    
87    C- Check for retired parameters:
88          IF ( tIter0 .NE. UNSET_I ) THEN
89    c      nRetired = nRetired+1
90           WRITE(msgBuf,'(A,A)')
91         &  'S/R GCHEM_READPARMS: Paramater "tIter0" is',
92         &  ' no longer allowed in file "data.gchem"'
93           CALL PRINT_ERROR( msgBuf , myThid)
94           WRITE(msgBuf,'(A,A)')
95         &  'S/R GCHEM_READPARMS: "tIter0" has been moved to',
96         &  ' PTRACERS_Iter0 in file "data.ptracers".'
97           CALL PRINT_ERROR( msgBuf , myThid)
98           STOP 'ABNORMAL END: S/R GCHEM_READPARMS'
99          ENDIF
100    
101        _END_MASTER(myThid)        _END_MASTER(myThid)
102    
103  C Everyone else must wait for the parameters to be loaded  C Everyone else must wait for the parameters to be loaded
104        _BARRIER        _BARRIER
105    
106  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_GCHEM */
107    
108        RETURN        RETURN
109        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22