outactfz.f

Go to the documentation of this file.

      subroutine outactfz(tmf1,tmf2)
c
c     subroutine outactfz outputs the actinic flux at all the levels
c     of the atmosphere
c***********************************************************************
c.....include the common blocks
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
      real*8 tma(nsz),tmb(nsz),tmc(nsz),tmfd(nsz),tmfu(nsz),tmf1(nsz,nlyr),
     1       tmf2(nlyr)
c***********************************************************************
c
      write(*,*)'welcome to outactfz'
      write(*,*)'bfr1(35)',bfr1(35)
      lsza=bfr1(35)+0.001
      do i=1,lsza,1
         m=msza(i)
         tma(i)=the0in(m)
         tmb(i)=fdirc(m)/pi
         tmc(i)=sbarz(m)
         tmfd(i)=fdown(m)/pi
         tmfu(i)=fup(m)/pi
         do j=1,nolyr+1
            tmf1(i,j)=acfg1(m,j)
            tmf2(j)=acfg2(j)
         enddo
      enddo
c
      ff0=1
      write(30,502)wvlth*1.0e+4
      write(30,503)ff0
      write(30,504)(tma(i),i=1,lsza)
      write(30,506)(tmb(i),i=1,lsza)
      write(30,508)(tmfd(i),i=1,lsza)
      write(30,510)(tmc(i),i=1,lsza)
      write(30,511)
      do k=1,nolyr+1
      write(30,512)k,ht(k),pl(k),taupl(k),tmf2(k),(tmf1(i,k),i=1,lsza)
      enddo
c
502   format('wavelength(um)',t20,1pe11.3)
503   format('solar flux(top)',t20,1pe11.3)
504   format('sza',t22,25(f4.1,7x))
506   format('f(dir)',t20,1p25e11.3)
508   format('f(diff)',t20,1p25e11.3)
510   format('sbar',t20,1p25e11.3)
511   format(t1,'lvl',t7,'height',t18,'pressure',t29,
     1       'tot. opt.',t40,'act_flx',t51,'act_flx'/
     2       t1,'no.',t8,'(km)',t19,'(atm.)',t29,'thickness',
     3       t40,'(ground=1)',t51,'(r=0)')
512   format(i3,1x,1p29e11.3)
      return
      end
c************************************************************************