[Commits] [svn:einsteintoolkit] NullExact/branches/tapir/src/ (Rev. 14)

[bela at caltech.edu]

User: szilagyi
Date: 2014/01/03 01:04 AM

Modified:
 /branches/tapir/src/
  NullExact_Error.F90

Log:
 remove references to some of the debugging variables

File Changes:

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

File [modified]: NullExact_Error.F90
Delta lines: +12 -10
===================================================================
--- branches/tapir/src/NullExact_Error.F90	2013-12-28 04:57:20 UTC (rev 13)
+++ branches/tapir/src/NullExact_Error.F90	2014-01-03 07:04:58 UTC (rev 14)
@@ -20,7 +20,9 @@
   ! local variables, as needed:
   CCTK_COMPLEX :: ii = (0,1.)
   CCTK_INT :: k, patch_ID
-  CCTK_COMPLEX, dimension(lsh(1), lsh(2), N_radial_pts) :: anajcn, anajcs, anaucn, anaucs, anadrucn, anadrucs, anaeth2jcn, anaeth2jcs, anaqcn, anaqcs
+  CCTK_COMPLEX, dimension(lsh(1), lsh(2), N_radial_pts) :: &
+    anajcn, anajcs, anaucn, anaucs, anadrucn, anadrucs, anaqcn, anaqcs
+  ! anaeth2jcn, anaeth2jcs,
 
 
   if (verbose.ne.0) call CCTK_INFO("checking for error in main complex evolution variables")
@@ -34,12 +36,12 @@
 !Compute analytic values in complex Characteristic Fields 
      call NullExact_Analytic_J_2D(patch_ID+1,null_lsh(1),null_lsh(2),stereo_q,stereo_p,null_xb(k),cctk_time,anajcn(:,:,k), Ylm_2)
      call NullExact_Analytic_U_2D(patch_ID+1,null_lsh(1),null_lsh(2),stereo_q,stereo_p,null_xbh(k),cctk_time,anaucn(:,:,k), Ylm_1)
-     call NullInterp_d2 (anaeth2jcn(:,:,k), anajcn(:,:,k), 2, 1, 1)
+     !call NullInterp_d2 (anaeth2jcn(:,:,k), anajcn(:,:,k), 2, 1, 1)
      call NullExact_Analytic_drU_2D(patch_ID+1,null_lsh(1),null_lsh(2),stereo_q,stereo_p,null_xb(k),cctk_time,anadrucn(:,:,k), Ylm_1)
 
      call NullExact_Analytic_J_2D(ip_s,null_lsh(1),null_lsh(2),stereo_q,stereo_p,null_xb(k),cctk_time,anajcs(:,:,k), Ylm_2)
      call NullExact_Analytic_U_2D(ip_s,null_lsh(1),null_lsh(2),stereo_q,stereo_p,null_xbh(k),cctk_time,anaucs(:,:,k), Ylm_1)
-     call NullInterp_d2 (anaeth2jcs(:,:,k), anajcs(:,:,k), 2, 1, 1)
+     !call NullInterp_d2 (anaeth2jcs(:,:,k), anajcs(:,:,k), 2, 1, 1)
      call NullExact_Analytic_drU_2D(ip_s,null_lsh(1),null_lsh(2),stereo_q,stereo_p,null_xb(k),cctk_time,anadrucs(:,:,k), Ylm_1)
 
 
@@ -57,8 +59,8 @@
      ucn_e(:,:,k) = EQ_mask * (ucn(:,:,k) - anaucn(:,:,k))
      ucs_e(:,:,k) = EQ_mask * (ucs(:,:,k) - anaucs(:,:,k))
 
-     eth2jcn_e(:,:,k) = EQ_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
-     eth2jcs_e(:,:,k) = EQ_mask * (eth2jcs(:,:,k) - anaeth2jcs(:,:,k))
+     ! eth2jcn_e(:,:,k) = EQ_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
+     ! eth2jcs_e(:,:,k) = EQ_mask * (eth2jcs(:,:,k) - anaeth2jcs(:,:,k))
 
      anaqcn(:,:,k) = (rwt*null_xb(k)/(1-null_xb(k)))**2&
                       *exp(-2.d0*bcn(:,:,k))*anadrucn(:,:,k)
@@ -74,31 +76,31 @@
 
        Mjcn_e(:,:,k) = (jcn(:,:,k) - anajcn(:,:,k))
        Mucn_e(:,:,k) = (ucn(:,:,k) - anaucn(:,:,k))
-       Meth2jcn_e(:,:,k) = (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
+       ! Meth2jcn_e(:,:,k) = (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
 
     else if (CCTK_EQUALS(error_mask_type, "EG_mask")) then
 
        Mjcn_e(:,:,k) = EG_mask * (jcn(:,:,k) - anajcn(:,:,k))
        Mucn_e(:,:,k) = EG_mask * (ucn(:,:,k) - anaucn(:,:,k))
-       Meth2jcn_e(:,:,k) = EG_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
+       ! Meth2jcn_e(:,:,k) = EG_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
 
     else if (CCTK_EQUALS(error_mask_type, "EQ_mask")) then
 
        Mjcn_e(:,:,k) = EQ_mask * (jcn(:,:,k) - anajcn(:,:,k))
        Mucn_e(:,:,k) = EQ_mask * (ucn(:,:,k) - anaucn(:,:,k))
-       Meth2jcn_e(:,:,k) = EQ_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
+       ! Meth2jcn_e(:,:,k) = EQ_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
 
     else if (CCTK_EQUALS(error_mask_type, "EV_mask")) then
 
        Mjcn_e(:,:,k) = EV_mask * (jcn(:,:,k) - anajcn(:,:,k))
        Mucn_e(:,:,k) = EV_mask * (ucn(:,:,k) - anaucn(:,:,k))
-       Meth2jcn_e(:,:,k) = EV_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
+       ! Meth2jcn_e(:,:,k) = EV_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
 
     else if (CCTK_EQUALS(error_mask_type, "guard_mask")) then
 
        Mjcn_e(:,:,k) = guard_mask * (jcn(:,:,k) - anajcn(:,:,k))
        Mucn_e(:,:,k) = guard_mask * (ucn(:,:,k) - anaucn(:,:,k))
-       Meth2jcn_e(:,:,k) = guard_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
+       ! Meth2jcn_e(:,:,k) = guard_mask * (eth2jcn(:,:,k) - anaeth2jcn(:,:,k))
 
     else
       call CCTK_WARN(0, "unsupported circular mask type")