Commits

Anonymous committed cba9e06 Merge

merge

Comments (0)

Files changed (2)

src/enzo/ClusterSMBHSumGasMass.C

 
 float ClusterSMBHColdGasMass;
 
-FILE *fptr=fopen("MT.out","w");
+FILE *fptr=fopen("MT.out","a");
 int ClusterSMBHSumGasMass(HierarchyEntry *Grids[], int NumberOfGrids, int level)
 {
 

src/enzo/Grid_ClusterInitializeGrid.C

 
   FILE *fptr = fopen("NFWProfile.out", "w");
   for (i = 0; i < NFW_POINTS; i++) {
-    NFWRadius[i] = 2.0*SphereRadius[sphere]*pow(10, -5*(float(i)/NFW_POINTS));
+    //NFWRadius[i] = 2.0*SphereRadius[sphere]*pow(10, -5*(float(i)/NFW_POINTS));
+    NFWRadius[i] = 0.5*pow(10, -5*(float(i)/NFW_POINTS));
     x1 = NFWRadius[i]/SphereCoreRadius[sphere];
   if (SphereType[sphere]!=6) {
      NFWDensity[i] = SphereDensity[sphere]/(x1*(1.0+x1)*(1.0+x1));
                    pow(fabs(z-SpherePosition[sphere][2]), 2) );
           r = max(r, 0.1*CellWidth[0][0]);
 
-          if (r < 2.0*SphereRadius[sphere]) {
-          //if (r < 0.45) {
+          //if (r < 2.0*SphereRadius[sphere]) {
+          if (r < 0.48) {
 
             /* 3) NFW profile (use look-up table for temperature and
                   velocity dispersion)*/