submodule(io:io_aurora) io_aurora_raw

use timeutils, only : date_filename

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