Time stepper routines

examples/coupled_heat_conduction/rod12.usr

@@ -14,10 +14,7 @@ c-----------------------------------------------------------------------
       include 'SIZE'
       include 'TOTAL'
       include 'NEKUSE'
-
       integer e,f,eg
-c     e = gllel(eg)
-
       udiff =0.
       utrans=0.
       return
@@ -27,17 +24,10 @@ c-----------------------------------------------------------------------
       include 'SIZE'
       include 'TOTAL'
       include 'NEKUSE'
-
       integer e,f,eg
-c     e = gllel(eg)
-
-
-c     Note: this is an acceleration term, NOT a force!
-c     Thus, ffx will subsequently be multiplied by rho(x,t).
       ffx=0.0
       ffy=0.0
       ffz=0.0     
-
       return
       end
 c-----------------------------------------------------------------------
@@ -95,8 +85,6 @@ c-----------------------------------------------------------------------
       temp_nek=sint1/sarea1
       
       if (nid.eq.0) then
-c         write(6,*)"*** Pressure: ",pmin," - ",pmax
-c         write(6,*)"*** Velocity: ",wmin," - ",wmax
          write(6,*)"*** Temperature: ",tmin," - ",tmax
          write(6,*)"*** Av. Temperature: ",temp_nek
       endif
@@ -175,3 +163,87 @@ c
       return
       end
 
+C=======================================================================
+      subroutine nek_init_step()
+      include 'SIZE'
+      include 'TSTEP'
+      include 'INPUT'
+      include 'CTIMER'
+      real*4 papi_mflops
+      integer*8 papi_flops
+      integer icall, kstep, i, pstep
+      common /cht_coupler/ pstep
+      save icall
+ 
+      if (icall.eq.0) then
+      call nekgsync()
+      if (instep.eq.0) then
+        if(nid.eq.0) write(6,'(/,A,/,A,/)')
+     &     ' nsteps=0 -> skip time loop',
+     &     ' running solver in post processing mode'
+      else
+        if(nio.eq.0) write(6,'(/,A,/)') 'Starting time loop ...'
+      endif
+      isyc  = 0
+      itime = 0
+      if(ifsync) isyc=1
+      itime = 1
+      call nek_comm_settings(isyc,itime)
+      call nek_comm_startstat()
+      istep  = 0
+      endif
+
+      istep=istep+1
+
+      if (lastep .eq. 1) then
+        pstep=2
+      else
+        call nek_advance
+        pstep=2
+      endif
+
+      return
+      end
+C=======================================================================
+      subroutine nek_step()
+
+      include 'SIZE'
+      include 'TSTEP'
+      include 'INPUT'
+      include 'CTIMER'
+      common /cht_coupler/ pstep
+      integer pstep
+
+      pstep=pstep+1
+      call heat(pstep)
+
+      return
+      end
+C=======================================================================
+      subroutine nek_finalize_step()
+
+      include 'SIZE'
+      include 'TSTEP'
+      include 'INPUT'
+      include 'CTIMER'
+      common /cht_coupler/ pstep
+      integer pstep 
+
+      real*4 papi_mflops
+      integer*8 papi_flops
+      integer icall, kstep, knstep, i
+
+      if (param(103).gt.0)   call q_filter      (param(103))
+      call setup_convect (pstep)  ! Save convective velocity _after_ filter
+
+      call userchk
+      call prepost (.false.,'his')
+      call in_situ_check()
+
+      if (mod(istep,nsteps).eq.0) lastep=1
+      call nek_comm_settings(isyc,0)
+      call comment
+
+      return
+      end
+c-----------------------------------------------------------------------