/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/source_MWtemp_WL.f90

Go to the documentation of this file.

subroutine source_MWtemp_WL(sou)
  use grid_parameter, only : nrg, ntg, npg
  use interface_sourceterm_MWtemp_CF
  use interface_sourceterm_MWtemp_current
  use interface_sourceterm_MWtemp_WL
  use make_array_3d
  implicit none
  real(8), pointer :: sou(:,:,:), sou1(:,:,:), sou2(:,:,:), sou3(:,:,:)
  call alloc_array3d(sou1,0,nrg,0,ntg,0,npg)
  call alloc_array3d(sou2,0,nrg,0,ntg,0,npg)
  call alloc_array3d(sou3,0,nrg,0,ntg,0,npg)
  call sourceterm_MWtemp_CF(sou1)
  call sourceterm_MWtemp_current(sou2)
  call sourceterm_MWtemp_WL(sou3)
  sou(0:nrg,0:ntg,0:npg) = sou1(0:nrg,0:ntg,0:npg) &
  &                      + sou2(0:nrg,0:ntg,0:npg) &
  &                      + sou3(0:nrg,0:ntg,0:npg)
  deallocate(sou1)
  deallocate(sou2)
  deallocate(sou3)
!
!  use def_velocity_potential, only : ramg
!  use def_metric, only : alph
!  use flusw, only : swflu
!        call interpo_linear_type0(ramgc,ramg,irg,itg,ipg)
!        call interpo_linear_type0(alpgc,alph,irg,itg,ipg)
!        utgc  = swflu *ramgc/(alpgc**2*hhgc) &
!           &      +(1.0d0-swflu)*hhgc/ber
!
end subroutine source_MWtemp_WL

---