16 real(kind=8):: d, e, ee, p, u, temp, grav, mu, kr, b
25 SUBROUTINE init (self)
27 logical,
save:: first_time=.true.
29 if (.not.first_time)
return 30 call trace%begin (
'scaling_t%init')
31 if (io%master) print *, &
32 '-------------------------------- scaling_mod -------------------------' 41 self%u = self%l/self%t
42 self%m = self%d*self%l**3
43 self%p = self%d*self%u**2
44 self%kr = 1.0/(self%d*self%l)
46 self%e = self%m*self%u**2
47 self%b = self%l**(-0.5)*self%m**0.5*self%t**(-1)
49 self%temp = self%mu*(cgs%m_u/self%m)/(cgs%k_b/self%e)
50 self%grav = cgs%grav*(self%m/self%l**3)*self%t**2
52 print 1,
' CODE UNITS: (cgs):' 53 print 1,
' length: ',self%l
54 print 1,
' mass: ',self%m
55 print 1,
' time: ',self%t
56 print 1,
' velocity: ',self%u
57 print 1,
' density: ',self%d
58 print 1,
' pressure: ',self%p
59 print 1,
'magnetic flux density: ',self%b
60 print 1,
' energy: ',self%e
61 print 1,
' gravity: ',self%grav
62 print 1,
' temperature: ',self%temp
63 1
format(1x,a,1p,e11.3)
Define code units, in terms of (here) CGS units.
Fundamental constants in CGS and SI units.
Fundamental constants in CGS and SI units.