-- FreeFem++ v4.6 (Thu Apr 2 15:47:38 CEST 2020 - git v4.6) Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // convection.edp 2 : // 3 : // Discussion: 4 : // 5 : // Solve a pure transport problem using the method of characteristics. 6 : // 7 : // Location: 8 : // 9 : // http://people.sc.fsu.edu/~jburkardt/freefem_src/convection/convection.edp 10 : // 11 : // Modified: 12 : // 13 : // 17 May 2015 14 : // 15 : // Author: 16 : // 17 : // Florian De Vuyst 18 : // 19 : // Reference: 20 : // 21 : // Florian De Vuyst, 22 : // Numerical modeling of transport problems using freefem++ software - 23 : // with examples in biology, CFD, traffic flow and energy transfer, 24 : // HAL id: cel-00842234 25 : // https://cel.archives-ouvertes.fr/cel-00842234 26 : // 27 : real Lx = 6; 28 : real Ly = 4; 29 : real dt = 3; 30 : real dtsnap = 100.0; 31 : real time = 0.0; *** Warning The identifier time hide a Global identifier 32 : real tsnap = dt; 33 : int itmax=300; 34 : // 35 : real[int] A(2); 36 : real[int] B(2); 37 : real[int] C(2); 38 : real[int] D(2); 39 : real[int] E(2); 40 : real[int] F(2); 41 : real[int] G(2); 42 : real[int] H(2); 43 : real[int] I(2); 44 : real[int] J(2); 45 : real[int] K(2); 46 : real[int] L(2); 47 : real[int] O(2); 48 : real[int] P(2) *** Warning The identifier P hide a Global identifier ; 49 : real[int] Q(2); 50 : real[int] R(2); 51 : real[int] itplot(itmax); 52 : real[int] mass(itmax); 53 : 54 : cout << "\n"; 55 : cout << "convection:\n"; 56 : cout << " A mass convection problem.\n"; 57 : 58 : A = [0,0]; 59 : B = [Lx,0]; 60 : C = [Lx,Ly]; 61 : D = [0,Ly]; 62 : E = [Lx/2,Ly/4]; 63 : F = [Lx/2,3*Ly/4]; 64 : G = [Lx/6,Ly/4]; 65 : H = [5*Lx/6,Ly/4]; 66 : I = [Lx/6,3*Ly/4]; 67 : J = [5*Lx/6,3*Ly/4]; 68 : O = [0, Ly/2]; 69 : P = [Lx/3,Ly/2]; 70 : Q = [2*Lx/3, Ly/2]; 71 : R = [Lx,Ly/2]; 72 : 73 : border c1(t=0,1) { x=(1-t)*A[0]+t*B[0]; y=(1-t)*A[1]+t*B[1]; } 74 : border c2(t=0,1) { x=(1-t)*B[0]+t*C[0]; y=(1-t)*B[1]+t*C[1]; } 75 : border c3(t=0,1) { x=(1-t)*C[0]+t*D[0]; y=(1-t)*C[1]+t*D[1]; } 76 : border c4(t=0,1) { x=(1-t)*D[0]+t*A[0]; y=(1-t)*D[1]+t*A[1]; } 77 : border c5(t=0,1) { x=(1-t)*E[0]+t*F[0]; y=(1-t)*E[1]+t*F[1]; } 78 : border c6(t=0,1) { x=(1-t)*G[0]+t*H[0]; y=(1-t)*G[1]+t*H[1]; } 79 : border c7(t=0,1) { x=(1-t)*I[0]+t*J[0]; y=(1-t)*I[1]+t*J[1]; } 80 : border c8(t=0,1) { x=(1-t)*O[0]+t*P[0]; y=(1-t)*O[1]+t*P[1]; } 81 : border c9(t=0,1) { x=(1-t)*Q[0]+t*R[0]; y=(1-t)*Q[1]+t*R[1]; } 82 : // 83 : // Build the mesh. 84 : // 85 : int nn = 20; 86 : 87 : mesh Th = buildmesh ( 88 : c1(8*nn) + c2(6*nn) + c3(8*nn) + c4(6*nn) 89 : + c5(4*nn) + c6(5*nn) + c7(5*nn) + c8(3*nn) + c9(3*nn) ); 90 : 91 : plot ( Th, wait = 0, ps = "convection_mesh.ps" ); 92 : // 93 : func fy = x - Lx / 2; 94 : // 95 : // Define the finite element spaces. 96 : // 97 : fespace Uh ( Th, P1b ); 98 : fespace Vh ( Th, P1 ); 99 : 100 : Uh u1; 101 : Uh u1h; 102 : Uh u2; 103 : Uh u2h; 104 : Vh p; 105 : Vh ph; 106 : Vh q; 107 : // 108 : // Define the Stokes problem: 109 : // 110 : problem Stokes ( [u1, u2, p], [u1h, u2h,ph] ) = 111 : int2d(Th) (dx(u1)*dx(u1h)+dy(u1)*dy(u1h)) 112 : + int2d(Th) (dx(u2)*dx(u2h)+dy(u2)*dy(u2h)) 113 : + int2d(Th) (dx(p)*u1h+dy(p)*u2h) 114 : + int2d(Th) (dx(u1)*ph+dy(u2)*ph) 115 : - int2d(Th) (fy*u2h) 116 : + on ( c1, c2, c3, c4, c5, c6, c7, c8, c9, u1=0, u2=0 ); 117 : // 118 : // Solve the Stokes problem. 119 : // 120 : Stokes; 121 : // 122 : // Plot the convection velocity field. 123 : // 124 : plot ( [u1, u2], ps = "convection_velocity.ps" ); 125 : // 126 : q = ( sqrt((x-Lx/2)^2+(y-Ly/8)^2)<0.2 ); 127 : // 128 : // Measure the initial total mass. 129 : // 130 : real mass0 = int2d ( Th ) ( q ); 131 : // 132 : plot ( q, nbiso=40, fill=1, wait = 0 ); 133 : // 134 : for ( int it = 0; it < itmax; it++ ) 135 : { 136 : itplot[it] = it; 137 : q = convect ( [u1,u2], -dt, q ); 138 : mass[it] = int2d(Th) (q); 139 : time += dt; 140 : plot ( q, nbiso=40, fill=1 ); 141 : cout << "Mass0 = " << mass0 << " Mass(" << time << ") = " << mass[it] << endl; 142 : if ( ( it + 1 ) % 100 == 0 ) 143 : { 144 : tsnap = tsnap + dtsnap; 145 : string filename = "q_time_"+time+".ps"; 146 : plot ( q, nbiso=60, fill=1, ps=filename ); 147 : } 148 : } 149 : 150 : plot ( [ itplot, mass ], cmm="Mass iteration", ps="mass_histogram.ps", value=1 ); 151 : // 152 : // Write a file. 153 : // 154 : ofstream of ( "mass_histogram.txt" ); 155 : for ( int it = 0; it < itmax; it++ ) 156 : { 157 : of << itplot[it] << " " << mass[it] << endl; 158 : } 159 : // 160 : // Terminate. 161 : // 162 : cout << "\n"; 163 : cout << "convection:\n"; 164 : cout << " Normal end of execution.\n"; 165 : 166 : exit ( 0 ); 167 : 168 : sizestack + 1024 =3744 ( 2720 ) convection: A mass convection problem. -- mesh: Nb of Triangles = 40618, Nb of Vertices 20590 -- Solve : min -0.0705034 max 0.070457 min -0.0601026 max 0.0601263 min -184.59 max -177.854 Mass0 = 0.125382 Mass(3) = 0.124936 Mass0 = 0.125382 Mass(6) = 0.124878 Mass0 = 0.125382 Mass(9) = 0.124851 Mass0 = 0.125382 Mass(12) = 0.125366 Mass0 = 0.125382 Mass(15) = 0.125021 Mass0 = 0.125382 Mass(18) = 0.125213 Mass0 = 0.125382 Mass(21) = 0.125001 Mass0 = 0.125382 Mass(24) = 0.12465 Mass0 = 0.125382 Mass(27) = 0.124465 Mass0 = 0.125382 Mass(30) = 0.124457 Mass0 = 0.125382 Mass(33) = 0.124764 Mass0 = 0.125382 Mass(36) = 0.124656 Mass0 = 0.125382 Mass(39) = 0.124724 Mass0 = 0.125382 Mass(42) = 0.12452 Mass0 = 0.125382 Mass(45) = 0.12445 Mass0 = 0.125382 Mass(48) = 0.123984 Mass0 = 0.125382 Mass(51) = 0.123419 Mass0 = 0.125382 Mass(54) = 0.123439 Mass0 = 0.125382 Mass(57) = 0.123378 Mass0 = 0.125382 Mass(60) = 0.123285 Mass0 = 0.125382 Mass(63) = 0.123349 Mass0 = 0.125382 Mass(66) = 0.123358 Mass0 = 0.125382 Mass(69) = 0.123324 Mass0 = 0.125382 Mass(72) = 0.123397 Mass0 = 0.125382 Mass(75) = 0.123592 Mass0 = 0.125382 Mass(78) = 0.123829 Mass0 = 0.125382 Mass(81) = 0.124042 Mass0 = 0.125382 Mass(84) = 0.124068 Mass0 = 0.125382 Mass(87) = 0.124181 Mass0 = 0.125382 Mass(90) = 0.124307 Mass0 = 0.125382 Mass(93) = 0.124201 Mass0 = 0.125382 Mass(96) = 0.123999 Mass0 = 0.125382 Mass(99) = 0.123631 Mass0 = 0.125382 Mass(102) = 0.123346 Mass0 = 0.125382 Mass(105) = 0.123209 Mass0 = 0.125382 Mass(108) = 0.123035 Mass0 = 0.125382 Mass(111) = 0.122875 Mass0 = 0.125382 Mass(114) = 0.122759 Mass0 = 0.125382 Mass(117) = 0.12267 Mass0 = 0.125382 Mass(120) = 0.122644 Mass0 = 0.125382 Mass(123) = 0.12262 Mass0 = 0.125382 Mass(126) = 0.122531 Mass0 = 0.125382 Mass(129) = 0.122407 Mass0 = 0.125382 Mass(132) = 0.122369 Mass0 = 0.125382 Mass(135) = 0.122422 Mass0 = 0.125382 Mass(138) = 0.12245 Mass0 = 0.125382 Mass(141) = 0.122392 Mass0 = 0.125382 Mass(144) = 0.122282 Mass0 = 0.125382 Mass(147) = 0.122157 Mass0 = 0.125382 Mass(150) = 0.122003 Mass0 = 0.125382 Mass(153) = 0.12182 Mass0 = 0.125382 Mass(156) = 0.121655 Mass0 = 0.125382 Mass(159) = 0.121541 Mass0 = 0.125382 Mass(162) = 0.121491 Mass0 = 0.125382 Mass(165) = 0.121465 Mass0 = 0.125382 Mass(168) = 0.121402 Mass0 = 0.125382 Mass(171) = 0.121298 Mass0 = 0.125382 Mass(174) = 0.121218 Mass0 = 0.125382 Mass(177) = 0.121201 Mass0 = 0.125382 Mass(180) = 0.121211 Mass0 = 0.125382 Mass(183) = 0.121206 Mass0 = 0.125382 Mass(186) = 0.121172 Mass0 = 0.125382 Mass(189) = 0.121114 Mass0 = 0.125382 Mass(192) = 0.121044 Mass0 = 0.125382 Mass(195) = 0.120932 Mass0 = 0.125382 Mass(198) = 0.12074 Mass0 = 0.125382 Mass(201) = 0.120531 Mass0 = 0.125382 Mass(204) = 0.120391 Mass0 = 0.125382 Mass(207) = 0.120297 Mass0 = 0.125382 Mass(210) = 0.120173 Mass0 = 0.125382 Mass(213) = 0.120033 Mass0 = 0.125382 Mass(216) = 0.119939 Mass0 = 0.125382 Mass(219) = 0.119872 Mass0 = 0.125382 Mass(222) = 0.119752 Mass0 = 0.125382 Mass(225) = 0.119572 Mass0 = 0.125382 Mass(228) = 0.119399 Mass0 = 0.125382 Mass(231) = 0.119249 Mass0 = 0.125382 Mass(234) = 0.119112 Mass0 = 0.125382 Mass(237) = 0.119019 Mass0 = 0.125382 Mass(240) = 0.118984 Mass0 = 0.125382 Mass(243) = 0.118986 Mass0 = 0.125382 Mass(246) = 0.119019 Mass0 = 0.125382 Mass(249) = 0.119096 Mass0 = 0.125382 Mass(252) = 0.119232 Mass0 = 0.125382 Mass(255) = 0.119417 Mass0 = 0.125382 Mass(258) = 0.119592 Mass0 = 0.125382 Mass(261) = 0.119707 Mass0 = 0.125382 Mass(264) = 0.119742 Mass0 = 0.125382 Mass(267) = 0.119686 Mass0 = 0.125382 Mass(270) = 0.119553 Mass0 = 0.125382 Mass(273) = 0.119382 Mass0 = 0.125382 Mass(276) = 0.119199 Mass0 = 0.125382 Mass(279) = 0.119009 Mass0 = 0.125382 Mass(282) = 0.118819 Mass0 = 0.125382 Mass(285) = 0.118639 Mass0 = 0.125382 Mass(288) = 0.118465 Mass0 = 0.125382 Mass(291) = 0.118283 Mass0 = 0.125382 Mass(294) = 0.118087 Mass0 = 0.125382 Mass(297) = 0.117874 Mass0 = 0.125382 Mass(300) = 0.117634 Mass0 = 0.125382 Mass(303) = 0.117371 Mass0 = 0.125382 Mass(306) = 0.117102 Mass0 = 0.125382 Mass(309) = 0.116849 Mass0 = 0.125382 Mass(312) = 0.116622 Mass0 = 0.125382 Mass(315) = 0.116419 Mass0 = 0.125382 Mass(318) = 0.116222 Mass0 = 0.125382 Mass(321) = 0.116015 Mass0 = 0.125382 Mass(324) = 0.115796 Mass0 = 0.125382 Mass(327) = 0.115569 Mass0 = 0.125382 Mass(330) = 0.115348 Mass0 = 0.125382 Mass(333) = 0.115155 Mass0 = 0.125382 Mass(336) = 0.115004 Mass0 = 0.125382 Mass(339) = 0.114885 Mass0 = 0.125382 Mass(342) = 0.114781 Mass0 = 0.125382 Mass(345) = 0.114688 Mass0 = 0.125382 Mass(348) = 0.114607 Mass0 = 0.125382 Mass(351) = 0.114537 Mass0 = 0.125382 Mass(354) = 0.114467 Mass0 = 0.125382 Mass(357) = 0.114386 Mass0 = 0.125382 Mass(360) = 0.114292 Mass0 = 0.125382 Mass(363) = 0.114186 Mass0 = 0.125382 Mass(366) = 0.114074 Mass0 = 0.125382 Mass(369) = 0.113967 Mass0 = 0.125382 Mass(372) = 0.113869 Mass0 = 0.125382 Mass(375) = 0.113768 Mass0 = 0.125382 Mass(378) = 0.113644 Mass0 = 0.125382 Mass(381) = 0.113483 Mass0 = 0.125382 Mass(384) = 0.113293 Mass0 = 0.125382 Mass(387) = 0.113088 Mass0 = 0.125382 Mass(390) = 0.112877 Mass0 = 0.125382 Mass(393) = 0.112666 Mass0 = 0.125382 Mass(396) = 0.112457 Mass0 = 0.125382 Mass(399) = 0.112257 Mass0 = 0.125382 Mass(402) = 0.112076 Mass0 = 0.125382 Mass(405) = 0.111927 Mass0 = 0.125382 Mass(408) = 0.111815 Mass0 = 0.125382 Mass(411) = 0.111737 Mass0 = 0.125382 Mass(414) = 0.111688 Mass0 = 0.125382 Mass(417) = 0.111661 Mass0 = 0.125382 Mass(420) = 0.111649 Mass0 = 0.125382 Mass(423) = 0.111642 Mass0 = 0.125382 Mass(426) = 0.111628 Mass0 = 0.125382 Mass(429) = 0.111598 Mass0 = 0.125382 Mass(432) = 0.111552 Mass0 = 0.125382 Mass(435) = 0.111494 Mass0 = 0.125382 Mass(438) = 0.111426 Mass0 = 0.125382 Mass(441) = 0.111351 Mass0 = 0.125382 Mass(444) = 0.111273 Mass0 = 0.125382 Mass(447) = 0.111193 Mass0 = 0.125382 Mass(450) = 0.111111 Mass0 = 0.125382 Mass(453) = 0.111021 Mass0 = 0.125382 Mass(456) = 0.110913 Mass0 = 0.125382 Mass(459) = 0.110782 Mass0 = 0.125382 Mass(462) = 0.11063 Mass0 = 0.125382 Mass(465) = 0.11046 Mass0 = 0.125382 Mass(468) = 0.110278 Mass0 = 0.125382 Mass(471) = 0.110087 Mass0 = 0.125382 Mass(474) = 0.109892 Mass0 = 0.125382 Mass(477) = 0.109694 Mass0 = 0.125382 Mass(480) = 0.109496 Mass0 = 0.125382 Mass(483) = 0.109296 Mass0 = 0.125382 Mass(486) = 0.109097 Mass0 = 0.125382 Mass(489) = 0.108898 Mass0 = 0.125382 Mass(492) = 0.108702 Mass0 = 0.125382 Mass(495) = 0.108512 Mass0 = 0.125382 Mass(498) = 0.108327 Mass0 = 0.125382 Mass(501) = 0.108147 Mass0 = 0.125382 Mass(504) = 0.107974 Mass0 = 0.125382 Mass(507) = 0.10781 Mass0 = 0.125382 Mass(510) = 0.107656 Mass0 = 0.125382 Mass(513) = 0.107509 Mass0 = 0.125382 Mass(516) = 0.107366 Mass0 = 0.125382 Mass(519) = 0.107223 Mass0 = 0.125382 Mass(522) = 0.107075 Mass0 = 0.125382 Mass(525) = 0.106919 Mass0 = 0.125382 Mass(528) = 0.106753 Mass0 = 0.125382 Mass(531) = 0.106577 Mass0 = 0.125382 Mass(534) = 0.106393 Mass0 = 0.125382 Mass(537) = 0.106207 Mass0 = 0.125382 Mass(540) = 0.106018 Mass0 = 0.125382 Mass(543) = 0.105825 Mass0 = 0.125382 Mass(546) = 0.105627 Mass0 = 0.125382 Mass(549) = 0.105425 Mass0 = 0.125382 Mass(552) = 0.105217 Mass0 = 0.125382 Mass(555) = 0.105003 Mass0 = 0.125382 Mass(558) = 0.104784 Mass0 = 0.125382 Mass(561) = 0.104564 Mass0 = 0.125382 Mass(564) = 0.104348 Mass0 = 0.125382 Mass(567) = 0.104139 Mass0 = 0.125382 Mass(570) = 0.103942 Mass0 = 0.125382 Mass(573) = 0.10376 Mass0 = 0.125382 Mass(576) = 0.103597 Mass0 = 0.125382 Mass(579) = 0.103454 Mass0 = 0.125382 Mass(582) = 0.103335 Mass0 = 0.125382 Mass(585) = 0.103238 Mass0 = 0.125382 Mass(588) = 0.103162 Mass0 = 0.125382 Mass(591) = 0.103099 Mass0 = 0.125382 Mass(594) = 0.103036 Mass0 = 0.125382 Mass(597) = 0.102966 Mass0 = 0.125382 Mass(600) = 0.10288 Mass0 = 0.125382 Mass(603) = 0.102776 Mass0 = 0.125382 Mass(606) = 0.102656 Mass0 = 0.125382 Mass(609) = 0.102522 Mass0 = 0.125382 Mass(612) = 0.102377 Mass0 = 0.125382 Mass(615) = 0.102226 Mass0 = 0.125382 Mass(618) = 0.10207 Mass0 = 0.125382 Mass(621) = 0.101911 Mass0 = 0.125382 Mass(624) = 0.101747 Mass0 = 0.125382 Mass(627) = 0.101577 Mass0 = 0.125382 Mass(630) = 0.101399 Mass0 = 0.125382 Mass(633) = 0.101212 Mass0 = 0.125382 Mass(636) = 0.101014 Mass0 = 0.125382 Mass(639) = 0.10081 Mass0 = 0.125382 Mass(642) = 0.100601 Mass0 = 0.125382 Mass(645) = 0.100389 Mass0 = 0.125382 Mass(648) = 0.100176 Mass0 = 0.125382 Mass(651) = 0.0999613 Mass0 = 0.125382 Mass(654) = 0.0997431 Mass0 = 0.125382 Mass(657) = 0.0995216 Mass0 = 0.125382 Mass(660) = 0.0992985 Mass0 = 0.125382 Mass(663) = 0.0990767 Mass0 = 0.125382 Mass(666) = 0.0988598 Mass0 = 0.125382 Mass(669) = 0.0986504 Mass0 = 0.125382 Mass(672) = 0.0984494 Mass0 = 0.125382 Mass(675) = 0.0982562 Mass0 = 0.125382 Mass(678) = 0.0980698 Mass0 = 0.125382 Mass(681) = 0.0978889 Mass0 = 0.125382 Mass(684) = 0.0977119 Mass0 = 0.125382 Mass(687) = 0.0975368 Mass0 = 0.125382 Mass(690) = 0.0973621 Mass0 = 0.125382 Mass(693) = 0.0971866 Mass0 = 0.125382 Mass(696) = 0.0970101 Mass0 = 0.125382 Mass(699) = 0.0968331 Mass0 = 0.125382 Mass(702) = 0.0966562 Mass0 = 0.125382 Mass(705) = 0.0964788 Mass0 = 0.125382 Mass(708) = 0.096299 Mass0 = 0.125382 Mass(711) = 0.0961137 Mass0 = 0.125382 Mass(714) = 0.0959205 Mass0 = 0.125382 Mass(717) = 0.0957182 Mass0 = 0.125382 Mass(720) = 0.0955074 Mass0 = 0.125382 Mass(723) = 0.0952894 Mass0 = 0.125382 Mass(726) = 0.0950662 Mass0 = 0.125382 Mass(729) = 0.0948404 Mass0 = 0.125382 Mass(732) = 0.094615 Mass0 = 0.125382 Mass(735) = 0.0943936 Mass0 = 0.125382 Mass(738) = 0.0941797 Mass0 = 0.125382 Mass(741) = 0.093976 Mass0 = 0.125382 Mass(744) = 0.0937842 Mass0 = 0.125382 Mass(747) = 0.0936047 Mass0 = 0.125382 Mass(750) = 0.0934371 Mass0 = 0.125382 Mass(753) = 0.0932797 Mass0 = 0.125382 Mass(756) = 0.0931305 Mass0 = 0.125382 Mass(759) = 0.0929868 Mass0 = 0.125382 Mass(762) = 0.0928464 Mass0 = 0.125382 Mass(765) = 0.0927076 Mass0 = 0.125382 Mass(768) = 0.0925696 Mass0 = 0.125382 Mass(771) = 0.092432 Mass0 = 0.125382 Mass(774) = 0.0922949 Mass0 = 0.125382 Mass(777) = 0.0921585 Mass0 = 0.125382 Mass(780) = 0.0920227 Mass0 = 0.125382 Mass(783) = 0.0918869 Mass0 = 0.125382 Mass(786) = 0.0917495 Mass0 = 0.125382 Mass(789) = 0.0916089 Mass0 = 0.125382 Mass(792) = 0.0914634 Mass0 = 0.125382 Mass(795) = 0.091312 Mass0 = 0.125382 Mass(798) = 0.0911543 Mass0 = 0.125382 Mass(801) = 0.0909904 Mass0 = 0.125382 Mass(804) = 0.0908207 Mass0 = 0.125382 Mass(807) = 0.090646 Mass0 = 0.125382 Mass(810) = 0.090467 Mass0 = 0.125382 Mass(813) = 0.0902843 Mass0 = 0.125382 Mass(816) = 0.0900986 Mass0 = 0.125382 Mass(819) = 0.0899107 Mass0 = 0.125382 Mass(822) = 0.0897213 Mass0 = 0.125382 Mass(825) = 0.0895312 Mass0 = 0.125382 Mass(828) = 0.0893412 Mass0 = 0.125382 Mass(831) = 0.0891519 Mass0 = 0.125382 Mass(834) = 0.088964 Mass0 = 0.125382 Mass(837) = 0.088778 Mass0 = 0.125382 Mass(840) = 0.0885945 Mass0 = 0.125382 Mass(843) = 0.0884136 Mass0 = 0.125382 Mass(846) = 0.0882353 Mass0 = 0.125382 Mass(849) = 0.0880589 Mass0 = 0.125382 Mass(852) = 0.0878834 Mass0 = 0.125382 Mass(855) = 0.0877077 Mass0 = 0.125382 Mass(858) = 0.0875307 Mass0 = 0.125382 Mass(861) = 0.0873515 Mass0 = 0.125382 Mass(864) = 0.0871696 Mass0 = 0.125382 Mass(867) = 0.0869851 Mass0 = 0.125382 Mass(870) = 0.0867979 Mass0 = 0.125382 Mass(873) = 0.0866082 Mass0 = 0.125382 Mass(876) = 0.0864159 Mass0 = 0.125382 Mass(879) = 0.0862209 Mass0 = 0.125382 Mass(882) = 0.086023 Mass0 = 0.125382 Mass(885) = 0.0858221 Mass0 = 0.125382 Mass(888) = 0.085618 Mass0 = 0.125382 Mass(891) = 0.0854111 Mass0 = 0.125382 Mass(894) = 0.0852022 Mass0 = 0.125382 Mass(897) = 0.0849925 Mass0 = 0.125382 Mass(900) = 0.0847833 convection: Normal end of execution. current line = 166 exit(0) err code 0 , mpirank 0 CodeAlloc : nb ptr 4400, size :505400 mpirank: 0 Ok: Normal End