/[MITgcm]/MITgcm/model/src/ini_parms.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_parms.F

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

revision 1.17 by cnh, Fri Jun 12 19:33:33 1998 UTC revision 1.18 by cnh, Mon Jun 15 05:13:56 1998 UTC
# Line 72  C--   Time stepping parammeters Line 72  C--   Time stepping parammeters
72        NAMELIST /PARM03/        NAMELIST /PARM03/
73       & nIter0, nTimeSteps, deltaT, deltaTmom, deltaTtracer, abEps, tauCD, rCD,       & nIter0, nTimeSteps, deltaT, deltaTmom, deltaTtracer, abEps, tauCD, rCD,
74       & startTime, endTime, chkPtFreq, dumpFreq, deltaTClock, pChkPtFreq,       & startTime, endTime, chkPtFreq, dumpFreq, deltaTClock, pChkPtFreq,
75       & cAdjFreq       & cAdjFreq, tauThetaClimRelax, tauSaltClimRelax
76    
77  C--   Gridding parameters  C--   Gridding parameters
78        NAMELIST /PARM04/        NAMELIST /PARM04/
# Line 83  C--   Gridding parameters Line 83  C--   Gridding parameters
83  C--   Input files  C--   Input files
84        NAMELIST /PARM05/        NAMELIST /PARM05/
85       & bathyFile, hydrogThetaFile, hydrogSaltFile,       & bathyFile, hydrogThetaFile, hydrogSaltFile,
86       & zonalWindFile, meridWindFile       & zonalWindFile, meridWindFile, thetaClimFile,
87         & saltClimFile
88    
89  C  C
90        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
# Line 170  C--   Set default "physical" parameters Line 171  C--   Set default "physical" parameters
171        GMdepth=1000.        GMdepth=1000.
172        GMkbackground=0.        GMkbackground=0.
173        tAlpha=2.d-4        tAlpha=2.d-4
174        sBeta=1.d-5        sBeta=7.4d-4
175        eosType='LINEAR'        eosType='LINEAR'
176        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
177        rigidLid            = .FALSE.        rigidLid            = .FALSE.
# Line 287  C--   Time stepping parameters Line 288  C--   Time stepping parameters
288        cAdjFreq       = -1. _d 0        cAdjFreq       = -1. _d 0
289        rCD            = -1. _d 0        rCD            = -1. _d 0
290        tauCD          = 0. _d 0        tauCD          = 0. _d 0
291          tauThetaClimRelax = 0. _d 0
292          doThetaClimRelax  = .FALSE.
293          tauSaltClimRelax  = 0. _d 0
294          doSaltClimRelax   = .FALSE.
295        READ(UNIT=iUnit,NML=PARM03,IOSTAT=errIO,err=7)        READ(UNIT=iUnit,NML=PARM03,IOSTAT=errIO,err=7)
296        IF ( errIO .GE. 0 ) GOTO 8        IF ( errIO .GE. 0 ) GOTO 8
297      7 CONTINUE      7 CONTINUE
# Line 323  C     o CD coupling Line 328  C     o CD coupling
328        IF ( rCD .LT. 0. ) THEN        IF ( rCD .LT. 0. ) THEN
329         rCD = 1. - deltaTMom/tauCD         rCD = 1. - deltaTMom/tauCD
330        ENDIF        ENDIF
331    C     o Temperature climatology relaxation time scale
332          IF ( tauThetaClimRelax .EQ. 0. _d 0 ) THEN
333           doThetaClimRelax     = .FALSE.
334           lambdaThetaClimRelax = 0. _d 0
335          ELSE
336           doThetaClimRelax     = .TRUE.
337           lambdaThetaClimRelax = 1./tauThetaClimRelax
338          ENDIF
339    C     o Salinity climatology relaxation time scale
340          IF ( tauSaltClimRelax .EQ. 0. _d 0 ) THEN
341           doSaltClimRelax     = .FALSE.
342           lambdaSaltClimRelax = 0. _d 0
343          ELSE
344           doSaltClimRelax     = .TRUE.
345           lambdaSaltClimRelax = 1./tauSaltClimRelax
346          ENDIF
347  C     o Time step count  C     o Time step count
348        IF ( endTime .NE. 0 ) THEN        IF ( endTime .NE. 0 ) THEN
349         IF ( deltaTClock .NE. 0 ) nTimeSteps =         IF ( deltaTClock .NE. 0 ) nTimeSteps =
# Line 419  C     Initialize EOS coefficients (3rd o Line 440  C     Initialize EOS coefficients (3rd o
440         usingSphericalPolarMterms = .FALSE.         usingSphericalPolarMterms = .FALSE.
441         metricTerms = .FALSE.         metricTerms = .FALSE.
442         mTFacMom = 0         mTFacMom = 0
443           useBetaPlaneF = .TRUE.
        useConstantF  = .FALSE.  
        useBetaPlaneF = .FALSE.  
        useSphereF    = .TRUE.  
        omega         = 2. _d 0 * PI / ( 3600. _d 0 * 24. _d 0 )  
        usingSphericalPolarMterms = .TRUE.  
        metricTerms = .TRUE.  
        mTFacMom = 1  
444        ENDIF        ENDIF
445        IF ( usingSphericalPolarGrid ) THEN        IF ( usingSphericalPolarGrid ) THEN
446         useConstantF  = .FALSE.         useConstantF  = .FALSE.
# Line 436  C     Initialize EOS coefficients (3rd o Line 450  C     Initialize EOS coefficients (3rd o
450         usingSphericalPolarMterms = .TRUE.         usingSphericalPolarMterms = .TRUE.
451         metricTerms = .TRUE.         metricTerms = .TRUE.
452         mTFacMom = 1         mTFacMom = 1
   
453        ENDIF        ENDIF
454    
455  C--   Input files  C--   Input files
# Line 445  C--   Input files Line 458  C--   Input files
458        hydrogThetaFile = ' '        hydrogThetaFile = ' '
459        zonalWindFile   = ' '        zonalWindFile   = ' '
460        meridWindFile   = ' '        meridWindFile   = ' '
461          thetaClimFile   = ' '
462        READ(UNIT=iUnit,NML=PARM05,IOSTAT=errIO,err=11)        READ(UNIT=iUnit,NML=PARM05,IOSTAT=errIO,err=11)
463        IF ( errIO .GE. 0 ) GOTO 12        IF ( errIO .GE. 0 ) GOTO 12
464     11 CONTINUE     11 CONTINUE
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22