“…, file=filename3) write(7,*) header write(7,1) numeos, zbar3, abar3, rho3, n3 write(7,2) (array3(i), i=1,n3) close (7) 1 format(1x,i5,4x,1p3e15.8,3x,i5) 2 format(1p5e15. 8) stop end c ********************************************************* c EXTERNAL FUNCTION XMAXVAL function xmaxval(a,na) real*4 a(100) xmaxval=0. do i=1,na if(a(i).gt.xmaxval) xmaxval=a(i) enddo return end c ********************************************************* c UNPACK a linear array into d(nd), t(nt), p1(i,j), p2(i,j) subroutine unpack(a,d,t,p1,p2,nd,nt,incr) real*4 a(10000), d(100), t(100), p1(100,100), p2(100,100) c (0.01 is to avoid roundoff errors in integers) nd = aint(a(1)+0.01) nt = aint(a(2)+0.01) incr = 2 do i=1,nd incr=incr+1 d(i) = a(incr) enddo do i=1,nt incr=incr+1 t(i) = a(incr) enddo do j=1,nt do i=1,nd incr=incr+1 p1(i,j) = a(incr) enddo enddo do j=1,nt do i=1,nd incr=incr+1 p2(i,j) = a(incr) enddo enddo return end c ********************************************************* c PACK 2d arrays d(nd), t(nt), p1(i,j), p2(i,j) into a subroutine pack(a,d,t,p1,p2,nd,nt,incr) real*4 a(10000), d(100), t(100), p1(100,100), p2(100,100) a(1) = nd a(2) = nt SiO 2 aerogel exposed to 4-ns laser irradiation UCSD-CER-05-01 59 incr = 2 do i=1,nd incr=incr+1 a(incr) = d(i) enddo do i=1,nt incr=incr+1 a(incr) = t(i) enddo do j=1,nt do i=1,nd incr=incr+1 a(incr) = p1(i,j) enddo enddo do j=1,nt do i=1,nd incr=incr+1 a(incr) = p2(i,j) enddo enddo return end…”