pro norh_hanaoka,rawfile,szfile,frame $ ,imgdir=imgdir $ ,outfile=outfile,rlphdiffile=rlphdiffile $ ,version=version,freq=freq $ ,cenfnl=cenfnl,cenunit=cenunit,size=size $ ,nfrcal=nfrcal,nfrint=nfrint $ ,clcra=clcra,clcrs=clcrs,clcrfac=clcrfac $ ,cendty=cendty,beamcor=beamcor,restoration=restoration if (n_elements(outfile) eq 0) then outfile=curdir()+'/'+'norh_synth_input' if (n_elements(imgdir) eq 0) then imgdir=curdir()+'/' if (n_elements(version) eq 0) then version='-1' if (version eq '-1') then version='6.1' if (n_elements(freq) eq 0) then freq=17 if (n_elements(cenunit) eq 0) then cenunit=1 if (n_elements(cenfnl) eq 0) then cenfnl=[0,0] if (n_elements(nfrcal) eq 0) then nfrcal=1 if (n_elements(nfrint) eq 0) then nfrint=1 if (n_elements(clcrfac) eq 0) then clcrfac=0.01 if (n_elements(cendty) eq 0) then cendty=[999,999] if (n_elements(beamcor) eq 0) then beamcor=1 if (n_elements(restoration) eq 0) then restoration=1 mfrx=n_elements(frame) case 1 of (version eq '5.0') : begin if (freq eq 17) then begin if (n_elements(clcra) eq 0) then clcra=-0.5 if (n_elements(clcrs) eq 0) then clcrs=-0.15 openw,unit,outfile,/get_lun printf,unit,rawfile printf,unit,szfile printf,unit,rlphdiffile printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,clcra printf,unit,clcrs printf,unit,cendty printf,unit,cenfnl printf,unit,beamcor printf,unit,restoration printf,unit,imgdir close,unit free_lun,unit endif else begin message,'Sorry, this version only supports 17GHz data' endelse end (version eq '5.1') or (version eq '6.0') : begin if (freq eq 17) then begin if (n_elements(clcra) eq 0) then clcra=-0.5 if (n_elements(clcrs) eq 0) then clcrs=-0.5 if (n_elements(size) eq 0) then size=512 openw,unit,outfile,/get_lun printf,unit,rawfile printf,unit,szfile printf,unit,rlphdiffile printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,clcra,clcrs printf,unit,clcrfac printf,unit,cendty printf,unit,cenunit,cenfnl printf,unit,beamcor printf,unit,restoration printf,unit,size printf,unit,imgdir close,unit free_lun,unit endif else begin if (n_elements(clcra) eq 0) then clcra=-1. if (n_elements(size) eq 0) then size=1024 openw,unit,outfile,/get_lun printf,unit,rawfile printf,unit,szfile printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,clcra printf,unit,clcrfac printf,unit,cendty printf,unit,cenunit,cenfnl printf,unit,beamcor printf,unit,restoration printf,unit,size printf,unit,imgdir close,unit free_lun,unit endelse end (version eq '6.1') : begin if (freq eq 17) then begin if (n_elements(clcra) eq 0) then clcra=-0.5 if (n_elements(clcrs) eq 0) then clcrs=-0.5 if (n_elements(size) eq 0) then size=512 openw,unit,outfile,/get_lun printf,unit,imgdir printf,unit,rawfile printf,unit,szfile printf,unit,size printf,unit,cenunit,cenfnl printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,clcra,clcrs printf,unit,clcrfac close,unit free_lun,unit endif else begin if (n_elements(clcra) eq 0) then clcra=-1. if (n_elements(size) eq 0) then size=1024 openw,unit,outfile,/get_lun printf,unit,imgdir printf,unit,rawfile printf,unit,szfile printf,unit,size printf,unit,cenunit,cenfnl printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,clcra printf,unit,clcrfac close,unit free_lun,unit endelse end endcase return end