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

tet mesh: writing out base level information to header, level written out is no w.r.t. this base level

Browse source
This commit is contained in:
Oliver Hahn 2013-10-24 18:48:35 +02:00
parent 91133a23a4
commit 2ef810611c
1 changed files with 12 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

20
plugins/output_gadget_tetmesh.cc
View file

@ -405,7 +405,7 @@ void delete_duplicates( void )
//////////////////////////////////////////////////////////////////////////////////
//////////////////////////////////////////////////////////////////////////////////
const int empty_fill_bytes = 60;
const int empty_fill_bytes = 54;
template< typename T_store=float >
class gadget_tetmesh_output_plugin : public output_plugin
@ -442,7 +442,10 @@ protected:
int flag_metals;
unsigned int npartTotalHighWord[6];
int flag_entropy_instead_u;
char fill[empty_fill_bytes];
//////////////////////////////////////
char fill[empty_fill_bytes];
//////////////////////////////////////
int tetgrid_baselevel;
}header;
@ -1179,8 +1182,6 @@ public:
break;
}
lidsum += lid;
}
@ -1227,6 +1228,11 @@ public:
rhoc *= 1e10;
header_.mass[1] = 0.0;
header_.mass[2] = 0.0;
header_.mass[5] = 0.0;
header_.tetgrid_baselevel = tetgrid_baselevel;
idmap.clear();
@ -1245,8 +1251,6 @@ public:
np_type2_ = num_p_t2;
np_type5_ = num_p_t5;
header_.npart[1] = num_p_t1;
header_.npartTotal[1] = num_p_t1;
@ -1390,13 +1394,13 @@ public:
for( size_t ip=0; ip<num_p; ++ip )
{
if( temp_dat.size() < block_buf_size_ )
temp_dat.push_back( P[ip].Level );
temp_dat.push_back( P[ip].Level - header_.tetgrid_baselevel );
else
{
ofs_temp.write( (char*)&temp_dat[0], sizeof(T_store)*block_buf_size_ );
nwritten += block_buf_size_;
temp_dat.clear();
temp_dat.push_back( P[ip].Level );
temp_dat.push_back( P[ip].Level - header_.tetgrid_baselevel );
}
}