Boundary routines:

reworked Robin to have conformal calling interfaces
fixed two bugs: the more interresting one: coordinates do not have timelevels.

Replaced most CCTK_WARNS with CCTK_VWarns to be more informative.

gerd


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@115 6a38eb6e-646e-4a02-a296-d141613ad6c4

1 files changed, 76 insertions, 6 deletions

diff --git a/src/RadiationBoundary.c b/src/RadiationBoundary.c
index 4ba041d..433634d 100644
--- a/src/RadiationBoundary.c
+++ b/src/RadiationBoundary.c
@@ -64,6 +64,19 @@ static int ApplyBndRadiative(cGH *GH,
  ********************    External Routines   ************************
  ********************************************************************/
 
+/*@@
+   @routine    BndRadiativeGI
+   @date       Tue Jul 18 18:04:07 2000
+   @author     Gerd Lanfermann
+   @desc 
+      Aply radiative BC's by group index
+   @enddesc 
+   @calls     
+   @calledby   
+   @history 
+   @endhistory 
+
+@@*/
 
 int BndRadiativeGI(cGH *GH, 
 		   int *sw,
@@ -89,7 +102,18 @@ void CCTK_FCALL CCTK_FNAME(BndRadiativeGI)
   *ierr = BndRadiativeGI(GH, sw, *var0, *v0, *gi, *gi_p);
 }
 
+/*@@
+   @routine    BndRadiativeGI
+   @date       Tue Jul 18 18:04:07 2000
+   @author     Gerd Lanfermann
+   @desc 
+      Aply radiative BC's by group name   @enddesc 
+   @calls     
+   @calledby   
+   @history 
+   @endhistory 
 
+@@*/
 
 int BndRadiativeGN(cGH *GH, 
 		   int *sw, 
@@ -126,7 +150,19 @@ void CCTK_FCALL CCTK_FNAME(BndRadiativeGN)
   return;
 }
 
+/*@@
+   @routine    BndRadiativeGI
+   @date       Tue Jul 18 18:04:07 2000
+   @author     Gerd Lanfermann
+   @desc 
+      Apply radiative BC's by variable index 
+   @enddesc 
+   @calls     
+   @calledby   
+   @history 
+   @endhistory 
 
+@@*/
 
 int BndRadiativeVI(cGH *GH, 
                     int *sw, 
@@ -146,6 +182,19 @@ void CCTK_FCALL CCTK_FNAME(BndRadiativeVI)
   return;
 }
 
+/*@@
+   @routine    BndRadiativeGI
+   @date       Tue Jul 18 18:04:07 2000
+   @author     Gerd Lanfermann
+   @desc 
+      Apply radiative BC's by variable name 
+   @enddesc 
+   @calls     
+   @calledby   
+   @history 
+   @endhistory 
+
+@@*/
 
 int BndRadiativeVN(cGH *GH, 
                    int *sw,
@@ -190,6 +239,19 @@ void CCTK_FCALL CCTK_FNAME(BndRadiativeVN)
 /********************************************************************
  *********************     Local Routines   *************************
  ********************************************************************/
+/*@@
+   @routine    BndRadiativeGI
+   @date       Tue Jul 18 18:04:07 2000
+   @author     Gerd Lanfermann
+   @desc 
+      Apply radiative BC's - internal routine
+   @enddesc 
+   @calls     
+   @calledby   
+   @history 
+   @endhistory 
+
+@@*/
 
 static int ApplyBndRadiative(cGH *GH, 
 			     int *sw,
@@ -245,7 +307,9 @@ static int ApplyBndRadiative(cGH *GH,
       ri = CCTK_CoordIndex(-1,"r","spher3d");
       break;
     default:
-      CCTK_WARN(2,"ApplyBndRadiative: Radiative boundary only for dim=3");
+      CCTK_VWarn(1,__LINE__,__FILE__,"Boundary",
+		 "Radiative boundary only for dim=3: grid varible '%s'",
+		 CCTK_VarName(var_current));
       break;
   }
 
@@ -313,13 +377,17 @@ static int ApplyBndRadiative(cGH *GH,
     {
       case 1: 
 	berr = -1; 
-	CCTK_WARN(0,"ApplyBndRadiative: Radiative boundary conditions" 
-                    "not implemented for dimension = 1");
+	CCTK_VWarn(1,__LINE__,__FILE__,"Boundary",
+		   "ApplyBndRadiative: Radiative boundary conditions" 
+		   "not implemented for dimension = 1, grid variable '%s'",
+		   CCTK_VarName(vi1));
 	break;
       case 2: 
 	berr = -1; 
-	CCTK_WARN(0,"ApplyBndRadiative: Radiative boundary conditions"
-                    "not implemented for dimension = 2");
+	CCTK_VWarn(1,__LINE__,__FILE__,"Boundary",
+		   "ApplyBndRadiative: Radiative boundary conditions"
+		   "not implemented for dimension = 2, grid variable '%s'",
+		   CCTK_VarName(vi1));
 	break;
       case 3:   
 	/* According to the Cactus spec, the true dx and dt values for a
@@ -348,7 +416,9 @@ static int ApplyBndRadiative(cGH *GH,
 	break;
       default : 
 	berr=-1; 
-	CCTK_WARN(1, "ApplyBndRadiative: No Radiative BC for dim>3");
+	CCTK_VWarn(1,__LINE__,__FILE__,"Boundary",
+		   "No Radiative BC for dim>3, grid variable '%s'",
+		   CCTK_VarName(vi1));
     }
     berr=(berr<0)?-1:0;
   }