;pro azam_prep_hinode,merlin_fitsfile,byte_ambig,afiles_dir=afiles_dir ; for now do not output byte_ambig pro azam_prep_hinode,merlin_fitsfile,afiles_dir=afiles_dir ; routine to generate the ASP-format a_* files needed for ; input to AZAM. This routine reads the FITS output from MERLIN, ; then prepares the needed a_files and writes them out ;INPUTS: merlin_fitsfile = path to MERLIN FITS inversion output ; ;OUTPUTS: byte_ambig = byte array with second bit altered ; for AZAM ambiguity resolution. First ; bit of this array is reserved for ; Metcalf auto ambig resolution. Set to ; zero means unchanged from original MERLIN ; observer-frame azimuth (a_azm), set to 1 means ; opposite observer-frame azimuth ;KEYWORDS: afiles_dir = output directory for a_* files. If not ; set, files output to the present working ; directory ; common area for inversion results read in rdmerlinftis.pro common com_merlinresults, $ fld, psi, azm, eta0, dmp, dop,$ bzero, b1mu, vmac, cen1, cen2, alpha, $ delta, cct, chisq, chisq_i, chisq_q, $ chisq_u, chisq_v, signed_v, niter, occt, $ rconv, netpol, netpol_full, $ fld_er, psi_er, azm_er, eta0_er, dmp_er, dop_er, $ bzero_er, b1mu_er, vmac_er, cen1_er, cen2_er, $ alpha_er, delta_er,xx,yy,slat,slong, $ yr,mnth,dy,hr,minit,sec outdir = './' if keyword_set(afiles_dir) then outdir = afiles_dir ; !!!check to see if directory exists, if not, create it!!! ; first read in the MERLIN inversion results rdmerlinfits,afiles_dir,merlin_fitsfile ; generate the a_files azam_hinode2afiles,outdir ;; prepare output files ;sz = size(cct) ;nxx = sz(1) & nyy = sz(2) ;byte_ambig = bytarr(nxx,nyy) ;; for present, set the byte to 0, 1. Later must change to ;; setting second bit ; ; ; ;stop end