From e7ac99eecabbe517b1abdd17bc2bd67636ecf804 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Wed, 25 Jan 2012 09:05:38 +0000 Subject: [PATCH] for spectral solver, --restart XX sets restartReadInc to XX-1, meaning restartReadInc is the step to read and XX the step at which the calculation begins --- code/FEsolving.f90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/code/FEsolving.f90 b/code/FEsolving.f90 index 7d1d744e5..4c60eafe3 100644 --- a/code/FEsolving.f90 +++ b/code/FEsolving.f90 @@ -25,13 +25,13 @@ use prec, only: pInt,pReal implicit none - integer(pInt) :: cycleCounter = 0_pInt, theInc = -1_pInt, restartReadInc = 1_pInt + integer(pInt) :: cycleCounter = 0_pInt, theInc = -1_pInt, restartReadInc = 0_pInt real(pReal) :: theTime = 0.0_pReal, theDelta = 0.0_pReal logical :: lastIncConverged = .false.,outdatedByNewInc = .false.,outdatedFFN1 = .false.,terminallyIll = .false. logical :: symmetricSolver = .false. logical :: parallelExecution = .true. logical :: restartWrite = .false. - logical :: restartRead = .false., restartReadSpectral = .false. + logical :: restartRead = .false. logical :: lastMode = .true., cutBack = .false. logical, dimension(:,:), allocatable :: calcMode integer(pInt), dimension(:,:), allocatable :: FEsolving_execIP @@ -75,6 +75,7 @@ if(start /= 0_pInt) then ! found something length = verify(commandLine(start:len(commandLine)),'0123456789',.false.) ! where is first non number after argument? read(commandLine(start:start+length),'(I12)') restartReadInc ! read argument + restartReadInc = restartReadInc - 1_pInt ! command line argument is inc to compute restartRead = max(0_pInt,restartReadInc) > 0_pInt if(restartReadInc < 0_pInt) call IO_warning(warning_ID=34_pInt) endif