diff options
-rw-r--r-- | interface.ccl | 4 | ||||
-rw-r--r-- | src/MoL_Init.F90 | 27 |
2 files changed, 19 insertions, 12 deletions
diff --git a/interface.ccl b/interface.ccl index 4eae7cc..09ad541 100644 --- a/interface.ccl +++ b/interface.ccl @@ -11,6 +11,10 @@ CCTK_INT FUNCTION MoLRegisterEvolved(CCTK_INT IN EvolvedIndex, \ CCTK_INT IN RHSIndex) USES FUNCTION MoLRegisterEvolved +CCTK_INT FUNCTION MoLRegisterEvolvedGroup(CCTK_INT IN EvolvedIndex, \ + CCTK_INT IN RHSIndex) +USES FUNCTION MoLRegisterEvolvedGroup + public: CCTK_REAL level_set TYPE=GF TIMELEVELS=2 diff --git a/src/MoL_Init.F90 b/src/MoL_Init.F90 index bccc870..6a43ddd 100644 --- a/src/MoL_Init.F90 +++ b/src/MoL_Init.F90 @@ -10,7 +10,7 @@ subroutine EHFinder_MoLRegister(CCTK_ARGUMENTS) implicit none CCTK_INT :: ierr, ierr_cum, varindex, rhsindex - CCTK_INT :: MoLRegisterEvolved + CCTK_INT :: MoLRegisterEvolved, MoLRegisterEvolvedGroup DECLARE_CCTK_PARAMETERS DECLARE_CCTK_ARGUMENTS @@ -25,20 +25,23 @@ subroutine EHFinder_MoLRegister(CCTK_ARGUMENTS) if ( evolve_generators .gt. 0 ) then -! call CCTK_GroupIndex (varindex, 'ehfinder::generators') -! call CCTK_GroupIndex(rhsindex, 'ehfinder::generators_rhs') + call CCTK_GroupIndex (varindex, 'ehfinder::generators') + call CCTK_GroupIndex(rhsindex, 'ehfinder::generators_rhs') - call CCTK_VarIndex(varindex, 'ehfinder::xg') - call CCTK_VarIndex(rhsindex, 'ehfinder::dxg') - ierr_cum = ierr_cum + MoLRegisterEvolved(varindex, rhsindex) +!!$ call CCTK_VarIndex(varindex, 'ehfinder::xg') +!!$ call CCTK_VarIndex(rhsindex, 'ehfinder::dxg') +!!$ ierr_cum = ierr_cum + MoLRegisterEvolved(varindex, rhsindex) +!!$ +!!$ call CCTK_VarIndex(varindex, 'ehfinder::yg') +!!$ call CCTK_VarIndex(rhsindex, 'ehfinder::dyg') +!!$ ierr_cum = ierr_cum + MoLRegisterEvolved(varindex, rhsindex) +!!$ +!!$ call CCTK_VarIndex(varindex, 'ehfinder::zg') +!!$ call CCTK_VarIndex(rhsindex, 'ehfinder::dzg') +!!$ ierr_cum = ierr_cum + MoLRegisterEvolved(varindex, rhsindex) - call CCTK_VarIndex(varindex, 'ehfinder::yg') - call CCTK_VarIndex(rhsindex, 'ehfinder::dyg') - ierr_cum = ierr_cum + MoLRegisterEvolved(varindex, rhsindex) + ierr_cum = ierr_cum + MoLRegisterEvolvedGroup(varindex, rhsindex) - call CCTK_VarIndex(varindex, 'ehfinder::zg') - call CCTK_VarIndex(rhsindex, 'ehfinder::dzg') - ierr_cum = ierr_cum + MoLRegisterEvolved(varindex, rhsindex) end if ! call CCTK_VarIndex(varindex, 'admbase::gxx') |