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