mirror of
https://github.com/cosmo-sims/MUSIC.git
synced 2024-09-19 17:03:46 +02:00
added explicit coarsening of all fields to increase stability when using new masking
This commit is contained in:
parent
cf9df8e6b7
commit
e9ffb9b953
1 changed files with 12 additions and 4 deletions
10
main.cc
10
main.cc
|
@ -640,6 +640,7 @@ int main (int argc, const char * argv[])
|
|||
//... displacement
|
||||
the_poisson_solver->gradient(icoord, u, data_forIO );
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing CDM displacements");
|
||||
the_output_plugin->write_dm_position(icoord, data_forIO );
|
||||
}
|
||||
|
@ -694,6 +695,7 @@ int main (int argc, const char * argv[])
|
|||
//... displacement
|
||||
the_poisson_solver->gradient(icoord, u, data_forIO );
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing baryon displacements");
|
||||
the_output_plugin->write_gas_position(icoord, data_forIO );
|
||||
|
||||
|
@ -781,6 +783,7 @@ int main (int argc, const char * argv[])
|
|||
double sigv = compute_finest_sigma( data_forIO );
|
||||
LOGINFO("sigma of %c-velocity of high-res particles is %f",'x'+icoord, sigv);
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing CDM velocities");
|
||||
the_output_plugin->write_dm_velocity(icoord, data_forIO);
|
||||
|
||||
|
@ -790,7 +793,6 @@ int main (int argc, const char * argv[])
|
|||
the_output_plugin->write_gas_velocity(icoord, data_forIO);
|
||||
}
|
||||
|
||||
|
||||
}
|
||||
|
||||
u.deallocate();
|
||||
|
@ -857,6 +859,7 @@ int main (int argc, const char * argv[])
|
|||
double sigv = compute_finest_sigma( data_forIO );
|
||||
LOGINFO("sigma of %c-velocity of high-res DM is %f",'x'+icoord, sigv);
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing CDM velocities");
|
||||
the_output_plugin->write_dm_velocity(icoord, data_forIO);
|
||||
}
|
||||
|
@ -908,6 +911,7 @@ int main (int argc, const char * argv[])
|
|||
double sigv = compute_finest_sigma( data_forIO );
|
||||
LOGINFO("sigma of %c-velocity of high-res baryons is %f",'x'+icoord, sigv);
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing baryon velocities");
|
||||
the_output_plugin->write_gas_velocity(icoord, data_forIO);
|
||||
}
|
||||
|
@ -1057,6 +1061,7 @@ int main (int argc, const char * argv[])
|
|||
double sigv = compute_finest_sigma( data_forIO );
|
||||
std::cerr << " - velocity component " << icoord << " : sigma = " << sigv << std::endl;
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing CDM velocities");
|
||||
the_output_plugin->write_dm_velocity(icoord, data_forIO);
|
||||
|
||||
|
@ -1144,6 +1149,7 @@ int main (int argc, const char * argv[])
|
|||
double sigv = compute_finest_sigma( data_forIO );
|
||||
std::cerr << " - velocity component " << icoord << " : sigma = " << sigv << std::endl;
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing baryon velocities");
|
||||
the_output_plugin->write_gas_velocity(icoord, data_forIO);
|
||||
}
|
||||
|
@ -1238,6 +1244,7 @@ int main (int argc, const char * argv[])
|
|||
else
|
||||
the_poisson_solver->gradient(icoord, u1, data_forIO );
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing CDM displacements");
|
||||
the_output_plugin->write_dm_position(icoord, data_forIO );
|
||||
}
|
||||
|
@ -1347,6 +1354,7 @@ int main (int argc, const char * argv[])
|
|||
else
|
||||
the_poisson_solver->gradient(icoord, u1, data_forIO );
|
||||
|
||||
coarsen_density( rh_Poisson, data_forIO, false );
|
||||
LOGUSER("Writing baryon displacements");
|
||||
the_output_plugin->write_gas_position(icoord, data_forIO );
|
||||
}
|
||||
|
|
Loading…
Reference in a new issue