pro task25,iout

; itask=25 Look at the PLT files you created, do the different azimuths give different radiances?

  printf,iout,'  '
  printf,iout,' task25: Look at the PLT files you created, do the different azimuths give different radiances?'

; ******************************
; Specify the ascii file listing
   flist='/ur/massie/libera/json/modtran6_sza15cld/zplt.list'
    ndat=1317
   nj=4
   jval=intarr(nj)
   jval=[0,1,2,3]

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

; ******************************
     jfiles=nj
    spectra=fltarr(ndat,jfiles)
    wavenm=fltarr(ndat,jfiles)

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

     jj=jval(j)
     print,j,jj

     idat=50
     openr,idat,files(jj)

     i1=0
     i2=ndat-1
     iskip=-1
    for i=i2,i1,iskip do begin
     readf,idat,a1,a2
     wavenm(i,j)=a1
; spectra values are in mW ster-1 cm-2 nm-1
     spectra(i,j)=a2
    endfor

    close,idat

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

   endfor

; ******************************
     printf,iout,'  '
     printf,iout,'  task25: spectra values are in mW ster-1 cm-2 nm-1'
     printf,iout,'  task25: i,wavenm(i,0),vec'
     vec=fltarr(nj)
    for i=i1,i2 do begin
     for j=0,nj-1 do begin
      vec(j)=spectra(i,j)
     endfor
     printf,iout,format='(2x,i4,1(1x,f10.4),2x,4(1x,e10.3))',$
       i,wavenm(i,0),vec
    endfor

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