Changes between Version 54 and Version 55 of u/bliu/MHDJets


Ignore:
Timestamp:
03/25/20 11:55:34 (5 years ago)
Author:
Baowei Liu
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/bliu/MHDJets

    v54 v55  
     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
    138= Runs following cooling curve 2.24.2020 =
    239