import numpy as np
import sys
import os.path
from astropy.io import fits
from scipy.io import FortranFile

#Read fits file name, output file name, map dimensions in npix*npix
args = sys.argv
fname_in = args[1]
fname_out = args[2]
#npix = int(args[3])

print('Reading', fname_in)

#Read in fits file:
hdul = fits.open(fname_in)
hdul[0].header
map = hdul[0].data
npix=np.shape(map)[0]
print('Got map, with size', npix)
print('Writing', fname_out)

#unformatted
f = FortranFile(fname_out, 'w')
f.write_record(np.array(map, dtype=np.float))
f.close()

#read these files in f90 as:
#real(8) :: input_map_d(npix,npix)
#character(*), parameter :: input_form=    'unformatted'
#character(*), parameter :: input_access= 'stream'
#open(unit=fu,file=fp,form=input_form,status='old',iostat=file_status, access=input_access)
#read(fu)  input_map_d