John Wise  committed fc40bd0 Merge

Merging in some more changes I forgot on the GT cluster.

  • Participants
  • Parent commits ccdc9c5, 599f469

Comments (0)

Files changed (2)

File src/enzo/Grid_CopyZonesFromGrid.C

   /* declarations */
-  int i,j,k,field,dim;
+  int i, j, k, dim, loop, iLoop, field, otherindexB;
+  float vx, vy, vz, v2, rho, bx, by, bz, b2;
   bool shiftPos, shiftNeg; float delta; FLOAT L;
   if (ShearingBoundaryDirection!=-1){
-    bool noMove=false;
     //printf("L: %"GSYM" Delta: %"GSYM" %"GSYM" (%"GSYM" %"GSYM")\n", L, delta, delta, AngularVelocity, VelocityGradient);
   //Shearing Boundary Variables
-  float rho, vx, vy, vz, v2, b2, bx, by, bz=0.0;
   int thisindex, otherindex=0;
   FLOAT a,b;  FLOAT val1=-9999;FLOAT val2=-9999;
   int addDim[3] = {1, OtherDim[0], OtherDim[0]*OtherDim[1]};
   int velocityTypes[3]={Velocity1, Velocity2, Velocity3};
-  int Zero[3] = {0,0,0};
   if (!isShearing)
-    //#pragma omp parallel for schedule(static)
+//#pragma omp parallel for schedule(static)
     for (field = 0; field < NumberOfBaryonFields; field++)
       FORTRAN_NAME(copy3drel)(OtherGrid->BaryonField[field], BaryonField[field],
 			      Dim, Dim+1, Dim+2,
 	    (k + StartOther[2])*OtherDim[0]*OtherDim[1];
 	  for (i = 0; i < Dim[0]; i++, thisindex++, otherindex++){
-	    int otherindexB=otherindex+ addDim[ShearingVelocityDirection];
+	    otherindexB=otherindex+ addDim[ShearingVelocityDirection];
 	    if (DualEnergyFormalism==0 && FieldType[field]==TotalEnergy) {
-	      for (int loop=0; loop<=1; loop++){
-		int iLoop=otherindex;
+	      for (loop=0; loop < 2; loop++){
+		iLoop=otherindex;
 		if (loop==1) iLoop= otherindexB;
-		float vx, vy, vz, v2, rho;
 		rho= OtherGrid->BaryonField[iden][iLoop];
 		vx= OtherGrid->BaryonField[ivx][iLoop];
 		else vz=0.0;
-		float bx, by, bz, b2;
 		if (useMHD) {
 		  bx= OtherGrid->BaryonField[iBx][iLoop];
 	thisindex = (0 + Start[0]) + (j + Start[1])*GridDimension[0] +
 	  (k + Start[2])*GridDimension[0]*GridDimension[1];
 	for (i = 0; i < Dim[0]; i++, thisindex++){
-	  float vx, vy, vz, v2, rho, bx, by, bz, b2;
 	  rho= BaryonField[iden][thisindex];  
 	  vx= BaryonField[ivx][thisindex];
 	  vy= BaryonField[ivy][thisindex];  

File src/enzo/Grid_WalkPhotonPackage.C

       return SUCCESS;
     // return in case we're out of photons
-    if ((*PP)->Photons < tiny_number || 
+    if ((*PP)->Photons < 1e-10*tiny_number || 
 	(*PP)->ColumnDensity > tau_delete) {
       if (DEBUG>1) {
 	fprintf(stderr, "PP-Photons: %"GSYM"  PP->Radius: %"GSYM