Conversion of standardized ReadMe file for
file /./ftp/cats/J/PASJ/62/273 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe C============================================================================= C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-20 C============================================================================= * This code was generated from the ReadMe file documenting a catalogue * according to the "Standard for Documentation of Astronomical Catalogues" * currently in use by the Astronomical Data Centers (CDS, ADC, A&A) * (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx) * Please report problems or questions to C============================================================================= implicit none * Unspecified or NULL values, generally corresponding to blank columns, * are assigned one of the following special values: * rNULL__ for unknown or NULL floating-point values * iNULL__ for unknown or NULL integer values real*4 rNULL__ integer*4 iNULL__ parameter (rNULL__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. J/PASJ/62/273 AKARI SMC photometry (Ita+, 2010) *================================================================================ *AKARI near- to mid-infrared imaging and spectroscopic observations of the *Small Magellanic Cloud. *I. Bright point source list. * Ita Y., Onaka T., Tanabe T., Matsunaga N., Matsuura M., Yamamura I., * Nakada Y., Izumiura H., Ueta T., Mito H., Fukushi H., Kato D. * <Publ. Astron. Soc. Jap., 62, 273-286 (2010)> * =2010PASJ...62..273I C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'n3.dat' ! Photometry in N3 (3.2um) band integer*4 nr__ parameter (nr__=12899) ! Number of records character*110 ar__ ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000) character*3 Band (nr__) ! Observed band (N3, N4, S7, S11, L15 or L24) real*4 mag (nr__) ! (mag) Magnitude in Band in Vega system real*4 e_mag (nr__) ! (mag) rms uncertainty on mag character*10 Obs_date (nr__) ! ("YYYY-MM-DD") Observation date character*8 Obs_time (nr__) ! ("h:m:s") Observation time (hh:mm:ss) integer*4 S0L1 (nr__) ! [0/1] S0L1 flag (photometry coming either * from short (0) or long (1) data) integer*4 A0F1 (nr__) ! [0/1] A0F1 flag (photometry method, * 0 for aperture photometry and * 1 for fit photometry,) integer*4 SAT (nr__) ! [0/1] SAT flag (set to 1 if a bright source * is in proximity) integer*4 COL (nr__) ! [0/1] COL flag (set to 1 if the star is * located in the column pulldown belt) integer*4 MUX (nr__) ! [0/1] MUX flag (set to 1 if the star is * located in the mux-bleed pulldown belt) integer*4 GHO (nr__) ! [0/1] GHO flag (set to 1 if the star is * suspected to be a ghost) real*8 RASdeg (nr__) ! (deg) Closest SAGE-SMC Archive/Full catalog source * right ascension in decimal degrees (J2000) real*8 DESdeg (nr__) ! (deg) Closest SAGE-SMC Archive/Full catalog source * declination in decimal degrees (J2000) real*4 Iband (nr__) ! (um) IRAC band (3.6, 4.5, 8.0 or 24) real*4 IRAC (nr__) ! (mag) IRAC magnitude in Iband in Vega system real*4 e_IRAC (nr__) ! (mag) rms uncertainty on IRAC c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'n4.dat' ! Photometry in N4 (4.1um) band integer*4 nr__1 parameter (nr__1=9748) ! Number of records character*110 ar__1 ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 RAdeg_1 (nr__1) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg_1 (nr__1) ! (deg) Declination in decimal degrees (J2000) character*3 Band_1 (nr__1) ! Observed band (N3, N4, S7, S11, L15 or L24) real*4 mag_1 (nr__1) ! (mag) Magnitude in Band in Vega system real*4 e_mag_1 (nr__1) ! (mag) rms uncertainty on mag character*10 Obs_date_1 (nr__1) ! ("YYYY-MM-DD") Observation date character*8 Obs_time_1 (nr__1) ! ("h:m:s") Observation time (hh:mm:ss) integer*4 S0L1_1 (nr__1) ! [0/1] S0L1 flag (photometry coming either * from short (0) or long (1) data) integer*4 A0F1_1 (nr__1) ! [0/1] A0F1 flag (photometry method, * 0 for aperture photometry and * 1 for fit photometry,) integer*4 SAT_1 (nr__1) ! [0/1] SAT flag (set to 1 if a bright source * is in proximity) integer*4 COL_1 (nr__1) ! [0/1] COL flag (set to 1 if the star is * located in the column pulldown belt) integer*4 MUX_1 (nr__1) ! [0/1] MUX flag (set to 1 if the star is * located in the mux-bleed pulldown belt) integer*4 GHO_1 (nr__1) ! [0/1] GHO flag (set to 1 if the star is * suspected to be a ghost) real*8 RASdeg_1 (nr__1) ! (deg) Closest SAGE-SMC Archive/Full catalog source * right ascension in decimal degrees (J2000) real*8 DESdeg_1 (nr__1) ! (deg) Closest SAGE-SMC Archive/Full catalog source * declination in decimal degrees (J2000) real*4 Iband_1 (nr__1) ! (um) IRAC band (3.6, 4.5, 8.0 or 24) real*4 IRAC_1 (nr__1) ! (mag) IRAC magnitude in Iband in Vega system real*4 e_IRAC_1 (nr__1) ! (mag) rms uncertainty on IRAC c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 's7.dat' ! Photometry in S7 (7.0um) band integer*4 nr__2 parameter (nr__2=1838) ! Number of records character*110 ar__2 ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 RAdeg_2 (nr__2) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg_2 (nr__2) ! (deg) Declination in decimal degrees (J2000) character*3 Band_2 (nr__2) ! Observed band (N3, N4, S7, S11, L15 or L24) real*4 mag_2 (nr__2) ! (mag) Magnitude in Band in Vega system real*4 e_mag_2 (nr__2) ! (mag) rms uncertainty on mag character*10 Obs_date_2 (nr__2) ! ("YYYY-MM-DD") Observation date character*8 Obs_time_2 (nr__2) ! ("h:m:s") Observation time (hh:mm:ss) integer*4 S0L1_2 (nr__2) ! [0/1] S0L1 flag (photometry coming either * from short (0) or long (1) data) integer*4 A0F1_2 (nr__2) ! [0/1] A0F1 flag (photometry method, * 0 for aperture photometry and * 1 for fit photometry,) integer*4 SAT_2 (nr__2) ! [0/1] SAT flag (set to 1 if a bright source * is in proximity) integer*4 COL_2 (nr__2) ! [0/1] COL flag (set to 1 if the star is * located in the column pulldown belt) integer*4 MUX_2 (nr__2) ! [0/1] MUX flag (set to 1 if the star is * located in the mux-bleed pulldown belt) integer*4 GHO_2 (nr__2) ! [0/1] GHO flag (set to 1 if the star is * suspected to be a ghost) real*8 RASdeg_2 (nr__2) ! (deg) Closest SAGE-SMC Archive/Full catalog source * right ascension in decimal degrees (J2000) real*8 DESdeg_2 (nr__2) ! (deg) Closest SAGE-SMC Archive/Full catalog source * declination in decimal degrees (J2000) real*4 Iband_2 (nr__2) ! (um) IRAC band (3.6, 4.5, 8.0 or 24) real*4 IRAC_2 (nr__2) ! (mag) IRAC magnitude in Iband in Vega system real*4 e_IRAC_2 (nr__2) ! (mag) rms uncertainty on IRAC c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 's11.dat' ! Photometry in S11 (11.0um) band integer*4 nr__3 parameter (nr__3=1045) ! Number of records character*110 ar__3 ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 RAdeg_3 (nr__3) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg_3 (nr__3) ! (deg) Declination in decimal degrees (J2000) character*3 Band_3 (nr__3) ! Observed band (N3, N4, S7, S11, L15 or L24) real*4 mag_3 (nr__3) ! (mag) Magnitude in Band in Vega system real*4 e_mag_3 (nr__3) ! (mag) rms uncertainty on mag character*10 Obs_date_3 (nr__3) ! ("YYYY-MM-DD") Observation date character*8 Obs_time_3 (nr__3) ! ("h:m:s") Observation time (hh:mm:ss) integer*4 S0L1_3 (nr__3) ! [0/1] S0L1 flag (photometry coming either * from short (0) or long (1) data) integer*4 A0F1_3 (nr__3) ! [0/1] A0F1 flag (photometry method, * 0 for aperture photometry and * 1 for fit photometry,) integer*4 SAT_3 (nr__3) ! [0/1] SAT flag (set to 1 if a bright source * is in proximity) integer*4 COL_3 (nr__3) ! [0/1] COL flag (set to 1 if the star is * located in the column pulldown belt) integer*4 MUX_3 (nr__3) ! [0/1] MUX flag (set to 1 if the star is * located in the mux-bleed pulldown belt) integer*4 GHO_3 (nr__3) ! [0/1] GHO flag (set to 1 if the star is * suspected to be a ghost) real*8 RASdeg_3 (nr__3) ! (deg) Closest SAGE-SMC Archive/Full catalog source * right ascension in decimal degrees (J2000) real*8 DESdeg_3 (nr__3) ! (deg) Closest SAGE-SMC Archive/Full catalog source * declination in decimal degrees (J2000) real*4 Iband_3 (nr__3) ! (um) IRAC band (3.6, 4.5, 8.0 or 24) real*4 IRAC_3 (nr__3) ! (mag) IRAC magnitude in Iband in Vega system real*4 e_IRAC_3 (nr__3) ! (mag) rms uncertainty on IRAC c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'l15.dat' ! Photometry in L15 (15.0um) band integer*4 nr__4 parameter (nr__4=479) ! Number of records character*110 ar__4 ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 RAdeg_4 (nr__4) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg_4 (nr__4) ! (deg) Declination in decimal degrees (J2000) character*3 Band_4 (nr__4) ! Observed band (N3, N4, S7, S11, L15 or L24) real*4 mag_4 (nr__4) ! (mag) Magnitude in Band in Vega system real*4 e_mag_4 (nr__4) ! (mag) rms uncertainty on mag character*10 Obs_date_4 (nr__4) ! ("YYYY-MM-DD") Observation date character*8 Obs_time_4 (nr__4) ! ("h:m:s") Observation time (hh:mm:ss) integer*4 S0L1_4 (nr__4) ! [0/1] S0L1 flag (photometry coming either * from short (0) or long (1) data) integer*4 A0F1_4 (nr__4) ! [0/1] A0F1 flag (photometry method, * 0 for aperture photometry and * 1 for fit photometry,) integer*4 SAT_4 (nr__4) ! [0/1] SAT flag (set to 1 if a bright source * is in proximity) integer*4 COL_4 (nr__4) ! [0/1] COL flag (set to 1 if the star is * located in the column pulldown belt) integer*4 MUX_4 (nr__4) ! [0/1] MUX flag (set to 1 if the star is * located in the mux-bleed pulldown belt) integer*4 GHO_4 (nr__4) ! [0/1] GHO flag (set to 1 if the star is * suspected to be a ghost) real*8 RASdeg_4 (nr__4) ! (deg) Closest SAGE-SMC Archive/Full catalog source * right ascension in decimal degrees (J2000) real*8 DESdeg_4 (nr__4) ! (deg) Closest SAGE-SMC Archive/Full catalog source * declination in decimal degrees (J2000) real*4 Iband_4 (nr__4) ! (um) IRAC band (3.6, 4.5, 8.0 or 24) real*4 IRAC_4 (nr__4) ! (mag) IRAC magnitude in Iband in Vega system real*4 e_IRAC_4 (nr__4) ! (mag) rms uncertainty on IRAC c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'l24.dat' ! Photometry in L24 (24.0um) band integer*4 nr__5 parameter (nr__5=356) ! Number of records character*110 ar__5 ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 RAdeg_5 (nr__5) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg_5 (nr__5) ! (deg) Declination in decimal degrees (J2000) character*3 Band_5 (nr__5) ! Observed band (N3, N4, S7, S11, L15 or L24) real*4 mag_5 (nr__5) ! (mag) Magnitude in Band in Vega system real*4 e_mag_5 (nr__5) ! (mag) rms uncertainty on mag character*10 Obs_date_5 (nr__5) ! ("YYYY-MM-DD") Observation date character*8 Obs_time_5 (nr__5) ! ("h:m:s") Observation time (hh:mm:ss) integer*4 S0L1_5 (nr__5) ! [0/1] S0L1 flag (photometry coming either * from short (0) or long (1) data) integer*4 A0F1_5 (nr__5) ! [0/1] A0F1 flag (photometry method, * 0 for aperture photometry and * 1 for fit photometry,) integer*4 SAT_5 (nr__5) ! [0/1] SAT flag (set to 1 if a bright source * is in proximity) integer*4 COL_5 (nr__5) ! [0/1] COL flag (set to 1 if the star is * located in the column pulldown belt) integer*4 MUX_5 (nr__5) ! [0/1] MUX flag (set to 1 if the star is * located in the mux-bleed pulldown belt) integer*4 GHO_5 (nr__5) ! [0/1] GHO flag (set to 1 if the star is * suspected to be a ghost) real*8 RASdeg_5 (nr__5) ! (deg) Closest SAGE-SMC Archive/Full catalog source * right ascension in decimal degrees (J2000) real*8 DESdeg_5 (nr__5) ! (deg) Closest SAGE-SMC Archive/Full catalog source * declination in decimal degrees (J2000) real*4 Iband_5 (nr__5) ! (um) IRAC band (3.6, 4.5, 8.0 or 24) real*4 IRAC_5 (nr__5) ! (mag) IRAC magnitude in Iband in Vega system real*4 e_IRAC_5 (nr__5) ! (mag) rms uncertainty on IRAC C============================================================================= C Loading file 'n3.dat' ! Photometry in N3 (3.2um) band C Format for file interpretation 1 format( + F9.6,1X,F10.6,1X,A3,1X,F6.3,3X,F5.3,1X,A10,1X,A8,1X,I1,1X,I1, + 1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.6,1X,F10.6,1X,F4.1,1X,F6.3,1X, + F5.3) C Effective file loading open(unit=1,status='old',file= +'n3.dat') write(6,*) '....Loading file: n3.dat' do i__=1,12899 read(1,'(A110)')ar__ read(ar__,1) + RAdeg(i__),DEdeg(i__),Band(i__),mag(i__),e_mag(i__), + Obs_date(i__),Obs_time(i__),S0L1(i__),A0F1(i__),SAT(i__), + COL(i__),MUX(i__),GHO(i__),RASdeg(i__),DESdeg(i__),Iband(i__), + IRAC(i__),e_IRAC(i__) c ..............Just test output........... write(6,1) + RAdeg(i__),DEdeg(i__),Band(i__),mag(i__),e_mag(i__), + Obs_date(i__),Obs_time(i__),S0L1(i__),A0F1(i__),SAT(i__), + COL(i__),MUX(i__),GHO(i__),RASdeg(i__),DESdeg(i__),Iband(i__), + IRAC(i__),e_IRAC(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'n4.dat' ! Photometry in N4 (4.1um) band C Format for file interpretation 2 format( + F9.6,1X,F10.6,1X,A3,1X,F6.3,3X,F5.3,1X,A10,1X,A8,1X,I1,1X,I1, + 1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.6,1X,F10.6,1X,F4.1,1X,F6.3,1X, + F5.3) C Effective file loading open(unit=1,status='old',file= +'n4.dat') write(6,*) '....Loading file: n4.dat' do i__=1,9748 read(1,'(A110)')ar__1 read(ar__1,2) + RAdeg_1(i__),DEdeg_1(i__),Band_1(i__),mag_1(i__),e_mag_1(i__), + Obs_date_1(i__),Obs_time_1(i__),S0L1_1(i__),A0F1_1(i__), + SAT_1(i__),COL_1(i__),MUX_1(i__),GHO_1(i__),RASdeg_1(i__), + DESdeg_1(i__),Iband_1(i__),IRAC_1(i__),e_IRAC_1(i__) c ..............Just test output........... write(6,2) + RAdeg_1(i__),DEdeg_1(i__),Band_1(i__),mag_1(i__),e_mag_1(i__), + Obs_date_1(i__),Obs_time_1(i__),S0L1_1(i__),A0F1_1(i__), + SAT_1(i__),COL_1(i__),MUX_1(i__),GHO_1(i__),RASdeg_1(i__), + DESdeg_1(i__),Iband_1(i__),IRAC_1(i__),e_IRAC_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 's7.dat' ! Photometry in S7 (7.0um) band C Format for file interpretation 3 format( + F9.6,1X,F10.6,1X,A3,1X,F6.3,3X,F5.3,1X,A10,1X,A8,1X,I1,1X,I1, + 1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.6,1X,F10.6,1X,F4.1,1X,F6.3,1X, + F5.3) C Effective file loading open(unit=1,status='old',file= +'s7.dat') write(6,*) '....Loading file: s7.dat' do i__=1,1838 read(1,'(A110)')ar__2 read(ar__2,3) + RAdeg_2(i__),DEdeg_2(i__),Band_2(i__),mag_2(i__),e_mag_2(i__), + Obs_date_2(i__),Obs_time_2(i__),S0L1_2(i__),A0F1_2(i__), + SAT_2(i__),COL_2(i__),MUX_2(i__),GHO_2(i__),RASdeg_2(i__), + DESdeg_2(i__),Iband_2(i__),IRAC_2(i__),e_IRAC_2(i__) c ..............Just test output........... write(6,3) + RAdeg_2(i__),DEdeg_2(i__),Band_2(i__),mag_2(i__),e_mag_2(i__), + Obs_date_2(i__),Obs_time_2(i__),S0L1_2(i__),A0F1_2(i__), + SAT_2(i__),COL_2(i__),MUX_2(i__),GHO_2(i__),RASdeg_2(i__), + DESdeg_2(i__),Iband_2(i__),IRAC_2(i__),e_IRAC_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 's11.dat' ! Photometry in S11 (11.0um) band C Format for file interpretation 4 format( + F9.6,1X,F10.6,1X,A3,1X,F6.3,3X,F5.3,1X,A10,1X,A8,1X,I1,1X,I1, + 1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.6,1X,F10.6,1X,F4.1,1X,F6.3,1X, + F5.3) C Effective file loading open(unit=1,status='old',file= +'s11.dat') write(6,*) '....Loading file: s11.dat' do i__=1,1045 read(1,'(A110)')ar__3 read(ar__3,4) + RAdeg_3(i__),DEdeg_3(i__),Band_3(i__),mag_3(i__),e_mag_3(i__), + Obs_date_3(i__),Obs_time_3(i__),S0L1_3(i__),A0F1_3(i__), + SAT_3(i__),COL_3(i__),MUX_3(i__),GHO_3(i__),RASdeg_3(i__), + DESdeg_3(i__),Iband_3(i__),IRAC_3(i__),e_IRAC_3(i__) c ..............Just test output........... write(6,4) + RAdeg_3(i__),DEdeg_3(i__),Band_3(i__),mag_3(i__),e_mag_3(i__), + Obs_date_3(i__),Obs_time_3(i__),S0L1_3(i__),A0F1_3(i__), + SAT_3(i__),COL_3(i__),MUX_3(i__),GHO_3(i__),RASdeg_3(i__), + DESdeg_3(i__),Iband_3(i__),IRAC_3(i__),e_IRAC_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'l15.dat' ! Photometry in L15 (15.0um) band C Format for file interpretation 5 format( + F9.6,1X,F10.6,1X,A3,1X,F6.3,3X,F5.3,1X,A10,1X,A8,1X,I1,1X,I1, + 1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.6,1X,F10.6,1X,F4.1,1X,F6.3,1X, + F5.3) C Effective file loading open(unit=1,status='old',file= +'l15.dat') write(6,*) '....Loading file: l15.dat' do i__=1,479 read(1,'(A110)')ar__4 read(ar__4,5) + RAdeg_4(i__),DEdeg_4(i__),Band_4(i__),mag_4(i__),e_mag_4(i__), + Obs_date_4(i__),Obs_time_4(i__),S0L1_4(i__),A0F1_4(i__), + SAT_4(i__),COL_4(i__),MUX_4(i__),GHO_4(i__),RASdeg_4(i__), + DESdeg_4(i__),Iband_4(i__),IRAC_4(i__),e_IRAC_4(i__) c ..............Just test output........... write(6,5) + RAdeg_4(i__),DEdeg_4(i__),Band_4(i__),mag_4(i__),e_mag_4(i__), + Obs_date_4(i__),Obs_time_4(i__),S0L1_4(i__),A0F1_4(i__), + SAT_4(i__),COL_4(i__),MUX_4(i__),GHO_4(i__),RASdeg_4(i__), + DESdeg_4(i__),Iband_4(i__),IRAC_4(i__),e_IRAC_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'l24.dat' ! Photometry in L24 (24.0um) band C Format for file interpretation 6 format( + F9.6,1X,F10.6,1X,A3,1X,F6.3,3X,F5.3,1X,A10,1X,A8,1X,I1,1X,I1, + 1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.6,1X,F10.6,1X,F4.1,1X,F6.3,1X, + F5.3) C Effective file loading open(unit=1,status='old',file= +'l24.dat') write(6,*) '....Loading file: l24.dat' do i__=1,356 read(1,'(A110)')ar__5 read(ar__5,6) + RAdeg_5(i__),DEdeg_5(i__),Band_5(i__),mag_5(i__),e_mag_5(i__), + Obs_date_5(i__),Obs_time_5(i__),S0L1_5(i__),A0F1_5(i__), + SAT_5(i__),COL_5(i__),MUX_5(i__),GHO_5(i__),RASdeg_5(i__), + DESdeg_5(i__),Iband_5(i__),IRAC_5(i__),e_IRAC_5(i__) c ..............Just test output........... write(6,6) + RAdeg_5(i__),DEdeg_5(i__),Band_5(i__),mag_5(i__),e_mag_5(i__), + Obs_date_5(i__),Obs_time_5(i__),S0L1_5(i__),A0F1_5(i__), + SAT_5(i__),COL_5(i__),MUX_5(i__),GHO_5(i__),RASdeg_5(i__), + DESdeg_5(i__),Iband_5(i__),IRAC_5(i__),e_IRAC_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end