/[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.75 by cnh, Thu Dec 13 17:27:53 2001 UTC revision 1.76 by jmc, Sun Feb 10 00:35:30 2002 UTC
# Line 127  C--   Continuous equation parameters Line 127  C--   Continuous equation parameters
127       & tempAdvScheme, saltAdvScheme, tracerAdvScheme,       & tempAdvScheme, saltAdvScheme, tracerAdvScheme,
128       & multiDimAdvection, useEnergyConservingCoriolis,       & multiDimAdvection, useEnergyConservingCoriolis,
129       & useJamartWetPoints,       & useJamartWetPoints,
130         & useRealFreshWaterFlux, convertFW2Salt,
131         & temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
132       & zonal_filt_lat       & zonal_filt_lat
133    
134  C--   Elliptic solver parameters  C--   Elliptic solver parameters
# Line 259  C--   Set default "physical" parameters Line 261  C--   Set default "physical" parameters
261        hFacMinDr           = UNSET_RL        hFacMinDr           = UNSET_RL
262        hFacMinDz           = UNSET_RL        hFacMinDz           = UNSET_RL
263        hFacMinDp           = UNSET_RL        hFacMinDp           = UNSET_RL
264          convertFW2Salt = UNSET_RL
265        READ(UNIT=iUnit,NML=PARM01) !,IOSTAT=errIO)        READ(UNIT=iUnit,NML=PARM01) !,IOSTAT=errIO)
266        IF ( errIO .LT. 0 ) THEN        IF ( errIO .LT. 0 ) THEN
267         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
# Line 349  C--   z,p,r coord input switching. Line 352  C--   z,p,r coord input switching.
352        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDp        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDp
353        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDrDefault        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDrDefault
354    
355          IF (convertFW2Salt.EQ.UNSET_RL) THEN
356            convertFW2Salt = 35.
357            IF (useRealFreshWaterFlux) convertFW2Salt=-1
358          ENDIF
359    
360        IF ( ivdc_kappa .NE. 0. .AND. .NOT. implicitDiffusion ) THEN        IF ( ivdc_kappa .NE. 0. .AND. .NOT. implicitDiffusion ) THEN
361         WRITE(msgBuf,'(A,A)')         WRITE(msgBuf,'(A,A)')
362       &  'S/R INI_PARMS: To use ivdc_kappa you must enable implicit',       &  'S/R INI_PARMS: To use ivdc_kappa you must enable implicit',
# Line 356  C--   z,p,r coord input switching. Line 364  C--   z,p,r coord input switching.
364         CALL PRINT_ERROR( msgBuf , myThid)         CALL PRINT_ERROR( msgBuf , myThid)
365         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
366        ENDIF        ENDIF
   
       IF ( implicitFreeSurface .AND.  rigidLid ) THEN  
        WRITE(msgBuf,'(A,A)')  
      & 'S/R INI_PARMS: Cannot select both implicitFreeSurface',  
      & ' and rigidLid.'  
        CALL PRINT_ERROR( msgBuf , myThid)  
        STOP 'ABNORMAL END: S/R INI_PARMS'  
       ENDIF  
       IF ( (implicSurfPress.NE.1. .OR. implicDiv2DFlow.NE.1.)  
      &    .AND. nonHydrostatic ) THEN  
        WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: nonHydrostatic',  
      & ' NOT SAFE with non-fully implicit Barotropic solver'  
        CALL PRINT_ERROR( msgBuf , myThid)  
        WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: To by-pass this',  
      &    'STOP, comment this test and re-compile ini_params'  
        CALL PRINT_ERROR( msgBuf , myThid)  
        STOP 'ABNORMAL END: S/R INI_PARMS'  
       ENDIF  
367    
368        coordsSet = 0        coordsSet = 0
369        IF ( zCoordInputData ) coordsSet = coordsSet + 1        IF ( zCoordInputData ) coordsSet = coordsSet + 1
Legend:
Removed from v.1.75  
changed lines
  Added in v.1.76
  ViewVC Help
Powered by ViewVC 1.1.22