CCC CONVECTION AND DIFFUSION OF A TEMPERATURE PULSE   3-1-15-3
CCC     PROGRAM THEOR
        SUBROUTINE ABQMAIN
C
        OPEN(UNIT=7,STATUS='UNKNOWN',FILE='3011503.DAT')
        XK=3.90625E-04
        XV=.25
        H=.03125
        TIME=2.0
        DO 25 L=1,3
        XX=0.
        IF(L.EQ.1)THEN
           DO 50 IIC=1,65
           T=EXP(-(XX-XV)**2./(4.*XK))
           WRITE(7,77)XX,T
           XX=XX+H
 50        CONTINUE
        ELSE
           DO 75 ITT=1,65
           P1=1./(SQRT(1.+TIME))
           PNUM=(XX-XV*(1.+TIME))**2.
           PDEN=4.*XK*(1.+TIME)
           T=P1*EXP(-(PNUM/PDEN))
           WRITE(7,77)XX,T
           XX=XX+H
 75        CONTINUE
           TIME=4.0
        ENDIF
 25     CONTINUE
 77     FORMAT(5X,F6.4,5X,F11.8)
        STOP
        END