Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/410/735 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-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. J/MNRAS/410/735 MAG waves in magnetic and rotating stars (Campos, 2011)
*================================================================================
*On magnetoacoustic-gravity-inertial (MAGI) waves.
*II. Application to magnetic and rotating stars.
* Campos L.M.B.C.
* <Mon. Not. R. Astron. Soc., 410, 735-761 (2011)>
* =2011MNRAS.410..735C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Magnetic fields (Gauss) of non-degenerate stars
integer*4 nr__
parameter (nr__=23) ! Number of records
character*41 ar__ ! Full-size record
integer*4 Seq (nr__) ! [1/23] Sequential number of the star
integer*4 HD (nr__) ! Star HD number
real*4 Vmag (nr__) ! (mag) V magnitude
integer*4 Bmin (nr__) ! (0.1mT) Minimum magnetic field strength (Gauss)
integer*4 Bmax (nr__) ! (0.1mT) Maximum magnetic field strength (Gauss)
real*4 log_Bmin (nr__) ! ([0.1mT]) log of minimum magnetic field strength
real*4 log_Bmax (nr__) ! ([0.1mT]) log of maximum magnetic field strength
real*4 G (nr__) ! ([Lsun]) ?=- Total brightness
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Magnetic field (Gauss) and rotation period (days)
of Ap stars
integer*4 nr__1
parameter (nr__1=51) ! Number of records
character*68 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! [24/75] Sequential number of the star
integer*4 HD_1 (nr__1) ! Star HD number
real*4 Vmag_1 (nr__1) ! (mag) V magnitude
integer*4 Bmin_1 (nr__1) ! (0.1mT) Minimum magnetic field strength (Gauss)
integer*4 Bmax_1 (nr__1) ! (0.1mT) Maximum magnetic field strength (Gauss)
real*8 tau (nr__1) ! (s) Rotation period {tau}
real*4 log_Bmin_1 (nr__1) ! ([0.1mT]) log of minimum magnetic field strength
real*4 log_Bmax_1 (nr__1) ! ([0.1mT]) log of maximum magnetic field strength
real*8 log_tau (nr__1) ! ([s]) log of rotation period [log({tau})]
real*4 Gm (nr__1) ! ([Lsun]) ?=- Magnetic brightness
real*4 Gr (nr__1) ! ([Lsun]) ?=- Rotational brightness
real*4 G_1 (nr__1) ! ([Lsun]) ?=- Total brightness
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Observed and calculated data for strong
non-rotating magnetic stars
integer*4 nr__2
parameter (nr__2=21) ! Number of records
character*144 ar__2 ! Full-size record
integer*4 Seq_2 (nr__2) ! [0/23] Star sequential number (0 for Sun)
real*4 M_ (nr__2) ! (Msun) Observed stellar mass
real*4 Bsup (nr__2) ! (0.1mT) Observed supremum of maximum and minimum (Gauss)
* magnetic field in modulus
real*4 I (nr__2) ! ([Lsun]) Stellar intensity (Eq. 30)
real*4 G_2 (nr__2) ! ([Lsun]) Stellar brightness (Eq. 31,34b)
real*4 T (nr__2) ! (K) Effective temperature (Eq. 38b)
real*4 L (nr__2) ! (10-7W) Luminostiy (Eq. 25,26b)
real*4 r (nr__2) ! (cm) Radius (Eq. 28)
real*4 p (nr__2) ! (dPa) Total pressure in dyn/cm2 (Eq. 51b)
real*4 rho (nr__2) ! (g/cm2) Mass density (Eq. 57b)
real*4 m (nr__2) ! (g) Stellar mass (Eq. 138)
real*4 MLR (nr__2) ! (g/s) Stellar mass-loss rate, -dm/dt (Eq. 83,78a)
real*4 Vr (nr__2) ! (cm/s) Radial velocity (Eq. 87)
real*4 c (nr__2) ! (cm/s) Sound speed (Eq. 108a)
real*4 Mr (nr__2) ! Radial Mach number (Eq. 123a)
real*4 a (nr__2) ! (cm/s) Alfven speed (Eq. 114a)
real*4 Nr (nr__2) ! Radial Alfven number (Eq. 126a)
real*4 u (nr__2) ! (cm/s) Magnetosonic wave speed (Eq. 121a)
real*4 Qr (nr__2) ! Radial hydromagnetic number (Eq. 130d)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Observed and calculated data for fast
rotating magnetic stars
integer*4 nr__3
parameter (nr__3=18) ! Number of records
character*225 ar__3 ! Full-size record
integer*4 Order (nr__3) ! [0/17] Order
integer*4 Seq_3 (nr__3) ! [28/73]?=- Star sequential number
real*4 M__1 (nr__3) ! (Msun) ?=- Observed stellar mass
real*4 Bsup_1 (nr__3) ! (0.1mT) Observed supremum of maximum and minimum of
* magnetic field in modulus (Gauss)
real*8 tau_1 (nr__3) ! (d) Observed period {tau}
real*4 I_1 (nr__3) ! ([Lsun]) Stellar intensity (Eq. 40)
real*4 G_3 (nr__3) ! ([Lsun]) Stellar brightness (Eq. 31)
real*4 T_1 (nr__3) ! (K) Effective temperature (Eq. 38b)
real*4 L_1 (nr__3) ! (mW/m2) Luminosity (Eq. 25,26b)
real*4 r_1 (nr__3) ! (cm) Radius (Eq. 28)
real*4 P_1 (nr__3) ! (dPa) Total pressure in dyn/cm2 (Eq. 51b)
real*4 rho_1 (nr__3) ! (g/cm3) Mass density (Eq. 57b)
real*4 m_1 (nr__3) ! (g) Stellar mass (Eq. 138)
real*4 MLR_1 (nr__3) ! (g/s) Stellar mass-loss rate, -dm/dt (Eq. 83,78a)
real*4 Vr_1 (nr__3) ! (cm/s) Radial velocity (Eq. 87)
real*4 Vt (nr__3) ! (cm/s) Tangential velocity (Eq. 71b)
real*4 V (nr__3) ! (cm/s) Total velocity (Eq. 93b)
real*4 alpha (nr__3) ! (deg) Spiral angle with the radial direction (Eq. 101b)
real*4 c_1 (nr__3) ! (cm/s) Sound speed (Eq. 125a)
real*4 Mr_1 (nr__3) ! Radial Mach number (Eq. 125a)
real*4 a_1 (nr__3) ! (cm/s) Alfven speed (Eq. 125a)
real*4 Nr_1 (nr__3) ! Radial Alfven number(Eq. 125a)
real*4 u_1 (nr__3) ! (cm/s) Magnetosonic wave speed (Eq. 121a)
real*4 Qr_1 (nr__3) ! Radial hydromagnetic number (Eq. 128a)
real*4 Mt (nr__3) ! Tangential Mach number (Eq. 122a)
real*4 Nt (nr__3) ! Tangential Alfven number (Eq. 125a)
real*4 Qt (nr__3) ! Tangential hydromagnetic number (Eq. 128a)
real*4 M_2 (nr__3) ! Total Mach number (Eq. 122a)
real*4 N (nr__3) ! Total Alfven number (Eq. 125a)
real*4 Q (nr__3) ! Total hydromagnetic number (Eq. 128a)
C=============================================================================
C Loading file 'table2.dat' ! Magnetic fields (Gauss) of non-degenerate stars
C Format for file interpretation
1 format(I2,1X,I6,1X,F3.1,2X,I5,1X,I5,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,23
read(1,'(A41)')ar__
read(ar__,1)
+ Seq(i__),HD(i__),Vmag(i__),Bmin(i__),Bmax(i__),log_Bmin(i__),
+ log_Bmax(i__),G(i__)
if (idig(ar__(38:41)).EQ.0) G(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),HD(i__),Vmag(i__),Bmin(i__),Bmax(i__),log_Bmin(i__),
+ log_Bmax(i__),G(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Magnetic field (Gauss) and rotation period (days)
* of Ap stars
C Format for file interpretation
2 format(
+ I2,1X,I6,1X,F3.1,1X,I6,1X,I5,1X,F9.5,1X,F4.2,1X,F4.2,1X,F7.4,
+ 1X,F3.1,1X,F5.2,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,51
read(1,'(A68)')ar__1
read(ar__1,2)
+ Seq_1(i__),HD_1(i__),Vmag_1(i__),Bmin_1(i__),Bmax_1(i__),
+ tau(i__),log_Bmin_1(i__),log_Bmax_1(i__),log_tau(i__),Gm(i__),
+ Gr(i__),G_1(i__)
if (idig(ar__1(56:58)).EQ.0) Gm(i__) = rNULL__
if (idig(ar__1(60:64)).EQ.0) Gr(i__) = rNULL__
if (idig(ar__1(65:68)).EQ.0) G_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),HD_1(i__),Vmag_1(i__),Bmin_1(i__),Bmax_1(i__),
+ tau(i__),log_Bmin_1(i__),log_Bmax_1(i__),log_tau(i__),Gm(i__),
+ Gr(i__),G_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Observed and calculated data for strong
* non-rotating magnetic stars
C Format for file interpretation
3 format(
+ I2,1X,F3.1,1X,E8.3,F5.2,1X,F6.3,1X,E7.3,1X,E8.4,1X,E8.3,1X,
+ E7.3,1X,E7.3,1X,E8.3,1X,E8.3,1X,E7.3,1X,E7.3,1X,E7.3,1X,E8.4,
+ 1X,E7.3,1X,E7.3,1X,E7.3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,21
read(1,'(A144)')ar__2
read(ar__2,3)
+ Seq_2(i__),M_(i__),Bsup(i__),I(i__),G_2(i__),T(i__),L(i__),
+ r(i__),p(i__),rho(i__),m(i__),MLR(i__),Vr(i__),c(i__),Mr(i__),
+ a(i__),Nr(i__),u(i__),Qr(i__)
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),M_(i__),Bsup(i__),I(i__),G_2(i__),T(i__),L(i__),
+ r(i__),p(i__),rho(i__),m(i__),MLR(i__),Vr(i__),c(i__),Mr(i__),
+ a(i__),Nr(i__),u(i__),Qr(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Observed and calculated data for fast
* rotating magnetic stars
C Format for file interpretation
4 format(
+ I2,1X,I2,1X,F3.1,1X,E7.3,1X,F8.5,1X,F5.2,1X,F6.3,1X,E7.3,1X,
+ E8.3,1X,E8.3,1X,E7.3,1X,E7.3,1X,E8.3,1X,E8.3,1X,E7.3,1X,E7.3,
+ 1X,E7.3,1X,F4.1,1X,E7.3,1X,E7.3,1X,E7.3,1X,E7.3,1X,E7.3,1X,
+ E7.3,1X,E7.3,1X,E7.3,1X,E7.3,1X,E7.3,1X,E7.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,18
read(1,'(A225)')ar__3
read(ar__3,4)
+ Order(i__),Seq_3(i__),M__1(i__),Bsup_1(i__),tau_1(i__),
+ I_1(i__),G_3(i__),T_1(i__),L_1(i__),r_1(i__),P_1(i__),
+ rho_1(i__),m_1(i__),MLR_1(i__),Vr_1(i__),Vt(i__),V(i__),
+ alpha(i__),c_1(i__),Mr_1(i__),a_1(i__),Nr_1(i__),u_1(i__),
+ Qr_1(i__),Mt(i__),Nt(i__),Qt(i__),M_2(i__),N(i__),Q(i__)
if (Seq_3(i__) .EQ. 45) Seq_3(i__) = iNULL__
if (idig(ar__3(7:9)).EQ.0) M__1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Order(i__),Seq_3(i__),M__1(i__),Bsup_1(i__),tau_1(i__),
+ I_1(i__),G_3(i__),T_1(i__),L_1(i__),r_1(i__),P_1(i__),
+ rho_1(i__),m_1(i__),MLR_1(i__),Vr_1(i__),Vt(i__),V(i__),
+ alpha(i__),c_1(i__),Mr_1(i__),a_1(i__),Nr_1(i__),u_1(i__),
+ Qr_1(i__),Mt(i__),Nt(i__),Qt(i__),M_2(i__),N(i__),Q(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end