Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/181/321 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-12
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/ApJS/181/321 Properties of Spitzer c2d dark clouds (Evans+, 2009)
*================================================================================
*The Spitzer c2d legacy results:
*star-formation rates and efficiencies: evolution and lifetimes.
* Evans N.J., Dunham M.M., Jorgensen J.K., Enoch M.L., Merin B.,
* Van Dishoeck E.F., Alcala J.M., Myers P.C., Stapelfeldt K.R., Huard T.L.,
* Allen L.E., Harvey P.M., Van Kempen T., Blake G.A., Koerner D.W.,
* Mundy L.G., Padgett D.L., Sargent A.I.
* <Astrophys. J. Suppl. Ser., 181, 321-350 (2009)>
* =2009ApJS..181..321E
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! 2MASS and Spitzer flux densities of YSOs in the
c2d clouds
integer*4 nr__
parameter (nr__=1024) ! Number of records
character*173 ar__ ! Full-size record
integer*4 v_EDJ2009_ (nr__) ! Running source index number
real*8 F1_25 (nr__) ! (mJy) ? 2MASS J band flux density
real*8 e_F1_25 (nr__) ! (mJy) ? Error in F1.25
real*8 F1_65 (nr__) ! (mJy) ? 2MASS H band flux density
real*8 e_F1_65 (nr__) ! (mJy) ? Error in F1.65
real*8 F2_17 (nr__) ! (mJy) ? 2MASS K_S_ band flux density
real*8 e_F2_17 (nr__) ! (mJy) ? Error in F2.17
real*8 F3_6 (nr__) ! (mJy) ? IRAC 3.6 micron flux density
real*8 e_F3_6 (nr__) ! (mJy) ? Error in F3.6
real*8 F4_5 (nr__) ! (mJy) ? IRAC 4.5 micron flux density
real*8 e_F4_5 (nr__) ! (mJy) ? Error in F4.5
real*8 F5_8 (nr__) ! (mJy) ? IRAC 5.8 micron flux density
real*8 e_F5_8 (nr__) ! (mJy) ? Error in F5.8
real*8 F8_0 (nr__) ! (mJy) ? IRAC 8.0 micron flux density
real*8 e_F8_0 (nr__) ! (mJy) ? Error in F8.0
real*8 F24 (nr__) ! (mJy) ? MIPS 24 micron flux density
real*8 e_F24 (nr__) ! (mJy) ? Error in F24
integer*4 F70 (nr__) ! (mJy) ? MIPS 70 micron flux density
real*8 e_F70 (nr__) ! (mJy) ? Error in F70
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! 0.36-0.96um flux densities of YSOs in the c2d clouds
integer*4 nr__1
parameter (nr__1=1024) ! Number of records
character*136 ar__1 ! Full-size record
integer*4 v_EDJ2009__1(nr__1) ! Running source index number
real*4 F0_36 (nr__1) ! (mJy) ? The U band flux density
real*8 e_F0_36 (nr__1) ! (mJy) ? Error in F0.36
real*8 F0_44 (nr__1) ! (mJy) ? The B band flux density
real*8 e_F0_44 (nr__1) ! (mJy) ? Error in F0.44
real*4 F0_55 (nr__1) ! (mJy) ? The V band flux density
real*4 e_F0_55 (nr__1) ! (mJy) ? Error in F0.55
real*8 F0_64 (nr__1) ! (mJy) ? The R band flux density
real*8 e_F0_64 (nr__1) ! (mJy) ? Error in F0.64
real*8 F0_665 (nr__1) ! (mJy) ? The H{alpha}12 flux density
real*8 e_F0_665 (nr__1) ! (mJy) ? Error in F0.665
real*8 F0_79 (nr__1) ! (mJy) ? The I band flux density
real*8 e_F0_79 (nr__1) ! (mJy) ? Error in F0.79
real*8 F0_915 (nr__1) ! (mJy) ? The m914 band flux density
real*8 e_F0_915 (nr__1) ! (mJy) ? Error in F0.915
real*8 F0_96 (nr__1) ! (mJy) ? The z band flux density
real*8 e_F0_96 (nr__1) ! (mJy) ? Error in F0.96
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table10.dat' ! 3.4-100um flux densities of YSOs in the c2d clouds
integer*4 nr__2
parameter (nr__2=1024) ! Number of records
character*104 ar__2 ! Full-size record
integer*4 v_EDJ2009__2(nr__2) ! Running source index number
integer*4 F3_4 (nr__2) ! (mJy) ? The L band flux density
real*4 e_F3_4 (nr__2) ! (mJy) ? Error in F3.4
integer*4 F5_0 (nr__2) ! (mJy) ? The M band flux density
integer*4 e_F5_0 (nr__2) ! (mJy) ? Error in F5.0
real*4 F6_7 (nr__2) ! (mJy) ? ISO 6.7 micron flux density
real*4 e_F6_7 (nr__2) ! (mJy) ? Error in F6.7
integer*4 F12 (nr__2) ! (mJy) ? IRAS 12 micron flux density
real*4 e_F12 (nr__2) ! (mJy) ? Error in F12
real*4 F14_3 (nr__2) ! (mJy) ? ISO 14.3 micron flux density
real*4 e_F14_3 (nr__2) ! (mJy) ? Error in F14.3
integer*4 F25 (nr__2) ! (mJy) ? IRAS 25 micron flux density
real*4 e_F25 (nr__2) ! (mJy) ? Error in F25
integer*4 F60 (nr__2) ! (mJy) ? IRAS 60 micron flux density
integer*4 e_F60 (nr__2) ! (mJy) ? Error in F60
integer*4 F100 (nr__2) ! (mJy) ? IRAS 100 micron flux density
integer*4 e_F100 (nr__2) ! (mJy) ? Error in F100
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! 160-1300um flux densities of YSOs in the c2d clouds
integer*4 nr__3
parameter (nr__3=1024) ! Number of records
character*79 ar__3 ! Full-size record
integer*4 v_EDJ2009__3(nr__3) ! Running source index number
integer*4 F160 (nr__3) ! (mJy) ? MIPS 160 micron flux density
integer*4 e_F160 (nr__3) ! (mJy) ? Error in F160
integer*4 F350 (nr__3) ! (mJy) ? The 350 micron flux density
integer*4 e_F350 (nr__3) ! (mJy) ? Error in F350
integer*4 F450 (nr__3) ! (mJy) ? The 450 micron flux density
integer*4 e_F450 (nr__3) ! (mJy) ? Error in F450
integer*4 F850 (nr__3) ! (mJy) ? The 850 micron flux density
integer*4 e_F850 (nr__3) ! (mJy) ? Error in F850
integer*4 F1100 (nr__3) ! (mJy) ? The 1100 micron flux density
real*4 e_F1100 (nr__3) ! (mJy) ? Error in F1100
integer*4 F1200 (nr__3) ! (mJy) ? The 1200 micron flux density
integer*4 e_F1200 (nr__3) ! (mJy) ? Error in F1200
integer*4 F1300 (nr__3) ! (mJy) ? The 1300 micron flux density
real*4 e_F1300 (nr__3) ! (mJy) ? Error in F1300
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! Properties of YSOs in the c2d clouds
integer*4 nr__4
parameter (nr__4=1024) ! Number of records
character*112 ar__4 ! Full-size record
integer*4 v_EDJ2009__4(nr__4) ! Running source index number
character*7 Cloud (nr__4) ! Cloud name containing the source
character*19 SSTc2d (nr__4) ! Spitzer source name (JHHMMSS.ss+DDMMSS.s)
character*19 Class (nr__4) ! c2d classification (1)
real*4 AV (nr__4) ! (mag) The V band extinction (2)
real*4 alpha (nr__4) ! ? Observed 2-24{mu}m spectral index
character*2 l_Tbol (nr__4) ! [<>= ] Limit flag on Tbol
integer*4 Tbol (nr__4) ! (K) ? Observed bolometric temperature
character*2 l_Lbol (nr__4) ! [<>= ] Limit flag on Lbol
real*8 Lbol (nr__4) ! (solLum) ? Observed bolometric luminosity
real*4 alphac (nr__4) ! ? Corrected 2-24{mu}m spectral index
character*2 l_Tbolc (nr__4) ! [<>= ] Limit flag on Tbolc
integer*4 Tbolc (nr__4) ! (K) ? Corrected bolometric temperature
character*2 l_Lbolc (nr__4) ! [<>= ] Limit flag on Lbolc
real*8 Lbolc (nr__4) ! (solLum) ? Corrected bolometric luminosity
character*1 Env (nr__4) ! Envelope evidence? (3)
character*4 Qual (nr__4) ! Quality of calculated Lbol and Tbol (4)
*Note (1): Only sources with detections in at least three bands could be
* classified at all. Sources are labelled as follows:
* Galc = Galaxy candidates
* YSOc = Candidates for YSOs. Detections in all four IRAC bands and MIPS-1;
* see selection criteria in section 3.1.
* red = Flux densities at 24um that are at least 3 times the flux density from
* the nearest available IRAC band.
* PAH_em = Have colors indicative of a peak in the 8um band.
* star+dust = The SED is consistent with that of a stellar photosphere for
* wavelengths shorter than a particular band but an excess in that band
* of at least 3{sigma}. This band was appended to make a full name such
* as "YSOc_star+dust(MP1)" where MP1 implies the first MIPS band,
* at 24um or "YSOc_star+dust(IR4)" where IR4 implies the fourth
* IRAC band, at 8.0{mu}m.
* rising = Lack photometry at enough wavelengths to be classified as YSOc, but
* were later added as YSOs based on other information. The available
* flux densities rise significantly to longer wavelengths (15 cases)
* red1 = Lack photometry at enough wavelengths to be classified as YSOc, but
* were later added as YSOs based on other information. Source is
* detected only in IRAC-4 or MIPS-1 (3 cases).
*Note (2): Value of A_V_ used for de-reddening.
*Note (3): Indicates whether or not the source is associated with an
* envelope as traced by extended millimeter emission (based primarily
* on the sample of Enoch et al., 2008ApJ...684.1240E).
*Note (4): Quality flag as follows:
* G = Existing photometry provides sufficient spectral coverage for
* reliable calculations.
* L_I = the source is associated with an envelope and thus likely an
* embedded source, but lacks spectral coverage between 24{mu}m
* and 1mm. The calculated values of Lbol and Tbol are thus
* considered lower and upper limits, respectively.
* L_II = the source is not associated with an envelope and has an observed
* {alpha}{<=}-0.3, thus is not likely an embedded source, but no
* extinction correction could be derived based on available data,
* usually because of a lack of sufficient coverage in the near-IR.
* The calculated values of Lbol and Tbol are thus considered
* lower limits.
* O = Other. These three sources were added to the list of YSOs based
* on the searches for embedded protostars presented by Jorgensen
* et al. (2007ApJ...656..293J; for Perseus) and Jorgensen et al.
* (2008ApJ...683..822J; for Ophiuchus), but they are not associated
* with envelopes according to Enoch et al. (2008ApJ...684.1240).
* They are only detected at one wavelength (24{mu}m), thus
* calculations of {alpha}, Tbol, and Lbol are not possible. They are
* counted as Class I objects when classifying by {alpha} and {alpha}c;
* they are not counted at all when classifying by Tbol and Tbolc.
C=============================================================================
C Loading file 'table8.dat' ! 2MASS and Spitzer flux densities of YSOs in the
* c2d clouds
C Format for file interpretation
1 format(
+ I4,1X,F9.3,1X,F11.5,1X,F9.3,1X,F9.4,1X,F8.2,1X,F8.3,1X,F9.4,
+ 1X,F9.5,1X,F8.3,1X,F9.5,1X,F9.3,1X,F8.3,1X,F8.2,1X,F8.3,1X,
+ F7.2,1X,F8.3,1X,I6,1X,F8.1)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,1024
read(1,'(A173)')ar__
read(ar__,1)
+ v_EDJ2009_(i__),F1_25(i__),e_F1_25(i__),F1_65(i__),
+ e_F1_65(i__),F2_17(i__),e_F2_17(i__),F3_6(i__),e_F3_6(i__),
+ F4_5(i__),e_F4_5(i__),F5_8(i__),e_F5_8(i__),F8_0(i__),
+ e_F8_0(i__),F24(i__),e_F24(i__),F70(i__),e_F70(i__)
if(ar__(6:14) .EQ. '') F1_25(i__) = rNULL__
if(ar__(16:26) .EQ. '') e_F1_25(i__) = rNULL__
if(ar__(28:36) .EQ. '') F1_65(i__) = rNULL__
if(ar__(38:46) .EQ. '') e_F1_65(i__) = rNULL__
if(ar__(48:55) .EQ. '') F2_17(i__) = rNULL__
if(ar__(57:64) .EQ. '') e_F2_17(i__) = rNULL__
if(ar__(66:74) .EQ. '') F3_6(i__) = rNULL__
if(ar__(76:84) .EQ. '') e_F3_6(i__) = rNULL__
if(ar__(86:93) .EQ. '') F4_5(i__) = rNULL__
if(ar__(95:103) .EQ. '') e_F4_5(i__) = rNULL__
if(ar__(105:113) .EQ. '') F5_8(i__) = rNULL__
if(ar__(115:122) .EQ. '') e_F5_8(i__) = rNULL__
if(ar__(124:131) .EQ. '') F8_0(i__) = rNULL__
if(ar__(133:140) .EQ. '') e_F8_0(i__) = rNULL__
if(ar__(142:148) .EQ. '') F24(i__) = rNULL__
if(ar__(150:157) .EQ. '') e_F24(i__) = rNULL__
if(ar__(159:164) .EQ. '') F70(i__) = iNULL__
if(ar__(166:173) .EQ. '') e_F70(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ v_EDJ2009_(i__),F1_25(i__),e_F1_25(i__),F1_65(i__),
+ e_F1_65(i__),F2_17(i__),e_F2_17(i__),F3_6(i__),e_F3_6(i__),
+ F4_5(i__),e_F4_5(i__),F5_8(i__),e_F5_8(i__),F8_0(i__),
+ e_F8_0(i__),F24(i__),e_F24(i__),F70(i__),e_F70(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! 0.36-0.96um flux densities of YSOs in the c2d clouds
C Format for file interpretation
2 format(
+ I4,1X,F4.2,1X,F7.5,1X,F8.5,1X,F7.5,1X,F5.2,1X,F6.4,1X,F9.5,1X,
+ F8.5,1X,F8.5,1X,F8.5,1X,F7.3,1X,F8.5,1X,F7.3,1X,F8.5,1X,F7.3,
+ 1X,F9.5)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,1024
read(1,'(A136)')ar__1
read(ar__1,2)
+ v_EDJ2009__1(i__),F0_36(i__),e_F0_36(i__),F0_44(i__),
+ e_F0_44(i__),F0_55(i__),e_F0_55(i__),F0_64(i__),e_F0_64(i__),
+ F0_665(i__),e_F0_665(i__),F0_79(i__),e_F0_79(i__),F0_915(i__),
+ e_F0_915(i__),F0_96(i__),e_F0_96(i__)
if(ar__1(6:9) .EQ. '') F0_36(i__) = rNULL__
if(ar__1(11:17) .EQ. '') e_F0_36(i__) = rNULL__
if(ar__1(19:26) .EQ. '') F0_44(i__) = rNULL__
if(ar__1(28:34) .EQ. '') e_F0_44(i__) = rNULL__
if(ar__1(36:40) .EQ. '') F0_55(i__) = rNULL__
if(ar__1(42:47) .EQ. '') e_F0_55(i__) = rNULL__
if(ar__1(49:57) .EQ. '') F0_64(i__) = rNULL__
if(ar__1(59:66) .EQ. '') e_F0_64(i__) = rNULL__
if(ar__1(68:75) .EQ. '') F0_665(i__) = rNULL__
if(ar__1(77:84) .EQ. '') e_F0_665(i__) = rNULL__
if(ar__1(86:92) .EQ. '') F0_79(i__) = rNULL__
if(ar__1(94:101) .EQ. '') e_F0_79(i__) = rNULL__
if(ar__1(103:109) .EQ. '') F0_915(i__) = rNULL__
if(ar__1(111:118) .EQ. '') e_F0_915(i__) = rNULL__
if(ar__1(120:126) .EQ. '') F0_96(i__) = rNULL__
if(ar__1(128:136) .EQ. '') e_F0_96(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ v_EDJ2009__1(i__),F0_36(i__),e_F0_36(i__),F0_44(i__),
+ e_F0_44(i__),F0_55(i__),e_F0_55(i__),F0_64(i__),e_F0_64(i__),
+ F0_665(i__),e_F0_665(i__),F0_79(i__),e_F0_79(i__),F0_915(i__),
+ e_F0_915(i__),F0_96(i__),e_F0_96(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table10.dat' ! 3.4-100um flux densities of YSOs in the c2d clouds
C Format for file interpretation
3 format(
+ I4,1X,I5,1X,F5.1,1X,I5,1X,I3,1X,F5.1,1X,F5.2,1X,I5,1X,F6.1,1X,
+ F6.1,1X,F4.1,1X,I6,1X,F6.1,1X,I6,1X,I5,1X,I7,1X,I5)
C Effective file loading
open(unit=1,status='old',file=
+'table10.dat')
write(6,*) '....Loading file: table10.dat'
do i__=1,1024
read(1,'(A104)')ar__2
read(ar__2,3)
+ v_EDJ2009__2(i__),F3_4(i__),e_F3_4(i__),F5_0(i__),e_F5_0(i__),
+ F6_7(i__),e_F6_7(i__),F12(i__),e_F12(i__),F14_3(i__),
+ e_F14_3(i__),F25(i__),e_F25(i__),F60(i__),e_F60(i__),
+ F100(i__),e_F100(i__)
if(ar__2(6:10) .EQ. '') F3_4(i__) = iNULL__
if(ar__2(12:16) .EQ. '') e_F3_4(i__) = rNULL__
if(ar__2(18:22) .EQ. '') F5_0(i__) = iNULL__
if(ar__2(24:26) .EQ. '') e_F5_0(i__) = iNULL__
if(ar__2(28:32) .EQ. '') F6_7(i__) = rNULL__
if(ar__2(34:38) .EQ. '') e_F6_7(i__) = rNULL__
if(ar__2(40:44) .EQ. '') F12(i__) = iNULL__
if(ar__2(46:51) .EQ. '') e_F12(i__) = rNULL__
if(ar__2(53:58) .EQ. '') F14_3(i__) = rNULL__
if(ar__2(60:63) .EQ. '') e_F14_3(i__) = rNULL__
if(ar__2(65:70) .EQ. '') F25(i__) = iNULL__
if(ar__2(72:77) .EQ. '') e_F25(i__) = rNULL__
if(ar__2(79:84) .EQ. '') F60(i__) = iNULL__
if(ar__2(86:90) .EQ. '') e_F60(i__) = iNULL__
if(ar__2(92:98) .EQ. '') F100(i__) = iNULL__
if(ar__2(100:104) .EQ. '') e_F100(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ v_EDJ2009__2(i__),F3_4(i__),e_F3_4(i__),F5_0(i__),e_F5_0(i__),
+ F6_7(i__),e_F6_7(i__),F12(i__),e_F12(i__),F14_3(i__),
+ e_F14_3(i__),F25(i__),e_F25(i__),F60(i__),e_F60(i__),
+ F100(i__),e_F100(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! 160-1300um flux densities of YSOs in the c2d clouds
C Format for file interpretation
4 format(
+ I4,1X,I5,1X,I5,1X,I5,1X,I4,1X,I5,1X,I5,1X,I4,1X,I4,1X,I4,1X,
+ F6.1,1X,I4,1X,I3,1X,I3,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,1024
read(1,'(A79)')ar__3
read(ar__3,4)
+ v_EDJ2009__3(i__),F160(i__),e_F160(i__),F350(i__),e_F350(i__),
+ F450(i__),e_F450(i__),F850(i__),e_F850(i__),F1100(i__),
+ e_F1100(i__),F1200(i__),e_F1200(i__),F1300(i__),e_F1300(i__)
if(ar__3(6:10) .EQ. '') F160(i__) = iNULL__
if(ar__3(12:16) .EQ. '') e_F160(i__) = iNULL__
if(ar__3(18:22) .EQ. '') F350(i__) = iNULL__
if(ar__3(24:27) .EQ. '') e_F350(i__) = iNULL__
if(ar__3(29:33) .EQ. '') F450(i__) = iNULL__
if(ar__3(35:39) .EQ. '') e_F450(i__) = iNULL__
if(ar__3(41:44) .EQ. '') F850(i__) = iNULL__
if(ar__3(46:49) .EQ. '') e_F850(i__) = iNULL__
if(ar__3(51:54) .EQ. '') F1100(i__) = iNULL__
if(ar__3(56:61) .EQ. '') e_F1100(i__) = rNULL__
if(ar__3(63:66) .EQ. '') F1200(i__) = iNULL__
if(ar__3(68:70) .EQ. '') e_F1200(i__) = iNULL__
if(ar__3(72:74) .EQ. '') F1300(i__) = iNULL__
if(ar__3(76:79) .EQ. '') e_F1300(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ v_EDJ2009__3(i__),F160(i__),e_F160(i__),F350(i__),e_F350(i__),
+ F450(i__),e_F450(i__),F850(i__),e_F850(i__),F1100(i__),
+ e_F1100(i__),F1200(i__),e_F1200(i__),F1300(i__),e_F1300(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table12.dat' ! Properties of YSOs in the c2d clouds
C Format for file interpretation
5 format(
+ I4,1X,A7,1X,A19,1X,A19,1X,F4.1,1X,F5.2,1X,A2,I4,1X,A2,F8.5,1X,
+ F5.2,1X,A2,I4,1X,A2,F8.4,1X,A1,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,1024
read(1,'(A112)')ar__4
read(ar__4,5)
+ v_EDJ2009__4(i__),Cloud(i__),SSTc2d(i__),Class(i__),AV(i__),
+ alpha(i__),l_Tbol(i__),Tbol(i__),l_Lbol(i__),Lbol(i__),
+ alphac(i__),l_Tbolc(i__),Tbolc(i__),l_Lbolc(i__),Lbolc(i__),
+ Env(i__),Qual(i__)
if(ar__4(59:63) .EQ. '') alpha(i__) = rNULL__
if(ar__4(67:70) .EQ. '') Tbol(i__) = iNULL__
if(ar__4(74:81) .EQ. '') Lbol(i__) = rNULL__
if(ar__4(83:87) .EQ. '') alphac(i__) = rNULL__
if(ar__4(91:94) .EQ. '') Tbolc(i__) = iNULL__
if(ar__4(98:105) .EQ. '') Lbolc(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ v_EDJ2009__4(i__),Cloud(i__),SSTc2d(i__),Class(i__),AV(i__),
+ alpha(i__),l_Tbol(i__),Tbol(i__),l_Lbol(i__),Lbol(i__),
+ alphac(i__),l_Tbolc(i__),Tbolc(i__),l_Lbolc(i__),Lbolc(i__),
+ Env(i__),Qual(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end