[Users] Debugging patch for SetMask_SphericalSurface

-- Erik Schnetter schnetter@cct.lsu.edu   http://www.cct.lsu.edu/~eschnett/ $ svn diff Index: param.ccl =================================================================== --- param.ccl (revision 107) +++ param.ccl (working copy) @@ -2,7 +2,7 @@ STRING SetMask_MaskName "Full variable name of the mask to be set" { - ".+" :: "Any full Cactus variable name, GF of type CCTK_BYTE" + ".+" :: "Any full Cactus variable name, GF of type CCTK_INT" } "HydroBase::hydro_excision_mask" #BOOLEAN SetMask_ResetAll "'yes' (default): Overwrite complete mask, 'no': Only set excision points" Index: src/SetMask.c =================================================================== --- src/SetMask.c (revision 107) +++ src/SetMask.c (working copy) @@ -10,9 +10,18 @@ DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; - CCTK_INT *mask = (CCTK_INT*) CCTK_VarDataPtr(cctkGH, 0, SetMask_MaskName); + int varindex = CCTK_VarIndex(SetMask_MaskName); + if (varindex<0) + CCTK_WARN(0, "No such variable"); + int grouptype = CCTK_GroupTypeFromVarI(varindex); + if (grouptype!=CCTK_GF) + CCTK_WARN(0, "Wrong group type, must be GF"); + int vartype = CCTK_VarTypeI(varindex); + if (vartype!=CCTK_VARIABLE_INT) + CCTK_WARN(0, "Wrong variable type, must be CCTK_INT"); + CCTK_INT *mask = (CCTK_INT*) CCTK_VarDataPtrI(cctkGH, 0, varindex); if (!mask) - CCTK_WARN(0, "No such variable, or no storage enabled"); + CCTK_WARN(0, "No storage enabled"); for (int smi = 0; smi < 10; smi++) {