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

Go to the documentation of this file.

subroutine source_trG_WL(sou)
  use grid_parameter, only : nrg, ntg, npg
  use interface_sourceterm_trG_CF
  use interface_sourceterm_trG_CF_corot
  use interface_sourceterm_trG_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_trG_CF(sou1)
  call sourceterm_trG_CF_corot(sou2)
  call sourceterm_trG_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)
!        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_trG_WL

---