00001 subroutine calc_ToverW 00002 use phys_constant, only : long 00003 use def_matter_parameter, only : ome, radi 00004 use def_quantities 00005 implicit none 00006 ! 00007 omega = ome/radi 00008 T_kinene = 0.5d0*omega*angmom 00009 W_gravene = propermass + T_kinene - admmass 00010 ToverW = T_kinene/dabs(W_gravene) 00011 I_inertia = angmom/omega 00012 ! 00013 end subroutine calc_ToverW