pro task30,iout

; itask=30 Look at the FLX files you created, do the different azimuths give different integrated fluxes?

   printf,iout,'  '
   printf,iout,'  task30: Look at the FLX files you created, do the different azimuths give different integrated fluxes?'

; ******************************
; Specify the ascii file listing
   iset=1
  if (iset eq 1) then begin
   flist='/ur/massie/libera/json/modtran6_sza15cld/zflx.list'
   nj=4
   jval=intarr(nj)
   jval=[0,1,2,3]

 ALTITUDES:          0.00000 KM                        100.00000 KM
          ----------------------------------  ----------------------------------
  WAVLEN      UPWARD    DOWNWARD      DIRECT      UPWARD    DOWNWARD      DIRECT
   (NM)      DIFFUSE     DIFFUSE       SOLAR     DIFFUSE     DIFFUSE       SOLAR

; W CM-2: 1.05865E-03 3.90264E-02 1.98837E-04 1.82129E-02 0.00000E+00 6.01955E-02     last line of the file

; For the same solar zenith angle and same sensor zenith angle, you get the same  upward diffuse 2 cm-2 flux
; of 1.82129e-2
;      0       0
; read in from /ur/massie/libera/json/modtran6_sza15cld/cld300700solarzen15sensorzen100reflaz0.flx
;      1       1
; read in from /ur/massie/libera/json/modtran6_sza15cld/cld300700solarzen15sensorzen100reflaz20.flx
;      2       2
; read in from /ur/massie/libera/json/modtran6_sza15cld/cld300700solarzen15sensorzen100reflaz40.flx
;      3       3
; read in from /ur/massie/libera/json/modtran6_sza15cld/cld300700solarzen15sensorzen100reflaz60.flx

   endif

; ******************************
; Read in the file pathnames
    noprr=1
   readflist,noprr,iout,flist,nfiles,files
   if (noprr eq 1) then begin
;   stop
   endif

; ****************
   nmax=8000
  words=strarr(1)
  laststr=strarr(1)
  datastr=strarr(nmax,nfiles)

; ******************************
; Loop over the nj files
   for j=0,nj-1 do begin

     jj=jval(j)
     print,j,jj

; ****************
   idat=30
  openr,idat,files(jj)

; ****************
; Read in all lines of the file
   nlines=-1
  for i=0,nmax-1 do begin
    icheck=EOF(idat)
   if (icheck eq 0) then begin
     nlines=nlines+1
    readf,idat,words
    datastr(nlines,j)=words
    laststr=words
   endif
   if (icheck eq 1) then begin
    goto,jump1
   endif
  endfor
  jump1: i=0
  nlines=nlines+1

  close,idat

; ****************
    print,'  read in from ',files(jj)
    printf,iout,'  '
    printf,iout,'  task25: read in from ',files(jj)
    printf,iout,'  task25: laststr ',laststr

   endfor

; ******************************
   return
   end