1
0
Fork 0
mirror of https://github.com/cosmo-sims/MUSIC.git synced 2024-09-19 17:03:46 +02:00

Removed option "shift_back" for Gadget2 output as it can cause problems with zoom sims overlapping the periodic boundaries.

This commit is contained in:
Oliver Hahn 2011-08-29 14:00:47 -07:00
parent 8ee0772b81
commit 29822726ec

View file

@ -918,7 +918,7 @@ public:
//... determine if we need to shift the coordinates back //... determine if we need to shift the coordinates back
double *shift = NULL; double *shift = NULL;
if( cf_.getValueSafe<bool>("output","shift_back",false ) ) /* if( cf_.getValueSafe<bool>("output","shift_back",false ) )
{ {
if( coord == 0 ) if( coord == 0 )
std::cout << " - gadget2 output plug-in will shift particle positions back...\n"; std::cout << " - gadget2 output plug-in will shift particle positions back...\n";
@ -929,7 +929,7 @@ public:
shift[1] = -(double)cf_.getValue<int>( "setup", "shift_y" )*h; shift[1] = -(double)cf_.getValue<int>( "setup", "shift_y" )*h;
shift[2] = -(double)cf_.getValue<int>( "setup", "shift_z" )*h; shift[2] = -(double)cf_.getValue<int>( "setup", "shift_z" )*h;
} }
*/
size_t npart = npfine+npcoarse; size_t npart = npfine+npcoarse;
size_t nwritten = 0; size_t nwritten = 0;
@ -1171,7 +1171,7 @@ public:
//... determine if we need to shift the coordinates back //... determine if we need to shift the coordinates back
double *shift = NULL; double *shift = NULL;
if( cf_.getValueSafe<bool>("output","shift_back",false ) ) /*if( cf_.getValueSafe<bool>("output","shift_back",false ) )
{ {
if( coord == 0 ) if( coord == 0 )
std::cout << " - gadget2 output plug-in will shift particle positions back...\n"; std::cout << " - gadget2 output plug-in will shift particle positions back...\n";
@ -1181,7 +1181,7 @@ public:
shift[0] = -(double)cf_.getValue<int>( "setup", "shift_x" )*h; shift[0] = -(double)cf_.getValue<int>( "setup", "shift_x" )*h;
shift[1] = -(double)cf_.getValue<int>( "setup", "shift_y" )*h; shift[1] = -(double)cf_.getValue<int>( "setup", "shift_y" )*h;
shift[2] = -(double)cf_.getValue<int>( "setup", "shift_z" )*h; shift[2] = -(double)cf_.getValue<int>( "setup", "shift_z" )*h;
} }*/
size_t npart = gh.count_leaf_cells(gh.levelmin(), gh.levelmax());; size_t npart = gh.count_leaf_cells(gh.levelmin(), gh.levelmax());;
size_t nwritten = 0; size_t nwritten = 0;