00001 subroutine source_virial_WL(sou_Tkin,sou_Pint,sou_Wgra) 00002 use phys_constant, only : long 00003 use interface_source_virial_matter 00004 use interface_source_virial_gravity_WL 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 ! --- sources for EM, and gravitational energies are defined on mid points. 00011 ! 00012 call source_virial_matter(sou_Tkin,sou_Pint) 00013 call source_virial_gravity_WL(sou_Wgra) 00014 ! 00015 end subroutine source_virial_WL