Conversion of standardized ReadMe file for
file /./ftp/cats/VII/299 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 2026-Oct-13
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/299 HELP DR1 photometry of extragalactic field sources (Shirley+, 2021)
*================================================================================
*HELP the Herschel Extragalactic Legacy Project.
* Shirley R., Duncan K., Campos-Varillas M.C., Hurley P.D., Malek K.,
* Roehlly Y., Smith M.W.L., Aussel H., Bakx T., Buat V., Burgarella D.,
* Christopher N., Duivenvoorden S., Eales S., Efstathiou A.,
* Gonzalez Solares E.A., Griffin M., Jarvis M., Faro B.L., Marchetti L.,
* McCheyne I., Papadopoulos A., Penner K., Pons E., Prescott M., Rigby E.,
* Rottgering H., Saxena A., Scudder J., Vaccari M., Wang L. and Oliver S.J.
* <Mon. Not. R. Astron. Soc. 507, 129-155 (2021)>
* =2021MNRAS.507..129S (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'help.sam' ! Herschel Extragalactic Legacy Project HELP
photometric catalogue of extragalactic field
sources
integer*4 nr__
parameter (nr__=1000) ! Number of records
character*1443 ar__ ! Full-size record
character*28 HELP (nr__) ! HELP identifier based on J2000 position
* as HELP_JHHMMSS.sss+DDMMSS.sss (help_id)
real*8 RAJ2000 (nr__) ! (deg) Right Ascension (J2000) (ra)
real*8 DEJ2000 (nr__) ! (deg) Declination (J2000) (dec)
character*18 Field (nr__) ! Name of the field in HELP survey as
* detailed in appendix A in our study
* (field)
integer*4 HEALPIX (nr__) ! HEALPix index of the source position at
* order 13 using the nested scheme (hp_idx)
real*8 E_B_V (nr__) ! (mag) Galactic extinction (Schlegel et al.
* 1998ApJ...500..525S) at the object
* position (ebv)
real*8 zphEAZY (nr__) ! ? Photometric redshift obtained with EAZY
* (redshift)
real*4 zsp (nr__) ! ? Spectroscopic redshift collated from
* the literature (zspec)
integer*4 q_zsp (nr__) ! [0/5]? Spectroscopic redshift quality
* flag FIXME code meaning (zspec_qual)
integer*4 f_zsp (nr__) ! [0/1] Boolean flag set to true when the
* association of a source to its redshift
* may be problematic possible
* mis-cross-association
* (zspec_association_flag)
real*8 stellar (nr__) ! ? Index evaluation is the source is
* point-source as 0 not point source to 1
* point source and generally name
* stellarity index (stellarity)
character*20 r_stellar (nr__) ! ? Name of Survey from which stellarity is
* calculated (stellarity_origin)
real*8 Fu (nr__) ! (uJy) ? Total flux density in the lowest error
* u band measurement (f_best_u)
real*8 e_Fu (nr__) ! (uJy) ? Uncertainty on f_best_u (ferr_best_u)
real*8 umag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* u band measurement (m_best_u)
real*4 e_umag (nr__) ! (mag) ? Uncertainty on m_best_u (merr_best_u)
integer*4 f_Fu (nr__) ! [0] Flag set to true for sources for
* spurious measurements in u (flag_best_u)
real*4 Fg (nr__) ! (uJy) ? Total flux density in the lowest error
* g band measurement (f_best_g)
real*4 e_Fg (nr__) ! (uJy) ? Uncertainty on f_best_g (ferr_best_g)
real*8 gmag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* g band measurement (m_best_g)
real*4 e_gmag (nr__) ! (mag) []? Uncertainty on m_best_g (merr_best_g)
integer*4 f_Fg (nr__) ! [0/1] Flag set to true for sources for
* spurious measurements in g (flag_best_g)
real*4 Fr (nr__) ! (uJy) ? Total flux density in the lowest error
* r band measurement (f_best_r)
real*4 e_Fr (nr__) ! (uJy) ? Uncertainty on f_best_r (ferr_best_r)
real*8 rmag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* g band measurement (m_best_r)
real*4 e_rmag (nr__) ! (mag) []? Uncertainty on m_best_r (merr_best_r)
integer*4 f_Fr (nr__) ! [0/1] Flag set to true for sources for
* spurious measurements in r (flag_best_r)
real*4 Fi (nr__) ! (uJy) ? Total flux density in the lowest error
* i band measurement (f_best_i)
real*4 e_Fi (nr__) ! (uJy) ? Uncertainty on f_best_i (ferr_best_i)
real*8 imag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* i band measurement (m_best_i)
real*4 e_imag (nr__) ! (mag) ? Uncertainty on m_best_i (merr_best_i)
integer*4 f_Fi (nr__) ! [0/1] Flag set to true for sources for
* spurious measurements in i (flag_best_i)
real*4 Fz (nr__) ! (uJy) ? Total flux density in the lowest error
* z band measurement (f_best_z)
real*4 e_Fz (nr__) ! (uJy) ? Uncertainty on f_best_z (ferr_best_z)
real*8 zmag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* z band measurement (m_best_z)
real*4 e_zmag (nr__) ! (mag) []? Uncertainty on m_best_z (merr_best_z)
integer*4 f_Fz (nr__) ! [0/1] Flag set to true for sources for
* spurious measurements in z (flag_best_z)
real*4 Fy (nr__) ! (uJy) ? Total flux density in the lowest error
* y band measurement (f_best_y)
real*4 e_Fy (nr__) ! (uJy) ? Uncertainty on f_best_y (ferr_best_y)
real*8 ymag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* y band measurement (m_best_y)
real*4 e_ymag (nr__) ! (mag) ? Uncertainty on m_best_y (merr_best_y)
integer*4 f_Fy (nr__) ! [0/1]? Flag set to true for sources for
* spurious measurements in y (flag_best_y)
real*4 Fj (nr__) ! (uJy) ? Total flux density in the lowest error
* j band measurement (f_best_j)
real*4 e_Fj (nr__) ! (uJy) ? Uncertainty on f_best_j (ferr_best_j)
real*8 jmag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* j band measurement (m_best_j)
real*4 e_jmag (nr__) ! (mag) ? Uncertainty on m_best_j (merr_best_j)
integer*4 f_Fj (nr__) ! [0/1]? Flag set to true for sources for
* spurious measurements in j (flag_best_j)
real*4 Fh (nr__) ! (uJy) ? Total flux density in the lowest error
* h band measurement (f_best_h)
real*4 e_Fh (nr__) ! (uJy) ? Uncertainty on f_best_h (ferr_best_h)
real*8 hmag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* h band measurement (m_best_h)
real*4 e_hmag (nr__) ! (mag) ? Uncertainty on m_best_h (merr_best_h)
integer*4 f_Fh (nr__) ! [0] Flag set to true for sources for
* spurious measurements in h (flag_best_h)
real*4 Fk (nr__) ! (uJy) ? Total flux density in the lowest error
* k band measurement (f_best_k)
real*4 e_Fk (nr__) ! (uJy) ? Uncertainty on f_best_k (ferr_best_k)
real*8 kmag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* k band measurement (m_best_k)
real*4 e_kmag (nr__) ! (mag) ? Uncertainty on m_best_k (merr_best_k)
integer*4 f_Fk (nr__) ! [0/1]? Flag set to true for sources for
* spurious measurements in k (flag_best_k)
real*4 Fks (nr__) ! (uJy) ? Total flux density in the lowest error
* ks band measurement (f_best_ks)
real*4 e_Fks (nr__) ! (uJy) ? Uncertainty on f_best_ks (ferr_best_ks)
real*8 ksmag (nr__) ! (mag) ? Total AB magnitude of the lowest error
* ks band measurement (m_best_ks)
real*4 e_ksmag (nr__) ! (mag) ? Uncertainty on m_best_ks (merr_best_ks)
integer*4 f_Fks (nr__) ! [0] Flag set to true for sources for
* spurious measurements in ks (flag_best_ks)
real*4 Fi1 (nr__) ! (uJy) ? Total flux density in the 3.6um band of
* The Infrared Array Camera IRAC on The
* Spitzer Space Telescope (f_irac_i1)
real*4 e_Fi1 (nr__) ! (uJy) ? Uncertainty on f_irac_i1 (ferr_irac_i1)
real*8 i1mag (nr__) ! (mag) ? Total AB magnitude in the i1 band
* 3.6{mu}m of the Infrared Array Camera
* IRAC on the Spitzer Space Telescope
* (m_irac_i1)
real*4 e_i1mag (nr__) ! (mag) ? Uncertainty on m_irac_i1 (merr_irac_i1)
integer*4 f_Fi1 (nr__) ! [0/1]? Flag set to true for sources for
* which niether the irac_i1 aperture nor
* total flux should be used for SED fitting
* (flag_irac_i1)
real*4 Fi2 (nr__) ! (uJy) ? Total flux density in the 4.5um band of
* The Infrared Array Camera (IRAC) on The
* Spitzer Space Telescope (f_irac_i2)
real*4 e_Fi2 (nr__) ! (uJy) ? Uncertainty on f_irac_i2 (ferr_irac_i2)
real*8 i2mag (nr__) ! (mag) ? Total AB magnitude in i2 band 4.5{mu}m
* of The Infrared Array Camera IRAC on The
* Spitzer Space Telescope (m_irac_i2)
real*4 e_i2mag (nr__) ! (mag) ? Uncertainty on m_irac_i2 (merr_irac_i2)
integer*4 f_Fi2 (nr__) ! [0]? Flag set to true for sources for
* which niether the irac_i2 aperture nor
* total flux should be used for SED fitting
* (flag_irac_i2)
real*4 Fi3 (nr__) ! (uJy) ? Total flux density in the 5.6um band of
* The Infrared Array Camera (IRAC) on The
* Spitzer Space Telescope (f_irac_i3)
real*4 e_Fi3 (nr__) ! (uJy) ? Uncertainty on f_irac_i3 (ferr_irac_i3)
real*8 i3mag (nr__) ! (mag) ? Total AB magnitude in the 5.6 {mu}m
* band of The Infrared Array Camera IRAC on
* The Spitzer Space Telescope (m_irac_i3)
real*4 e_i3mag (nr__) ! (mag) ? Uncertainty on m_irac_i3 (merr_irac_i3)
integer*4 f_Fi3 (nr__) ! [0]? Flag set to true for sources for
* which niether the irac_i3 aperture nor
* total flux should be used for SED fitting
* (flag_irac_i3)
real*4 Fi4 (nr__) ! (uJy) ? Total flux density in the 8.0um band of
* The Infrared Array Camera (IRAC) on The
* Spitzer Space Telescope (f_irac_i4)
real*4 e_Fi4 (nr__) ! (uJy) ? Uncertainty on f_irac_i4 (ferr_irac_i4)
real*8 i4mag (nr__) ! (mag) ? Total AB magnitude in the 8.0 um band
* of The Infrared Array Camera (IRAC) on
* The Spitzer Space Telescope (m_irac_i4)
real*4 e_i4mag (nr__) ! (mag) ? Uncertainty on m_irac_i4 (merr_irac_i4)
integer*4 f_Fi4 (nr__) ! [0]? Flag set to true for sources for
* which niether the irac_i4 aperture nor
* total flux should be used for SED fitting
* (flag_irac_i4)
real*4 Fm24 (nr__) ! (uJy) ? Total flux density in the 15-30um band
* of The Multiband Imaging Photometer MIPS
* on The Spitzer Space Telescope (f_mips_24)
real*4 e_Fm24 (nr__) ! (uJy) ? Uncertainty on f_mips_24 (ferr_mips_24)
integer*4 f_Fm24 (nr__) ! [0/1]? Flag set to true for sources for
* which niether the mips_24 aperture nor
* total flux should be used for SED fitting
* (flag_mips_24)
real*4 Fpacsg (nr__) ! (uJy) ? Total flux density in the 100um green
* band of Photoconductor Array Camera and
* Spectrometer PACS on The Herschel Space
* Telescope (f_pacs_green)
real*4 e_Fpacsg (nr__) ! (uJy) ? Uncertainty on f_pacs_green
* (ferr_pacs_green)
integer*4 f_Fpacsg (nr__) ! [0/1]? Flag set to true for sources for
* which niether the pacs_green aperture nor
* total flux should be used for SED fitting
* (see documentation) (flag_pacs_green)
real*4 Fpacsr (nr__) ! (uJy) ? Total flux density in the 160um red
* band of Photoconductor Array Camera and
* Spectrometer PACS on The Herschel Space
* Telescope (f_pacs_red)
real*4 e_Fpacsr (nr__) ! (uJy) ? Uncertainty on f_pacs_red
* (ferr_pacs_red)
integer*4 f_Fpacsr (nr__) ! [0/1]? Flag set to true for sources for
* which niether the pacs_red aperture nor
* total flux should be used for SED fitting
* (flag_pacs_red)
real*4 F250sp (nr__) ! (uJy) ? Total flux density in the 750-1500GHz
* band of The Spectral and Photometric
* Imaging Reciever SPIRE on The Herschel
* Space Telescope (f_spire_250)
real*4 e_F250sp (nr__) ! (uJy) ? Uncertainty on f_spire_250
* (ferr_spire_250)
integer*4 f_F250sp (nr__) ! [0/1]? Flag set to true for sources for
* which niether the spire_250 aperture nor
* total flux should be used for SED fitting
* (flag_spire_250)
real*4 F350sp (nr__) ! (uJy) ? Total flux density in the 750-1500GHz
* band of The Spectral and Photometric
* Imaging Reciever (SPIRE) on The Herschel
* Space Telescope (f_spire_350)
real*4 e_F350sp (nr__) ! (uJy) ? Uncertainty on f_spire_350
* (ferr_spire_350)
integer*4 f_F350sp (nr__) ! [0/1]? Flag set to true for sources for
* which niether the spire_350 aperture nor
* total flux should be used for SED fitting
* (flag_spire_350)
real*4 F500sp (nr__) ! (uJy) ? Total flux density in the 400-750GHz
* band of The Spectral and Photometric
* Imaging Reciever SPIRE on The Herschel
* Space Telescope (f_spire_500)
real*4 e_F500sp (nr__) ! (uJy) ? Uncertainty on f_spire_500
* (ferr_spire_500)
integer*4 f_F500sp (nr__) ! [0/1]? Flag set to true for sources for
* which niether the spire_500 aperture nor
* total flux should be used for SED fitting
* (see documentation) (flag_spire_500)
real*4 SFR (nr__) ! (Msun) ? Star Formation Rate estimated by SED
* fitting with CIGALE (cigale_sfr)
real*4 e_SFR (nr__) ! (Msun) ? Uncertainty on CIGALE_SFR
* (cigale_sfr_err)
real*4 M_ (nr__) ! (Msun) ? Stellar Mass estimated by SED fitting
* with CIGALE (cigale_mstar)
real*4 e_M_ (nr__) ! (Msun) ? Uncertainty on CIGALE_Mstar
* (cigale_mstar_err)
real*4 Ldust (nr__) ! (W) ? Dust luminosity estimated by SED
* fitting with CIGALE using both optical
* and infra-red information
* (cigale_dustlumin)
real*4 e_Ldust (nr__) ! (W) ? Uncertainty on CIGALE_DustLumin
* (cigale_dustlumin_err)
real*4 LdustIR (nr__) ! (W) ? Dust luminosity estimated by SED
* fitting with CIGALE using only the
* infra-red information
* (cigale_dustlumin_ironly)
real*4 e_LdustIR (nr__) ! (W) ? Uncertainty on CIGALE_DustLumin_IRonly
* (cigale_dustlumin_ironly_err)
real*4 Chi2Opt (nr__) ! ? Reduced optical chi squared of best
* CIGALE fit to all data (cigale_chi2_opt)
real*4 Chi2IRbest (nr__) ! ? Reduced infrared chi squared of best
* CIGALE fit to all data (cigale_chi2_ir)
real*4 Chi2best (nr__) ! ? Reduced chi squared of best CIGALE fit
* (cigale_chi2_red)
real*4 Chi2IR (nr__) ! ? Reduced chi squared of best CIGALE fit
* to infrared data (cigale_chi2_ironly)
integer*4 fGaia (nr__) ! [0/3]? Flag indicating the increasing
* probability of the source being a Gaia
* object (flag_gaia)
integer*4 fClean (nr__) ! [0/1]? Boolean flag denoting that the
* source was associated to really near
* object that were removed during the
* cleaning procedure (flag_cleaned)
integer*4 fMerg (nr__) ! [0/1]? Boolean flag denoting a possible
* mis-association during the cross-matching
* of the various catalogues (flag_merged)
integer*4 fOptNIRo (nr__) ! [0/7]? Flag indicating in which
* wavelength regimes the source was
* observed (flag_optnir_obs)
integer*4 fOptNIRd (nr__) ! [0/7]? Flag indicating in which
* wavelength regimes the source was
* detected (flag_optnir_det)
C=============================================================================
C Loading file 'help.sam' ! Herschel Extragalactic Legacy Project HELP
* photometric catalogue of extragalactic field
* sources
C Format for file interpretation
1 format(
+ A28,1X,F16.12,1X,F16.12,1X,A18,1X,I9,1X,F11.9,1X,F11.9,1X,
+ E15.9,1X,I1,1X,I1,1X,F11.9,1X,A20,1X,F18.9,1X,F15.9,1X,F12.9,
+ 1X,E15.9,1X,I1,1X,E16.9,1X,E16.9,1X,F12.9,1X,E15.9,1X,I1,1X,
+ E16.9,1X,E15.9,1X,F12.9,1X,E15.9,1X,I1,1X,E16.9,1X,E16.9,1X,
+ F13.9,1X,E15.9,1X,I1,1X,E16.9,1X,E16.9,1X,F12.9,1X,E15.9,1X,
+ I1,1X,E16.9,1X,E14.9,1X,F12.9,1X,E15.9,1X,I1,1X,E15.9,1X,
+ E15.9,1X,F12.9,1X,E15.9,1X,I1,1X,E15.9,1X,E14.9,1X,F12.9,1X,
+ E15.9,1X,I1,1X,E16.9,1X,E15.9,1X,F13.9,1X,E15.9,1X,I1,1X,
+ E15.9,1X,E15.9,1X,F12.9,1X,E15.9,1X,I1,1X,E15.9,1X,E14.9,1X,
+ F12.9,1X,E15.9,1X,I1,1X,E15.9,1X,E14.9,1X,F12.9,1X,E15.9,1X,
+ I1,1X,E15.9,1X,E15.9,1X,F12.9,1X,E15.9,1X,I1,1X,E15.9,1X,
+ E15.9,1X,F12.9,1X,E15.9,1X,I1,1X,E15.9,1X,E14.9,1X,I1,1X,
+ E13.9,1X,E13.9,1X,I1,1X,E13.9,1X,E13.9,1X,I1,1X,E15.9,1X,
+ E14.9,1X,I1,1X,E15.9,1X,E14.9,1X,I1,1X,E15.9,1X,E14.9,1X,I1,
+ 1X,E14.9,1X,E14.9,1X,E14.9,1X,E14.9,1X,E14.9,1X,E14.9,1X,
+ E14.9,1X,E14.9,1X,E14.9,1X,E14.9,1X,E14.9,1X,E15.9,1X,I1,1X,
+ I1,1X,I1,1X,I1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'help.sam')
write(6,*) '....Loading file: help.sam'
do i__=1,1000
read(1,'(A1443)')ar__
read(ar__,1)
+ HELP(i__),RAJ2000(i__),DEJ2000(i__),Field(i__),HEALPIX(i__),
+ E_B_V(i__),zphEAZY(i__),zsp(i__),q_zsp(i__),f_zsp(i__),
+ stellar(i__),r_stellar(i__),Fu(i__),e_Fu(i__),umag(i__),
+ e_umag(i__),f_Fu(i__),Fg(i__),e_Fg(i__),gmag(i__),e_gmag(i__),
+ f_Fg(i__),Fr(i__),e_Fr(i__),rmag(i__),e_rmag(i__),f_Fr(i__),
+ Fi(i__),e_Fi(i__),imag(i__),e_imag(i__),f_Fi(i__),Fz(i__),
+ e_Fz(i__),zmag(i__),e_zmag(i__),f_Fz(i__),Fy(i__),e_Fy(i__),
+ ymag(i__),e_ymag(i__),f_Fy(i__),Fj(i__),e_Fj(i__),jmag(i__),
+ e_jmag(i__),f_Fj(i__),Fh(i__),e_Fh(i__),hmag(i__),e_hmag(i__),
+ f_Fh(i__),Fk(i__),e_Fk(i__),kmag(i__),e_kmag(i__),f_Fk(i__),
+ Fks(i__),e_Fks(i__),ksmag(i__),e_ksmag(i__),f_Fks(i__),
+ Fi1(i__),e_Fi1(i__),i1mag(i__),e_i1mag(i__),f_Fi1(i__),
+ Fi2(i__),e_Fi2(i__),i2mag(i__),e_i2mag(i__),f_Fi2(i__),
+ Fi3(i__),e_Fi3(i__),i3mag(i__),e_i3mag(i__),f_Fi3(i__),
+ Fi4(i__),e_Fi4(i__),i4mag(i__),e_i4mag(i__),f_Fi4(i__),
+ Fm24(i__),e_Fm24(i__),f_Fm24(i__),Fpacsg(i__),e_Fpacsg(i__),
+ f_Fpacsg(i__),Fpacsr(i__),e_Fpacsr(i__),f_Fpacsr(i__),
+ F250sp(i__),e_F250sp(i__),f_F250sp(i__),F350sp(i__),
+ e_F350sp(i__),f_F350sp(i__),F500sp(i__),e_F500sp(i__),
+ f_F500sp(i__),SFR(i__),e_SFR(i__),M_(i__),e_M_(i__),
+ Ldust(i__),e_Ldust(i__),LdustIR(i__),e_LdustIR(i__),
+ Chi2Opt(i__),Chi2IRbest(i__),Chi2best(i__),Chi2IR(i__),
+ fGaia(i__),fClean(i__),fMerg(i__),fOptNIRo(i__),fOptNIRd(i__)
if(ar__(105:115) .EQ. '') zphEAZY(i__) = rNULL__
if(ar__(117:131) .EQ. '') zsp(i__) = rNULL__
if(ar__(133:133) .EQ. '') q_zsp(i__) = iNULL__
if(ar__(137:147) .EQ. '') stellar(i__) = rNULL__
if(ar__(170:187) .EQ. '') Fu(i__) = rNULL__
if(ar__(189:203) .EQ. '') e_Fu(i__) = rNULL__
if(ar__(205:216) .EQ. '') umag(i__) = rNULL__
if(ar__(218:232) .EQ. '') e_umag(i__) = rNULL__
if(ar__(236:251) .EQ. '') Fg(i__) = rNULL__
if(ar__(253:268) .EQ. '') e_Fg(i__) = rNULL__
if(ar__(270:281) .EQ. '') gmag(i__) = rNULL__
if(ar__(283:297) .EQ. '') e_gmag(i__) = rNULL__
if(ar__(301:316) .EQ. '') Fr(i__) = rNULL__
if(ar__(318:332) .EQ. '') e_Fr(i__) = rNULL__
if(ar__(334:345) .EQ. '') rmag(i__) = rNULL__
if(ar__(347:361) .EQ. '') e_rmag(i__) = rNULL__
if(ar__(365:380) .EQ. '') Fi(i__) = rNULL__
if(ar__(382:397) .EQ. '') e_Fi(i__) = rNULL__
if(ar__(399:411) .EQ. '') imag(i__) = rNULL__
if(ar__(413:427) .EQ. '') e_imag(i__) = rNULL__
if(ar__(431:446) .EQ. '') Fz(i__) = rNULL__
if(ar__(448:463) .EQ. '') e_Fz(i__) = rNULL__
if(ar__(465:476) .EQ. '') zmag(i__) = rNULL__
if(ar__(478:492) .EQ. '') e_zmag(i__) = rNULL__
if(ar__(496:511) .EQ. '') Fy(i__) = rNULL__
if(ar__(513:526) .EQ. '') e_Fy(i__) = rNULL__
if(ar__(528:539) .EQ. '') ymag(i__) = rNULL__
if(ar__(541:555) .EQ. '') e_ymag(i__) = rNULL__
if(ar__(557:557) .EQ. '') f_Fy(i__) = iNULL__
if(ar__(559:573) .EQ. '') Fj(i__) = rNULL__
if(ar__(575:589) .EQ. '') e_Fj(i__) = rNULL__
if(ar__(591:602) .EQ. '') jmag(i__) = rNULL__
if(ar__(604:618) .EQ. '') e_jmag(i__) = rNULL__
if(ar__(620:620) .EQ. '') f_Fj(i__) = iNULL__
if(ar__(622:636) .EQ. '') Fh(i__) = rNULL__
if(ar__(638:651) .EQ. '') e_Fh(i__) = rNULL__
if(ar__(653:664) .EQ. '') hmag(i__) = rNULL__
if(ar__(666:680) .EQ. '') e_hmag(i__) = rNULL__
if(ar__(684:699) .EQ. '') Fk(i__) = rNULL__
if(ar__(701:715) .EQ. '') e_Fk(i__) = rNULL__
if(ar__(717:729) .EQ. '') kmag(i__) = rNULL__
if(ar__(731:745) .EQ. '') e_kmag(i__) = rNULL__
if(ar__(747:747) .EQ. '') f_Fk(i__) = iNULL__
if(ar__(749:763) .EQ. '') Fks(i__) = rNULL__
if(ar__(765:779) .EQ. '') e_Fks(i__) = rNULL__
if(ar__(781:792) .EQ. '') ksmag(i__) = rNULL__
if(ar__(794:808) .EQ. '') e_ksmag(i__) = rNULL__
if(ar__(812:826) .EQ. '') Fi1(i__) = rNULL__
if(ar__(828:841) .EQ. '') e_Fi1(i__) = rNULL__
if(ar__(843:854) .EQ. '') i1mag(i__) = rNULL__
if(ar__(856:870) .EQ. '') e_i1mag(i__) = rNULL__
if(ar__(872:872) .EQ. '') f_Fi1(i__) = iNULL__
if(ar__(874:888) .EQ. '') Fi2(i__) = rNULL__
if(ar__(890:903) .EQ. '') e_Fi2(i__) = rNULL__
if(ar__(905:916) .EQ. '') i2mag(i__) = rNULL__
if(ar__(918:932) .EQ. '') e_i2mag(i__) = rNULL__
if(ar__(934:934) .EQ. '') f_Fi2(i__) = iNULL__
if(ar__(936:950) .EQ. '') Fi3(i__) = rNULL__
if(ar__(952:966) .EQ. '') e_Fi3(i__) = rNULL__
if(ar__(968:979) .EQ. '') i3mag(i__) = rNULL__
if(ar__(981:995) .EQ. '') e_i3mag(i__) = rNULL__
if(ar__(997:997) .EQ. '') f_Fi3(i__) = iNULL__
if(ar__(999:1013) .EQ. '') Fi4(i__) = rNULL__
if(ar__(1015:1029) .EQ. '') e_Fi4(i__) = rNULL__
if(ar__(1031:1042) .EQ. '') i4mag(i__) = rNULL__
if(ar__(1044:1058) .EQ. '') e_i4mag(i__) = rNULL__
if(ar__(1060:1060) .EQ. '') f_Fi4(i__) = iNULL__
if(ar__(1062:1076) .EQ. '') Fm24(i__) = rNULL__
if(ar__(1078:1091) .EQ. '') e_Fm24(i__) = rNULL__
if(ar__(1093:1093) .EQ. '') f_Fm24(i__) = iNULL__
if(ar__(1095:1107) .EQ. '') Fpacsg(i__) = rNULL__
if(ar__(1109:1121) .EQ. '') e_Fpacsg(i__) = rNULL__
if(ar__(1123:1123) .EQ. '') f_Fpacsg(i__) = iNULL__
if(ar__(1125:1137) .EQ. '') Fpacsr(i__) = rNULL__
if(ar__(1139:1151) .EQ. '') e_Fpacsr(i__) = rNULL__
if(ar__(1153:1153) .EQ. '') f_Fpacsr(i__) = iNULL__
if(ar__(1155:1169) .EQ. '') F250sp(i__) = rNULL__
if(ar__(1171:1184) .EQ. '') e_F250sp(i__) = rNULL__
if(ar__(1186:1186) .EQ. '') f_F250sp(i__) = iNULL__
if(ar__(1188:1202) .EQ. '') F350sp(i__) = rNULL__
if(ar__(1204:1217) .EQ. '') e_F350sp(i__) = rNULL__
if(ar__(1219:1219) .EQ. '') f_F350sp(i__) = iNULL__
if(ar__(1221:1235) .EQ. '') F500sp(i__) = rNULL__
if(ar__(1237:1250) .EQ. '') e_F500sp(i__) = rNULL__
if(ar__(1252:1252) .EQ. '') f_F500sp(i__) = iNULL__
if(ar__(1254:1267) .EQ. '') SFR(i__) = rNULL__
if(ar__(1269:1282) .EQ. '') e_SFR(i__) = rNULL__
if(ar__(1284:1297) .EQ. '') M_(i__) = rNULL__
if(ar__(1299:1312) .EQ. '') e_M_(i__) = rNULL__
if(ar__(1314:1327) .EQ. '') Ldust(i__) = rNULL__
if(ar__(1329:1342) .EQ. '') e_Ldust(i__) = rNULL__
if(ar__(1344:1357) .EQ. '') LdustIR(i__) = rNULL__
if(ar__(1359:1372) .EQ. '') e_LdustIR(i__) = rNULL__
if(ar__(1374:1387) .EQ. '') Chi2Opt(i__) = rNULL__
if(ar__(1389:1402) .EQ. '') Chi2IRbest(i__) = rNULL__
if(ar__(1404:1417) .EQ. '') Chi2best(i__) = rNULL__
if(ar__(1419:1433) .EQ. '') Chi2IR(i__) = rNULL__
if(ar__(1435:1435) .EQ. '') fGaia(i__) = iNULL__
if(ar__(1437:1437) .EQ. '') fClean(i__) = iNULL__
if(ar__(1439:1439) .EQ. '') fMerg(i__) = iNULL__
if(ar__(1441:1441) .EQ. '') fOptNIRo(i__) = iNULL__
if(ar__(1443:1443) .EQ. '') fOptNIRd(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ HELP(i__),RAJ2000(i__),DEJ2000(i__),Field(i__),HEALPIX(i__),
+ E_B_V(i__),zphEAZY(i__),zsp(i__),q_zsp(i__),f_zsp(i__),
+ stellar(i__),r_stellar(i__),Fu(i__),e_Fu(i__),umag(i__),
+ e_umag(i__),f_Fu(i__),Fg(i__),e_Fg(i__),gmag(i__),e_gmag(i__),
+ f_Fg(i__),Fr(i__),e_Fr(i__),rmag(i__),e_rmag(i__),f_Fr(i__),
+ Fi(i__),e_Fi(i__),imag(i__),e_imag(i__),f_Fi(i__),Fz(i__),
+ e_Fz(i__),zmag(i__),e_zmag(i__),f_Fz(i__),Fy(i__),e_Fy(i__),
+ ymag(i__),e_ymag(i__),f_Fy(i__),Fj(i__),e_Fj(i__),jmag(i__),
+ e_jmag(i__),f_Fj(i__),Fh(i__),e_Fh(i__),hmag(i__),e_hmag(i__),
+ f_Fh(i__),Fk(i__),e_Fk(i__),kmag(i__),e_kmag(i__),f_Fk(i__),
+ Fks(i__),e_Fks(i__),ksmag(i__),e_ksmag(i__),f_Fks(i__),
+ Fi1(i__),e_Fi1(i__),i1mag(i__),e_i1mag(i__),f_Fi1(i__),
+ Fi2(i__),e_Fi2(i__),i2mag(i__),e_i2mag(i__),f_Fi2(i__),
+ Fi3(i__),e_Fi3(i__),i3mag(i__),e_i3mag(i__),f_Fi3(i__),
+ Fi4(i__),e_Fi4(i__),i4mag(i__),e_i4mag(i__),f_Fi4(i__),
+ Fm24(i__),e_Fm24(i__),f_Fm24(i__),Fpacsg(i__),e_Fpacsg(i__),
+ f_Fpacsg(i__),Fpacsr(i__),e_Fpacsr(i__),f_Fpacsr(i__),
+ F250sp(i__),e_F250sp(i__),f_F250sp(i__),F350sp(i__),
+ e_F350sp(i__),f_F350sp(i__),F500sp(i__),e_F500sp(i__),
+ f_F500sp(i__),SFR(i__),e_SFR(i__),M_(i__),e_M_(i__),
+ Ldust(i__),e_Ldust(i__),LdustIR(i__),e_LdustIR(i__),
+ Chi2Opt(i__),Chi2IRbest(i__),Chi2best(i__),Chi2IR(i__),
+ fGaia(i__),fClean(i__),fMerg(i__),fOptNIRo(i__),fOptNIRd(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end