ew 1027310400000 ew 7/22/02 8/14/02 10000 0 Value_MIN double -1135038327 0.0 0 Value_MAX double -1135038327 102400.0 0 Value double -1135038327 re 0 MULTIPLIER double 0 100.0 1014872400000 /home/wender/eclpssmodels/laplace/sim2/Temp.esv Temp 1014923858148 0 Execution Group 1 ew 1014872400000 ew 2/28/02 7/22/02 7/22/02 7/22/02 7/22/02 7/22/02 7/22/02 "silly" simultation to show heat diffusion. 0 Value_MIN double -1135038327 0.0 0 Value_MAX double -1135038327 102400.0 0 Value double -1135038327 re 0 MULTIPLIER double 0 100.0 1014872400000 /home/wender/eclpssmodels/laplace/sim2/Temp.esv Temp 1014923858148 /home/tsmart/public_html/classes/eclpssSaves/HeatDiffusion/Components/CoolDown.cpt /** *<UL><BR> *<LI><BIG>This is an <font color = "mediumblue"><i> ECLPSS Component</i></font></BIG></LI><BR> *<LI><b>Author: </b>ew</LI><BR> *<LI><b>Date Created: </b>2/28/02</LI><BR> *<LI><b>Date Saved: </b>7/22/02</LI><BR> *<LI><b>Version: </b>7/22/02</LI><BR><br> *<LI><b>Purpose: </b>7/22/02</LI><BR> *<LI><b>Algorithm Citation: </b>7/22/02</LI><BR> *<LI><b>Modification History: </b>7/22/02</LI><BR> *<LI><b>Assumptions: </b>7/22/02</LI><BR> *<LI><b>Limitations: </b></LI><BR> *<LI><b>Known Problems: </b>"silly" simultation to show heat diffusion.</LI><BR> *<LI><b>Suggested Improvements: </b></LI><BR> *<LI><b>Testing Statistics: </b></LI><BR> *<LI><b>Models used in: </b></LI><BR> *<LI><b>General Comments: </b></LI><BR> */ // make an array of hot cell Points (for current and previous time // frames. Point [] hotCellsCurrent = new Point [2]; Point [] hotCellsPrevious = new Point [2]; static int xLeft, xRight, yLeft, yRight; // ECLPSSgrid is a "generic" grid - the framework replaces // it with the correct name at model run time. void getBounds(ECLPSSGrid grid, Point current) { xLeft = grid.getAbsXLB()-1; xRight = grid.getAbsXUB()-1; yLeft = grid.getAbsYLB()-1; yRight = grid.getAbsYUB()-1; } // relative offsets from current static final Point N = new Point( 0, -1, BACK1 ); static final Point S = new Point( 0, 1, BACK1 ); static final Point E = new Point( 1, 0, BACK1 ); static final Point W = new Point( -1, 0, BACK1 ); /** checks to see if the 2-D grid cell at absolute address (x,y) is a corner cell. @param x - absolute address of x @param y - absolute address of y @return true if a corner cell, false otherwise */ public boolean isCorner(int x, int y) { return ( (x==xLeft && y==yLeft) || (x==xRight && y==yLeft) || (x==xLeft && y==yRight) || (x==xRight && y==yRight) ); } /** checks to see if the 2-D grid cell at absolute address (x,y) is on a side of the grid. @param x - absolute address of x @param y - absolute address of y @return true if (x,y) the side of grid, false otherwise */ public boolean isSide(int x, int y) { return (!isCorner(x,y) && (x==xLeft || x==xRight || y==yLeft || y==yRight) ); } CoolDown 2 // initialize the hot cells to be in the upper left and lower right int xlbindex = grid.getAbsXLB()-1; int ylbindex = grid.getAbsYLB()-1; int yub= grid.getAbsYUB(); int xub= grid.getAbsXUB(); // iterate over the entire grid for (int i=xlbindex; i< xub; i++) { for (int j = ylbindex; j<yub; j++) { // the designated cells hot // and make the rest of the grid cold if (i == xlbindex && j == ylbindex) { // initialize grid grid.writeAbsTemp(current,new Point(i,j,NOW)).add(new Temp(Temp.Value_MAX)); // initialize arrays hotCellsCurrent[0] = new Point(i,j,NOW); hotCellsPrevious[0] = new Point(i,j,BACK1); } if (i == xub-1 && j == yub-1 ) { // initialize grid grid.writeAbsTemp(current,new Point(i,j,0)).add(new Temp(Temp.Value_MAX)); // initialize arrays hotCellsCurrent[1] = new Point(i,j,NOW); hotCellsPrevious[1] = new Point(i,j,BACK1); } else { // just a regular cell - initialize grid grid.writeAbsTemp(current,new Point(i,j,NOW)).add(new Temp(Temp.Value_MIN)); } } } //write out the hot cell addresses for (int i=0; i< hotCellsCurrent.length; i++) { System.out.println("hotCellsCurrent["+i+"] = "+hotCellsCurrent[i]); System.out.println("hotCellsPrevious["+i+"] = "+hotCellsPrevious[i]); } int x = current.getX(); int y = current.getY(); int it = model.getTimeIteration(); // print the iteration number to the console modulo displayRate //if (x==0 && y==0 && ((it-1) % model.displayRate == 0) ) //System.out.println("Iteration "+it); double temp = 0.0; double north = 0.0; double south = 0.0; double east = 0.0; double west = 0.0; try { south = grid.readTemp(current,S).getValue(); } catch (Exception e) {} try { north = grid.readTemp(current,N).getValue(); } catch (Exception e) {} try { east = grid.readTemp(current,E).getValue(); } catch (Exception e) {} try { west = grid.readTemp(current,W).getValue(); } catch (Exception e) {} getBounds(grid, current); // since all points are either correct (within the grid boundaries) // or 0.0 (all points are initialized to 0, and any out-of-bound gets // are caught) //we can add them all up, and then just divide by the appropriate amount. //I'm multiplying (where there is no roundoff error) since it is faster than divide if (isCorner(x,y)) { grid.writeTemp(current,CURRENT).add(new Temp( (north+south+east+west)*0.5)); } else if (isSide(x,y)) { grid.writeTemp(current,CURRENT).add(new Temp( (north+south+east+west)/3.0)); } else { grid.writeTemp(current,CURRENT).add(new Temp( (north+south+east+west)*0.25)); } for (int i=0; i< hotCellsCurrent.length; i++) { if (x==hotCellsCurrent[i].getX() && y == hotCellsCurrent[i].getY()) { // get previous temperature temp = grid.readAbsTemp(hotCellsPrevious[i]).getValue(); //System.out.println("read temp "+temp+" from "+hotCellsPrevious[i]); //System.out.println(model.getTimeIteration()); // write it here and now grid.writeAbsTemp(hotCellsCurrent[i]).setValue(temp); //System.out.println("write temp to "+hotCellsCurrent[i]); } } 1027349505733 ew 1014872400000 ew 2/28/02 7/23/02 7/22/02 2/28/02 Display the grid temperatures as ranging from "red" (hottest) to "blue" (coldest). 0 Value_MIN double -1135038327 0.0 0 Value_MAX double -1135038327 102400.0 0 Value double -1135038327 re 0 MULTIPLIER double 0 100.0 1014872400000 /home/wender/eclpssmodels/laplace/sim2/Temp.esv Temp 1014923858148 /home/tsmart/public_html/classes/eclpssSaves/HeatDiffusion/Components/Display.cpt true 600 /** *<UL><BR> *<LI><BIG>This is an <font color = "mediumblue"><i> ECLPSS Component</i></font></BIG></LI><BR> *<LI><b>Author: </b>ew</LI><BR> *<LI><b>Date Created: </b>2/28/02</LI><BR> *<LI><b>Date Saved: </b>7/23/02</LI><BR> *<LI><b>Version: </b>7/22/02</LI><BR><br> *<LI><b>Purpose: </b>2/28/02</LI><BR> *<LI><b>Algorithm Citation: </b></LI><BR> *<LI><b>Modification History: </b>Display the grid temperatures as ranging from "red" (hottest) to "blue" (coldest).</LI><BR> *<LI><b>Assumptions: </b></LI><BR> *<LI><b>Limitations: </b></LI><BR> *<LI><b>Known Problems: </b></LI><BR> *<LI><b>Suggested Improvements: </b></LI><BR> *<LI><b>Testing Statistics: </b></LI><BR> *<LI><b>Models used in: </b></LI><BR> *<LI><b>General Comments: </b></LI><BR> */ public void paint(Graphics g) { //java.awt.Container con = getContentPane(); g.drawImage(image,15,30,Color.blue,this); } JPanel canvas = null; Graphics2D offScreen = null; int frameWidth = 500; int frameHeight = 500; java.awt.image.BufferedImage image = new java.awt.image.BufferedImage (frameWidth, frameHeight, java.awt.image.BufferedImage.TYPE_INT_RGB); Display 2 int x = current.getX(); int y = current.getY(); if( ((model.getTimeIteration()-1) % model.displayRate) == 0 ) { if( x == grid.getXUB() - 1 && y == grid.getYUB() - 1) { // at the end of the iteration if(offScreen == null) { offScreen = image.createGraphics(); //canvas = new JPanel(); //getContentPane().add(canvas); //validate(); } int xUB = grid.getAbsXUB(); int yUB = grid.getAbsYUB(); int blockWidth = frameWidth / xUB; int blockHeight = frameHeight / yUB; int temp = 0; int rgb = 0; ChangeablePoint p = new ChangeablePoint(0,0,-1); boolean raised = true; //if(model.getTimeIteration() % 2 == 0){ //raised = !raised; //offScreen.rotate(Math.PI/2,250,250); //} // if(model.getTimeIteration() % 2 == 0){ // raised = !raised; // for(int xx=0; xx<=xUB; xx++) { // p.setX(xUB-xx); // for(int yy=0; yy<=yUB; yy++) { // p.setY(yUB-yy-1); // try{ // temp = (int) grid.readAbsTemp(current,p).getValue(); // } // catch(Exception e){} // temp /= Temp.MULTIPLIER; // // hold blue at 255, increase amount of green // if(temp < 255) { // rgb = (temp << 8) + 255; // } // // hold green at 255, decrease amount of blue // else if(temp < 511) { // rgb = (255 << 8) + (255 - (temp - 255)); // } // // hold green at 255 increase amount of red // else if(temp < 767) { // rgb = ((temp - 511)<<16) + (255 << 8); // } // // hold red at 255 decrease amount of green // else if(temp < 1023) { // rgb = (255 << 16) + ((255 - (temp - 767)) << 8); // } // else { // rgb = 0xff0000; // } // offScreen.setColor(new java.awt.Color(rgb)); // // //offScreen.fill3DRect(xx*blockWidth, yy*blockWidth, // // blockWidth, blockHeight, true); // //boolean raised = (xx%2 == 0)&&(yy%2 == 0); //cool display // raised = !raised; // offScreen.fill3DRect((xUB-xx)*blockWidth,(yy)*blockWidth, // blockWidth, blockHeight, raised); // } // } // } // else{ for(int xx=0; xx<xUB; xx++) { p.setX(xx); // raised = !raised; for(int yy=0; yy<yUB; yy++) { p.setY(yy); try{ temp = (int) grid.readAbsTemp(current,p).getValue(); } catch(Exception e){} temp /= Temp.MULTIPLIER; // hold blue at 255, increase amount of green if(temp < 255) { rgb = (temp << 8) + 255; } // hold green at 255, decrease amount of blue else if(temp < 511) { rgb = (255 << 8) + (255 - (temp - 255)); } // hold green at 255 increase amount of red else if(temp < 767) { rgb = ((temp - 511)<<16) + (255 << 8); } // hold red at 255 decrease amount of green else if(temp < 1023) { rgb = (255 << 16) + ((255 - (temp - 767)) << 8); } else { rgb = 0xff0000; } offScreen.setColor(new java.awt.Color(rgb)); //offScreen.fill3DRect(xx*blockWidth, yy*blockWidth, // blockWidth, blockHeight, true); //boolean raised = (xx%2 == 0)&&(yy%2 == 0); //cool display // raised = !raised; offScreen.fill3DRect(xx*blockWidth, yy*blockWidth, blockWidth, blockHeight, true); } } repaint(); } } 1027444357647 600 /home/tsmart/public_html/classes/eclpssSaves/HeatDiffusion/HeatDiffusion.mdl ew Grid units aren't necessary. 1002427200000 /home/tsmart/public_html/classes/eclpssSaves/HeatDiffusion/Grid/SimGrid.grd SimGrid 2 2 1029372789259 4 1.0 100 -1 4 1.0 100 -1 true 1 /** *<UL><BR> *<LI><b>Author: </b>ew</LI><BR> *<LI><b>Date Created: </b>7/22/02</LI><BR> *<LI><b>Date Saved: </b>8/14/02</LI><BR><br> *<LI><b>Version: </b></LI><BR> *<LI><b>Purpose: </b></LI><BR> *<LI><b>Type of Grid: </b></LI><BR> *<LI><b>Time step, time scale: </b></LI><BR> *<LI><b>Usage Notes: </b></LI><BR> *<LI><b>Miscellaneous Notes: </b></LI><BR> */ HeatDiffusion 9999 ew 1014872400000 ew 2/28/02 7/22/02 7/22/02 7/22/02 7/22/02 7/22/02 7/22/02 "silly" simultation to show heat diffusion. /home/tsmart/public_html/classes/eclpssSaves/HeatDiffusion/Components/CoolDown.cpt /** *<UL><BR> *<LI><BIG>This is an <font color = "mediumblue"><i> ECLPSS Component</i></font></BIG></LI><BR> *<LI><b>Author: </b>ew</LI><BR> *<LI><b>Date Created: </b>2/28/02</LI><BR> *<LI><b>Date Saved: </b>7/22/02</LI><BR> *<LI><b>Version: </b>7/22/02</LI><BR><br> *<LI><b>Purpose: </b>7/22/02</LI><BR> *<LI><b>Algorithm Citation: </b>7/22/02</LI><BR> *<LI><b>Modification History: </b>7/22/02</LI><BR> *<LI><b>Assumptions: </b>7/22/02</LI><BR> *<LI><b>Limitations: </b></LI><BR> *<LI><b>Known Problems: </b>"silly" simultation to show heat diffusion.</LI><BR> *<LI><b>Suggested Improvements: </b></LI><BR> *<LI><b>Testing Statistics: </b></LI><BR> *<LI><b>Models used in: </b></LI><BR> *<LI><b>General Comments: </b></LI><BR> */ // make an array of hot cell Points (for current and previous time // frames. Point [] hotCellsCurrent = new Point [2]; Point [] hotCellsPrevious = new Point [2]; static int xLeft, xRight, yLeft, yRight; // ECLPSSgrid is a "generic" grid - the framework replaces // it with the correct name at model run time. void getBounds(ECLPSSGrid grid, Point current) { xLeft = grid.getAbsXLB()-1; xRight = grid.getAbsXUB()-1; yLeft = grid.getAbsYLB()-1; yRight = grid.getAbsYUB()-1; } // relative offsets from current static final Point N = new Point( 0, -1, BACK1 ); static final Point S = new Point( 0, 1, BACK1 ); static final Point E = new Point( 1, 0, BACK1 ); static final Point W = new Point( -1, 0, BACK1 ); /** checks to see if the 2-D grid cell at absolute address (x,y) is a corner cell. @param x - absolute address of x @param y - absolute address of y @return true if a corner cell, false otherwise */ public boolean isCorner(int x, int y) { return ( (x==xLeft && y==yLeft) || (x==xRight && y==yLeft) || (x==xLeft && y==yRight) || (x==xRight && y==yRight) ); } /** checks to see if the 2-D grid cell at absolute address (x,y) is on a side of the grid. @param x - absolute address of x @param y - absolute address of y @return true if (x,y) the side of grid, false otherwise */ public boolean isSide(int x, int y) { return (!isCorner(x,y) && (x==xLeft || x==xRight || y==yLeft || y==yRight) ); } CoolDown 2 // initialize the hot cells to be in the upper left and lower right int xlbindex = grid.getAbsXLB()-1; int ylbindex = grid.getAbsYLB()-1; int yub= grid.getAbsYUB(); int xub= grid.getAbsXUB(); // iterate over the entire grid for (int i=xlbindex; i< xub; i++) { for (int j = ylbindex; j<yub; j++) { // the designated cells hot // and make the rest of the grid cold if (i == xlbindex && j == ylbindex) { // initialize grid grid.writeAbsTemp(current,new Point(i,j,NOW)).add(new Temp(Temp.Value_MAX)); // initialize arrays hotCellsCurrent[0] = new Point(i,j,NOW); hotCellsPrevious[0] = new Point(i,j,BACK1); } if (i == xub-1 && j == yub-1 ) { // initialize grid grid.writeAbsTemp(current,new Point(i,j,0)).add(new Temp(Temp.Value_MAX)); // initialize arrays hotCellsCurrent[1] = new Point(i,j,NOW); hotCellsPrevious[1] = new Point(i,j,BACK1); } else { // just a regular cell - initialize grid grid.writeAbsTemp(current,new Point(i,j,NOW)).add(new Temp(Temp.Value_MIN)); } } } //write out the hot cell addresses for (int i=0; i< hotCellsCurrent.length; i++) { System.out.println("hotCellsCurrent["+i+"] = "+hotCellsCurrent[i]); System.out.println("hotCellsPrevious["+i+"] = "+hotCellsPrevious[i]); } int x = current.getX(); int y = current.getY(); int it = model.getTimeIteration(); // print the iteration number to the console modulo displayRate //if (x==0 && y==0 && ((it-1) % model.displayRate == 0) ) //System.out.println("Iteration "+it); double temp = 0.0; double north = 0.0; double south = 0.0; double east = 0.0; double west = 0.0; try { south = grid.readTemp(current,S).getValue(); } catch (Exception e) {} try { north = grid.readTemp(current,N).getValue(); } catch (Exception e) {} try { east = grid.readTemp(current,E).getValue(); } catch (Exception e) {} try { west = grid.readTemp(current,W).getValue(); } catch (Exception e) {} getBounds(grid, current); // since all points are either correct (within the grid boundaries) // or 0.0 (all points are initialized to 0, and any out-of-bound gets // are caught) //we can add them all up, and then just divide by the appropriate amount. //I'm multiplying (where there is no roundoff error) since it is faster than divide if (isCorner(x,y)) { grid.writeTemp(current,CURRENT).add(new Temp( (north+south+east+west)*0.5)); } else if (isSide(x,y)) { grid.writeTemp(current,CURRENT).add(new Temp( (north+south+east+west)/3.0)); } else { grid.writeTemp(current,CURRENT).add(new Temp( (north+south+east+west)*0.25)); } for (int i=0; i< hotCellsCurrent.length; i++) { if (x==hotCellsCurrent[i].getX() && y == hotCellsCurrent[i].getY()) { // get previous temperature temp = grid.readAbsTemp(hotCellsPrevious[i]).getValue(); //System.out.println("read temp "+temp+" from "+hotCellsPrevious[i]); //System.out.println(model.getTimeIteration()); // write it here and now grid.writeAbsTemp(hotCellsCurrent[i]).setValue(temp); //System.out.println("write temp to "+hotCellsCurrent[i]); } } 1027349505733 1 true displayRate public static final int displayRate = 1; int 1 1029372947592