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

Go to the documentation of this file.

subroutine source_MWspatial_WL(souvec)
  use grid_parameter, only : nrg, ntg, npg
  use interface_sourceterm_MWspatial_current_CF
  use interface_sourceterm_MWspatial_current_WL
  use interface_sourceterm_MWspatial_CF
  use interface_sourceterm_MWspatial_WL
  use make_array_4d
  implicit none
  real(8), pointer :: souvec(:,:,:,:)
  real(8), pointer :: souvec1(:,:,:,:), souvec2(:,:,:,:), 
                     souvec3(:,:,:,:), souvec4(:,:,:,:)
!
  call alloc_array4d(souvec1,0,nrg,0,ntg,0,npg,1,3)
  call alloc_array4d(souvec2,0,nrg,0,ntg,0,npg,1,3)
  call alloc_array4d(souvec3,0,nrg,0,ntg,0,npg,1,3)
  call alloc_array4d(souvec4,0,nrg,0,ntg,0,npg,1,3)
!
  call sourceterm_MWspatial_current_CF(souvec1)
  call sourceterm_MWspatial_current_WL(souvec2)
  call sourceterm_MWspatial_CF(souvec3)
  call sourceterm_MWspatial_WL(souvec4)
  souvec(0:nrg,0:ntg,0:npg,1:3) = souvec1(0:nrg,0:ntg,0:npg,1:3) &
  &                             + souvec2(0:nrg,0:ntg,0:npg,1:3) &
  &                             + souvec3(0:nrg,0:ntg,0:npg,1:3) &
  &                             + souvec4(0:nrg,0:ntg,0:npg,1:3)
!
  deallocate(souvec1)
  deallocate(souvec2)
  deallocate(souvec3)
  deallocate(souvec4)
!
!          call interpo_linear_type0(ramgc,ramg,irg,itg,ipg)
!          call interpo_linear_type0(psigc, psi,irg,itg,ipg)
!          call interpo_linear_type0(alpgc,alph,irg,itg,ipg)
!          call interpo_linear_type0(fnc2gc,fnc2,irg,itg,ipg)
!          afnc2inv = alpgc/fnc2gc
!          call interpo_linear_type0(emdgc,emdg,irg,itg,ipg)
!          call peos_q2hprho(emdgc,hhgc,pregc,rhogc,ene)
!
!          utgc  =       swflu *ramgc/(alpgc**2*hhgc) &
!             &      +(1.0d0-swflu)*hhgc/ber
end subroutine source_MWspatial_WL

---