| 1 | = Debugging Cooling Length Field 3.23.2020 = |
| 2 | |
| 3 | {{{ |
| 4 | ... |
| 5 | ne = (6d0/7d0)*q(1)*nScale ! number densities based on Z = 6 |
| 6 | nAl = (1d0/7d0)*q(1)*nScale |
| 7 | logne = LOG10(ne) |
| 8 | Temp = Temp/1.16045d4 ! K to eV |
| 9 | logtemp = LOG10(Temp) |
| 10 | AlCoolingStrength=ne*nAl*GetTableValue(AlTab, (/logtemp, logne/))/(nScale**2d0) |
| 11 | AlCoolingStrength=AlCoolingStrength*ScaleCool |
| 12 | |
| 13 | ... |
| 14 | CASE(CoolingLength_Field) |
| 15 | CALL cons_to_source(q) |
| 16 | vg = SQRT(SUM(q(2:m_high)**2d0)) |
| 17 | !cs = SQRT(gamma*(gamma-1d0)*q(iE)/q(1)) |
| 18 | !v_rel = ABS(-vg/3d0 + SQRT(4d0/9d0*vg**2d0-cs**2d0)) |
| 19 | !cs_floor = SQRT(gamma*FloorTemp/TempScale) |
| 20 | !v = MAX(v_rel,cs_floor) |
| 21 | GetField = vg*q(iE)/GetCoolingStrength(q) |
| 22 | CALL source_to_cons(q) |
| 23 | |
| 24 | ... |
| 25 | |
| 26 | SUBROUTINE Cons_To_Source(q) |
| 27 | REAL(KIND=qpREC), DIMENSION(:) :: q |
| 28 | CALL ConvertTotalToInternalEnergy(q) |
| 29 | !DIR$ NOUNROLL |
| 30 | q(m_low:m_high)=q(m_low:m_high)/q(1) |
| 31 | |
| 32 | |
| 33 | }}} |
| 34 | |
| 35 | |
| 36 | ---- |
| 37 | |