CCC CONVECTION AND DIFFUSION OF A TEMPERATURE PULSE   
CCC     PROGRAM ICS
        SUBROUTINE ABQMAIN
C
        OPEN(UNIT=17,STATUS='UNKNOWN',FILE='CONVECT.TEM')
        DH=0.025
        X=0.0
        Y=0.0
        INCC=0
        INC=200
        DO 30 IV=1,21
        DO 50 II=1,81,2
        T=EXP(-((X-0.175)**2+(Y-0.175)**2)/0.00125)
        IF(T.LT.1.E-10)T=0.0D0
        NODE=II+INCC
        IF(T.GT.1.E-8)WRITE(17,888)NODE,T
        X=X+DH
 50     CONTINUE
        X=0.0
        Y=Y+DH
        INCC=INCC+INC
 30     CONTINUE
 888    FORMAT(I5,2H, ,F11.9)
        STOP
        END