--- MITgcm/model/src/ini_model_io.F 2007/10/17 22:08:13 1.23 +++ MITgcm/model/src/ini_model_io.F 2007/10/19 14:36:47 1.24 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_model_io.F,v 1.23 2007/10/17 22:08:13 jmc Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_model_io.F,v 1.24 2007/10/19 14:36:47 jmc Exp $ C $Name: $ #include "PACKAGES_CONFIG.h" @@ -10,7 +10,7 @@ C !INTERFACE: SUBROUTINE INI_MODEL_IO( myThid ) -C !DESCRIPTION: +C !DESCRIPTION: C Pass specific setup data for mdsio/rw. C !USES: @@ -19,6 +19,7 @@ #include "EEPARAMS.h" #include "EESUPPORT.h" #include "PARAMS.h" +#include "RESTART.h" C !INPUT/OUTPUT PARAMETERS: C myThid :: my Thread Id number @@ -49,12 +50,12 @@ C------ C GlobalFiles option with Multi-processors execution (with MPI) is not C safe: dependending on the platform & compiler, it may produce: -C - incomplete output files (wrong size) +C - incomplete output files (wrong size) C - wrong isolated values in some output files C - missing tiles (all zeros) in some output files. -C A safe alternative is to set "useSingleCpuIO=.TRUE." in file "data", +C A safe alternative is to set "useSingleCpuIO=.TRUE." in file "data", C namelist PARAM01 (and to keep the default value of globalFiles=FALSE) -C or if you are really sure that the globalFile works well on our platform +C or if you are really sure that the globalFile works well on our platform C & compiler, comment out the above "stop" C----- WRITE(msgBuf,'(2A)') @@ -71,11 +72,27 @@ C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| -C- Flags specific to RW and MDSIO - C- Only Master-thread updates IO-parameter in Common blocks: _BEGIN_MASTER( myThid ) +C- Initialise AB starting level +C notes: those could be modified when reading a pickup that does +C correspond to what is actually needed. + tempStartAB = nIter0 + saltStartAB = nIter0 + mom_StartAB = nIter0 + nHydStartAB = nIter0 + IF ( startFromPickupAB2 ) tempStartAB = MIN( nIter0 , 1 ) + saltStartAB = tempStartAB + mom_StartAB = tempStartAB + +C- Initialise Alternating pickup-suffix + nCheckLev = 1 + checkPtSuff(1) = 'ckptA' + checkPtSuff(2) = 'ckptB' + +C- Flags specific to RW and MDSIO + C- now we make local directories with myProcessStr appended IF ( mdsioLocalDir .NE. ' ' ) THEN pIL = ILNBLNK( mdsioLocalDir ) @@ -87,7 +104,7 @@ CALL SET_WRITE_GLOBAL_FLD( globalFiles ) C Set globalFiles flag for READ_WRITE_REC package CALL SET_WRITE_GLOBAL_REC( globalFiles ) -C Set globalFiles flag for READ_WRITE_PICKUP +C Set globalFiles flag for READ_WRITE_PICKUP CALL SET_WRITE_GLOBAL_PICKUP( globalFiles ) _END_MASTER( myThid )