FORTRAN Generation
(/./ftp/cats/II/307/)

Conversion of standardized ReadMe file for file /./ftp/cats/II/307/ 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-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. II/307              WISE Preliminary Data Release         (Cutri+ 2011)
*================================================================================
*WISE Preliminary Data Release
*     Cutri R.M., et al.
*    <IPAC/Caltech (2011) >
*   =2011yCat.2307....0C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'wise.sam'	! Sample of the "slim" version of the WISE
                              preliminary release (among 257,310,278 sources)

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

C  J2000 position composed of: RAdeg DEdeg
      character*19  WISEP       ! WISE Preliminary Release Catalog name,
*                                 based on J2000 position (designation)
      real*8        RAdeg       ! (deg) Right ascension (J2000)
      real*8        DEdeg       ! (deg) Declination (J2000)
      real*4        eeMaj       ! (arcsec) Semi-major axis of the error ellipse (6)
      real*4        eeMin       ! (arcsec) Semi-minor axis of the error ellipse (6)
      integer*4     eePA        ! (deg) Position angle of the error ellipse (6)
      real*4        W1mag       ! (mag) ? W1 magnitude (3.35um), Vega scale
      real*4        W2mag       ! (mag) ? W2 magnitude (4.6um), Vega scale
      real*4        W3mag       ! (mag) ? W3 magnitude (11.6um), Vega scale
      real*4        W4mag       ! (mag) ? W4 magnitude (22.1um), Vega scale
      real*4        Jmag        ! (mag) ? 2MASS J magnitude (1.2um)
      real*4        Hmag        ! (mag) ? 2MASS H magnitude (1.6um)
      real*4        Kmag        ! (mag) ? 2MASS Ks magnitude (2.2um)
      real*4        e_W1mag     ! (mag) ? Mean error on W1 magnitude (3.35um)
      real*4        e_W2mag     ! (mag) ? Mean error on W2 magnitude (4.6um)
      real*4        e_W3mag     ! (mag) ? Mean error on W3 magnitude (11.6um)
      real*4        e_W4mag     ! (mag) ? Mean error on W4 magnitude (22.1um)
      real*4        e_Jmag      ! (mag) ? Mean error on J magnitude
      real*4        e_Hmag      ! (mag) ? Mean error on H magnitude
      real*4        e_Kmag      ! (mag) ? Mean error on Ks magnitude
      character*19  cntr        ! Unique identification number for the source
      real*4        snr1        ! ? Signal to noise ratio for W1 filter
      real*4        chi2W1      ! ? Reduced {chi}^2^ of the W1 profile fit
      real*4        snr2        ! ? Signal to noise ratio for W2 filter
      real*4        chi2W2      ! ? Reduced {chi}^2^ of the W2 profile fit
      real*4        snr3        ! ? Signal to noise ratio for W3 filter
      real*4        chi2W3      ! ? Reduced {chi}^2^ of the W3 profile fit
      real*4        snr4        ! ? Signal to noise ratio for W4 filter
      real*4        chi2W4      ! ? Reduced {chi}^2^ of the W4 profile fit
      integer*4     nb          ! Number of PSF components used simultaneously
*                                 in the profile-fitting
      integer*4     na          ! [0/1] Active deblending flag (1 if a detection
*                                 was split into multiple sources)
      real*4        sat1        ! [0/1] Saturated pixel fraction in W1 (w1sat)
      real*4        sat2        ! [0/1] Saturated pixel fraction in W2 (w2sat)
      real*4        sat3        ! [0/1] Saturated pixel fraction in W3 (w3sat)
      real*4        sat4        ! [0/1] Saturated pixel fraction in W4 (w4sat)
      character*4   ccf         ! [DdPpHhOo0] Contamination and confusion flag,
*                                       one per band (cc_flags) (1)
      integer*4     ex          ! [0/1] Extended source flag (ext_flg)
*                                       (0 = consistent with point source)
      character*4   var         ! [0-9] Variability flag, one per band (2)
      integer*4     nW1         ! Frame detection count in W1 (w1nm) (4)
      integer*4     mW1         ! Integer frame coverage in W1 (w1m) (5)
      integer*4     nW2         ! Frame detection count in W2 (w2nm) (4)
      integer*4     mW2         ! Integer frame coverage in W2 (w2m) (5)
      integer*4     nW3         ! Frame detection count in W3 (w3nm) (4)
      integer*4     mW3         ! Integer frame coverage in W3 (w3m) (5)
      integer*4     nW4         ! Frame detection count in W4 (w4nm) (4)
      integer*4     mW4         ! Integer frame coverage in W4 (w4m) (5)
      integer*8     v2Mkey      ! ? 2MASS PSC association (not identification!)
      real*4        d2M         ! (arcsec) [0,3]? Distance separating the positions of the
*                                 WISE source and associated 2MASS PSC source
*Note (1): One character per band (W1/W2/W3/W4) that indicates that the
*     photometry and/or position measurements of a source may be
*     contaminated or biased due to proximity to an image artifact.
* D,d =  Diffraction spike. Source may be a spurious detection of or
*     contaminated by a diffraction spike from a nearby bright star on
*     the same image
* P,p = Persistence. Source may be a spurious detection of or contaminated
*       by a latent image left by a bright star
* H,h = Halo. Source photometry may be a spurious detection of or
*       contaminated by the scattered light halo surrounding a nearby
*       bright source
* O,o = (letter "o") Optical ghost. Source may be a spurious detection
*     of or contaminated by an optical ghost image caused by a nearby bright
*     source
* 0  = (number zero) Source is unaffected by known artifacts.
*Note (2): The variability flag is a 4-character string, one character per
*     band, containing a value related to the probability that the source
*     flux measured on the individual WISE exposures is variable.
*  * value "0" indicates insufficient or inadequate data to make a
*    determination  (<11 exposures)
*  * values 1 thru 9 indicate increasing probabilities of variation;
*    - values 1-4 are most likely not variable,
*    - values 5-7 are likely variables (but susceptible of false-positive
*      variability)
*    - values >7 have the highest probability of being true variables
*Note (4): number of individual exposures on which this source was detected
*     with SNR>3 in the profile-fit measurement. This number can be zero for
*     sources that are well-detected on the coadded Atlas Image, but too
*     faint for detection on the single exposures.
*Note (5): number of individual exposures on which a profile-fit measurement
*     of this source was possible.
*Note (6): the parameters of the error ellipse are computed from the
*     1-{sigma} in RA and Dec (sigra, sigdec) and the co-{sigma} sigradec
*     with the formulae:
*     {Delta} = (sigra^2^-sigdec^2^)^2^ + 4*sigradec^2^
*     eeMaj^2^ = 0.5*(sigra^2^+sigdec^2^+sqrt({Delta}))
*     eeMin^2^ = 0.5*(sigra^2^+sigdec^2^-sqrt({Delta}))
*     tan(eePA) = (eeMaj^2^-sigdec^2^)/(sigradec*|sigradec|)
*               = (sigradec*|sigradec|)/(eeMaj^2^-sigra^2^)
*
*     Conversely, the sigra/sigdec are given by:
*     sigra^2^  = eeMaj^2^sin^2^(eePA) + eeMin^2^cos^2^(eePA)
*     sigdec^2^ = eeMaj^2^cos^2^(eePA) + eeMin^2^sin^2^(eePA)

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

C  Loading file 'wise.sam'	! Sample of the "slim" version of the WISE
*                              preliminary release (among 257,310,278 sources)

C  Format for file interpretation

    1 format(
     +  A19,1X,F9.5,1X,F9.5,1X,F4.2,1X,F4.2,1X,I3,2X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,2X,F5.3,1X,F5.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,A19,2X,F6.1,1X,E9.3,1X,
     +  F6.1,1X,E9.3,1X,F6.1,1X,E9.3,1X,F6.1,1X,E9.3,2X,I1,1X,I1,2X,
     +  F5.3,1X,F5.3,1X,F5.3,1X,F5.3,2X,A4,1X,I1,1X,A4,7X,I4,1X,I4,1X,
     +  I4,1X,I4,1X,I4,1X,I4,1X,I4,1X,I4,1X,I10,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'wise.sam')
      write(6,*) '....Loading file: wise.sam'
      do i__=1,150
        read(1,'(A342)')ar__
        read(ar__,1)
     +  WISEP,RAdeg,DEdeg,eeMaj,eeMin,eePA,W1mag,W2mag,W3mag,W4mag,
     +  Jmag,Hmag,Kmag,e_W1mag,e_W2mag,e_W3mag,e_W4mag,e_Jmag,e_Hmag,
     +  e_Kmag,cntr,snr1,chi2W1,snr2,chi2W2,snr3,chi2W3,snr4,chi2W4,
     +  nb,na,sat1,sat2,sat3,sat4,ccf,ex,var,nW1,mW1,nW2,mW2,nW3,mW3,
     +  nW4,mW4,v2Mkey,d2M
        if(ar__(56:61) .EQ. '') W1mag = rNULL__
        if(ar__(63:68) .EQ. '') W2mag = rNULL__
        if(ar__(70:75) .EQ. '') W3mag = rNULL__
        if(ar__(77:82) .EQ. '') W4mag = rNULL__
        if(ar__(84:89) .EQ. '') Jmag = rNULL__
        if(ar__(91:96) .EQ. '') Hmag = rNULL__
        if(ar__(98:103) .EQ. '') Kmag = rNULL__
        if(ar__(106:110) .EQ. '') e_W1mag = rNULL__
        if(ar__(112:116) .EQ. '') e_W2mag = rNULL__
        if(ar__(118:122) .EQ. '') e_W3mag = rNULL__
        if(ar__(124:128) .EQ. '') e_W4mag = rNULL__
        if(ar__(130:134) .EQ. '') e_Jmag = rNULL__
        if(ar__(136:140) .EQ. '') e_Hmag = rNULL__
        if(ar__(142:146) .EQ. '') e_Kmag = rNULL__
        if(ar__(169:174) .EQ. '') snr1 = rNULL__
        if(ar__(176:184) .EQ. '') chi2W1 = rNULL__
        if(ar__(186:191) .EQ. '') snr2 = rNULL__
        if(ar__(193:201) .EQ. '') chi2W2 = rNULL__
        if(ar__(203:208) .EQ. '') snr3 = rNULL__
        if(ar__(210:218) .EQ. '') chi2W3 = rNULL__
        if(ar__(220:225) .EQ. '') snr4 = rNULL__
        if(ar__(227:235) .EQ. '') chi2W4 = rNULL__
        if(ar__(326:335) .EQ. '') v2Mkey = iNULL__
        if(ar__(337:342) .EQ. '') d2M = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  WISEP,RAdeg,DEdeg,eeMaj,eeMin,eePA,W1mag,W2mag,W3mag,W4mag,
     +  Jmag,Hmag,Kmag,e_W1mag,e_W2mag,e_W3mag,e_W4mag,e_Jmag,e_Hmag,
     +  e_Kmag,cntr,snr1,chi2W1,snr2,chi2W2,snr3,chi2W3,snr4,chi2W4,
     +  nb,na,sat1,sat2,sat3,sat4,ccf,ex,var,nW1,mW1,nW2,mW2,nW3,mW3,
     +  nW4,mW4,v2Mkey,d2M
c    .......End.of.Just test output...........
      end do
      close(1)

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