Conversion of standardized ReadMe file for
file /./ftp/cats/VII/299 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 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 ! HELP identifier based on J2000 position
* as HELP_JHHMMSS.sss+DDMMSS.sss (help_id)
real*8 RAJ2000 ! (deg) Right Ascension (J2000) (ra)
real*8 DEJ2000 ! (deg) Declination (J2000) (dec)
character*18 Field ! Name of the field in HELP survey as
* detailed in appendix A in our study
* (field)
integer*4 HEALPIX ! HEALPix index of the source position at
* order 13 using the nested scheme (hp_idx)
real*8 E_B_V ! (mag) Galactic extinction (Schlegel et al.
* 1998ApJ...500..525S) at the object
* position (ebv)
real*8 zphEAZY ! ? Photometric redshift obtained with EAZY
* (redshift)
real*4 zsp ! ? Spectroscopic redshift collated from
* the literature (zspec)
integer*4 q_zsp ! [0/5]? Spectroscopic redshift quality
* flag FIXME code meaning (zspec_qual)
integer*4 f_zsp ! [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 ! ? 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 ! ? Name of Survey from which stellarity is
* calculated (stellarity_origin)
real*8 Fu ! (uJy) ? Total flux density in the lowest error
* u band measurement (f_best_u)
real*8 e_Fu ! (uJy) ? Uncertainty on f_best_u (ferr_best_u)
real*8 umag ! (mag) ? Total AB magnitude of the lowest error
* u band measurement (m_best_u)
real*4 e_umag ! (mag) ? Uncertainty on m_best_u (merr_best_u)
integer*4 f_Fu ! [0] Flag set to true for sources for
* spurious measurements in u (flag_best_u)
real*4 Fg ! (uJy) ? Total flux density in the lowest error
* g band measurement (f_best_g)
real*4 e_Fg ! (uJy) ? Uncertainty on f_best_g (ferr_best_g)
real*8 gmag ! (mag) ? Total AB magnitude of the lowest error
* g band measurement (m_best_g)
real*4 e_gmag ! (mag) []? Uncertainty on m_best_g (merr_best_g)
integer*4 f_Fg ! [0/1] Flag set to true for sources for
* spurious measurements in g (flag_best_g)
real*4 Fr ! (uJy) ? Total flux density in the lowest error
* r band measurement (f_best_r)
real*4 e_Fr ! (uJy) ? Uncertainty on f_best_r (ferr_best_r)
real*8 rmag ! (mag) ? Total AB magnitude of the lowest error
* g band measurement (m_best_r)
real*4 e_rmag ! (mag) []? Uncertainty on m_best_r (merr_best_r)
integer*4 f_Fr ! [0/1] Flag set to true for sources for
* spurious measurements in r (flag_best_r)
real*4 Fi ! (uJy) ? Total flux density in the lowest error
* i band measurement (f_best_i)
real*4 e_Fi ! (uJy) ? Uncertainty on f_best_i (ferr_best_i)
real*8 imag ! (mag) ? Total AB magnitude of the lowest error
* i band measurement (m_best_i)
real*4 e_imag ! (mag) ? Uncertainty on m_best_i (merr_best_i)
integer*4 f_Fi ! [0/1] Flag set to true for sources for
* spurious measurements in i (flag_best_i)
real*4 Fz ! (uJy) ? Total flux density in the lowest error
* z band measurement (f_best_z)
real*4 e_Fz ! (uJy) ? Uncertainty on f_best_z (ferr_best_z)
real*8 zmag ! (mag) ? Total AB magnitude of the lowest error
* z band measurement (m_best_z)
real*4 e_zmag ! (mag) []? Uncertainty on m_best_z (merr_best_z)
integer*4 f_Fz ! [0/1] Flag set to true for sources for
* spurious measurements in z (flag_best_z)
real*4 Fy ! (uJy) ? Total flux density in the lowest error
* y band measurement (f_best_y)
real*4 e_Fy ! (uJy) ? Uncertainty on f_best_y (ferr_best_y)
real*8 ymag ! (mag) ? Total AB magnitude of the lowest error
* y band measurement (m_best_y)
real*4 e_ymag ! (mag) ? Uncertainty on m_best_y (merr_best_y)
integer*4 f_Fy ! [0/1]? Flag set to true for sources for
* spurious measurements in y (flag_best_y)
real*4 Fj ! (uJy) ? Total flux density in the lowest error
* j band measurement (f_best_j)
real*4 e_Fj ! (uJy) ? Uncertainty on f_best_j (ferr_best_j)
real*8 jmag ! (mag) ? Total AB magnitude of the lowest error
* j band measurement (m_best_j)
real*4 e_jmag ! (mag) ? Uncertainty on m_best_j (merr_best_j)
integer*4 f_Fj ! [0/1]? Flag set to true for sources for
* spurious measurements in j (flag_best_j)
real*4 Fh ! (uJy) ? Total flux density in the lowest error
* h band measurement (f_best_h)
real*4 e_Fh ! (uJy) ? Uncertainty on f_best_h (ferr_best_h)
real*8 hmag ! (mag) ? Total AB magnitude of the lowest error
* h band measurement (m_best_h)
real*4 e_hmag ! (mag) ? Uncertainty on m_best_h (merr_best_h)
integer*4 f_Fh ! [0] Flag set to true for sources for
* spurious measurements in h (flag_best_h)
real*4 Fk ! (uJy) ? Total flux density in the lowest error
* k band measurement (f_best_k)
real*4 e_Fk ! (uJy) ? Uncertainty on f_best_k (ferr_best_k)
real*8 kmag ! (mag) ? Total AB magnitude of the lowest error
* k band measurement (m_best_k)
real*4 e_kmag ! (mag) ? Uncertainty on m_best_k (merr_best_k)
integer*4 f_Fk ! [0/1]? Flag set to true for sources for
* spurious measurements in k (flag_best_k)
real*4 Fks ! (uJy) ? Total flux density in the lowest error
* ks band measurement (f_best_ks)
real*4 e_Fks ! (uJy) ? Uncertainty on f_best_ks (ferr_best_ks)
real*8 ksmag ! (mag) ? Total AB magnitude of the lowest error
* ks band measurement (m_best_ks)
real*4 e_ksmag ! (mag) ? Uncertainty on m_best_ks (merr_best_ks)
integer*4 f_Fks ! [0] Flag set to true for sources for
* spurious measurements in ks (flag_best_ks)
real*4 Fi1 ! (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 ! (uJy) ? Uncertainty on f_irac_i1 (ferr_irac_i1)
real*8 i1mag ! (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 ! (mag) ? Uncertainty on m_irac_i1 (merr_irac_i1)
integer*4 f_Fi1 ! [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 ! (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 ! (uJy) ? Uncertainty on f_irac_i2 (ferr_irac_i2)
real*8 i2mag ! (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 ! (mag) ? Uncertainty on m_irac_i2 (merr_irac_i2)
integer*4 f_Fi2 ! [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 ! (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 ! (uJy) ? Uncertainty on f_irac_i3 (ferr_irac_i3)
real*8 i3mag ! (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 ! (mag) ? Uncertainty on m_irac_i3 (merr_irac_i3)
integer*4 f_Fi3 ! [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 ! (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 ! (uJy) ? Uncertainty on f_irac_i4 (ferr_irac_i4)
real*8 i4mag ! (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 ! (mag) ? Uncertainty on m_irac_i4 (merr_irac_i4)
integer*4 f_Fi4 ! [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 ! (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 ! (uJy) ? Uncertainty on f_mips_24 (ferr_mips_24)
integer*4 f_Fm24 ! [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 ! (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 ! (uJy) ? Uncertainty on f_pacs_green
* (ferr_pacs_green)
integer*4 f_Fpacsg ! [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 ! (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 ! (uJy) ? Uncertainty on f_pacs_red
* (ferr_pacs_red)
integer*4 f_Fpacsr ! [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 ! (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 ! (uJy) ? Uncertainty on f_spire_250
* (ferr_spire_250)
integer*4 f_F250sp ! [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 ! (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 ! (uJy) ? Uncertainty on f_spire_350
* (ferr_spire_350)
integer*4 f_F350sp ! [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 ! (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 ! (uJy) ? Uncertainty on f_spire_500
* (ferr_spire_500)
integer*4 f_F500sp ! [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 ! (Msun) ? Star Formation Rate estimated by SED
* fitting with CIGALE (cigale_sfr)
real*4 e_SFR ! (Msun) ? Uncertainty on CIGALE_SFR
* (cigale_sfr_err)
real*4 M_ ! (Msun) ? Stellar Mass estimated by SED fitting
* with CIGALE (cigale_mstar)
real*4 e_M_ ! (Msun) ? Uncertainty on CIGALE_Mstar
* (cigale_mstar_err)
real*4 Ldust ! (W) ? Dust luminosity estimated by SED
* fitting with CIGALE using both optical
* and infra-red information
* (cigale_dustlumin)
real*4 e_Ldust ! (W) ? Uncertainty on CIGALE_DustLumin
* (cigale_dustlumin_err)
real*4 LdustIR ! (W) ? Dust luminosity estimated by SED
* fitting with CIGALE using only the
* infra-red information
* (cigale_dustlumin_ironly)
real*4 e_LdustIR ! (W) ? Uncertainty on CIGALE_DustLumin_IRonly
* (cigale_dustlumin_ironly_err)
real*4 Chi2Opt ! ? Reduced optical chi squared of best
* CIGALE fit to all data (cigale_chi2_opt)
real*4 Chi2IRbest ! ? Reduced infrared chi squared of best
* CIGALE fit to all data (cigale_chi2_ir)
real*4 Chi2best ! ? Reduced chi squared of best CIGALE fit
* (cigale_chi2_red)
real*4 Chi2IR ! ? Reduced chi squared of best CIGALE fit
* to infrared data (cigale_chi2_ironly)
integer*4 fGaia ! [0/3]? Flag indicating the increasing
* probability of the source being a Gaia
* object (flag_gaia)
integer*4 fClean ! [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 ! [0/1]? Boolean flag denoting a possible
* mis-association during the cross-matching
* of the various catalogues (flag_merged)
integer*4 fOptNIRo ! [0/7]? Flag indicating in which
* wavelength regimes the source was
* observed (flag_optnir_obs)
integer*4 fOptNIRd ! [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,RAJ2000,DEJ2000,Field,HEALPIX,E_B_V,zphEAZY,zsp,q_zsp,
+ f_zsp,stellar,r_stellar,Fu,e_Fu,umag,e_umag,f_Fu,Fg,e_Fg,gmag,
+ e_gmag,f_Fg,Fr,e_Fr,rmag,e_rmag,f_Fr,Fi,e_Fi,imag,e_imag,f_Fi,
+ Fz,e_Fz,zmag,e_zmag,f_Fz,Fy,e_Fy,ymag,e_ymag,f_Fy,Fj,e_Fj,
+ jmag,e_jmag,f_Fj,Fh,e_Fh,hmag,e_hmag,f_Fh,Fk,e_Fk,kmag,e_kmag,
+ f_Fk,Fks,e_Fks,ksmag,e_ksmag,f_Fks,Fi1,e_Fi1,i1mag,e_i1mag,
+ f_Fi1,Fi2,e_Fi2,i2mag,e_i2mag,f_Fi2,Fi3,e_Fi3,i3mag,e_i3mag,
+ f_Fi3,Fi4,e_Fi4,i4mag,e_i4mag,f_Fi4,Fm24,e_Fm24,f_Fm24,Fpacsg,
+ e_Fpacsg,f_Fpacsg,Fpacsr,e_Fpacsr,f_Fpacsr,F250sp,e_F250sp,
+ f_F250sp,F350sp,e_F350sp,f_F350sp,F500sp,e_F500sp,f_F500sp,
+ SFR,e_SFR,M_,e_M_,Ldust,e_Ldust,LdustIR,e_LdustIR,Chi2Opt,
+ Chi2IRbest,Chi2best,Chi2IR,fGaia,fClean,fMerg,fOptNIRo,
+ fOptNIRd
if(ar__(105:115) .EQ. '') zphEAZY = rNULL__
if(ar__(117:131) .EQ. '') zsp = rNULL__
if(ar__(133:133) .EQ. '') q_zsp = iNULL__
if(ar__(137:147) .EQ. '') stellar = rNULL__
if(ar__(170:187) .EQ. '') Fu = rNULL__
if(ar__(189:203) .EQ. '') e_Fu = rNULL__
if(ar__(205:216) .EQ. '') umag = rNULL__
if(ar__(218:232) .EQ. '') e_umag = rNULL__
if(ar__(236:251) .EQ. '') Fg = rNULL__
if(ar__(253:268) .EQ. '') e_Fg = rNULL__
if(ar__(270:281) .EQ. '') gmag = rNULL__
if(ar__(283:297) .EQ. '') e_gmag = rNULL__
if(ar__(301:316) .EQ. '') Fr = rNULL__
if(ar__(318:332) .EQ. '') e_Fr = rNULL__
if(ar__(334:345) .EQ. '') rmag = rNULL__
if(ar__(347:361) .EQ. '') e_rmag = rNULL__
if(ar__(365:380) .EQ. '') Fi = rNULL__
if(ar__(382:397) .EQ. '') e_Fi = rNULL__
if(ar__(399:411) .EQ. '') imag = rNULL__
if(ar__(413:427) .EQ. '') e_imag = rNULL__
if(ar__(431:446) .EQ. '') Fz = rNULL__
if(ar__(448:463) .EQ. '') e_Fz = rNULL__
if(ar__(465:476) .EQ. '') zmag = rNULL__
if(ar__(478:492) .EQ. '') e_zmag = rNULL__
if(ar__(496:511) .EQ. '') Fy = rNULL__
if(ar__(513:526) .EQ. '') e_Fy = rNULL__
if(ar__(528:539) .EQ. '') ymag = rNULL__
if(ar__(541:555) .EQ. '') e_ymag = rNULL__
if(ar__(557:557) .EQ. '') f_Fy = iNULL__
if(ar__(559:573) .EQ. '') Fj = rNULL__
if(ar__(575:589) .EQ. '') e_Fj = rNULL__
if(ar__(591:602) .EQ. '') jmag = rNULL__
if(ar__(604:618) .EQ. '') e_jmag = rNULL__
if(ar__(620:620) .EQ. '') f_Fj = iNULL__
if(ar__(622:636) .EQ. '') Fh = rNULL__
if(ar__(638:651) .EQ. '') e_Fh = rNULL__
if(ar__(653:664) .EQ. '') hmag = rNULL__
if(ar__(666:680) .EQ. '') e_hmag = rNULL__
if(ar__(684:699) .EQ. '') Fk = rNULL__
if(ar__(701:715) .EQ. '') e_Fk = rNULL__
if(ar__(717:729) .EQ. '') kmag = rNULL__
if(ar__(731:745) .EQ. '') e_kmag = rNULL__
if(ar__(747:747) .EQ. '') f_Fk = iNULL__
if(ar__(749:763) .EQ. '') Fks = rNULL__
if(ar__(765:779) .EQ. '') e_Fks = rNULL__
if(ar__(781:792) .EQ. '') ksmag = rNULL__
if(ar__(794:808) .EQ. '') e_ksmag = rNULL__
if(ar__(812:826) .EQ. '') Fi1 = rNULL__
if(ar__(828:841) .EQ. '') e_Fi1 = rNULL__
if(ar__(843:854) .EQ. '') i1mag = rNULL__
if(ar__(856:870) .EQ. '') e_i1mag = rNULL__
if(ar__(872:872) .EQ. '') f_Fi1 = iNULL__
if(ar__(874:888) .EQ. '') Fi2 = rNULL__
if(ar__(890:903) .EQ. '') e_Fi2 = rNULL__
if(ar__(905:916) .EQ. '') i2mag = rNULL__
if(ar__(918:932) .EQ. '') e_i2mag = rNULL__
if(ar__(934:934) .EQ. '') f_Fi2 = iNULL__
if(ar__(936:950) .EQ. '') Fi3 = rNULL__
if(ar__(952:966) .EQ. '') e_Fi3 = rNULL__
if(ar__(968:979) .EQ. '') i3mag = rNULL__
if(ar__(981:995) .EQ. '') e_i3mag = rNULL__
if(ar__(997:997) .EQ. '') f_Fi3 = iNULL__
if(ar__(999:1013) .EQ. '') Fi4 = rNULL__
if(ar__(1015:1029) .EQ. '') e_Fi4 = rNULL__
if(ar__(1031:1042) .EQ. '') i4mag = rNULL__
if(ar__(1044:1058) .EQ. '') e_i4mag = rNULL__
if(ar__(1060:1060) .EQ. '') f_Fi4 = iNULL__
if(ar__(1062:1076) .EQ. '') Fm24 = rNULL__
if(ar__(1078:1091) .EQ. '') e_Fm24 = rNULL__
if(ar__(1093:1093) .EQ. '') f_Fm24 = iNULL__
if(ar__(1095:1107) .EQ. '') Fpacsg = rNULL__
if(ar__(1109:1121) .EQ. '') e_Fpacsg = rNULL__
if(ar__(1123:1123) .EQ. '') f_Fpacsg = iNULL__
if(ar__(1125:1137) .EQ. '') Fpacsr = rNULL__
if(ar__(1139:1151) .EQ. '') e_Fpacsr = rNULL__
if(ar__(1153:1153) .EQ. '') f_Fpacsr = iNULL__
if(ar__(1155:1169) .EQ. '') F250sp = rNULL__
if(ar__(1171:1184) .EQ. '') e_F250sp = rNULL__
if(ar__(1186:1186) .EQ. '') f_F250sp = iNULL__
if(ar__(1188:1202) .EQ. '') F350sp = rNULL__
if(ar__(1204:1217) .EQ. '') e_F350sp = rNULL__
if(ar__(1219:1219) .EQ. '') f_F350sp = iNULL__
if(ar__(1221:1235) .EQ. '') F500sp = rNULL__
if(ar__(1237:1250) .EQ. '') e_F500sp = rNULL__
if(ar__(1252:1252) .EQ. '') f_F500sp = iNULL__
if(ar__(1254:1267) .EQ. '') SFR = rNULL__
if(ar__(1269:1282) .EQ. '') e_SFR = rNULL__
if(ar__(1284:1297) .EQ. '') M_ = rNULL__
if(ar__(1299:1312) .EQ. '') e_M_ = rNULL__
if(ar__(1314:1327) .EQ. '') Ldust = rNULL__
if(ar__(1329:1342) .EQ. '') e_Ldust = rNULL__
if(ar__(1344:1357) .EQ. '') LdustIR = rNULL__
if(ar__(1359:1372) .EQ. '') e_LdustIR = rNULL__
if(ar__(1374:1387) .EQ. '') Chi2Opt = rNULL__
if(ar__(1389:1402) .EQ. '') Chi2IRbest = rNULL__
if(ar__(1404:1417) .EQ. '') Chi2best = rNULL__
if(ar__(1419:1433) .EQ. '') Chi2IR = rNULL__
if(ar__(1435:1435) .EQ. '') fGaia = iNULL__
if(ar__(1437:1437) .EQ. '') fClean = iNULL__
if(ar__(1439:1439) .EQ. '') fMerg = iNULL__
if(ar__(1441:1441) .EQ. '') fOptNIRo = iNULL__
if(ar__(1443:1443) .EQ. '') fOptNIRd = iNULL__
c ..............Just test output...........
write(6,1)
+ HELP,RAJ2000,DEJ2000,Field,HEALPIX,E_B_V,zphEAZY,zsp,q_zsp,
+ f_zsp,stellar,r_stellar,Fu,e_Fu,umag,e_umag,f_Fu,Fg,e_Fg,gmag,
+ e_gmag,f_Fg,Fr,e_Fr,rmag,e_rmag,f_Fr,Fi,e_Fi,imag,e_imag,f_Fi,
+ Fz,e_Fz,zmag,e_zmag,f_Fz,Fy,e_Fy,ymag,e_ymag,f_Fy,Fj,e_Fj,
+ jmag,e_jmag,f_Fj,Fh,e_Fh,hmag,e_hmag,f_Fh,Fk,e_Fk,kmag,e_kmag,
+ f_Fk,Fks,e_Fks,ksmag,e_ksmag,f_Fks,Fi1,e_Fi1,i1mag,e_i1mag,
+ f_Fi1,Fi2,e_Fi2,i2mag,e_i2mag,f_Fi2,Fi3,e_Fi3,i3mag,e_i3mag,
+ f_Fi3,Fi4,e_Fi4,i4mag,e_i4mag,f_Fi4,Fm24,e_Fm24,f_Fm24,Fpacsg,
+ e_Fpacsg,f_Fpacsg,Fpacsr,e_Fpacsr,f_Fpacsr,F250sp,e_F250sp,
+ f_F250sp,F350sp,e_F350sp,f_F350sp,F500sp,e_F500sp,f_F500sp,
+ SFR,e_SFR,M_,e_M_,Ldust,e_Ldust,LdustIR,e_LdustIR,Chi2Opt,
+ Chi2IRbest,Chi2best,Chi2IR,fGaia,fClean,fMerg,fOptNIRo,
+ fOptNIRd
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end