/[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.7 by cnh, Thu May 21 18:25:49 1998 UTC revision 1.8 by cnh, Mon May 25 16:17:36 1998 UTC
# Line 51  C     K, I, J - Loop counters Line 51  C     K, I, J - Loop counters
51    
52  C--   Continuous equation parameters  C--   Continuous equation parameters
53        NAMELIST /PARM01/        NAMELIST /PARM01/
54       & gravity, rhonil, tAlpha, sBeta, f0, beta,       & gravity, gBaro, rhonil, tAlpha, sBeta, f0, beta,
55       &   viscAh,  viscAz,  viscA4,       &   viscAh,  viscAz,  viscA4,
56       &  diffKhT, diffKzT, diffK4T,       &  diffKhT, diffKzT, diffK4T,
57       &  diffKhS, diffKzS, diffK4S,       &  diffKhS, diffKzS, diffK4S,
# Line 59  C--   Continuous equation parameters Line 59  C--   Continuous equation parameters
59       &  tRef, sRef,       &  tRef, sRef,
60       & momViscosity,  momAdvection, momForcing, useCoriolis,       & momViscosity,  momAdvection, momForcing, useCoriolis,
61       & tempDiffusion, tempAdvection, tempForcing,       & tempDiffusion, tempAdvection, tempForcing,
62       & saltDiffusion, saltAdvection, saltForcing       & saltDiffusion, saltAdvection, saltForcing,
63         & implicitFreeSurface, rigidLid, freeSurfFac
64    
65  C--   Elliptic solver parameters  C--   Elliptic solver parameters
66        NAMELIST /PARM02/        NAMELIST /PARM02/
# Line 142  C--   Set default "physical" parameters Line 143  C--   Set default "physical" parameters
143         tRef(K) = 30.D0 - FLOAT(K)         tRef(K) = 30.D0 - FLOAT(K)
144        ENDDO        ENDDO
145        gravity =   9.81 D0        gravity =   9.81 D0
146          gBaro   = gravity
147        rhoNil   = 999.8 D0        rhoNil   = 999.8 D0
148        tAlpha=2.d-4        tAlpha=2.d-4
149        sBeta=1.d-5        sBeta=1.d-5
# Line 161  C--   Set default "physical" parameters Line 163  C--   Set default "physical" parameters
163        GMalpha=0.        GMalpha=0.
164        GMdepth=1000.        GMdepth=1000.
165        GMkbackground=0.        GMkbackground=0.
166          implicitFreeSurface = .TRUE.
167          rigidLid            = .FALSE.
168          freeSurfFac         = 1. _d 0
169        READ(UNIT=iUnit,NML=PARM01,IOSTAT=errIO,err=3)        READ(UNIT=iUnit,NML=PARM01,IOSTAT=errIO,err=3)
170        IF ( errIO .GE. 0 ) GOTO 4        IF ( errIO .GE. 0 ) GOTO 4
171      3 CONTINUE      3 CONTINUE
# Line 179  C--   Set default "physical" parameters Line 184  C--   Set default "physical" parameters
184         CALL MODELDATA_EXAMPLE( myThid )         CALL MODELDATA_EXAMPLE( myThid )
185         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
186     4  CONTINUE     4  CONTINUE
187          IF ( implicitFreeSurface ) freeSurfFac = 1. _d 0
188          IF ( rigidLid            ) freeSurfFac = 0. _d 0
189    
190          IF ( implicitFreeSurface .AND.  rigidLid ) THEN
191           WRITE(msgBuf,'(A)')
192         &  'S/R INI_PARMS: Cannot select implicitFreeSurface and rigidLid.'
193           CALL PRINT_ERROR( msgBuf , myThid)
194           STOP 'ABNORMAL END: S/R INI_PARMS'
195          ENDIF
196    
197  C--   Elliptic solver parameters  C--   Elliptic solver parameters
198        cg2dMaxIters   = 150        cg2dMaxIters   = 150
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22