| 221 | | Until now we silently omitted one important aspect of the PIERNIK simulaitons: the delta timestep is not constant and depends on the current state of both modules: MHD and MC. In the original code the minimum of the dt,,mc,, and dt,,mhd,, values were always chosen. |
| 222 | | {{{ |
| 223 | | |
| | 221 | Until now we silently omitted one important aspect of the PIERNIK simulaitons: the delta timestep is not constant and depends on the current state of both modules: MHD and MC. In the original code the minimum of the dt,,mc,, and dt,,mhd,, values were always chosen. In order to do the same in MUSCLE variant of FORTRAN code, we must define two new conduits: one for sending dt,,mc,, from MC to MHD and second one for sending back the final dt. The conduits definition looks now as follows: |
| | 222 | {{{ |
| | 223 | cs.attach('mhd' => 'mc') { |
| | 224 | tie('rho_gas', 'rho_gas') |
| | 225 | tie('m_gas', 'm_gas') |
| | 226 | tie('dt_final', 'dt_final') |
| | 227 | } |
| | 228 | |
| | 229 | cs.attach('mc' => 'mhd') { |
| | 230 | tie('dt_mc', 'dt_mc') |
| | 231 | } |
