00001 subroutine source_virial_CF(sou_Tkin,sou_Pint,sou_Wgra) 00002 use phys_constant, only : long 00003 use interface_source_virial_matter 00004 use interface_source_virial_gravity_CF 00005 implicit none 00006 real(long), pointer :: sou_Tkin(:,:,:), sou_Pint(:,:,:), sou_Wgra(:,:,:) 00007 ! 00008 ! --- Compute source terms for virial relations. 00009 ! --- sources for kinetic and internal energies are defined on SCF grid points. 00010 ! --- source for gravitational energies is defined on mid points. 00011 ! 00012 call source_virial_matter(sou_Tkin,sou_Pint) 00013 call source_virial_gravity_CF(sou_Wgra) 00014 ! 00015 end subroutine source_virial_CF