Added spinup phase variables

src/DOUAR.f90

@@ -246,7 +246,7 @@ call compute_plastic_params (params,mat)
 call show_time (total,step,inc,1,'define velocity octree$')
 call define_ov (ov,params,threadinfo)
 
-istep=1+params%irestart 
+istep=1+params%irestart
 
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
@@ -254,7 +254,7 @@ istep=1+params%irestart
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
 
-usecourant = (params%dt .lt. 0.d0)
+!usecourant = (params%dt .lt. 0.d0)
 
 !------------------------------------------------------------------------------|
 !------------------------------------------------------------------------------|
@@ -263,7 +263,29 @@ usecourant = (params%dt .lt. 0.d0)
 !------------------------------------------------------------------------------|
 increase_current_level=.false.
 do while (istep.le.params%nstep) 
-   
+
+   !---------------------------------------------------------------------------|
+   !---------------------------------------------------------------------------|
+   !     Are we in the main calculation phase or the spinup phase?
+   !---------------------------------------------------------------------------|
+   !---------------------------------------------------------------------------|
+   if (istep.le.params%nstep_spinup) then
+      params%dt=params%dt_spinup
+      params%griditer=params%griditer_spinup
+      params%nonlinear_iterations=params%nonlinear_iterations_spinup
+   else
+      params%dt=params%dt_main
+      params%griditer=params%griditer_main
+      params%nonlinear_iterations=params%nonlinear_iterations_main
+   endif
+
+   !---------------------------------------------------------------------------|
+   !---------------------------------------------------------------------------|
+   !     is the CFL condition used for the timestep ?
+   !---------------------------------------------------------------------------|
+   !---------------------------------------------------------------------------|
+   usecourant = (params%dt .lt. 0.d0)
+
    call write_streepjes(6,2)
    call write_streepjes(8,2)
    call show_time (total,step,inc,-istep,'Start of Step $')
@@ -967,7 +989,6 @@ do while (istep.le.params%nstep)
       ! finds courant condition time step in case this is the first iteration
       !--------------------------------------------------------------------------------------
       !--------------------------------------------------------------------------------------
-
       if (iter.eq.1) then
          umax=0.d0
          do i=1,ov%nnode