Team:Edinburgh/modelling(reallifemodelling)

if (index==1){
Bact[xbact][ybact]=1;
}

if (index==2){
Bact[xbact-1][ybact]=1;
xbact=xbact-1;
}

if (index==3){
Bact[xbact+1][ybact]=1;
xbact=xbact+1;
}

if (index==4){
Bact[xbact][ybact-1]=1;
ybact=ybact-1;
}

if (index==5){
Bact[xbact][ybact+1]=1;
ybact=ybact+1;
}

display ('***********************************************************')
%Finite difference model for the diffusion of TNT through soil

xsize=100;
ysize=100;
MaxSteps=8760;
TimeStep=0.01;
SpaceStep=1;
diffusivityC=1.18*(10-6)*3600/100/100;
%calculates s and s2 for TNT based on the diffusion coefficient of TNT
sC=(diffusivityC*TimeStep/SpaceStep/SpaceStep);
s2C=(1-4*sC);
diffusivityN=1.18*(10-6)*3600/100/100;
%Calculates s and s2 for Nitrites based on the diffusion coefficient of
%Nitrites
sN=(diffusivityN*TimeStep/SpaceStep/SpaceStep);
s2N=(1-4*sC);


% Sets up 2 matricies for the calculation
C=zeros (xsize,ysize);
CNew = zeros (xsize,ysize);
N= zeros (xsize,ysize);
NNew= zeros (xsize,ysize);


%Set up the concentration and position of the landmine
MineX=5;
MineY=5;
C(MineX,MineY)=107800;
N (MineX,MineY)=1;

React=0.99;
Produced=1/React;

for k=1:MaxSteps,
for i=2:xsize-1,
for j=2 :ysize-1,
%Calculates concentrations due to finite difference equation
CNew(i,j)= React*(sC*(C((i-1),j)+ C((i+1),j)+C(i,(j-1))+C(i,(j+1)))+s2C*C(i,j));
NNew(i,j)=Produced*(sN*(N((i-1),j)+ N((i+1),j)+N(i,(j-1))+N(i,(j+1)))+s2N*N(i,j));

end
end

C=CNew;
N=NNew;
display (C)
display(N)
>br /> surfc (C)

refreshdata

end

import java.applet.*;
import java.awt.*;
import java.awt.event.*;
import java.util.Random;

import javax.swing.JOptionPane;

public class ChemotaxisModel extends Applet implements Runnable {

/**
*
*/
private static final long serialVersionUID = 1L;

int h,j,k1,score,x1,y1,xmove,ymove,cloud,xsize1,ysize1,mx,my,cx,cy, c1;
Thread t;
Graphics offscreen;
Image workspace;
Color vdkgreen;
int MineX, MineY;
String UserMineX, UserMineY;{

UserMineX=JOptionPane.showInputDialog( "Enter x-location for the landmine" );
UserMineY=JOptionPane.showInputDialog("Enter y-location for the landmine");

MineX=Integer.parseInt(UserMineX);
MineY=Integer.parseInt(UserMineY); }

public ChemotaxisModel () {}

public void init() {
xsize1=getSize().width;
ysize1=getSize().height;
workspace=createImage(xsize1,ysize1);
offscreen=workspace.getGraphics();
x1=100;y1=100;xmove=0;ymove=0;

vdkgreen=new Color(0,102,0);

}

public void start (){
if (t==null){
t=new Thread (this);
t.start();

}
}
public void run () {}

public void stop (){
if (t != null){
t.stop ();
t=null;
}
}

public void paint (Graphics screen){

offscreen.setColor(vdkgreen);
offscreen.fillRect(0,0,xsize1,ysize1);

drawBacteria (offscreen);
drawLandmine (offscreen);
screen.drawImage(workspace,0,0,this);
}

public void update(Graphics screen){
paint(screen);
}

void drawBacteria (Graphics g){

int xsize, ysize, c1, c2, c3, indexx,indexy,xbact,ybact;
long c4,c5,c6,c7,c8;
double C[][],N[][],Bact[][],BactSurr[];
double NewC[][],NewN[][], NewBact[][];
double TimeStep, SpaceStep, Alpha,s, s2, k,x,maxC,max,index;

double TNTLoss,N2Prod,TotalTNT, TotalNitrite;

int MaxSteps=1000;

xsize=10;
ysize=10;
C=new double[xsize][ysize];
NewC=new double[xsize][ysize];
N=new double [xsize][ysize];
NewN=new double [xsize][ysize];
Bact=new double [xsize][ysize];
NewBact=new double [xsize][ysize];
BactSurr=new double[6];
TimeStep=0.01;
SpaceStep=1;
Alpha=0.2;
s=(Alpha*TimeStep/SpaceStep/SpaceStep);
s2=(1-4*s);

for (c1=0; c1 for (c2=0; c2 C[c1][c2]=0;
N[c1][c2]=0;
Bact [c1][c2]=0;
}
}
C[MineX][MineY]=1000;
N[MineX][MineY]=10;
k=0.06;
TNTLoss=1;
N2Prod=1/TNTLoss;
TotalTNT=0;
TotalNitrite=0;
x=0;
index=0;
indexx=0;
indexy=0;
max=0;

max=0;
Random generator = new Random();
Random generator2 = new Random();
xbact = generator.nextInt( xsize-5 )+2;
ybact = generator2.nextInt( ysize-5 )+2;
Bact[xbact][ybact]=1;
for (c3=0;c3 for (c1=1; c1< (xsize-1); c1++) {
for (c2=1; c2< (ysize-1); c2++) {
NewC[c1][c2]=TNTLoss*(s*(C[c1-1][c2]+C[c1+1][c2]+C[c1][c2-1]+C[c1][c2+1])+s2*C[c1][c2]);
NewN[c1][c2]=N2Prod*(s*(N[c1-1][c2]+N[c1+1][c2]+N[c1][c2-1]+N[c1][c2+1])+s2*N[c1][c2]);
TotalTNT +=(C[c1][c2]);
TotalNitrite +=N[c1][c2];
}
}

System.out.println("TNT");
for (c1=1; c1< (xsize-1); c1++) {
for (c2=1; c2< (ysize-1); c2++) {

c4=Math.round(C[c1][c2]);
System.out.print(c4);
System.out.print("-");
C[c1][c2]=NewC[c1][c2];
}
System.out.println(".");
}
c6=Math.round(TotalTNT);
System.out.println("*");
System.out.println(c6);
System.out.println("*");
TotalTNT=0;
}

for (c3=94;c3
System.out.println("Bacteria");
for (c1=1; c1< (xsize-1); c1++) {
for (c2=1; c2< (ysize-1); c2++) {

c8=Math.round(Bact[c1][c2]);
System.out.print(c8);
System.out.print("-");
Bact[c1][c2]=NewBact[c1][c2];
}
System.out.println(".");
}

BactSurr[1]=C[xbact][ybact];
BactSurr[2]=C[xbact-1][ybact];
BactSurr[3]=C[xbact+1][ybact];
BactSurr[4]=C[xbact][ybact-1];
BactSurr[5]=C[xbact][ybact+1];

System.out.println(BactSurr[1]);
System.out.println(BactSurr[2]);
System.out.println(BactSurr[3]);
System.out.println(BactSurr[4]);
System.out.println(BactSurr[5]);

max=0;
index=0;

for (c1=1; c1<6;c1++){
if (max max=BactSurr[c1];
index=c1;
System.out.println(index);
}
} System.out.println(index);

if (index==1){
Bact[xbact][ybact]=1;
}

if (index==2){
Bact[xbact-1][ybact]=1;
xbact=xbact-1;
}

if (index==3){
Bact[xbact+1][ybact]=1;
xbact=xbact+1;
}

if (index==4){
Bact[xbact][ybact-1]=1;
ybact=ybact-1;
}

if (index==5){
Bact[xbact][ybact+1]=1;
ybact=ybact+1;
}

g.setColor(Color.yellow);
g.fillOval(xbact*20, ybact*20, 5, 5);
}
}

void drawLandmine (Graphics g){
g.setColor (Color.black);
g.fillOval(MineX*20,MineY*20, 10,10);
}
}