cwpearson/astaroth
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Apparently my edit this time was unsuccessful, tons of error messeages showed up when I ran "make -j" command in working directory. However this commit is mainly for education purposes, it's so that Miikka (main) and others can see what I changed easily and can use as teaching reference.

Browse Source
This commit is contained in:
JackHsu
2019-07-26 16:22:04 +08:00
parent 67d9f19006
commit 9d625688ac
2 changed files with 14 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
acc/mhd_solver/stencil_defines.h
View File

@@ -72,6 +72,12 @@
FUNC(AC_star_pos_y),\
FUNC(AC_star_pos_z),\
FUNC(AC_M_star),\
/* properties of sink particle*/\
FUNC(AC_sink_pos_x),\
FUNC(AC_sink_pos_y),\
FUNC(AC_sink_pos_z),\
FUNC(AC_M_sink),\
//The above is function for sink particle, not sure if it's supposed to be here but just a try
/* Run params */\
FUNC(AC_cdt), \
FUNC(AC_cdtv), \

8
src/standalone/simulation.cc
View File

@@ -88,6 +88,14 @@ write_mesh_info(const AcMeshInfo* config)
fprintf(infotxt, "real AC_cs2_sound %e \n", (double)config->real_params[AC_cs2_sound]);
fprintf(infotxt, "real AC_cv_sound %e \n", (double)config->real_params[AC_cv_sound]);
//Here I'm still trying to copy the structure of the code above, and see if this will work for sink particle.
//I haven't fully undertand what these lines do but I'll read up on them soon. This is still yet experimental.
// Sink particle
fprintf(infotext, "real AC_sink_pos_x %i \n", config->real_params[AC_sink_pos_x]);
fprintf(infotext, "real AC_sink_pos_y %i \n", config->real_params[AC_sink_pos_y]);
fprintf(infotext, "real AC_sink_pos_z %i \n", config->real_params[AC_sink_pos_z]);
fprintf(infotext, "real AC_M_sink %e \n", (double)config->real_params[AC_M_sink]);
fclose(infotxt);
}