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 2024-Apr-19
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