/[MITgcm]/MITgcm/pkg/ptracers/ptracers_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_readparms.F

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

revision 1.34 by jmc, Tue Apr 28 22:11:12 2009 UTC revision 1.35 by jahn, Fri Jun 26 23:10:10 2009 UTC
# Line 17  C     !USES: Line 17  C     !USES:
17        IMPLICIT NONE        IMPLICIT NONE
18  #include "SIZE.h"  #include "SIZE.h"
19  #include "EEPARAMS.h"  #include "EEPARAMS.h"
20    #ifdef ALLOW_LONGSTEP
21    #include "LONGSTEP_PARAMS.h"
22    #endif
23  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
24  #include "PTRACERS_PARAMS.h"  #include "PTRACERS_PARAMS.h"
25  #include "PARAMS.h"  #include "PARAMS.h"
# Line 129  C     Set defaults values for parameters Line 132  C     Set defaults values for parameters
132        PTRACERS_pickup_read_mnc  = .FALSE.        PTRACERS_pickup_read_mnc  = .FALSE.
133  #endif  #endif
134        tauTr1ClimRelax = 0.        tauTr1ClimRelax = 0.
135          DO k = 1,Nr
136    #ifdef ALLOW_LONGSTEP
137            PTRACERS_dTLev(k) = LS_nIter*dTtracerLev(k)
138    #else
139            PTRACERS_dTLev(k) = dTtracerLev(k)
140    #endif
141          ENDDO
142    
143  C     Open and read the data.ptracers file  C     Open and read the data.ptracers file
144        WRITE(msgBuf,'(A)') ' PTRACERS_READPARMS: opening data.ptracers'        WRITE(msgBuf,'(A)') ' PTRACERS_READPARMS: opening data.ptracers'
# Line 194  C     Check that enough parameters were Line 204  C     Check that enough parameters were
204          STOP 'ABNORMAL END: S/R PTRACERS_READPARMS'          STOP 'ABNORMAL END: S/R PTRACERS_READPARMS'
205         ENDIF         ENDIF
206        ENDDO        ENDDO
207          IF ( PTRACERS_dTLev(1).NE.PTRACERS_dTLev(Nr)
208         &     .AND. implicitDiffusion ) THEN
209            WRITE(msgBuf,'(A)')
210         &   'PTRACERS_READPARMS: #undef INCLUDE_IMPLVERTADV_CODE'
211            CALL PRINT_ERROR( msgBuf , myThid)
212            WRITE(msgBuf,'(2A)') 'PTRACERS_READPARMS: ',
213         &   'but implicitDiffusion=T with non-uniform PTRACERS_dTLev'
214            CALL PRINT_ERROR( msgBuf , myThid)
215            STOP 'ABNORMAL END: S/R PTRACERS_READPARMS'
216          ENDIF
217  #endif  #endif
218        DO iTracer=1,PTRACERS_numInUse        DO iTracer=1,PTRACERS_numInUse
219          IF ( PTRACERS_useGMRedi(iTracer) .AND. .NOT.useGMRedi ) THEN          IF ( PTRACERS_useGMRedi(iTracer) .AND. .NOT.useGMRedi ) THEN
# Line 274  C--   Print a summary of pTracer paramet Line 294  C--   Print a summary of pTracer paramet
294        CALL WRITE_0D_I( PTRACERS_Iter0, INDEX_NONE,        CALL WRITE_0D_I( PTRACERS_Iter0, INDEX_NONE,
295       &   'PTRACERS_Iter0 =',       &   'PTRACERS_Iter0 =',
296       &   ' /* timestep number when tracers are initialized */')       &   ' /* timestep number when tracers are initialized */')
297          CALL WRITE_1D_RL( PTRACERS_dTLev, Nr, INDEX_K,
298         &   'PTRACERS_dTLev =',
299         &'   /* Ptracer timestep ( s ) */')
300        CALL WRITE_0D_RL(PTRACERS_dumpFreq, INDEX_NONE,        CALL WRITE_0D_RL(PTRACERS_dumpFreq, INDEX_NONE,
301       &   'PTRACERS_dumpFreq =',       &   'PTRACERS_dumpFreq =',
302       &   ' /* Frequency^-1 for snapshot output (s) */')       &   ' /* Frequency^-1 for snapshot output (s) */')
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22