aurora_raw.f90 Source File

submodule(io:io_aurora) io_aurora_raw

use timeutils, only : date_filename
implicit none

contains

module procedure output_aur_root
! subroutine output_aur_root(outdir,flagglow,ymd,UTsec,iver)
!! COLLECT COMPLETE DATA FROM WORKERS AND PROCESS FOR OUTPUT.
!! NO GHOST CELLS (I HOPE)

!real(wp), dimension(1:lx2,1:lwave,1:lx3) :: ivertmp
!real(wp), dimension(1:lx2all,1:lwave,1:lx3all) :: iverall
real(wp), dimension(1:lwave,1:lx2,1:lx3) :: ivertmp
real(wp), dimension(1:lwave,1:lx2all,1:lx3all) :: iverall

real(wp), dimension(1:lx2,1:lx3) :: emistmp                !single emission subgrid
real(wp), dimension(1:lx2all,1:lx3all) :: emisall          !single emission total grid
real(wp), dimension(1:lx2all,1:lx3all,1:lwave) :: iverout  !output array in the order scripts expect


character(:), allocatable :: outdir_composite, filenamefull
integer :: iwave

!!ivertmp=reshape(iver,[lx2,lwave,lx3],order=[1,3,2])
!ivertmp=reshape(iver,[lwave,lx2,lx3],order=[3,1,2])
!call gather_recv(ivertmp,tagAur,iverall)
do iwave=1,lwave
  emistmp=iver(:,:,iwave)
  call gather_recv(emistmp,tagAur,emisall)
  iverout(:,:,iwave)=emisall
end do


!FORM THE INPUT FILE NAME
outdir_composite=outdir//'/aurmaps/'

filenamefull = date_filename(outdir_composite,ymd,UTsec) // '.dat'

print *, '  Output file name (auroral maps):  ',filenamefull
block
  integer :: u
  open(newunit=u,file=filenamefull,status='replace',form='unformatted',access='stream',action='write')

  if(flagswap/=1) then
  !  write(u) reshape(iverall,[lx2all,lwave,lx3all],order=[2,1,3])
  !  write(u) reshape(iverall,[lx2all,lx3all,lwave],order=[2,3,1])
    write(u) iverout
  else
  !  write(u) reshape(iverall,[lx3all,lx2all,lwave],order=[3,2,1])
    write(u) reshape(iverout,[lx3all,lx2all,lwave],order=[2,1,3])
  end if

  close(u)
end block

if(.not. all(ieee_is_finite(iverout))) error stop 'iverout: non-finite value(s)'

end procedure output_aur_root

end submodule io_aurora_raw