Conversion of standardized ReadMe file for
file /./ftp/cats/VII/92A/ 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 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. 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 ! (deg) Right Ascension 1950 real*8 DEdeg ! (deg) Declination 1950 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ClSeq ! [1/64] Cluster Sequence code (1) integer*4 ClNum ! Number inside Cluster Sequence integer*4 RA2000h ! (h) Right Ascension 2000 (hours) real*4 RA2000m ! (min) Right Ascension 2000 (minutes) character*1 DE2000_ ! Declination 2000 (sign) integer*4 DE2000d ! (deg) Declination 2000 (degrees) integer*4 DE2000m ! (arcmin) Declination 2000 (minutes) integer*4 RAh ! (h) Right Ascension 1950 (hours) real*4 RAm ! (min) Right Ascension 1950 (minutes) character*1 DE_ ! Declination 1950 (sign) integer*4 DEd ! (deg) Declination 1950 (degrees) integer*4 DEm ! (arcmin) Declination 1950 (minutes) real*4 precRA ! (0.1min/yr) Precession in RA (for 10 years) real*4 precDE ! (0.1arcmin/yr) Precession in DE (for 10 years) real*4 GLON ! (deg) Galactic longitude real*4 GLAT ! (deg) Galactic latitude integer*4 LundRec ! [1/1212]+ Record number on Lund disc integer*4 OCL ! [1/1039]? OCL number (ref.19) real*4 Diam ! (arcmin) ? Selected Angular diameter integer*4 r_Diam ! [-1/527]? Reference for angular * diameter (7) integer*4 Dist ! (pc) ? Distance integer*4 r_Dist ! [1/527]? reference for Distance real*4 log_Age ! ([yr]) [6/10]? log(age), years integer*4 r_log_Age ! [1/527]? reference for log(age) real*4 v_Fe_H_ ! ([Sun]) [-.8/.3]? Metallicity, logarithmic scale integer*4 r__Fe_H_ ! [1/527]? reference for metallicity real*4 E_B_V ! (mag) [0/2]? Reddening character*2 r_E_B_V ! [-*0-9] reference for Reddening (5) character*2 ClTyp ! [DO] "DO" indicates at doubtful cluster integer*4 r_ClTyp ! [19]? reference for ClTyp integer*4 N_Ap ! Number of Ap stars integer*4 N_Am ! Number of Am stars integer*4 N_mAp ! Number of marginal Ap stars integer*4 N_mAm ! Number of marginal Am stars integer*4 N_WR ! Number of WR stars integer*4 N_Of ! Number of Of stars integer*4 N_Be ! Number of Be stars integer*4 N_sh ! Number of Shell stars integer*4 N_Be_sh ! Number of Be+Shell stars integer*4 N_He_weak ! Number of He weak stars integer*4 N_He_rich ! Number of He rich stars integer*4 N_Cep ! Number of Cepheids integer*4 N_PN ! Number of Planetary nebulae integer*4 N_C_ ! Number of Carbon stars integer*4 N_LPV ! Number of Long period variable stars integer*4 N_dSct ! Number of Delta Scuti stars integer*4 N_bCep ! Number of Beta Cephei stars integer*4 N_RCrB ! Number of R Cor Bor stars integer*4 N_BaII ! Number of Ba II stars integer*4 N_UGem ! Number of ? U Geminorum stars character*20 refN ! References for numbers of * peculiar stars (format 5I4) *-------------------------------------------------------------------------------- * Part of Dixon's catalogue of non-stellar objects (ref.128) *-------------------------------------------------------------------------------- character*2 POSS_type ! ? POSS chart type (So/O/E=1/2/3) (2) integer*4 POSS ! ? POSS chart number (2) integer*4 RAPOSSh ! (h) ? Right Ascension POSS (hours) (2) integer*4 RAPOSSm ! (min) ? Right Ascension POSS (minutes) (2) integer*4 DEPOSSd ! (deg) ? Declination POSS (degrees) (2) integer*4 xPOSS ! (mm) ? Position position on POSS chart, mm (2) integer*4 yPOSS ! (mm) ? Position position on POSS chart, mm (2) integer*4 maxDiam ! (arcmin) ? Max angular diam, min of arc (2) real*4 maxBr ! (mag) ? Brightest star (2) integer*4 TrConc ! [0/4]? Trumpler concentration class integer*4 TrRange ! [0/3]? Trumpler Range class character*2 TrRich ! [mpr ] Trumpler Richness class character*1 TrNeb ! [nue ] Trumpler nebulosity (n/u/e=1/2/3) integer*4 inRef_128 ! [0/1]? '1' if cluster is in (ref.128) (2) integer*4 POSS_128 ! [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 ! (arcmin) ? Selected angular diameter (2) integer*4 Dist_50 ! (pc) ? Distance (2) real*4 maxBr_50 ! (mag) ? Brightest star integer*4 Sp_50 ! [10005/10308]? Spectral class code real*4 totMag_50 ! (mag) ? Total magnitude real*4 E_B_V__50 ! (mag) ? Colour excess (2) integer*4 inRef_50 ! [0/1]? '1' if cluster is in (ref.50) (2) *-------------------------------------------------------------------------------- * Data received from B. Skiff (ref.422): *-------------------------------------------------------------------------------- real*4 totMag_422 ! (mag) ? Total magnitude real*4 i_B_V__422 ! (mag) ? Integrated color (B-V) integer*4 N_422 ! ? Number of stars considered *-------------------------------------------------------------------------------- * Data received from K. Janes (ref.265): *-------------------------------------------------------------------------------- integer*4 Num_265 ! [0/938]? Janes-Adler star number integer*4 Class_265 ! [0/5]? Janes-Adler class integer*4 Ref_265 ! ? Janes-Adler references (2) integer*4 maxClass_265 ! [0/5]? Janes-Adler maximum class integer*4 Rich_265 ! [0/5]? Janes-Adler richness class real*4 E_B_V__265 ! ? Janes-Adler Color Excess E(B-V) integer*4 r_E_B_V__265 ! [1/527]? Janes-Adler ref. for E(B-V) real*4 turn_265 ! (mag) ? Janes-Adler (B-V) turn-off point integer*4 r_turn_265 ! [1/527]? Janes-Adler ref for (B-V) * turn-off point *-------------------------------------------------------------------------------- * Results of 1982 examination (ref.320) *-------------------------------------------------------------------------------- integer*4 TrConc_320 ! [0/4]? Trumpler concentration class integer*4 TrRange_320 ! [0/3]? Trumpler range class character*2 TrRich_320 ! [mpr ] Trumpler richness class character*1 TrNeb_320 ! [n-] Trumpler nebulosity class integer*4 Members_320 ! ? Number of member stars real*4 Diam_320 ! (arcmin) ? Angular diameter *-------------------------------------------------------------------------------- * Chart used by Lynga (ref.320) or, if not in that survey, *-------------------------------------------------------------------------------- character*2 POSS_type_320 ! [sorjeir ] Survey chart type (3) integer*4 POSS_320 ! [0/1619]? chart number integer*4 RAPOSS_320h ! (h) ? Right Ascension POSS (hours) integer*4 RAPOSS_320m ! (min) ? Right Ascension POSS (minutes) integer*4 DEPOSS_320d ! (deg) ? Declination POSS (degrees) integer*4 xPOSS_320 ! (mm) [0/350]? distance from left edge of chart integer*4 yPOSS_320 ! (mm) [0/350]? distance from bottom edge of chart integer*4 hasClass_320 ! [0/1]? '1' if cluster classified * in ref.320 integer*4 ESO_SERC ! [0/1]? '1' if cluster is on ESO/SERC * Survey real*4 w_RVel ! [0/1.0]? Weights for radial velocity (6) integer*4 RVel ! (km/s) ? Radial Velocity integer*4 RVelClass ! [0/5]? Weight class for radial velocity character*20 r_RVel ! references for radial velocities * (format 5I4) *-------------------------------------------------------------------------------- * Buscombe Compilation (ref.64): *-------------------------------------------------------------------------------- integer*4 RVel_64 ! (km/s) ? Radial Velocity (2) integer*4 Dist_64 ! (pc) ? Distance (2) real*4 E_B_V__64 ! (mag) ? E(B-V) (2) real*4 iMag_64 ! (mag) ? Integrated apparent pg. magnitude (2) real*4 iAbsMag_64 ! (mag) ? Integrated absolute Bmag (2) integer*4 Diam_64 ! (pc) ? Absolute Diameter (2) real*4 log_age__64 ! ([yr]) ? log(age), years (2) real*4 log_Mass_64 ! ([solMass]) ? log (Mass), solar masses (2) real*4 inRef_64 ! [0/1]? '1' if cluster is in (ref.64) (2) *-------------------------------------------------------------------------------- * Basel data (ref.37): *-------------------------------------------------------------------------------- character*2 phSystem_37 ! Photometric system for following data (2) real*4 E_B_V__37 ! (mag) ? Colour excess (2) real*4 Av_37 ! (mag) ? Interstellar extinction (2) integer*4 Dist_37 ! (pc) ? Distance (2) integer*4 Sp_37 ! [0/11000]? Spectral type of hottest star integer*4 ColBlue_37 ! [0/10320]? Colour type of bluest star real*4 Diam_37 ! (pc) ? Cluster diameter (2) integer*4 ObsType_37 ! ? Type of observation (2) integer*4 inRef_37 ! [0/1]? '1' if cluster is in (ref.37) (2) *-------------------------------------------------------------------------------- * van den Bergh-Hagen data (ref.487): *-------------------------------------------------------------------------------- real*4 Diam_487 ! (arcmin) ? Apparent Diameter (2) character*2 Rich_487 ! Richness: vp, p, m, r, vr (2) integer*4 inRef_487 ! [0/1]? '1' if cluster is in (ref.487) (2) character*9 Cname ! [C0-9+-] Cluster "C" designation (4) *-------------------------------------------------------------------------------- * Various tracers: *-------------------------------------------------------------------------------- integer*4 tr_RA ! [0/1212]? tracer for RA, 1950.0 (8) integer*4 tr_GLON ! [0/1212]? tracer for Galactic longitude (8) integer*4 tr_ClNum ! [0/1212]? tracer for number in sequence (8) integer*4 tr_RA_1 ! [0/1212]? tracer for RA, 1950.0 (8) integer*4 tr_GLON_1 ! [0/1212]? tracer for Galactic longitude (8) integer*4 tr_ClNum_1 ! [0/1212]? tracer for number in sequence (8) *-------------------------------------------------------------------------------- * Data from joint project between Janes, Tilley and Lynga (1988): *-------------------------------------------------------------------------------- integer*4 Dist_jdl ! (pc) ? jdl distance modulus integer*4 w_Dist_jdl ! [0/5]? weight class real*4 turn_jdl ! (mag) [-0.4/0.7]? jdl turn-off point colour real*4 Age_jdl ! (Myr) [0/10000]? jdl derived age integer*4 w_Age_jdl ! [0/5]? weight class real*4 E_B_V__jdl ! (mag) [0/1.75]? jdl reddening character*1 n_E_B_V__jdl ! [v] 'v' if reddening varies integer*4 w_E_B_V__jdl ! [0/5]? weight class real*4 v_Fe_H__jdl ! (Sun) [0]? jdl abundance value character*18 refs_jdl ! 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 ! [1/527]+= Reference number character*126 Text ! 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 ! (deg) Galactic longitude real*4 GLAT_1 ! (deg) Galactic latitude integer*4 OCL_1 ! [1/1039]? OCL number (ref.19) character*9 Cname_1 ! Cluster name character*15 Name ! Cluster Sequence name and number integer*4 LundRec_1 ! [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,ClNum,RA2000h,RA2000m,DE2000_,DE2000d,DE2000m,RAh,RAm, + DE_,DEd,DEm,precRA,precDE,GLON,GLAT,LundRec,OCL,Diam,r_Diam, + Dist,r_Dist,log_Age,r_log_Age,v_Fe_H_,r__Fe_H_,E_B_V,r_E_B_V, + ClTyp,r_ClTyp,N_Ap,N_Am,N_mAp,N_mAm,N_WR,N_Of,N_Be,N_sh, + N_Be_sh,N_He_weak,N_He_rich,N_Cep,N_PN,N_C_,N_LPV,N_dSct, + N_bCep,N_RCrB,N_BaII,N_UGem,refN,POSS_type,POSS,RAPOSSh, + RAPOSSm,DEPOSSd,xPOSS,yPOSS,maxDiam,maxBr,TrConc,TrRange, + TrRich,TrNeb,inRef_128,POSS_128,Diam_50,Dist_50,maxBr_50, + Sp_50,totMag_50,E_B_V__50,inRef_50,totMag_422,i_B_V__422, + N_422,Num_265,Class_265,Ref_265,maxClass_265,Rich_265, + E_B_V__265,r_E_B_V__265,turn_265,r_turn_265,TrConc_320, + TrRange_320,TrRich_320,TrNeb_320,Members_320,Diam_320, + POSS_type_320,POSS_320,RAPOSS_320h,RAPOSS_320m,DEPOSS_320d, + xPOSS_320,yPOSS_320,hasClass_320,ESO_SERC,w_RVel,RVel, + RVelClass,r_RVel,RVel_64,Dist_64,E_B_V__64,iMag_64,iAbsMag_64, + Diam_64,log_age__64,log_Mass_64,inRef_64,phSystem_37, + E_B_V__37,Av_37,Dist_37,Sp_37,ColBlue_37,Diam_37,ObsType_37, + inRef_37,Diam_487,Rich_487,inRef_487,Cname,tr_RA,tr_GLON, + tr_ClNum,tr_RA_1,tr_GLON_1,tr_ClNum_1,Dist_jdl,w_Dist_jdl, + turn_jdl,Age_jdl,w_Age_jdl,E_B_V__jdl,n_E_B_V__jdl, + w_E_B_V__jdl,v_Fe_H__jdl,refs_jdl if(ar__(55:58) .EQ. '') OCL = iNULL__ if(ar__(59:64) .EQ. '') Diam = rNULL__ if(ar__(65:68) .EQ. '') r_Diam = iNULL__ if(ar__(69:72) .EQ. '') Dist = iNULL__ if(ar__(74:76) .EQ. '') r_Dist = iNULL__ if(ar__(77:81) .EQ. '') log_Age = rNULL__ if(ar__(82:84) .EQ. '') r_log_Age = iNULL__ if(ar__(85:89) .EQ. '') v_Fe_H_ = rNULL__ if(ar__(90:92) .EQ. '') r__Fe_H_ = iNULL__ if(ar__(93:98) .EQ. '') E_B_V = rNULL__ if(ar__(103:104) .EQ. '') r_ClTyp = iNULL__ if(ar__(167:170) .EQ. '') POSS = iNULL__ if(ar__(172:173) .EQ. '') RAPOSSh = iNULL__ if(ar__(175:176) .EQ. '') RAPOSSm = iNULL__ if(ar__(178:180) .EQ. '') DEPOSSd = iNULL__ if(ar__(181:183) .EQ. '') xPOSS = iNULL__ if(ar__(184:186) .EQ. '') yPOSS = iNULL__ if(ar__(187:190) .EQ. '') maxDiam = iNULL__ if(ar__(191:194) .EQ. '') maxBr = rNULL__ if(ar__(195:196) .EQ. '') TrConc = iNULL__ if(ar__(197:198) .EQ. '') TrRange = iNULL__ if(ar__(204:204) .EQ. '') inRef_128 = iNULL__ if(ar__(206:206) .EQ. '') POSS_128 = iNULL__ if(ar__(207:210) .EQ. '') Diam_50 = iNULL__ if(ar__(211:214) .EQ. '') Dist_50 = iNULL__ if(ar__(215:218) .EQ. '') maxBr_50 = rNULL__ if(ar__(220:224) .EQ. '') Sp_50 = iNULL__ if(ar__(225:228) .EQ. '') totMag_50 = rNULL__ if(ar__(229:232) .EQ. '') E_B_V__50 = rNULL__ if(ar__(233:234) .EQ. '') inRef_50 = iNULL__ if(ar__(235:238) .EQ. '') totMag_422 = rNULL__ if(ar__(239:242) .EQ. '') i_B_V__422 = rNULL__ if(ar__(243:246) .EQ. '') N_422 = iNULL__ if(ar__(247:250) .EQ. '') Num_265 = iNULL__ if(ar__(252:252) .EQ. '') Class_265 = iNULL__ if(ar__(254:256) .EQ. '') Ref_265 = iNULL__ if(ar__(258:258) .EQ. '') maxClass_265 = iNULL__ if(ar__(260:260) .EQ. '') Rich_265 = iNULL__ if(ar__(261:264) .EQ. '') E_B_V__265 = rNULL__ if(ar__(266:268) .EQ. '') r_E_B_V__265 = iNULL__ if(ar__(269:272) .EQ. '') turn_265 = rNULL__ if(ar__(274:276) .EQ. '') r_turn_265 = iNULL__ if(ar__(278:278) .EQ. '') TrConc_320 = iNULL__ if(ar__(280:280) .EQ. '') TrRange_320 = iNULL__ if(ar__(286:288) .EQ. '') Members_320 = iNULL__ if(ar__(290:294) .EQ. '') Diam_320 = rNULL__ if(ar__(297:300) .EQ. '') POSS_320 = iNULL__ if(ar__(302:303) .EQ. '') RAPOSS_320h = iNULL__ if(ar__(305:306) .EQ. '') RAPOSS_320m = iNULL__ if(ar__(308:310) .EQ. '') DEPOSS_320d = iNULL__ if(ar__(311:313) .EQ. '') xPOSS_320 = iNULL__ if(ar__(314:316) .EQ. '') yPOSS_320 = iNULL__ if(ar__(318:318) .EQ. '') hasClass_320 = iNULL__ if(ar__(320:320) .EQ. '') ESO_SERC = iNULL__ if(ar__(322:324) .EQ. '') w_RVel = rNULL__ if(ar__(325:328) .EQ. '') RVel = iNULL__ if(ar__(330:330) .EQ. '') RVelClass = iNULL__ if(ar__(351:354) .EQ. '') RVel_64 = iNULL__ if(ar__(355:358) .EQ. '') Dist_64 = iNULL__ if(ar__(359:362) .EQ. '') E_B_V__64 = rNULL__ if(ar__(363:366) .EQ. '') iMag_64 = rNULL__ if(ar__(367:370) .EQ. '') iAbsMag_64 = rNULL__ if(ar__(371:374) .EQ. '') Diam_64 = iNULL__ if(ar__(375:378) .EQ. '') log_age__64 = rNULL__ if(ar__(379:382) .EQ. '') log_Mass_64 = rNULL__ if(ar__(383:384) .EQ. '') inRef_64 = rNULL__ if(ar__(387:390) .EQ. '') E_B_V__37 = rNULL__ if(ar__(391:394) .EQ. '') Av_37 = rNULL__ if(ar__(395:398) .EQ. '') Dist_37 = iNULL__ if(ar__(399:404) .EQ. '') Sp_37 = iNULL__ if(ar__(405:410) .EQ. '') ColBlue_37 = iNULL__ if(ar__(411:416) .EQ. '') Diam_37 = rNULL__ if(ar__(417:418) .EQ. '') ObsType_37 = iNULL__ if(ar__(419:420) .EQ. '') inRef_37 = iNULL__ if(ar__(421:424) .EQ. '') Diam_487 = rNULL__ if(ar__(428:428) .EQ. '') inRef_487 = iNULL__ if(ar__(439:442) .EQ. '') tr_RA = iNULL__ if(ar__(443:446) .EQ. '') tr_GLON = iNULL__ if(ar__(447:450) .EQ. '') tr_ClNum = iNULL__ if(ar__(451:454) .EQ. '') tr_RA_1 = iNULL__ if(ar__(455:458) .EQ. '') tr_GLON_1 = iNULL__ if(ar__(459:462) .EQ. '') tr_ClNum_1 = iNULL__ if(ar__(464:468) .EQ. '') Dist_jdl = iNULL__ if(ar__(470:470) .EQ. '') w_Dist_jdl = iNULL__ if(ar__(471:475) .EQ. '') turn_jdl = rNULL__ if(ar__(476:481) .EQ. '') Age_jdl = rNULL__ if(ar__(483:483) .EQ. '') w_Age_jdl = iNULL__ if(ar__(484:488) .EQ. '') E_B_V__jdl = rNULL__ if(ar__(491:491) .EQ. '') w_E_B_V__jdl = iNULL__ if(ar__(492:496) .EQ. '') v_Fe_H__jdl = rNULL__ RAdeg = rNULL__ DEdeg = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh .GT. -180) RAdeg=RAh*15. if(RAm .GT. -180) RAdeg=RAdeg+RAm/4. if(DEd .GE. 0) DEdeg=DEd if(DEm .GE. 0) DEdeg=DEdeg+DEm/60. if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg c ..............Just test output........... write(6,1) + ClSeq,ClNum,RA2000h,RA2000m,DE2000_,DE2000d,DE2000m,RAh,RAm, + DE_,DEd,DEm,precRA,precDE,GLON,GLAT,LundRec,OCL,Diam,r_Diam, + Dist,r_Dist,log_Age,r_log_Age,v_Fe_H_,r__Fe_H_,E_B_V,r_E_B_V, + ClTyp,r_ClTyp,N_Ap,N_Am,N_mAp,N_mAm,N_WR,N_Of,N_Be,N_sh, + N_Be_sh,N_He_weak,N_He_rich,N_Cep,N_PN,N_C_,N_LPV,N_dSct, + N_bCep,N_RCrB,N_BaII,N_UGem,refN,POSS_type,POSS,RAPOSSh, + RAPOSSm,DEPOSSd,xPOSS,yPOSS,maxDiam,maxBr,TrConc,TrRange, + TrRich,TrNeb,inRef_128,POSS_128,Diam_50,Dist_50,maxBr_50, + Sp_50,totMag_50,E_B_V__50,inRef_50,totMag_422,i_B_V__422, + N_422,Num_265,Class_265,Ref_265,maxClass_265,Rich_265, + E_B_V__265,r_E_B_V__265,turn_265,r_turn_265,TrConc_320, + TrRange_320,TrRich_320,TrNeb_320,Members_320,Diam_320, + POSS_type_320,POSS_320,RAPOSS_320h,RAPOSS_320m,DEPOSS_320d, + xPOSS_320,yPOSS_320,hasClass_320,ESO_SERC,w_RVel,RVel, + RVelClass,r_RVel,RVel_64,Dist_64,E_B_V__64,iMag_64,iAbsMag_64, + Diam_64,log_age__64,log_Mass_64,inRef_64,phSystem_37, + E_B_V__37,Av_37,Dist_37,Sp_37,ColBlue_37,Diam_37,ObsType_37, + inRef_37,Diam_487,Rich_487,inRef_487,Cname,tr_RA,tr_GLON, + tr_ClNum,tr_RA_1,tr_GLON_1,tr_ClNum_1,Dist_jdl,w_Dist_jdl, + turn_jdl,Age_jdl,w_Age_jdl,E_B_V__jdl,n_E_B_V__jdl, + w_E_B_V__jdl,v_Fe_H__jdl,refs_jdl write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg 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,Text c ..............Just test output........... write(6,2)RefNum,Text 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,GLAT_1,OCL_1,Cname_1,Name,LundRec_1 if(ar__2(18:21) .EQ. '') OCL_1 = iNULL__ c ..............Just test output........... write(6,3)GLON_1,GLAT_1,OCL_1,Cname_1,Name,LundRec_1 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end