[Commits] [svn:einsteintoolkit] NullGrid/branches/tapir/ (Rev. 17)

bela at caltech.edu bela at caltech.edu
Tue Jan 7 16:38:56 CST 2014


User: szilagyi
Date: 2014/01/07 04:38 PM

Modified:
 /branches/tapir/
  schedule.ccl
 /branches/tapir/src/
  NullGrid_InitCoord.F90

Log:
 work on adaptive time-stepping algorithm

File Changes:

Directory: /branches/tapir/src/
===============================

File [modified]: NullGrid_InitCoord.F90
Delta lines: +8 -5
===================================================================
--- branches/tapir/src/NullGrid_InitCoord.F90	2014-01-03 07:23:01 UTC (rev 16)
+++ branches/tapir/src/NullGrid_InitCoord.F90	2014-01-07 22:38:56 UTC (rev 17)
@@ -32,7 +32,6 @@
   return
 end function NullGrid_RegisterCoords
 
-
 subroutine NullGrid_InitTime(CCTK_ARGUMENTS)
 
   use cctk
@@ -40,9 +39,9 @@
 
   DECLARE_CCTK_ARGUMENTS
 
-  null_time = cctk_time
-  null_time_p   = cctk_time -     cctk_delta_time
-  null_time_p_p = cctk_time - 2.0*cctk_delta_time
+  null_time     = 0
+  null_time_p   = 0
+  null_time_p_p = 0
 
 end subroutine NullGrid_InitTime
 
@@ -54,7 +53,11 @@
 
   DECLARE_CCTK_ARGUMENTS
 
-  null_time = cctk_time
+  null_time     = cctk_time
+  if(abs(null_time_p-null_time_p_p)<1.e-10*cctk_delta_time) then
+    null_time_p   = cctk_time -     cctk_delta_time
+    null_time_p_p = cctk_time - 2.0*cctk_delta_time
+  end if
 
 end subroutine NullGrid_CopyTime
 

Directory: /branches/tapir/
===========================

File [modified]: schedule.ccl
Delta lines: +9 -1
===================================================================
--- branches/tapir/schedule.ccl	2014-01-03 07:23:01 UTC (rev 16)
+++ branches/tapir/schedule.ccl	2014-01-07 22:38:56 UTC (rev 17)
@@ -28,12 +28,20 @@
   OPTIONS: global
 } "setup variables"
 
-schedule NullGrid_InitTime at CCTK_INITIAL
+schedule NullGrid_InitTime at CCTK_BASEGRID
 {
   LANG: Fortran
   OPTIONS: global
 } "initialize null time variable"
 
+schedule NullGrid_CopyTime at CCTK_INITIAL \
+  BEFORE NullEvol_BoundaryInit \
+  BEFORE NullEvol_Initial
+{
+  LANG: Fortran
+  OPTIONS: global
+} "copy value of time to local var"
+
 schedule NullGrid_CopyTime at CCTK_PRESTEP
 {
   LANG: Fortran



More information about the Commits mailing list