reading in data... ATMS 502/CSE 566 - Program 4, Fall 2013 Case 1, domain (106,106), cone type=1, radius=0.08, origin ( 0.00, 0.30), (*****,-0.30) Nest refinement 3:1, check every 5 steps Time step 0.008727, # of steps -360 Time steps to be plotted: -1 Method: 2nd-order Crowley Nest feedback to coarse grid is ON. Time interpolation of nest borders is OFF. Step Time Max at I , J Min at I , J --------------------------------------------------- 0 0.00000 9.90086 53 85 0.00000 1 1 nsteps reset to 40 INITIAL NEST at time step 1 Max trunc error = 2.3406 bounded by ( 45: 62 , 80: 90); center = ( 53, 85) PLACING nest; N= 1 TruncErr center 53 85 nest bounds: X = ( 36: 71), Y = ( 68:103) Interpolating > coarse to nest; step 1, nestX1,Y1 36 68, ratio 3 Interpolating > coarse to nest; step 1, nestX1,Y1 36 68, ratio 3 time_interp F Feedback from nest to coarse grid. Interpolating > nest to coarse; step 1, nestX1,Y1 36 68, ratio 3 1 0.00873 9.92718 53 85 -0.01542 61 83 Plotting coarse + nested grids at n= 1 Interpolating > nest to coarse; step 2, nestX1,Y1 36 68, ratio 3 2 0.01745 9.87591 52 85 -0.00820 48 92 Plotting coarse + nested grids at n= 2 Interpolating > nest to coarse; step 3, nestX1,Y1 36 68, ratio 3 3 0.02618 9.93417 52 85 -0.02325 60 84 Plotting coarse + nested grids at n= 3 Interpolating > nest to coarse; step 4, nestX1,Y1 36 68, ratio 3 4 0.03491 9.89707 51 85 -0.01658 47 92 Plotting coarse + nested grids at n= 4 NEST CHECK at time step 5 Max trunc error = 2.2867 bounded by ( 43: 60 , 80: 90); center = ( 51, 85) Old nest saved: X1,X2,Y1,Y2= 36 71 68 103 MOVING nest; N= 5 TruncErr center 51 85 nest bounds: X = ( 34: 69), Y = ( 68:103) Interpolating > coarse to nest w/nest copy; step 5, nestX1,Y1 34 68, ratio 3 Dointerp: copying from old nest ( 1:100, 1:106) to new ( 7:106, 1:106); 94.3% overlap Interpolating > coarse to nest w/nest copy; step 5, nestX1,Y1 34 68, ratio 3 Dointerp: copying from old nest ( 1:100, 1:106) to new ( 7:106, 1:106); 94.3% overlap Interpolating > nest to coarse; step 5, nestX1,Y1 34 68, ratio 3 5 0.04363 9.90939 51 85 -0.02636 59 84 Plotting coarse + nested grids at n= 5 Interpolating > nest to coarse; step 6, nestX1,Y1 34 68, ratio 3 6 0.05236 9.90018 50 85 -0.02234 58 87 Plotting coarse + nested grids at n= 6 Interpolating > nest to coarse; step 7, nestX1,Y1 34 68, ratio 3 7 0.06109 9.86447 50 85 -0.03351 58 85 Plotting coarse + nested grids at n= 7 Interpolating > nest to coarse; step 8, nestX1,Y1 34 68, ratio 3 8 0.06981 9.87158 49 85 -0.03426 57 87 Plotting coarse + nested grids at n= 8 Interpolating > nest to coarse; step 9, nestX1,Y1 34 68, ratio 3 9 0.07854 9.78646 49 85 -0.03960 57 85 Plotting coarse + nested grids at n= 9 NEST CHECK at time step 10 Max trunc error = 2.3474 bounded by ( 40: 57 , 80: 90); center = ( 48, 85) Old nest saved: X1,X2,Y1,Y2= 34 69 68 103 MOVING nest; N= 10 TruncErr center 48 85 nest bounds: X = ( 31: 66), Y = ( 68:103) Interpolating > coarse to nest w/nest copy; step 10, nestX1,Y1 31 68, ratio 3 Dointerp: copying from old nest ( 1: 97, 1:106) to new ( 10:106, 1:106); 91.5% overlap Interpolating > coarse to nest w/nest copy; step 10, nestX1,Y1 31 68, ratio 3 Dointerp: copying from old nest ( 1: 97, 1:106) to new ( 10:106, 1:106); 91.5% overlap Interpolating > nest to coarse; step 10, nestX1,Y1 31 68, ratio 3 10 0.08727 9.80016 48 85 -0.04184 56 87 Plotting coarse + nested grids at n= 10 Interpolating > nest to coarse; step 11, nestX1,Y1 31 68, ratio 3 Plotting coarse + nested grids at n= 11 Interpolating > nest to coarse; step 12, nestX1,Y1 31 68, ratio 3 Plotting coarse + nested grids at n= 12 Interpolating > nest to coarse; step 13, nestX1,Y1 31 68, ratio 3 Plotting coarse + nested grids at n= 13 Interpolating > nest to coarse; step 14, nestX1,Y1 31 68, ratio 3 Plotting coarse + nested grids at n= 14 NEST CHECK at time step 15 Max trunc error = 2.3526 bounded by ( 38: 55 , 79: 89); center = ( 46, 84) Old nest saved: X1,X2,Y1,Y2= 31 66 68 103 MOVING nest; N= 15 TruncErr center 46 84 nest bounds: X = ( 29: 64), Y = ( 67:102) Interpolating > coarse to nest w/nest copy; step 15, nestX1,Y1 29 67, ratio 3 Dointerp: copying from old nest ( 1:100, 1:103) to new ( 7:106, 4:106); 91.7% overlap Interpolating > coarse to nest w/nest copy; step 15, nestX1,Y1 29 67, ratio 3 Dointerp: copying from old nest ( 1:100, 1:103) to new ( 7:106, 4:106); 91.7% overlap Interpolating > nest to coarse; step 15, nestX1,Y1 29 67, ratio 3 Plotting coarse + nested grids at n= 15 Interpolating > nest to coarse; step 16, nestX1,Y1 29 67, ratio 3 Plotting coarse + nested grids at n= 16 Interpolating > nest to coarse; step 17, nestX1,Y1 29 67, ratio 3 Plotting coarse + nested grids at n= 17 Interpolating > nest to coarse; step 18, nestX1,Y1 29 67, ratio 3 Plotting coarse + nested grids at n= 18 Interpolating > nest to coarse; step 19, nestX1,Y1 29 67, ratio 3 Plotting coarse + nested grids at n= 19 NEST CHECK at time step 20 Max trunc error = 2.1572 bounded by ( 35: 52 , 76: 91); center = ( 43, 83) Old nest saved: X1,X2,Y1,Y2= 29 64 67 102 MOVING nest; N= 20 TruncErr center 43 83 nest bounds: X = ( 26: 61), Y = ( 66:101) Interpolating > coarse to nest w/nest copy; step 20, nestX1,Y1 26 66, ratio 3 Dointerp: copying from old nest ( 1: 97, 1:103) to new ( 10:106, 4:106); 88.9% overlap Interpolating > coarse to nest w/nest copy; step 20, nestX1,Y1 26 66, ratio 3 Dointerp: copying from old nest ( 1: 97, 1:103) to new ( 10:106, 4:106); 88.9% overlap Interpolating > nest to coarse; step 20, nestX1,Y1 26 66, ratio 3 20 0.17453 9.86183 43 83 -0.07591 51 84 Plotting coarse + nested grids at n= 20 Interpolating > nest to coarse; step 21, nestX1,Y1 26 66, ratio 3 Plotting coarse + nested grids at n= 21 Interpolating > nest to coarse; step 22, nestX1,Y1 26 66, ratio 3 Plotting coarse + nested grids at n= 22 Interpolating > nest to coarse; step 23, nestX1,Y1 26 66, ratio 3 Plotting coarse + nested grids at n= 23 Interpolating > nest to coarse; step 24, nestX1,Y1 26 66, ratio 3 Plotting coarse + nested grids at n= 24 NEST CHECK at time step 25 Max trunc error = 2.2867 bounded by ( 32: 50 , 74: 90); center = ( 41, 82) Old nest saved: X1,X2,Y1,Y2= 26 61 66 101 MOVING nest; N= 25 TruncErr center 41 82 nest bounds: X = ( 24: 59), Y = ( 65:100) Interpolating > coarse to nest w/nest copy; step 25, nestX1,Y1 24 65, ratio 3 Dointerp: copying from old nest ( 1:100, 1:103) to new ( 7:106, 4:106); 91.7% overlap Interpolating > coarse to nest w/nest copy; step 25, nestX1,Y1 24 65, ratio 3 Dointerp: copying from old nest ( 1:100, 1:103) to new ( 7:106, 4:106); 91.7% overlap Interpolating > nest to coarse; step 25, nestX1,Y1 24 65, ratio 3 Plotting coarse + nested grids at n= 25 Interpolating > nest to coarse; step 26, nestX1,Y1 24 65, ratio 3 Plotting coarse + nested grids at n= 26 Interpolating > nest to coarse; step 27, nestX1,Y1 24 65, ratio 3 Plotting coarse + nested grids at n= 27 Interpolating > nest to coarse; step 28, nestX1,Y1 24 65, ratio 3 Plotting coarse + nested grids at n= 28 Interpolating > nest to coarse; step 29, nestX1,Y1 24 65, ratio 3 Plotting coarse + nested grids at n= 29 NEST CHECK at time step 30 Max trunc error = 2.2244 bounded by ( 30: 47 , 73: 89); center = ( 38, 81) Old nest saved: X1,X2,Y1,Y2= 24 59 65 100 MOVING nest; N= 30 TruncErr center 38 81 nest bounds: X = ( 21: 56), Y = ( 64: 99) Interpolating > coarse to nest w/nest copy; step 30, nestX1,Y1 21 64, ratio 3 Dointerp: copying from old nest ( 1: 97, 1:103) to new ( 10:106, 4:106); 88.9% overlap Interpolating > coarse to nest w/nest copy; step 30, nestX1,Y1 21 64, ratio 3 Dointerp: copying from old nest ( 1: 97, 1:103) to new ( 10:106, 4:106); 88.9% overlap Interpolating > nest to coarse; step 30, nestX1,Y1 21 64, ratio 3 30 0.26180 9.83384 38 81 -0.10039 46 83 Plotting coarse + nested grids at n= 30 Interpolating > nest to coarse; step 31, nestX1,Y1 21 64, ratio 3 Plotting coarse + nested grids at n= 31 Interpolating > nest to coarse; step 32, nestX1,Y1 21 64, ratio 3 Plotting coarse + nested grids at n= 32 Interpolating > nest to coarse; step 33, nestX1,Y1 21 64, ratio 3 Plotting coarse + nested grids at n= 33 Interpolating > nest to coarse; step 34, nestX1,Y1 21 64, ratio 3 Plotting coarse + nested grids at n= 34 NEST CHECK at time step 35 Max trunc error = 2.2888 bounded by ( 28: 45 , 72: 88); center = ( 36, 80) Old nest saved: X1,X2,Y1,Y2= 21 56 64 99 MOVING nest; N= 35 TruncErr center 36 80 nest bounds: X = ( 19: 54), Y = ( 63: 98) Interpolating > coarse to nest w/nest copy; step 35, nestX1,Y1 19 63, ratio 3 Dointerp: copying from old nest ( 1:100, 1:103) to new ( 7:106, 4:106); 91.7% overlap Interpolating > coarse to nest w/nest copy; step 35, nestX1,Y1 19 63, ratio 3 Dointerp: copying from old nest ( 1:100, 1:103) to new ( 7:106, 4:106); 91.7% overlap Interpolating > nest to coarse; step 35, nestX1,Y1 19 63, ratio 3 Plotting coarse + nested grids at n= 35 Interpolating > nest to coarse; step 36, nestX1,Y1 19 63, ratio 3 Plotting coarse + nested grids at n= 36 Interpolating > nest to coarse; step 37, nestX1,Y1 19 63, ratio 3 Plotting coarse + nested grids at n= 37 Interpolating > nest to coarse; step 38, nestX1,Y1 19 63, ratio 3 Plotting coarse + nested grids at n= 38 Interpolating > nest to coarse; step 39, nestX1,Y1 19 63, ratio 3 Plotting coarse + nested grids at n= 39 NEST CHECK at time step 40 Max trunc error = 2.2775 bounded by ( 26: 43 , 70: 87); center = ( 34, 78) Old nest saved: X1,X2,Y1,Y2= 19 54 63 98 MOVING nest; N= 40 TruncErr center 34 78 nest bounds: X = ( 17: 52), Y = ( 61: 96) Interpolating > coarse to nest w/nest copy; step 40, nestX1,Y1 17 61, ratio 3 Dointerp: copying from old nest ( 1:100, 1:100) to new ( 7:106, 7:106); 89.0% overlap Interpolating > coarse to nest w/nest copy; step 40, nestX1,Y1 17 61, ratio 3 Dointerp: copying from old nest ( 1:100, 1:100) to new ( 7:106, 7:106); 89.0% overlap Interpolating > nest to coarse; step 40, nestX1,Y1 17 61, ratio 3 40 0.34907 9.78443 33 78 -0.12665 41 81 Plotting coarse + nested grids at n= 40 Plotting smax Plotting error field Error plotted. Grid 1 errors: max corr -.0091 error: tot 0.59246 dissip 0.00002 disper 0.59244 Interpolating > coarse to nest; step 41, nestX1,Y1 17 61, ratio 3 Grid 2 errors: max corr -.0532 error: tot 3.37743 dissip 0.60778 disper 2.76965