Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/789/L38 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-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. J/ApJ/789/L38 LRLL54361 protostar Herschel/PACS fluxes (Balog+, 2014)
*================================================================================
*The extraordinary far-infrared variation of a protostar: Herschel/PACS
*observations of LRLL54361.
* Balog Z., Muzerolle J., Flaherty K., Detre O.H., Bouwmann J., Furlan E.,
* Gutermuth R., Juhasz A., Bally J., Nielbock M., Klaas U., Krause O.,
* Henning T., Marton G.
* <Astrophys. J., 789, L38 (2014)>
* =2014ApJ...789L..38B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Observing log of the observations and measured
fluxes
integer*4 nr__
parameter (nr__=33) ! Number of records
character*80 ar__ ! Full-size record
integer*4 Cycle (nr__) ! ? Cycle number (1)
integer*4 Obs_Y (nr__) ! (yr) Year of the observation
integer*4 Obs_M (nr__) ! ("month") Month of the observation
integer*4 Obs_D (nr__) ! (d) Day of the observation
integer*4 Obs_h (nr__) ! (h) Hour of the observation
integer*4 Obs_m_1 (nr__) ! (min) Minute of the observation
real*4 Obs_s (nr__) ! (s) Second of the observation
real*8 JD (nr__) ! (d) Julian Date of the observation
real*4 Phase (nr__) ! Phase (2)
real*4 F70 (nr__) ! (Jy) ? Herschel/PACS 70 micron band flux density
real*4 e_F70 (nr__) ! (Jy) ? Error in F70
real*4 F100 (nr__) ! (Jy) ? Herschel/PACS 100 micron band flux density
real*4 e_F100 (nr__) ! (Jy) ? Error in F100
real*4 F160 (nr__) ! (Jy) ? Herschel/PACS 160 micron band flux density
real*4 e_F160 (nr__) ! (Jy) ? Error in F160
*Note (1): The last four lines without Cycle values give the spectrocopy
* information. The first 29 lines with Cycle values are the photometry.
*Note (2): We calculated the phase of the observations based on the period of
* 25.34 days and zero point JD 2455121.203 reported by
* Muzerolle et al. (2013, J/other/Nat/493.378). Phase=0 corresponds to
* the peak flux, as measured from the Spitzer/IRAC observations.
C=============================================================================
C Loading file 'table1.dat' ! Observing log of the observations and measured
* fluxes
C Format for file interpretation
1 format(
+ I1,1X,I4,1X,I2,1X,I2,1X,I2,1X,I2,1X,F4.1,1X,F11.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,33
read(1,'(A80)')ar__
read(ar__,1)
+ Cycle(i__),Obs_Y(i__),Obs_M(i__),Obs_D(i__),Obs_h(i__),
+ Obs_m_1(i__),Obs_s(i__),JD(i__),Phase(i__),F70(i__),
+ e_F70(i__),F100(i__),e_F100(i__),F160(i__),e_F160(i__)
if(ar__(1:1) .EQ. '') Cycle(i__) = iNULL__
if(ar__(43:48) .EQ. '') F70(i__) = rNULL__
if(ar__(50:54) .EQ. '') e_F70(i__) = rNULL__
if(ar__(56:61) .EQ. '') F100(i__) = rNULL__
if(ar__(63:67) .EQ. '') e_F100(i__) = rNULL__
if(ar__(69:74) .EQ. '') F160(i__) = rNULL__
if(ar__(76:80) .EQ. '') e_F160(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Cycle(i__),Obs_Y(i__),Obs_M(i__),Obs_D(i__),Obs_h(i__),
+ Obs_m_1(i__),Obs_s(i__),JD(i__),Phase(i__),F70(i__),
+ e_F70(i__),F100(i__),e_F100(i__),F160(i__),e_F160(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end