FORTRAN Generation
(/./ftp/cats/VII/92A/)

Conversion of standardized ReadMe file for file /./ftp/cats/VII/92A/ 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 2017-Sep-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. VII/92A             Open Cluster Data 5th Edition    (Lynga 1987)
*================================================================================
*Catalogue of Open Cluster Data (5th Ed.)
*     LYNGA G.: 1987
*     <Lund Observatory>
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'catalog'	! The Catalogue Data File

      integer*4 nr__
      parameter (nr__=1151)	! Number of records
      character*514 ar__   	! Full-size record

C  1950 position composed of: RAh RAm DE- DEd DEm
      real*8        RAdeg      (nr__) ! (deg) Right Ascension 1950
      real*8        DEdeg      (nr__) ! (deg)     Declination 1950
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     ClSeq      (nr__) ! [1/64] Cluster Sequence code (1)
      integer*4     ClNum      (nr__) ! Number inside Cluster Sequence
      integer*4     RA2000h    (nr__) ! (h) Right Ascension 2000 (hours)
      real*4        RA2000m    (nr__) ! (min) Right Ascension 2000 (minutes)
      character*1   DE2000_    (nr__) ! Declination 2000 (sign)
      integer*4     DE2000d    (nr__) ! (deg) Declination 2000 (degrees)
      integer*4     DE2000m    (nr__) ! (arcmin) Declination 2000 (minutes)
      integer*4     RAh        (nr__) ! (h) Right Ascension 1950 (hours)
      real*4        RAm        (nr__) ! (min) Right Ascension 1950 (minutes)
      character*1   DE_        (nr__) ! Declination 1950 (sign)
      integer*4     DEd        (nr__) ! (deg) Declination 1950 (degrees)
      integer*4     DEm        (nr__) ! (arcmin) Declination 1950 (minutes)
      real*4        precRA     (nr__) ! (0.1min/yr) Precession in RA (for 10 years)
      real*4        precDE     (nr__) ! (0.1arcmin/yr) Precession in DE (for 10 years)
      real*4        GLON       (nr__) ! (deg) Galactic longitude
      real*4        GLAT       (nr__) ! (deg) Galactic latitude
      integer*4     LundRec    (nr__) ! [1/1212]+ Record number on Lund disc
      integer*4     OCL        (nr__) ! [1/1039]? OCL number (ref.19)
      real*4        Diam       (nr__) ! (arcmin) ? Selected Angular diameter
      integer*4     r_Diam     (nr__) ! [-1/527]? Reference for angular
*                                               diameter (7)
      integer*4     Dist       (nr__) ! (pc) ? Distance
      integer*4     r_Dist     (nr__) ! [1/527]? reference for Distance
      real*4        log_Age    (nr__) ! ([yr]) [6/10]? log(age), years
      integer*4     r_log_Age  (nr__) ! [1/527]? reference for log(age)
      real*4        v_Fe_H_    (nr__) ! ([Sun]) [-.8/.3]? Metallicity, logarithmic scale
      integer*4     r__Fe_H_   (nr__) ! [1/527]? reference for metallicity
      real*4        E_B_V      (nr__) ! (mag) [0/2]? Reddening
      character*2   r_E_B_V    (nr__) ! [-*0-9] reference for Reddening (5)
      character*2   ClTyp      (nr__) ! [DO] "DO" indicates at doubtful cluster
      integer*4     r_ClTyp    (nr__) ! [19]? reference for ClTyp
      integer*4     N_Ap       (nr__) ! Number of Ap stars
      integer*4     N_Am       (nr__) ! Number of Am stars
      integer*4     N_mAp      (nr__) ! Number of marginal Ap stars
      integer*4     N_mAm      (nr__) ! Number of marginal Am stars
      integer*4     N_WR       (nr__) ! Number of WR stars
      integer*4     N_Of       (nr__) ! Number of Of stars
      integer*4     N_Be       (nr__) ! Number of Be stars
      integer*4     N_sh       (nr__) ! Number of Shell stars
      integer*4     N_Be_sh    (nr__) ! Number of Be+Shell stars
      integer*4     N_He_weak  (nr__) ! Number of He weak stars
      integer*4     N_He_rich  (nr__) ! Number of He rich stars
      integer*4     N_Cep      (nr__) ! Number of Cepheids
      integer*4     N_PN       (nr__) ! Number of Planetary nebulae
      integer*4     N_C_       (nr__) ! Number of Carbon stars
      integer*4     N_LPV      (nr__) ! Number of Long period variable stars
      integer*4     N_dSct     (nr__) ! Number of Delta Scuti stars
      integer*4     N_bCep     (nr__) ! Number of Beta Cephei stars
      integer*4     N_RCrB     (nr__) ! Number of R Cor Bor stars
      integer*4     N_BaII     (nr__) ! Number of Ba II stars
      integer*4     N_UGem     (nr__) ! Number of ? U Geminorum stars
      character*20  refN       (nr__) ! References for numbers of
*                                              peculiar stars (format 5I4)
*--------------------------------------------------------------------------------
*        Part of Dixon's catalogue of non-stellar objects (ref.128)
*--------------------------------------------------------------------------------
      character*2   POSS_type  (nr__) ! ? POSS chart type (So/O/E=1/2/3) (2)
      integer*4     POSS       (nr__) ! ? POSS chart number (2)
      integer*4     RAPOSSh    (nr__) ! (h) ? Right Ascension POSS (hours) (2)
      integer*4     RAPOSSm    (nr__) ! (min) ? Right Ascension POSS (minutes) (2)
      integer*4     DEPOSSd    (nr__) ! (deg) ? Declination POSS (degrees) (2)
      integer*4     xPOSS      (nr__) ! (mm) ? Position position on POSS chart, mm (2)
      integer*4     yPOSS      (nr__) ! (mm) ? Position position on POSS chart, mm (2)
      integer*4     maxDiam    (nr__) ! (arcmin) ? Max angular diam, min of arc (2)
      real*4        maxBr      (nr__) ! (mag) ? Brightest star (2)
      integer*4     TrConc     (nr__) ! [0/4]? Trumpler concentration class
      integer*4     TrRange    (nr__) ! [0/3]? Trumpler Range class
      character*2   TrRich     (nr__) ! [mpr ] Trumpler Richness class
      character*1   TrNeb      (nr__) ! [nue ] Trumpler nebulosity (n/u/e=1/2/3)
      integer*4     inRef_128  (nr__) ! [0/1]? '1' if cluster is in (ref.128) (2)
      integer*4     POSS_128   (nr__) ! [0/1]? '1' if cluster has POSS data
*                                            in (ref.128) (2)
*--------------------------------------------------------------------------------
*        Part of Strasbourg-Bochum file (ref.50, Cat. <VII/5>):
*--------------------------------------------------------------------------------
      integer*4     Diam_50    (nr__) ! (arcmin) ? Selected angular diameter (2)
      integer*4     Dist_50    (nr__) ! (pc) ? Distance (2)
      real*4        maxBr_50   (nr__) ! (mag) ? Brightest star
      integer*4     Sp_50      (nr__) ! [10005/10308]? Spectral class code
      real*4        totMag_50  (nr__) ! (mag) ? Total magnitude
      real*4        E_B_V__50  (nr__) ! (mag) ? Colour excess (2)
      integer*4     inRef_50   (nr__) ! [0/1]?  '1' if cluster is in (ref.50) (2)
*--------------------------------------------------------------------------------
*        Data received from B. Skiff (ref.422):
*--------------------------------------------------------------------------------
      real*4        totMag_422 (nr__) ! (mag) ? Total magnitude
      real*4        i_B_V__422 (nr__) ! (mag) ? Integrated color (B-V)
      integer*4     N_422      (nr__) ! ? Number of stars considered
*--------------------------------------------------------------------------------
*        Data received from K. Janes (ref.265):
*--------------------------------------------------------------------------------
      integer*4     Num_265    (nr__) ! [0/938]? Janes-Adler star number
      integer*4     Class_265  (nr__) ! [0/5]?   Janes-Adler class
      integer*4     Ref_265    (nr__) ! ? Janes-Adler references (2)
      integer*4     maxClass_265(nr__) ! [0/5]? Janes-Adler maximum class
      integer*4     Rich_265   (nr__) ! [0/5]? Janes-Adler richness class
      real*4        E_B_V__265 (nr__) ! ? Janes-Adler Color Excess E(B-V)
      integer*4     r_E_B_V__265(nr__) ! [1/527]? Janes-Adler ref. for E(B-V)
      real*4        turn_265   (nr__) ! (mag) ? Janes-Adler (B-V) turn-off point
      integer*4     r_turn_265 (nr__) ! [1/527]? Janes-Adler ref for (B-V)
*                                              turn-off point
*--------------------------------------------------------------------------------
*        Results of 1982 examination (ref.320)
*--------------------------------------------------------------------------------
      integer*4     TrConc_320 (nr__) ! [0/4]? Trumpler concentration class
      integer*4     TrRange_320(nr__) ! [0/3]? Trumpler range class
      character*2   TrRich_320 (nr__) ! [mpr ] Trumpler richness class
      character*1   TrNeb_320  (nr__) ! [n-]   Trumpler nebulosity class
      integer*4     Members_320(nr__) ! ? Number of member stars
      real*4        Diam_320   (nr__) ! (arcmin) ? Angular diameter
*--------------------------------------------------------------------------------
*        Chart used by Lynga (ref.320) or, if not in that survey,
*--------------------------------------------------------------------------------
      character*2   POSS_type_320(nr__) ! [sorjeir ] Survey chart type (3)
      integer*4     POSS_320   (nr__) ! [0/1619]? chart number
      integer*4     RAPOSS_320h(nr__) ! (h) ? Right Ascension POSS (hours)
      integer*4     RAPOSS_320m(nr__) ! (min) ? Right Ascension POSS (minutes)
      integer*4     DEPOSS_320d(nr__) ! (deg) ? Declination POSS (degrees)
      integer*4     xPOSS_320  (nr__) ! (mm) [0/350]? distance from left edge of chart
      integer*4     yPOSS_320  (nr__) ! (mm) [0/350]? distance from bottom edge of chart
      integer*4     hasClass_320(nr__) ! [0/1]?  '1' if cluster classified
*                                             in ref.320
      integer*4     ESO_SERC   (nr__) ! [0/1]?  '1' if cluster is on ESO/SERC
*                                             Survey
      real*4        w_RVel     (nr__) ! [0/1.0]? Weights for radial velocity (6)
      integer*4     RVel       (nr__) ! (km/s) ? Radial Velocity
      integer*4     RVelClass  (nr__) ! [0/5]? Weight class for radial velocity
      character*20  r_RVel     (nr__) ! references for radial velocities
*                                                 (format 5I4)
*--------------------------------------------------------------------------------
*        Buscombe Compilation (ref.64):
*--------------------------------------------------------------------------------
      integer*4     RVel_64    (nr__) ! (km/s) ? Radial Velocity (2)
      integer*4     Dist_64    (nr__) ! (pc) ? Distance (2)
      real*4        E_B_V__64  (nr__) ! (mag) ? E(B-V) (2)
      real*4        iMag_64    (nr__) ! (mag) ? Integrated apparent pg. magnitude (2)
      real*4        iAbsMag_64 (nr__) ! (mag) ? Integrated absolute Bmag (2)
      integer*4     Diam_64    (nr__) ! (pc) ? Absolute Diameter (2)
      real*4        log_age__64(nr__) ! ([yr]) ? log(age), years (2)
      real*4        log_Mass_64(nr__) ! ([solMass]) ? log (Mass), solar masses (2)
      real*4        inRef_64   (nr__) ! [0/1]? '1' if cluster is in (ref.64) (2)
*--------------------------------------------------------------------------------
*        Basel data (ref.37):
*--------------------------------------------------------------------------------
      character*2   phSystem_37(nr__) ! Photometric system for following data (2)
      real*4        E_B_V__37  (nr__) ! (mag) ? Colour excess (2)
      real*4        Av_37      (nr__) ! (mag) ? Interstellar extinction (2)
      integer*4     Dist_37    (nr__) ! (pc) ? Distance (2)
      integer*4     Sp_37      (nr__) ! [0/11000]? Spectral type of hottest star
      integer*4     ColBlue_37 (nr__) ! [0/10320]? Colour type of bluest star
      real*4        Diam_37    (nr__) ! (pc) ? Cluster diameter (2)
      integer*4     ObsType_37 (nr__) ! ? Type of observation (2)
      integer*4     inRef_37   (nr__) ! [0/1]? '1' if cluster is in (ref.37) (2)
*--------------------------------------------------------------------------------
*        van den Bergh-Hagen data (ref.487):
*--------------------------------------------------------------------------------
      real*4        Diam_487   (nr__) ! (arcmin) ? Apparent Diameter (2)
      character*2   Rich_487   (nr__) ! Richness: vp, p, m, r, vr (2)
      integer*4     inRef_487  (nr__) ! [0/1]? '1' if cluster is in (ref.487) (2)
      character*9   Cname      (nr__) ! [C0-9+-] Cluster "C" designation (4)
*--------------------------------------------------------------------------------
*        Various tracers:
*--------------------------------------------------------------------------------
      integer*4     tr_RA      (nr__) ! [0/1212]? tracer for RA, 1950.0 (8)
      integer*4     tr_GLON    (nr__) ! [0/1212]? tracer for Galactic longitude (8)
      integer*4     tr_ClNum   (nr__) ! [0/1212]? tracer for number in sequence (8)
      integer*4     tr_RA_1    (nr__) ! [0/1212]? tracer for RA, 1950.0 (8)
      integer*4     tr_GLON_1  (nr__) ! [0/1212]? tracer for Galactic longitude (8)
      integer*4     tr_ClNum_1 (nr__) ! [0/1212]? tracer for number in sequence (8)
*--------------------------------------------------------------------------------
*        Data from joint project between Janes, Tilley and Lynga (1988):
*--------------------------------------------------------------------------------
      integer*4     Dist_jdl   (nr__) ! (pc) ? jdl distance modulus
      integer*4     w_Dist_jdl (nr__) ! [0/5]? weight class
      real*4        turn_jdl   (nr__) ! (mag) [-0.4/0.7]? jdl turn-off point colour
      real*4        Age_jdl    (nr__) ! (Myr) [0/10000]? jdl derived age
      integer*4     w_Age_jdl  (nr__) ! [0/5]? weight class
      real*4        E_B_V__jdl (nr__) ! (mag) [0/1.75]? jdl reddening
      character*1   n_E_B_V__jdl(nr__) ! [v] 'v' if reddening varies
      integer*4     w_E_B_V__jdl(nr__) ! [0/5]? weight class
      real*4        v_Fe_H__jdl(nr__) ! (Sun) [0]? jdl abundance value
      character*18  refs_jdl   (nr__) ! references for jdl input
*                                              (format 6I3)
*Note (1): The Cluster Sequence Code has the following values:
*     1 NGC            2 IC              3 Berkeley          4 Czernik
*     5 Dolidze        6 Collinder       7 Upgren            8 Tombaugh
*     9 Ruprecht      10 King           11 Stock
*    13 Trumpler      14 Markarian                          16 Haffner
*    17 Hogg          18 Sher           19 Feinstein        20 Harvard
*    21 Lynga         22 Westerlund     23 Basel            24 Blanco
*    25 Baractova     26 Biurakan       27 Melotte          28 Pismis
*                     30 Trapezium                          32 Pleiades
*    33 Graff         34 Iskudarian     35 Stephenson       36 Roslund
*    37 Hyades                          41 van den Bergh-Hagen
*                     42 Bochum         43 Dolidze-Dzimselejsvili
*    45 Antalova      46 Moffat         47 Havlen-Moffat    48 Frolov
*                     50 van den Bergh  51 Mayer            52 Latysev
*    53 Sigma Ori     54 Graham         55 Aveni-Hunter     56 Loden
*    57 Grasdalen     58 Waterloo       59 Auner
*    61 Schuster      62 Danks          63 Muzzio           64 =ref.378
*Note (2): this field is always blank.
*Note (3): the survey may be SO, O, E, J, R, I or SR
*Note (4): the "C" cluster designation has the general "format"
*        CHHMM+DDd, i.e.  starts with the letter 'C',
*        followed by the right ascension (1950) in hours and minutes.
*        followed by the declination sign, degree and tenth of degree.
*Note (5): this field is apparently in error, it contains frequently
*        "-1" or "**"
*Note (6): w_RVel has only the values 0.0 or 1.0
*Note (7): the unexplained reference (-1) exists for the clusters
*          C0728-168, C0728-169, C1726-324
*Note (8): tracers are indexes of the previous or next entry according
*          to the parameter; it has been omitted for C2357+606

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'refs'	! Reference File

      integer*4 nr__1
      parameter (nr__1=670)	! Number of records
      character*132 ar__1  	! Full-size record

      integer*4     RefNum     (nr__1) ! [1/527]+= Reference number
      character*126 Text       (nr__1) ! Text of reference (may be on several lines)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'galactic'	! Galactic Coordinate File

      integer*4 nr__2
      parameter (nr__2=1111)	! Number of records
      character*54 ar__2  	! Full-size record

C  Position composed of: GLON GLAT
      real*4        GLON_1     (nr__2) ! (deg) Galactic longitude
      real*4        GLAT_1     (nr__2) ! (deg) Galactic latitude
      integer*4     OCL_1      (nr__2) ! [1/1039]? OCL number (ref.19)
      character*9   Cname_1    (nr__2) ! Cluster name
      character*15  Name       (nr__2) ! Cluster Sequence name and number
      integer*4     LundRec_1  (nr__2) ! [1/1212] Record number of Lund disc

C=============================================================================

C  Loading file 'catalog'	! The Catalogue Data File

C  Format for file interpretation

    1 format(
     +  I2,I4,1X,I2,F4.1,A1,I2,I2,1X,I2,F4.1,A1,I2,I2,F4.2,F4.1,F6.2,
     +  F6.2,I4,I4,F6.1,I4,I4,1X,I3,F5.2,I3,F5.2,I3,F6.2,A2,A2,I2,I2,
     +  I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,I2,A20,
     +  A2,I4,1X,I2,1X,I2,1X,I3,I3,I3,I4,F4.1,I2,I2,A2,1X,A1,1X,I1,1X,
     +  I1,I4,I4,F4.1,1X,I5,F4.1,F4.2,I2,F4.1,F4.2,I4,I4,1X,I1,1X,I3,
     +  1X,I1,1X,I1,F4.2,1X,I3,F4.2,1X,I3,1X,I1,1X,I1,A2,1X,A1,1X,I3,
     +  1X,F5.1,A2,I4,1X,I2,1X,I2,1X,I3,I3,I3,1X,I1,1X,I1,1X,F3.1,I4,
     +  1X,I1,A20,I4,I4,F4.2,F4.1,F4.1,I4,F4.2,F4.2,F2.1,A2,F4.2,F4.2,
     +  I4,I6,I6,F6.1,I2,I2,F4.1,A2,1X,I1,1X,A9,I4,I4,I4,I4,I4,I4,1X,
     +  I5,1X,I1,F5.2,F6.0,1X,I1,F5.2,A1,1X,I1,F5.2,A18)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog')
      write(6,*) '....Loading file: catalog'
      do i__=1,1151
        read(1,'(A514)')ar__
        read(ar__,1)
     +  ClSeq(i__),ClNum(i__),RA2000h(i__),RA2000m(i__),DE2000_(i__),
     +  DE2000d(i__),DE2000m(i__),RAh(i__),RAm(i__),DE_(i__),DEd(i__),
     +  DEm(i__),precRA(i__),precDE(i__),GLON(i__),GLAT(i__),
     +  LundRec(i__),OCL(i__),Diam(i__),r_Diam(i__),Dist(i__),
     +  r_Dist(i__),log_Age(i__),r_log_Age(i__),v_Fe_H_(i__),
     +  r__Fe_H_(i__),E_B_V(i__),r_E_B_V(i__),ClTyp(i__),r_ClTyp(i__),
     +  N_Ap(i__),N_Am(i__),N_mAp(i__),N_mAm(i__),N_WR(i__),N_Of(i__),
     +  N_Be(i__),N_sh(i__),N_Be_sh(i__),N_He_weak(i__),
     +  N_He_rich(i__),N_Cep(i__),N_PN(i__),N_C_(i__),N_LPV(i__),
     +  N_dSct(i__),N_bCep(i__),N_RCrB(i__),N_BaII(i__),N_UGem(i__),
     +  refN(i__),POSS_type(i__),POSS(i__),RAPOSSh(i__),RAPOSSm(i__),
     +  DEPOSSd(i__),xPOSS(i__),yPOSS(i__),maxDiam(i__),maxBr(i__),
     +  TrConc(i__),TrRange(i__),TrRich(i__),TrNeb(i__),
     +  inRef_128(i__),POSS_128(i__),Diam_50(i__),Dist_50(i__),
     +  maxBr_50(i__),Sp_50(i__),totMag_50(i__),E_B_V__50(i__),
     +  inRef_50(i__),totMag_422(i__),i_B_V__422(i__),N_422(i__),
     +  Num_265(i__),Class_265(i__),Ref_265(i__),maxClass_265(i__),
     +  Rich_265(i__),E_B_V__265(i__),r_E_B_V__265(i__),turn_265(i__),
     +  r_turn_265(i__),TrConc_320(i__),TrRange_320(i__),
     +  TrRich_320(i__),TrNeb_320(i__),Members_320(i__),Diam_320(i__),
     +  POSS_type_320(i__),POSS_320(i__),RAPOSS_320h(i__),
     +  RAPOSS_320m(i__),DEPOSS_320d(i__),xPOSS_320(i__),
     +  yPOSS_320(i__),hasClass_320(i__),ESO_SERC(i__),w_RVel(i__),
     +  RVel(i__),RVelClass(i__),r_RVel(i__),RVel_64(i__),
     +  Dist_64(i__),E_B_V__64(i__),iMag_64(i__),iAbsMag_64(i__),
     +  Diam_64(i__),log_age__64(i__),log_Mass_64(i__),inRef_64(i__),
     +  phSystem_37(i__),E_B_V__37(i__),Av_37(i__),Dist_37(i__),
     +  Sp_37(i__),ColBlue_37(i__),Diam_37(i__),ObsType_37(i__),
     +  inRef_37(i__),Diam_487(i__),Rich_487(i__),inRef_487(i__),
     +  Cname(i__),tr_RA(i__),tr_GLON(i__),tr_ClNum(i__),tr_RA_1(i__),
     +  tr_GLON_1(i__),tr_ClNum_1(i__),Dist_jdl(i__),w_Dist_jdl(i__),
     +  turn_jdl(i__),Age_jdl(i__),w_Age_jdl(i__),E_B_V__jdl(i__),
     +  n_E_B_V__jdl(i__),w_E_B_V__jdl(i__),v_Fe_H__jdl(i__),
     +  refs_jdl(i__)
        if(ar__(55:58) .EQ. '') OCL(i__) = iNULL__
        if(ar__(59:64) .EQ. '') Diam(i__) = rNULL__
        if(ar__(65:68) .EQ. '') r_Diam(i__) = iNULL__
        if(ar__(69:72) .EQ. '') Dist(i__) = iNULL__
        if(ar__(74:76) .EQ. '') r_Dist(i__) = iNULL__
        if(ar__(77:81) .EQ. '') log_Age(i__) = rNULL__
        if(ar__(82:84) .EQ. '') r_log_Age(i__) = iNULL__
        if(ar__(85:89) .EQ. '') v_Fe_H_(i__) = rNULL__
        if(ar__(90:92) .EQ. '') r__Fe_H_(i__) = iNULL__
        if(ar__(93:98) .EQ. '') E_B_V(i__) = rNULL__
        if(ar__(103:104) .EQ. '') r_ClTyp(i__) = iNULL__
        if(ar__(167:170) .EQ. '') POSS(i__) = iNULL__
        if(ar__(172:173) .EQ. '') RAPOSSh(i__) = iNULL__
        if(ar__(175:176) .EQ. '') RAPOSSm(i__) = iNULL__
        if(ar__(178:180) .EQ. '') DEPOSSd(i__) = iNULL__
        if(ar__(181:183) .EQ. '') xPOSS(i__) = iNULL__
        if(ar__(184:186) .EQ. '') yPOSS(i__) = iNULL__
        if(ar__(187:190) .EQ. '') maxDiam(i__) = iNULL__
        if(ar__(191:194) .EQ. '') maxBr(i__) = rNULL__
        if(ar__(195:196) .EQ. '') TrConc(i__) = iNULL__
        if(ar__(197:198) .EQ. '') TrRange(i__) = iNULL__
        if(ar__(204:204) .EQ. '') inRef_128(i__) = iNULL__
        if(ar__(206:206) .EQ. '') POSS_128(i__) = iNULL__
        if(ar__(207:210) .EQ. '') Diam_50(i__) = iNULL__
        if(ar__(211:214) .EQ. '') Dist_50(i__) = iNULL__
        if(ar__(215:218) .EQ. '') maxBr_50(i__) = rNULL__
        if(ar__(220:224) .EQ. '') Sp_50(i__) = iNULL__
        if(ar__(225:228) .EQ. '') totMag_50(i__) = rNULL__
        if(ar__(229:232) .EQ. '') E_B_V__50(i__) = rNULL__
        if(ar__(233:234) .EQ. '') inRef_50(i__) = iNULL__
        if(ar__(235:238) .EQ. '') totMag_422(i__) = rNULL__
        if(ar__(239:242) .EQ. '') i_B_V__422(i__) = rNULL__
        if(ar__(243:246) .EQ. '') N_422(i__) = iNULL__
        if(ar__(247:250) .EQ. '') Num_265(i__) = iNULL__
        if(ar__(252:252) .EQ. '') Class_265(i__) = iNULL__
        if(ar__(254:256) .EQ. '') Ref_265(i__) = iNULL__
        if(ar__(258:258) .EQ. '') maxClass_265(i__) = iNULL__
        if(ar__(260:260) .EQ. '') Rich_265(i__) = iNULL__
        if(ar__(261:264) .EQ. '') E_B_V__265(i__) = rNULL__
        if(ar__(266:268) .EQ. '') r_E_B_V__265(i__) = iNULL__
        if(ar__(269:272) .EQ. '') turn_265(i__) = rNULL__
        if(ar__(274:276) .EQ. '') r_turn_265(i__) = iNULL__
        if(ar__(278:278) .EQ. '') TrConc_320(i__) = iNULL__
        if(ar__(280:280) .EQ. '') TrRange_320(i__) = iNULL__
        if(ar__(286:288) .EQ. '') Members_320(i__) = iNULL__
        if(ar__(290:294) .EQ. '') Diam_320(i__) = rNULL__
        if(ar__(297:300) .EQ. '') POSS_320(i__) = iNULL__
        if(ar__(302:303) .EQ. '') RAPOSS_320h(i__) = iNULL__
        if(ar__(305:306) .EQ. '') RAPOSS_320m(i__) = iNULL__
        if(ar__(308:310) .EQ. '') DEPOSS_320d(i__) = iNULL__
        if(ar__(311:313) .EQ. '') xPOSS_320(i__) = iNULL__
        if(ar__(314:316) .EQ. '') yPOSS_320(i__) = iNULL__
        if(ar__(318:318) .EQ. '') hasClass_320(i__) = iNULL__
        if(ar__(320:320) .EQ. '') ESO_SERC(i__) = iNULL__
        if(ar__(322:324) .EQ. '') w_RVel(i__) = rNULL__
        if(ar__(325:328) .EQ. '') RVel(i__) = iNULL__
        if(ar__(330:330) .EQ. '') RVelClass(i__) = iNULL__
        if(ar__(351:354) .EQ. '') RVel_64(i__) = iNULL__
        if(ar__(355:358) .EQ. '') Dist_64(i__) = iNULL__
        if(ar__(359:362) .EQ. '') E_B_V__64(i__) = rNULL__
        if(ar__(363:366) .EQ. '') iMag_64(i__) = rNULL__
        if(ar__(367:370) .EQ. '') iAbsMag_64(i__) = rNULL__
        if(ar__(371:374) .EQ. '') Diam_64(i__) = iNULL__
        if(ar__(375:378) .EQ. '') log_age__64(i__) = rNULL__
        if(ar__(379:382) .EQ. '') log_Mass_64(i__) = rNULL__
        if(ar__(383:384) .EQ. '') inRef_64(i__) = rNULL__
        if(ar__(387:390) .EQ. '') E_B_V__37(i__) = rNULL__
        if(ar__(391:394) .EQ. '') Av_37(i__) = rNULL__
        if(ar__(395:398) .EQ. '') Dist_37(i__) = iNULL__
        if(ar__(399:404) .EQ. '') Sp_37(i__) = iNULL__
        if(ar__(405:410) .EQ. '') ColBlue_37(i__) = iNULL__
        if(ar__(411:416) .EQ. '') Diam_37(i__) = rNULL__
        if(ar__(417:418) .EQ. '') ObsType_37(i__) = iNULL__
        if(ar__(419:420) .EQ. '') inRef_37(i__) = iNULL__
        if(ar__(421:424) .EQ. '') Diam_487(i__) = rNULL__
        if(ar__(428:428) .EQ. '') inRef_487(i__) = iNULL__
        if(ar__(439:442) .EQ. '') tr_RA(i__) = iNULL__
        if(ar__(443:446) .EQ. '') tr_GLON(i__) = iNULL__
        if(ar__(447:450) .EQ. '') tr_ClNum(i__) = iNULL__
        if(ar__(451:454) .EQ. '') tr_RA_1(i__) = iNULL__
        if(ar__(455:458) .EQ. '') tr_GLON_1(i__) = iNULL__
        if(ar__(459:462) .EQ. '') tr_ClNum_1(i__) = iNULL__
        if(ar__(464:468) .EQ. '') Dist_jdl(i__) = iNULL__
        if(ar__(470:470) .EQ. '') w_Dist_jdl(i__) = iNULL__
        if(ar__(471:475) .EQ. '') turn_jdl(i__) = rNULL__
        if(ar__(476:481) .EQ. '') Age_jdl(i__) = rNULL__
        if(ar__(483:483) .EQ. '') w_Age_jdl(i__) = iNULL__
        if(ar__(484:488) .EQ. '') E_B_V__jdl(i__) = rNULL__
        if(ar__(491:491) .EQ. '') w_E_B_V__jdl(i__) = iNULL__
        if(ar__(492:496) .EQ. '') v_Fe_H__jdl(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  ClSeq(i__),ClNum(i__),RA2000h(i__),RA2000m(i__),DE2000_(i__),
     +  DE2000d(i__),DE2000m(i__),RAh(i__),RAm(i__),DE_(i__),DEd(i__),
     +  DEm(i__),precRA(i__),precDE(i__),GLON(i__),GLAT(i__),
     +  LundRec(i__),OCL(i__),Diam(i__),r_Diam(i__),Dist(i__),
     +  r_Dist(i__),log_Age(i__),r_log_Age(i__),v_Fe_H_(i__),
     +  r__Fe_H_(i__),E_B_V(i__),r_E_B_V(i__),ClTyp(i__),r_ClTyp(i__),
     +  N_Ap(i__),N_Am(i__),N_mAp(i__),N_mAm(i__),N_WR(i__),N_Of(i__),
     +  N_Be(i__),N_sh(i__),N_Be_sh(i__),N_He_weak(i__),
     +  N_He_rich(i__),N_Cep(i__),N_PN(i__),N_C_(i__),N_LPV(i__),
     +  N_dSct(i__),N_bCep(i__),N_RCrB(i__),N_BaII(i__),N_UGem(i__),
     +  refN(i__),POSS_type(i__),POSS(i__),RAPOSSh(i__),RAPOSSm(i__),
     +  DEPOSSd(i__),xPOSS(i__),yPOSS(i__),maxDiam(i__),maxBr(i__),
     +  TrConc(i__),TrRange(i__),TrRich(i__),TrNeb(i__),
     +  inRef_128(i__),POSS_128(i__),Diam_50(i__),Dist_50(i__),
     +  maxBr_50(i__),Sp_50(i__),totMag_50(i__),E_B_V__50(i__),
     +  inRef_50(i__),totMag_422(i__),i_B_V__422(i__),N_422(i__),
     +  Num_265(i__),Class_265(i__),Ref_265(i__),maxClass_265(i__),
     +  Rich_265(i__),E_B_V__265(i__),r_E_B_V__265(i__),turn_265(i__),
     +  r_turn_265(i__),TrConc_320(i__),TrRange_320(i__),
     +  TrRich_320(i__),TrNeb_320(i__),Members_320(i__),Diam_320(i__),
     +  POSS_type_320(i__),POSS_320(i__),RAPOSS_320h(i__),
     +  RAPOSS_320m(i__),DEPOSS_320d(i__),xPOSS_320(i__),
     +  yPOSS_320(i__),hasClass_320(i__),ESO_SERC(i__),w_RVel(i__),
     +  RVel(i__),RVelClass(i__),r_RVel(i__),RVel_64(i__),
     +  Dist_64(i__),E_B_V__64(i__),iMag_64(i__),iAbsMag_64(i__),
     +  Diam_64(i__),log_age__64(i__),log_Mass_64(i__),inRef_64(i__),
     +  phSystem_37(i__),E_B_V__37(i__),Av_37(i__),Dist_37(i__),
     +  Sp_37(i__),ColBlue_37(i__),Diam_37(i__),ObsType_37(i__),
     +  inRef_37(i__),Diam_487(i__),Rich_487(i__),inRef_487(i__),
     +  Cname(i__),tr_RA(i__),tr_GLON(i__),tr_ClNum(i__),tr_RA_1(i__),
     +  tr_GLON_1(i__),tr_ClNum_1(i__),Dist_jdl(i__),w_Dist_jdl(i__),
     +  turn_jdl(i__),Age_jdl(i__),w_Age_jdl(i__),E_B_V__jdl(i__),
     +  n_E_B_V__jdl(i__),w_E_B_V__jdl(i__),v_Fe_H__jdl(i__),
     +  refs_jdl(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'refs'	! Reference File

C  Format for file interpretation

    2 format(1X,I3,2X,A126)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs')
      write(6,*) '....Loading file: refs'
      do i__=1,670
        read(1,'(A132)')ar__1
        read(ar__1,2)RefNum(i__),Text(i__)
c    ..............Just test output...........
        write(6,2)RefNum(i__),Text(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'galactic'	! Galactic Coordinate File

C  Format for file interpretation

    3 format(2X,F6.2,2X,F6.2,1X,I4,2X,A9,2X,A15,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'galactic')
      write(6,*) '....Loading file: galactic'
      do i__=1,1111
        read(1,'(A54)')ar__2
        read(ar__2,3)
     +  GLON_1(i__),GLAT_1(i__),OCL_1(i__),Cname_1(i__),Name(i__),
     +  LundRec_1(i__)
        if(ar__2(18:21) .EQ. '') OCL_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  GLON_1(i__),GLAT_1(i__),OCL_1(i__),Cname_1(i__),Name(i__),
     +  LundRec_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end