Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/193/25 into FORTRAN code for reading data files line by line.
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 2026-Sep-07
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/ApJS/193/25 The NAN complex. II. MIPS observations (Rebull+, 2011)
*================================================================================
*The North American and Pelican nebulae.
*II. MIPS observations and analysis.
* Rebull L.M., Guieu S., Stauffer J.R., Hillenbrand L.A., Noriega-Crespo A.,
* Stapelfeldt K.R., Carey S.J., Carpenter J.M., Cole D.M., Padgett D.L.,
* Strom S.E., Wolff S.C.
* <Astrophys. J. Suppl. Ser., 193, 25 (2011)>
* =2011ApJS..193...25R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Complete list of YSOs I. Spitzer measurements
integer*4 nr__
parameter (nr__=1286) ! Number of records
character*187 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*18 v_RGS2011_ ! Source name (HHMMSS.ss+DDMMSS.s on J2000)
real*8 RAdeg ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000)
character*16 v2MASS ! 2MASS name (HHMMSSss+DDMMSSs)
character*23 AName ! Alternate name
character*12 SpType ! MK spectral type
character*22 r_SpType ! Reference for SpType (1)
character*1 l__3_6_ ! Limit flag on 3.6]
real*4 v_3_6_ ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude
real*4 e__3_6_ ! (mag) ? Uncertainty in 3.6mag
character*1 l__4_5_ ! Limit flag on 4.5]
real*4 v_4_5_ ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
real*4 e__4_5_ ! (mag) ? Uncertainty in 4.5mag
character*1 l__5_8_ ! Limit flag on 5.8]
real*4 v_5_8_ ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude
real*4 e__5_8_ ! (mag) ? Uncertainty in 5.8mag
character*1 l__8_0_ ! Limit flag on 8.0]
real*4 v_8_0_ ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
real*4 e__8_0_ ! (mag) ? Uncertainty in 8.0mag
character*1 l__24_ ! Limit flag on 24]
real*4 v_24_ ! (mag) ? Spitzer/MIPS 24 micron band magnitude
real*4 e__24_ ! (mag) ? Uncertainty in 24mag
real*4 v_70_ ! (mag) ? Spitzer/MIPS 70 micron band magnitude
real*4 e__70_ ! (mag) ? Uncertainty in 70mag
*Note (1): References as follows:
* BallyScoville = 1980ApJ...239..121B
* CohenKuhi79 = 1979ApJS...41..743C
* Comeron05 = 2005A&A...430..541C; <[CP2005] NN> in Simbad
* Corbally09 = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad
* Laugalys02 = 2002BaltA..11..205L
* Laugalys06a = 2006BaltA..15..483L (Cat. J/BaltA/15/483;
* <[LSV2006] Area II NNN> in Simbad)
* Laugalys06b = 2006BaltA..15..327L
* MerrillBurwell49 = 1949ApJ...110..387M
* Ogura2002 = 2002AJ....123.2597O (Cat. J/AJ/123/2597;
* <[OSP2002] BRC 31 NN> in Simbad)
* Shevchenko93 = 1993Ap&SS.202..121S
* SIMBAD = http://simbad.u-strasbg.fr/simbad/
* Straizys08a = 2008BaltA..17....1S (Cat. J/BaltA/17/1)
* Straizys08b = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad
* Straizys09 = 2009BaltA..18..141S (Cat. J/BaltA/18/141)
* Straizys89 = 1989VilOB..83....3S
* Straizys93 = 1993BaltA...2..171S (Cat. II/196)
* Straizys99 = 1999BaltA...8..355S
* Terranegra94 = 1994A&AS..104..557T (Cat. J/A+AS/104/557)
* Vill04 = 2004A&A...419..149V (Cat. J/A+A/419/149;
* <Cl* NGC 6996 VBC NNNN> in Simbad)
* Welin73 = 1973A&AS....9..183W
* Zdan90 = 1990Ap&SS.173..309Z
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Complete list of YSOs II. YSO criteria
integer*4 nr__1
parameter (nr__1=1286) ! Number of records
character*82 ar__1 ! Full-size record
character*18 v_RGS2011__1 ! Source name (HHMMSS.ss+DDMMSS.s)
character*8 Clu ! Cluster Name (1)
character*1 Q ! [ABC] YSO quality (A for good confidence)
character*4 cl ! YSO class (2)
character*3 cr1 ! [yes/no] Selected from [3.6] vs [3.6]-[24] CMD (3)
character*3 cr2 ! [yes/no] Not too faint in [3.6] vs [3.6]-[24]
* CMD (3)
character*3 cr3 ! [yes/no] Selected from Ks vs Ks-[24] CMD (3)
character*3 cr4 ! [yes/no] Detected at 70 microns (3)
character*3 cr5 ! [yes/no] Found in Guieu et al. 2009
* Cat. J/ApJ/697/787 ([GRS2009] in Simbad)
character*3 cr6 ! [yes/no] Close to the middle of the nebula (3)
character*3 cr7 ! [yes/no] In a cluster (3)
character*3 cr8 ! [yes/no] Previously identified YSO (3)
character*3 cr9 ! [yes/no] Selected from Ic vs. V-Ic CMD (3)
character*3 cr10 ! [yes/no] Selected from IPHAS r' vs. r'-i' CMD (3)
character*3 cr11 ! [yes/no] A bright IPHAS H{alpha} source (3)
character*3 cr12 ! [yes/no] Selected from IPHAS r'-H{alpha} vs
* r'-i' diagram (3)
*Note (1): Clusters as follows:
* Pelican = IC 5070. See section 6.3 for further details.
* PelHat = Pelican's Hat. North of the Pelican's head, in roughly the same
* location as would be found a hat on its head, is a previously
* unidentified dark lane at 24um, similar to but much smaller than
* the Gulf of Mexico. See section 6.4 for further details.
* Gulf = The "Gulf of Mexico" region is the area of highest extinction in
* the NAN (The North American (NGC 7000) and Pelican (IC 5070)
* Nebulae) complex. See section 6.2 for further details.
*Note (2): In the spirit of Wilking et al. (2001ApJ...551..357W), we define the
* near- to mid-IR slope of the SED,
* {alpha}=dlog{lambda}F_{lambda}_/dlog{lambda}, where {alpha}>0.3 for a
* Class I, 0.3 to -0.3 for a flat-spectrum source, -0.3 to -1.6 for a
* Class II, and <-1.6 for a Class III.
*Note (3): For each individual criterion, a "yes" means the object met that
* given criterion for YSO selection, "no" means it did not, and blank
* means that it was unable to be assessed for that criterion.
* For example, column S3624 indicates whether each object was selected
* as a YSO candidate in the [3.6] vs. [3.6]-[24] diagram ("yes"), or
* it was not ("no"), or, it was missing a 3.6 or 24 micron
* measurement, and thus could not be assessed in the [3.6] vs.
* [3.6]-[24] diagram (blank).
* See section 4.2 for further details on each criterion.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! IRAC-selected YSOs not recovered by our
MIPS-based search
integer*4 nr__2
parameter (nr__2=796) ! Number of records
character*187 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*18 v_RGS2011__2 ! Source name (HHMMSS.ss+DDMMSS.s on J2000)
real*8 RAdeg_1 ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_1 ! (deg) Declination in decimal degrees (J2000)
character*16 v2MASS_1 ! 2MASS name (HHMMSSss+DDMMSSs)
character*23 AName_1 ! Alternate name
character*12 SpType_1 ! MK spectral type
character*22 r_SpType_1 ! Reference for SpType (1)
character*1 l__3_6__1 ! Limit flag on 3.6]
real*4 v_3_6__1 ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude
real*4 e__3_6__1 ! (mag) ? Uncertainty in 3.6mag
character*1 l__4_5__1 ! Limit flag on 4.5]
real*4 v_4_5__1 ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
real*4 e__4_5__1 ! (mag) ? Uncertainty in 4.5mag
character*1 l__5_8__1 ! Limit flag on 5.8]
real*4 v_5_8__1 ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude
real*4 e__5_8__1 ! (mag) ? Uncertainty in 5.8mag
character*1 l__8_0__1 ! Limit flag on 8.0]
real*4 v_8_0__1 ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
real*4 e__8_0__1 ! (mag) ? Uncertainty in 8.0mag
character*1 l__24__1 ! Limit flag on 24]
real*4 v_24__1 ! (mag) ? Spitzer/MIPS 24 micron band magnitude
real*4 e__24__1 ! (mag) ? Uncertainty in 24mag
real*4 v_70__1 ! (mag) ? Spitzer/MIPS 70 micron band magnitude
real*4 e__70__1 ! (mag) ? Uncertainty in 70mag
*Note (1): References as follows:
* BallyScoville = 1980ApJ...239..121B
* CohenKuhi79 = 1979ApJS...41..743C
* Comeron05 = 2005A&A...430..541C; <[CP2005] NN> in Simbad
* Corbally09 = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad
* Laugalys02 = 2002BaltA..11..205L
* Laugalys06a = 2006BaltA..15..483L (Cat. J/BaltA/15/483;
* <[LSV2006] Area II NNN> in Simbad)
* Laugalys06b = 2006BaltA..15..327L
* MerrillBurwell49 = 1949ApJ...110..387M
* Ogura2002 = 2002AJ....123.2597O (Cat. J/AJ/123/2597;
* <[OSP2002] BRC 31 NN> in Simbad)
* Shevchenko93 = 1993Ap&SS.202..121S
* SIMBAD = http://simbad.u-strasbg.fr/simbad/
* Straizys08a = 2008BaltA..17....1S (Cat. J/BaltA/17/1)
* Straizys08b = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad
* Straizys09 = 2009BaltA..18..141S (Cat. J/BaltA/18/141)
* Straizys89 = 1989VilOB..83....3S
* Straizys93 = 1993BaltA...2..171S (Cat. II/196)
* Straizys99 = 1999BaltA...8..355S
* Terranegra94 = 1994A&AS..104..557T (Cat. J/A+AS/104/557)
* Vill04 = 2004A&A...419..149V (Cat. J/A+A/419/149;
* <Cl* NGC 6996 VBC NNNN> in Simbad)
* Welin73 = 1973A&AS....9..183W
* Zdan90 = 1990Ap&SS.173..309Z
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Previously suggested YSOs not recovered by our
MIPS or IRAC searches
integer*4 nr__3
parameter (nr__3=114) ! Number of records
character*187 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*18 v_RGS2011__3 ! Source name (HHMMSS.ss+DDMMSS.s on J2000)
real*8 RAdeg_2 ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_2 ! (deg) Declination in decimal degrees (J2000)
character*16 v2MASS_2 ! 2MASS name (HHMMSSss+DDMMSSs)
character*23 AName_2 ! Alternate name
character*12 SpType_2 ! MK spectral type
character*22 r_SpType_2 ! Reference for SpType (1)
character*1 l__3_6__2 ! Limit flag on 3.6]
real*4 v_3_6__2 ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude
real*4 e__3_6__2 ! (mag) ? Uncertainty in 3.6mag
character*1 l__4_5__2 ! Limit flag on 4.5]
real*4 v_4_5__2 ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
real*4 e__4_5__2 ! (mag) ? Uncertainty in 4.5mag
character*1 l__5_8__2 ! Limit flag on 5.8]
real*4 v_5_8__2 ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude
real*4 e__5_8__2 ! (mag) ? Uncertainty in 5.8mag
character*1 l__8_0__2 ! Limit flag on 8.0]
real*4 v_8_0__2 ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
real*4 e__8_0__2 ! (mag) ? Uncertainty in 8.0mag
character*1 l__24__2 ! Limit flag on 24]
real*4 v_24__2 ! (mag) ? Spitzer/MIPS 24 micron band magnitude
real*4 e__24__2 ! (mag) ? Uncertainty in 24mag
real*4 v_70__2 ! (mag) ? Spitzer/MIPS 70 micron band magnitude
real*4 e__70__2 ! (mag) ? Uncertainty in 70mag
*Note (1): References as follows:
* BallyScoville = 1980ApJ...239..121B
* CohenKuhi79 = 1979ApJS...41..743C
* Comeron05 = 2005A&A...430..541C; <[CP2005] NN> in Simbad
* Corbally09 = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad
* Laugalys02 = 2002BaltA..11..205L
* Laugalys06a = 2006BaltA..15..483L (Cat. J/BaltA/15/483;
* <[LSV2006] Area II NNN> in Simbad)
* Laugalys06b = 2006BaltA..15..327L
* MerrillBurwell49 = 1949ApJ...110..387M
* Ogura2002 = 2002AJ....123.2597O (Cat. J/AJ/123/2597;
* <[OSP2002] BRC 31 NN> in Simbad)
* Shevchenko93 = 1993Ap&SS.202..121S
* SIMBAD = http://simbad.u-strasbg.fr/simbad/
* Straizys08a = 2008BaltA..17....1S (Cat. J/BaltA/17/1)
* Straizys08b = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad
* Straizys09 = 2009BaltA..18..141S (Cat. J/BaltA/18/141)
* Straizys89 = 1989VilOB..83....3S
* Straizys93 = 1993BaltA...2..171S (Cat. II/196)
* Straizys99 = 1999BaltA...8..355S
* Terranegra94 = 1994A&AS..104..557T (Cat. J/A+AS/104/557)
* Vill04 = 2004A&A...419..149V (Cat. J/A+A/419/149;
* <Cl* NGC 6996 VBC NNNN> in Simbad)
* Welin73 = 1973A&AS....9..183W
* Zdan90 = 1990Ap&SS.173..309Z
C=============================================================================
C Loading file 'table3.dat' ! Complete list of YSOs I. Spitzer measurements
C Format for file interpretation
1 format(
+ A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2,
+ 1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
+ F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,1286
read(1,'(A187)')ar__
read(ar__,1)
+ v_RGS2011_,RAdeg,DEdeg,v2MASS,AName,SpType,r_SpType,l__3_6_,
+ v_3_6_,e__3_6_,l__4_5_,v_4_5_,e__4_5_,l__5_8_,v_5_8_,e__5_8_,
+ l__8_0_,v_8_0_,e__8_0_,l__24_,v_24_,e__24_,v_70_,e__70_
if(ar__(119:123) .EQ. '') v_3_6_ = rNULL__
if(ar__(125:128) .EQ. '') e__3_6_ = rNULL__
if(ar__(131:135) .EQ. '') v_4_5_ = rNULL__
if(ar__(137:140) .EQ. '') e__4_5_ = rNULL__
if(ar__(143:147) .EQ. '') v_5_8_ = rNULL__
if(ar__(149:152) .EQ. '') e__5_8_ = rNULL__
if(ar__(155:159) .EQ. '') v_8_0_ = rNULL__
if(ar__(161:164) .EQ. '') e__8_0_ = rNULL__
if(ar__(167:171) .EQ. '') v_24_ = rNULL__
if(ar__(173:176) .EQ. '') e__24_ = rNULL__
if(ar__(178:182) .EQ. '') v_70_ = rNULL__
if(ar__(184:187) .EQ. '') e__70_ = rNULL__
c ..............Just test output...........
write(6,1)
+ v_RGS2011_,RAdeg,DEdeg,v2MASS,AName,SpType,r_SpType,l__3_6_,
+ v_3_6_,e__3_6_,l__4_5_,v_4_5_,e__4_5_,l__5_8_,v_5_8_,e__5_8_,
+ l__8_0_,v_8_0_,e__8_0_,l__24_,v_24_,e__24_,v_70_,e__70_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Complete list of YSOs II. YSO criteria
C Format for file interpretation
2 format(
+ A18,1X,A8,1X,A1,1X,A4,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,
+ A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,1286
read(1,'(A82)')ar__1
read(ar__1,2)
+ v_RGS2011__1,Clu,Q,cl,cr1,cr2,cr3,cr4,cr5,cr6,cr7,cr8,cr9,
+ cr10,cr11,cr12
c ..............Just test output...........
write(6,2)
+ v_RGS2011__1,Clu,Q,cl,cr1,cr2,cr3,cr4,cr5,cr6,cr7,cr8,cr9,
+ cr10,cr11,cr12
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! IRAC-selected YSOs not recovered by our
* MIPS-based search
C Format for file interpretation
3 format(
+ A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2,
+ 1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
+ F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,796
read(1,'(A187)')ar__2
read(ar__2,3)
+ v_RGS2011__2,RAdeg_1,DEdeg_1,v2MASS_1,AName_1,SpType_1,
+ r_SpType_1,l__3_6__1,v_3_6__1,e__3_6__1,l__4_5__1,v_4_5__1,
+ e__4_5__1,l__5_8__1,v_5_8__1,e__5_8__1,l__8_0__1,v_8_0__1,
+ e__8_0__1,l__24__1,v_24__1,e__24__1,v_70__1,e__70__1
if(ar__2(119:123) .EQ. '') v_3_6__1 = rNULL__
if(ar__2(125:128) .EQ. '') e__3_6__1 = rNULL__
if(ar__2(131:135) .EQ. '') v_4_5__1 = rNULL__
if(ar__2(137:140) .EQ. '') e__4_5__1 = rNULL__
if(ar__2(143:147) .EQ. '') v_5_8__1 = rNULL__
if(ar__2(149:152) .EQ. '') e__5_8__1 = rNULL__
if(ar__2(155:159) .EQ. '') v_8_0__1 = rNULL__
if(ar__2(161:164) .EQ. '') e__8_0__1 = rNULL__
if(ar__2(167:171) .EQ. '') v_24__1 = rNULL__
if(ar__2(173:176) .EQ. '') e__24__1 = rNULL__
if(ar__2(178:182) .EQ. '') v_70__1 = rNULL__
if(ar__2(184:187) .EQ. '') e__70__1 = rNULL__
c ..............Just test output...........
write(6,3)
+ v_RGS2011__2,RAdeg_1,DEdeg_1,v2MASS_1,AName_1,SpType_1,
+ r_SpType_1,l__3_6__1,v_3_6__1,e__3_6__1,l__4_5__1,v_4_5__1,
+ e__4_5__1,l__5_8__1,v_5_8__1,e__5_8__1,l__8_0__1,v_8_0__1,
+ e__8_0__1,l__24__1,v_24__1,e__24__1,v_70__1,e__70__1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Previously suggested YSOs not recovered by our
* MIPS or IRAC searches
C Format for file interpretation
4 format(
+ A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2,
+ 1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
+ F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,114
read(1,'(A187)')ar__3
read(ar__3,4)
+ v_RGS2011__3,RAdeg_2,DEdeg_2,v2MASS_2,AName_2,SpType_2,
+ r_SpType_2,l__3_6__2,v_3_6__2,e__3_6__2,l__4_5__2,v_4_5__2,
+ e__4_5__2,l__5_8__2,v_5_8__2,e__5_8__2,l__8_0__2,v_8_0__2,
+ e__8_0__2,l__24__2,v_24__2,e__24__2,v_70__2,e__70__2
if(ar__3(119:123) .EQ. '') v_3_6__2 = rNULL__
if(ar__3(125:128) .EQ. '') e__3_6__2 = rNULL__
if(ar__3(131:135) .EQ. '') v_4_5__2 = rNULL__
if(ar__3(137:140) .EQ. '') e__4_5__2 = rNULL__
if(ar__3(143:147) .EQ. '') v_5_8__2 = rNULL__
if(ar__3(149:152) .EQ. '') e__5_8__2 = rNULL__
if(ar__3(155:159) .EQ. '') v_8_0__2 = rNULL__
if(ar__3(161:164) .EQ. '') e__8_0__2 = rNULL__
if(ar__3(167:171) .EQ. '') v_24__2 = rNULL__
if(ar__3(173:176) .EQ. '') e__24__2 = rNULL__
if(ar__3(178:182) .EQ. '') v_70__2 = rNULL__
if(ar__3(184:187) .EQ. '') e__70__2 = rNULL__
c ..............Just test output...........
write(6,4)
+ v_RGS2011__3,RAdeg_2,DEdeg_2,v2MASS_2,AName_2,SpType_2,
+ r_SpType_2,l__3_6__2,v_3_6__2,e__3_6__2,l__4_5__2,v_4_5__2,
+ e__4_5__2,l__5_8__2,v_5_8__2,e__5_8__2,l__8_0__2,v_8_0__2,
+ e__8_0__2,l__24__2,v_24__2,e__24__2,v_70__2,e__70__2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end