00001 subroutine source_AHarea_CF(sousf)
00002 use phys_constant, only : long, pi
00003 use grid_parameter, only : ntg, npg
00004 use def_metric, only : psi
00005 use interface_interpo_linear_surface_type0
00006 implicit none
00007 real(long), pointer :: sousf(:,:)
00008 real(long) :: val
00009 integer :: irg, itg, ipg
00010
00011 irg = 0
00012 do ipg = 1, npg
00013 do itg = 1, ntg
00014 call interpo_linear_surface_type0(val,psi,irg,itg,ipg)
00015 sousf(itg,ipg) = val**4
00016 end do
00017 end do
00018
00019 end subroutine source_AHarea_CF