%flow

clear all;
close all;

cnn_setenv;

mCNN.TemGroup='temlib_diff';

Input=loadCNN('flow_input.bmp');


InitState=zeros(size(Input));
State=1.*InitState;
for repeat=1:100
    mCNN.INPUT1=1*Input;
    mCNN.STATE=1*State;
    mCNN.Boundary=2;
    mCNN.TimeStep=0.001;
    mCNN.IterNum=1;

    loadtem('HEATING')
    runtem;
    
    State=1.*mCNN.STATE;
    surf(mCNN.STATE);
    
    drawnow
end