/[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.10 by cnh, Mon May 25 21:29:45 1998 UTC revision 1.11 by adcroft, Thu May 28 16:19:50 1998 UTC
# Line 56  C--   Continuous equation parameters Line 56  C--   Continuous equation parameters
56       &  diffKhT, diffKzT, diffK4T,       &  diffKhT, diffKzT, diffK4T,
57       &  diffKhS, diffKzS, diffK4S,       &  diffKhS, diffKzS, diffK4S,
58       &  GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,       &  GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,
59       &  tRef, sRef,       &  tRef, sRef, eosType,
60       & momViscosity,  momAdvection, momForcing, useCoriolis,       & momViscosity,  momAdvection, momForcing, useCoriolis,
61       & momPressureForcing,       & momPressureForcing,
62       & tempDiffusion, tempAdvection, tempForcing,       & tempDiffusion, tempAdvection, tempForcing,
# Line 148  C--   Set default "physical" parameters Line 148  C--   Set default "physical" parameters
148        gravity =   9.81 D0        gravity =   9.81 D0
149        gBaro   = gravity        gBaro   = gravity
150        rhoNil   = 999.8 D0        rhoNil   = 999.8 D0
       tAlpha=2.d-4  
       sBeta=1.d-5  
151        f0=1.D-4        f0=1.D-4
152        beta = 1. _d -11        beta = 1. _d -11
153        viscAh=1.d3        viscAh=1.d3
# Line 166  C--   Set default "physical" parameters Line 164  C--   Set default "physical" parameters
164        GMalpha=0.        GMalpha=0.
165        GMdepth=1000.        GMdepth=1000.
166        GMkbackground=0.        GMkbackground=0.
167          tAlpha=2.d-4
168          sBeta=1.d-5
169          eosType='POLY3'
170        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
171        rigidLid            = .FALSE.        rigidLid            = .FALSE.
172        freeSurfFac         = 1. _d 0        freeSurfFac         = 1. _d 0
# Line 357  C     In spherical polar distances are i Line 358  C     In spherical polar distances are i
358         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
359    10  CONTINUE    10  CONTINUE
360    
361    C     Initialize EOS coefficients (3rd order polynomial)
362          IF (eostype.eq.'POLY3') THEN
363           OPEN(37,FILE='POLY3.COEFFS',STATUS='OLD',FORM='FORMATTED')
364           READ(37,*) I
365           IF (I.NE.Nz) THEN
366            WRITE(0,*) 'ini_parms: attempt to read POLY3.COEFFS failed'
367            WRITE(0,*) '           because bad # of levels in data'
368            STOP 'Bad data in POLY3.COEFFS'
369           ENDIF
370           READ(37,*) (eosRefT(K),eosRefS(K),eosSig0(K),K=1,Nz)
371           DO K=1,Nz
372            READ(37,*) (eosC(I,K),I=1,9)
373            write(0,'(i3,13f8.3)') K,eosRefT(K),eosRefS(K),eosSig0(K),
374         &                (eosC(I,K),I=1,9)
375           ENDDO
376           CLOSE(37)
377          ENDIF
378    
379        goptCount = 0        goptCount = 0
380        IF ( usingCartesianGrid )      goptCount = goptCount+1        IF ( usingCartesianGrid )      goptCount = goptCount+1
381        IF ( usingSphericalPolarGrid ) goptCount = goptCount+1        IF ( usingSphericalPolarGrid ) goptCount = goptCount+1
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22