c cformat.f c c
      program format
c
c    Examples of various uses of the Format Statement, and a small
c    warning about precision of constants.
c
c    John Mahaffy  2/26/96
c
      integer n
      parameter (n=5)
      real r,cv,t(n),p(n),rho(n),u(n),tc(n),tf(n)
c
      double precision rd,cvd,td,pd,rhod
      double precision rd1
c
      parameter (r=287.0478,cv=1004.832-r)
c
c    Watch the output to see what happens when I casually move
c    single precision numbers into double precision parameters
c    (or constants).  The computer is pretty causual about what
c    it fills into the extra digits of precision.
c
      parameter (rd1=r)
c
c    Here I carefully define constants as DOUBLE PRECISION
c
      parameter (rd=287.0478d0,cvd=1.004832d3-r)
c
c    In Fortran90 I could have used a different notation in conjuction
c    with KIND.  Comment out old code and uncomment the following to try.
c
c     integer, parameter :: r8=selected_real_kind(15,30)
c     real (r8) rd,cvd,td,pd,rhod,ud
c     parameter (rd=287.0478_r8, cvd=1.004832e3_r8)
c
      data p,t /n*1.e5,n*300./
      data pd,td/1.d5,3.d2/
      open(11,file='format.out')
      print *, ' Casual  Definition of Double Precision Gas Constant ='
     &  ,rd1
      print *, ' Careful Definition of Double Precision Gas Constant ='
     &  ,rd
c
      print *
c
      print 2222
 2222 format ('Look what I''ve got to do to include a single quote')
c
      print *
c
      rho(1)= p(1)/(r*t(1))
      rhod=pd/(rd*td)
c
c    Start by writing a header
c
      write(*,1000)
 1000 format(14x,' Pressure',13x,'Temperature',17x,'Density')
      write(*,2000) p(1),t(1),rho(1)
      write(*,2001) pd,td,rhod
      write(*,2002) pd,td,rhod
c
c  Note how numbers in the following are right justified in the output
c
 2000 format(1p,3(1x,e23.7))
 2001 format(1p,3(1x,d23.15))
c
c   You can get away with "e" edit descriptor on DOUBLE PRECISION
c
 2002 format(1p,3(1x,e23.15))
      tc(1)=t(1)-273.16
      tf(1)=32.+1.8*tc(1)
      u(1)=cv*t(1)
      rho(1)=p(1)/(r*t(1))
      do 100 i=2,n
         t(i)=t(i-1)+100.
         tc(i)=t(i)-273.16
         tf=32.+1.8*tc(i)
         p(i)=p(i-1)+1.e4
         u(i)=cv*t(i)
 100     rho(i)=p(i)/(r*t(i))
      write(11,2020)    2005
      write(11,2004)
 2004 format(//,63x,'Internal',/,
     &  1x,'Cell',3(5x,'Temp '), 7x,'P',8x,'Density',6x,'Energy',
     &  /,2x,'Num',6x,'(K)',7x,'(C)',7x,'(F)',6x,'(Pa)',6x,
     &  '(kg/m**3)',5x,'(J/kg)' )
      write(11,2005) (i,t(i),tc(i),tf(i),p(i),rho(i),u(i),i=1,n)
 2005 format(1x,i4,0p,3f10.1,1p,3e12.3)
c
c   Watch What gets reused in the 2006 FORMAT
c
      write(11,2020)    2006
      write(11,2004)
      write(11,2006) (i,t(i),tc(i),tf(i),p(i),rho(i),u(i),i=1,n)
 2006 format(71('-'),/,(1x,i4,0p,3f10.1,1p,3e12.3))
c
c   Watch What gets reused in the 2007 FORMAT, containing an extra
c   bounding parentheses, and a repeat factor on the main interior
c   block of the FORMAT
c
      write(11,2020)    2007
      write(11,2004)
      write(11,2007) (i,t(i),tc(i),tf(i),p(i),rho(i),u(i),i=1,n)
 2007 format((71('-'),/,1(1x,i4,0p,3f10.1,1p,3e12.3)))
c
c   Next I get fancy with Repeat Operators and some vertical lines
c
      write(11,2020)    2008
      write(11,2004)
      write(11,2008) (i,t(i),tc(i),tf(i),p(i),rho(i),u(i),i=1,n)
 2008 format((1x,i4,0p,3(2x,'|',2x,f5.1),1p, 3(1x,'|',1x,e9.3)))
c
      write(11,2020)    2009
      write(11,2004)
c
c   Look how the seemingly harmless drop of one pair of parentheses
c   can totally wreck the pattern of reuse in the FORMAT
c
      write(11,2009) (i,t(i),tc(i),tf(i),p(i),rho(i),u(i),i=1,n)
 2009 format(1x,i4,0p,3(2x,'|',2x,f5.1),1p, 3(1x,'|',1x,e9.3))
c
c   Try replacing the above failed format with one that produces a
c   table containing both vertical and horizontal bounding lines.
c
 2020 format(//,'  Format Number ',i4)
      stop
      end
c
c c