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")