%flow

clear all;
close all;

cnn_setenv;

mCNN.TemGroup='temlib_diff';


Input=loadCNN('flow_input.bmp');
Gates=loadCNN('gates.bmp');
mCNN.UseMask=1;
mCNN.MASK=Gates;


InitState=zeros(size(Input));
State=1.*InitState;
for repeat=1:500
    mCNN.INPUT1=1*Input;
    mCNN.STATE=1*State;
    mCNN.Boundary=2;
    mCNN.TimeStep=0.1;
    mCNN.IterNum=5;
    loadtem('FLOW')
    runtem;  
    State=1.*mCNN.STATE;
    surf(mCNN.OUTPUT);   
    drawnow
end