c subroutine sumry(jmuz) c************************************************************************ c subroutine sumry prints out sb, t, i0 pol etc c************************************************************************ c implicit real*8(a-h,o-z),integer*4(i-n) c include "contrl.inc" include "out.inc" include "totals.inc" include "emm.inc" include "atmos.inc" include "in.inc" include "buff1.inc" include "buff2.inc" include "buff3.inc" include "buff4.inc" include "depolt.inc" c maxpp=itmax+2 c c write(9,50) do i=1,maxpp do j=1,imu do k=1,iazmth c write(9,100)i,thta(j),azmth(k),e0za(i,j),e1za(i,j,k), c 1 e2za(i,j,k),ttz(i,j),sbz(i),ggz(i) enddo enddo enddo c write(9,150)thnot(jmuz),wavlen write(9,200) do i=1,nalb c do j=1,imu c do k=1,iazmth do k=1,iazmth do j=1,imu write(9,225)thta(j),azmth(k),eittl(j,k,i), 1 eittr(j,k,i),eitotz(j,k,i),polz(j,k,i),alb(i) enddo enddo enddo c return c 50 format(t2,'itr',t10,'the',t15,'phi',t26,'ei0',t37,'ei1', 1 t48,'ei2',t59,'t ',t67,'sb',t77,'gg') 100 format(i5,1x,f7.1,1x,f7.1,1p4e11.3,0pf8.5,1pe11.3) 150 format('solar zenith angle=',f6.1,' wavelength=',f7.1) 200 format(t4,'the',t12,'phi',t24,'eil',t35,'eir',t44,'eitot', 1 t55,'pol',t62,'alb') 225 format(f7.1,1x,f7.1,1x,1p3e11.3,0pf11.3,0pf7.3) c end