diff --git a/Thermics/Q.pos b/Thermics/Q.pos
new file mode 100644
index 0000000000000000000000000000000000000000..4e60d13e3e6fe9498ffdeeac637a25d2f69361e5
--- /dev/null
+++ b/Thermics/Q.pos
@@ -0,0 +1,252 @@
+View "New view_MathEval" {
+ST(0.02634916017313363,0.02577951485829298,0,0.02646397875083406,0.02390163971325145,0,0.02832346729453371,0.0238260250095,0){6266831.340785229,7578556.193563664,8398833.273041854};
+ST(0.035,0.015,0,0.035,0.01699999999999357,0,0.03300243786989122,0.01699756213010876,0){6065306.597126331,7117703.227623346,8350257.942833206};
+ST(0.015,0.015,0,0.01699999999999357,0.015,0,0.01699756213010875,0.01699756213010876,0){820849.986238988,1437039.497774626,1685066.208253873};
+ST(0.035,0.035,0,0.0330000000000061,0.035,0,0.03300243786989124,0.03300243786989122,0){820849.9862389868,963276.3823045766,1685066.208253887};
+ST(0.015,0.035,0,0.015,0.0330000000000061,0,0.01699756213010876,0.03300243786989126,0){111089.965382423,194482.1474535455,340043.1634134258};
+ST(0.035,0.015,0,0.03300243786989122,0.01699756213010876,0,0.03299999999999222,0.015,0){6065306.597126331,8350257.942833206,7117703.227629418};
+ST(0.035,0.035,0,0.03300243786989124,0.03300243786989122,0,0.035,0.03299999999999222,0){820849.9862389868,1685066.208253887,1437039.497779934};
+ST(0.015,0.015,0,0.01699756213010875,0.01699756213010876,0,0.015,0.01700000000000744,0){820849.986238988,1685066.208253873,963276.3823053603};
+ST(0.015,0.035,0,0.01699756213010876,0.03300243786989126,0,0.01700000000000744,0.035,0){111089.965382423,340043.1634134258,194482.14745423};
+ST(0.02787134264717452,0.02636461712411311,0,0.02634916017313363,0.02577951485829298,0,0.02832346729453371,0.0238260250095,0){6373761.262323579,6266831.340785229,8398833.273041854};
+ST(0.0218094427360229,0.0263682556397269,0,0.02004288562623822,0.02710123186687632,0,0.02102282424207957,0.02483178028678816,0){3408205.856754326,2240880.53857153,3536817.246116393};
+ST(0.02630937339024055,0.02808524626481607,0,0.02709387158569049,0.02994395049886663,0,0.02481685402985528,0.02896618798686593,0){4538814.13250844,3418862.023887427,3421263.069148063};
+ST(0.0238260250094999,0.02167653270546622,0,0.02292305321427705,0.02002627187241877,0,0.0251688693519739,0.02080806878318098,0){6641234.082622026,6060221.998352256,7866845.121381707};
+ST(0.02832346729453371,0.0238260250095,0,0.02997372812758119,0.02292305321427708,0,0.02919889109015292,0.02517140974684725,0){8398833.273041854,9180996.867303077,7604442.001547791};
+ST(0.02633972125435513,0.02219686411149813,0,0.02646397875083406,0.02390163971325145,0,0.02417751621551469,0.023393463480524,0){8334025.300326461,7578556.193563664,6349753.383152546};
+ST(0.02355978257831163,0.02772127893890192,0,0.02321297569192563,0.02525629976722429,0,0.0251772925710883,0.02680645346909019,0){3638747.773110651,4785838.434520846,4986438.53885695};
+ST(0.02355978257831163,0.02772127893890192,0,0.0218094427360229,0.0263682556397269,0,0.02321297569192563,0.02525629976722429,0){3638747.773110651,3408205.856754326,4785838.434520846};
+ST(0.02633972125435513,0.02219686411149813,0,0.02417751621551469,0.023393463480524,0,0.0238260250094999,0.02167653270546622,0){8334025.300326461,6349753.383152546,6641234.082622026};
+ST(0.02102282424207957,0.02483178028678816,0,0.02004288562623822,0.02710123186687632,0,0.0192147217946296,0.02486333911250368,0){3536817.246116393,2240880.53857153,2466603.640555887};
+ST(0.02481685402985528,0.02896618798686593,0,0.02709387158569049,0.02994395049886663,0,0.0248633391125036,0.03078527820537037,0){3421263.069148063,3418862.023887427,2400098.990608844};
+ST(0.0251688693519739,0.02080806878318098,0,0.02292305321427705,0.02002627187241877,0,0.0251366608874963,0.01921472179462946,0){7866845.121381707,6060221.998352256,7845168.002306522};
+ST(0.02919889109015292,0.02517140974684725,0,0.02997372812758119,0.02292305321427708,0,0.03078527820537049,0.02513666088749638,0){7604442.001547791,9180996.867303077,7633646.319945103};
+ST(0.02787134264717452,0.02636461712411311,0,0.02630937339024055,0.02808524626481607,0,0.02634916017313363,0.02577951485829298,0){6373761.262323579,4538814.13250844,6266831.340785229};
+ST(0.015,0.0250000000000237,0,0.015,0.02300000000001896,0,0.01718800098674322,0.02423878267049852,0){820849.9862370425,963276.3823038344,1618401.802394528};
+ST(0.0250000000000237,0.035,0,0.02300000000001896,0.035,0,0.02423878267049847,0.03281199901325677,0){820849.9862409331,645703.4689333984,1389849.759568753};
+ST(0.02499999999997461,0.015,0,0.02699999999997969,0.015,0,0.02576121732950113,0.01718800098674297,0){6065306.597110934,7117703.227617423,7720189.599155798};
+ST(0.035,0.02499999999997461,0,0.035,0.02699999999997969,0,0.03281199901325701,0.02576121732950117,0){6065306.597141732,4771139.155223909,6629938.030430425};
+ST(0.015,0.02300000000001896,0,0.01643969820243331,0.02212372169562202,0,0.01718800098674322,0.02423878267049852,0){963276.3823038344,1519931.557947947,1618401.802394528};
+ST(0.02300000000001896,0.035,0,0.02212372169562198,0.03356030179756669,0,0.02423878267049847,0.03281199901325677,0){645703.4689333984,855054.256044248,1389849.759568753};
+ST(0.02699999999997969,0.015,0,0.02787627830437758,0.0164396982024332,0,0.02576121732950113,0.01718800098674297,0){7117703.227617423,8420975.407718396,7720189.599155798};
+ST(0.035,0.02699999999997969,0,0.0335603017975668,0.02787627830437762,0,0.03281199901325701,0.02576121732950117,0){4771139.155223909,4737312.561715177,6629938.030430425};
+ST(0.02633972125435513,0.02219686411149813,0,0.0238260250094999,0.02167653270546622,0,0.0251688693519739,0.02080806878318098,0){8334025.300326461,6641234.082622026,7866845.121381707};
+ST(0.02787134264717452,0.02636461712411311,0,0.02832346729453371,0.0238260250095,0,0.02919889109015292,0.02517140974684725,0){6373761.262323579,8398833.273041854,7604442.001547791};
+ST(0.02355978257831163,0.02772127893890192,0,0.02192345518181211,0.02793885538913906,0,0.0218094427360229,0.0263682556397269,0){3638747.773110651,2773264.32237239,3408205.856754326};
+ST(0.02787134264717452,0.02636461712411311,0,0.02794496480441602,0.02806239939751506,0,0.02630937339024055,0.02808524626481607,0){6373761.262323579,5004465.210416885,4538814.13250844};
+ST(0.02630937339024055,0.02808524626481607,0,0.0251772925710883,0.02680645346909019,0,0.02634916017313363,0.02577951485829298,0){4538814.13250844,4986438.53885695,6266831.340785229};
+ST(0.025,0.025,0,0.02417751621551469,0.023393463480524,0,0.02646397875083406,0.02390163971325145,0){6065306.597126335,6349753.383152546,7578556.193563664};
+ST(0.025,0.025,0,0.0251772925710883,0.02680645346909019,0,0.02321297569192563,0.02525629976722429,0){6065306.597126335,4986438.53885695,4785838.434520846};
+ST(0.01802267861831516,0.02218579774721222,0,0.01718800098674322,0.02423878267049852,0,0.01643969820243331,0.02212372169562202,0){2271057.768506458,1618401.802394528,1519931.557947947};
+ST(0.02218579774721212,0.03197732138168483,0,0.02423878267049847,0.03281199901325677,0,0.02212372169562198,0.03356030179756669,0){1293569.239122899,1389849.759568753,855054.256044248};
+ST(0.02781420225278781,0.01802267861831496,0,0.02576121732950113,0.01718800098674297,0,0.02787627830437758,0.0164396982024332,0){9167915.807340225,7720189.599155798,8420975.407718396};
+ST(0.03197732138168502,0.0278142022527879,0,0.03281199901325701,0.02576121732950117,0,0.0335603017975668,0.02787627830437762,0){5221942.849583705,6629938.030430425,4737312.561715177};
+ST(0.035,0.0189999999999931,0,0.035,0.02099999999998355,0,0.03304669728426029,0.02021852793711807,0){7710515.858034595,7710515.858038196,9109191.0106201};
+ST(0.0189999999999931,0.015,0,0.02099999999998355,0.015,0,0.02021852793711804,0.01695330271573968,0){2322362.747294063,3464558.103290316,3500792.923489964};
+ST(0.03100000000000718,0.035,0,0.02900000000000975,0.035,0,0.02978147206288093,0.033046697284261,0){1043504.847547499,1043504.847547853,1822017.581094539};
+ST(0.015,0.03100000000000718,0,0.015,0.02900000000000975,0,0.01695330271573897,0.02978147206288096,0){314297.6201831806,468876.952199062,700227.5226125802};
+ST(0.02899999999998503,0.015,0,0.03099999999999255,0.015,0,0.02980597193963426,0.01698781140710807,0){7710515.858033354,7710515.858036811,9129179.508373139};
+ST(0.015,0.02100000000001463,0,0.015,0.01900000000000812,0,0.01698781140710817,0.02019402806036576,0){1043504.847547345,1043504.847547819,1838661.946960948};
+ST(0.035,0.02899999999998503,0,0.035,0.03099999999999255,0,0.03301218859289196,0.02980597193963428,0){3464558.103309908,2322362.747301393,3491325.400350199};
+ST(0.02100000000001463,0.035,0,0.01900000000000812,0.035,0,0.02019402806036574,0.03301218859289186,0){468876.9522011194,314297.6201842385,703170.2194260015};
+ST(0.01802267861832114,0.02781420225278695,0,0.02004288562623822,0.02710123186687632,0,0.01897469406756322,0.02955587648676389,0){1293569.239124938,2240880.53857153,1189893.288594959};
+ST(0.02781420225278688,0.03197732138167883,0,0.02709387158569049,0.02994395049886663,0,0.029555876486764,0.03102530593243676,0){2271057.768509638,3418862.023887427,2959555.02311246};
+ST(0.0221857977472131,0.0180226786183216,0,0.02292305321427705,0.02002627187241877,0,0.02044412351323598,0.01897469406756427,0){5221942.849585891,6060221.998352256,3970630.174131208};
+ST(0.03197732138167837,0.02218579774721317,0,0.02997372812758119,0.02292305321427708,0,0.03102530593243569,0.0204441235132358,0){9167915.807342248,9180996.867303077,9875926.345166853};
+ST(0.0209606099800624,0.02937278878959456,0,0.01897469406756322,0.02955587648676389,0,0.02004288562623822,0.02710123186687632,0){1834895.766186544,1189893.288594959,2240880.53857153};
+ST(0.02937278878959463,0.02903939001993757,0,0.029555876486764,0.03102530593243676,0,0.02709387158569049,0.02994395049886663,0){4399739.940052692,2959555.02311246,3418862.023887427};
+ST(0.02062721121040539,0.02096060998006265,0,0.02044412351323598,0.01897469406756427,0,0.02292305321427705,0.02002627187241877,0){4115933.579576971,3970630.174131208,6060221.998352256};
+ST(0.0290393900199373,0.0206272112104053,0,0.03102530593243569,0.0204441235132358,0,0.02997372812758119,0.02292305321427708,0){9869245.814602781,9875926.345166853,9180996.867303077};
+ST(0.015,0.0250000000000237,0,0.01718800098674322,0.02423878267049852,0,0.01666349631902295,0.0263865043232969,0){820849.9862370425,1618401.802394528,1123094.041931162};
+ST(0.0250000000000237,0.035,0,0.02423878267049847,0.03281199901325677,0,0.02638650432329687,0.03333650368097703,0){820849.9862409331,1389849.759568753,1481993.704088334};
+ST(0.02499999999997461,0.015,0,0.02576121732950113,0.01718800098674297,0,0.02361349567670159,0.01666349631902287,0){6065306.597110934,7720189.599155798,5949981.737036595};
+ST(0.035,0.02499999999997461,0,0.03281199901325701,0.02576121732950117,0,0.0333365036809771,0.02361349567670162,0){6065306.597141732,6629938.030430425,7851377.66252339};
+ST(0.01802267861831516,0.02218579774721222,0,0.01698781140710817,0.02019402806036576,0,0.01900766867126057,0.02025952709593554,0){2271057.768506458,1838661.946960948,2984994.992984212};
+ST(0.02781420225278781,0.01802267861831496,0,0.02980597193963426,0.01698781140710807,0,0.02974047290406433,0.01900766867126047,0){9167915.807340225,9129179.508373139,9895343.918024337};
+ST(0.02218579774721212,0.03197732138168483,0,0.02019402806036574,0.03301218859289186,0,0.02025952709593563,0.03099233132873949,0){1293569.239122899,703170.2194260015,1156621.769609863};
+ST(0.03197732138168502,0.0278142022527879,0,0.03301218859289196,0.02980597193963428,0,0.03099233132873955,0.02974047290406432,0){5221942.849583705,3491325.400350199,3834234.301988425};
+ST(0.0192147217946296,0.02486333911250368,0,0.01836860548693934,0.02621782439543604,0,0.01718800098674322,0.02423878267049852,0){2466603.640555887,1756069.912092069,1618401.802394528};
+ST(0.0248633391125036,0.03078527820537037,0,0.02621744120711934,0.0316307091801319,0,0.02423878267049847,0.03281199901325677,0){2400098.990608844,2240667.074053982,1389849.759568753};
+ST(0.0251366608874963,0.01921472179462946,0,0.02380323230261221,0.01834737080814979,0,0.02576121732950113,0.01718800098674297,0){7845168.002306522,6627806.876270836,7720189.599155798};
+ST(0.03078527820537049,0.02513666088749638,0,0.03165262919185018,0.02380323230261227,0,0.03281199901325701,0.02576121732950117,0){7633646.319945103,8420148.806340775,6629938.030430425};
+ST(0.02219686411149813,0.02366027874564487,0,0.0200943405656333,0.02291758750152122,0,0.02215001361588054,0.02196095145324617,0){4757486.064843449,3442648.322572604,5196103.016175205};
+ST(0.02633972125435513,0.02219686411149813,0,0.02708241249847881,0.02009434056563312,0,0.02803904854675385,0.0221500136158806,0){8334025.300326461,9183174.245855754,9188072.738583822};
+ST(0.02355978257831163,0.02772127893890192,0,0.02291758750152112,0.02990565943436665,0,0.02192345518181211,0.02793885538913906,0){3638747.773110651,2269951.661359691,2773264.32237239};
+ST(0.02787134264717452,0.02636461712411311,0,0.02990565943436682,0.0270824124984789,0,0.02794496480441602,0.02806239939751506,0){6373761.262323579,6055036.612150463,5004465.210416885};
+ST(0.01841005291005176,0.01841005291005398,0,0.02021852793711804,0.01695330271573968,0,0.02044412351323598,0.01897469406756427,0){2544761.773764699,3500792.923489964,3970630.174131208};
+ST(0.03158994708994604,0.01841005291005184,0,0.03304669728426029,0.02021852793711807,0,0.03102530593243569,0.0204441235132358,0){9506981.828722859,9109191.0106201,9875926.345166853};
+ST(0.01841005291005399,0.03158994708994826,0,0.01695330271573897,0.02978147206288096,0,0.01897469406756322,0.02955587648676389,0){681163.8648190071,700227.5226125802,1189893.288594959};
+ST(0.0315899470899482,0.03158994708994603,0,0.02978147206288093,0.033046697284261,0,0.029555876486764,0.03102530593243676,0){2544761.773765826,1822017.581094539,2959555.02311246};
+ST(0.02062721121040539,0.02096060998006265,0,0.02215001361588054,0.02196095145324617,0,0.0200943405656333,0.02291758750152122,0){4115933.579576971,5196103.016175205,3442648.322572604};
+ST(0.0290393900199373,0.0206272112104053,0,0.02803904854675385,0.0221500136158806,0,0.02708241249847881,0.02009434056563312,0){9869245.814602781,9188072.738583822,9183174.245855754};
+ST(0.0209606099800624,0.02937278878959456,0,0.02192345518181211,0.02793885538913906,0,0.02291758750152112,0.02990565943436665,0){1834895.766186544,2773264.32237239,2269951.661359691};
+ST(0.02937278878959463,0.02903939001993757,0,0.02794496480441602,0.02806239939751506,0,0.02990565943436682,0.0270824124984789,0){4399739.940052692,5004465.210416885,6055036.612150463};
+ST(0.02219686411149813,0.02366027874564487,0,0.02215001361588054,0.02196095145324617,0,0.02417751621551469,0.023393463480524,0){4757486.064843449,5196103.016175205,6349753.383152546};
+ST(0.0238260250094999,0.02167653270546622,0,0.02417751621551469,0.023393463480524,0,0.02215001361588054,0.02196095145324617,0){6641234.082622026,6349753.383152546,5196103.016175205};
+ST(0.02633972125435513,0.02219686411149813,0,0.02803904854675385,0.0221500136158806,0,0.02646397875083406,0.02390163971325145,0){8334025.300326461,9188072.738583822,7578556.193563664};
+ST(0.02832346729453371,0.0238260250095,0,0.02646397875083406,0.02390163971325145,0,0.02803904854675385,0.0221500136158806,0){8398833.273041854,7578556.193563664,9188072.738583822};
+ST(0.015,0.02700000000001855,0,0.015,0.0250000000000237,0,0.01666349631902295,0.0263865043232969,0){645703.4689300079,820849.9862370425,1123094.041931162};
+ST(0.02700000000001855,0.035,0,0.0250000000000237,0.035,0,0.02638650432329687,0.03333650368097703,0){963276.382306002,820849.9862409331,1481993.704088334};
+ST(0.02299999999997874,0.015,0,0.02499999999997461,0.015,0,0.02361349567670159,0.01666349631902287,0){4771139.155196142,6065306.597110934,5949981.737036595};
+ST(0.035,0.02299999999997874,0,0.035,0.02499999999997461,0,0.0333365036809771,0.02361349567670162,0){7117703.227635174,6065306.597141732,7851377.66252339};
+ST(0.0192147217946296,0.02486333911250368,0,0.02004288562623822,0.02710123186687632,0,0.01836860548693934,0.02621782439543604,0){2466603.640555887,2240880.53857153,1756069.912092069};
+ST(0.0248633391125036,0.03078527820537037,0,0.02709387158569049,0.02994395049886663,0,0.02621744120711934,0.0316307091801319,0){2400098.990608844,3418862.023887427,2240667.074053982};
+ST(0.0251366608874963,0.01921472179462946,0,0.02292305321427705,0.02002627187241877,0,0.02380323230261221,0.01834737080814979,0){7845168.002306522,6060221.998352256,6627806.876270836};
+ST(0.03078527820537049,0.02513666088749638,0,0.02997372812758119,0.02292305321427708,0,0.03165262919185018,0.02380323230261227,0){7633646.319945103,9180996.867303077,8420148.806340775};
+ST(0.035,0.0189999999999931,0,0.03304669728426029,0.02021852793711807,0,0.0334005718855488,0.01859796559780545,0){7710515.858034595,9109191.0106201,8734574.207698569};
+ST(0.0189999999999931,0.015,0,0.02021852793711804,0.01695330271573968,0,0.01859796559780539,0.01659942811445112,0){2322362.747294063,3500792.923489964,2427549.509936762};
+ST(0.03100000000000718,0.035,0,0.02978147206288093,0.033046697284261,0,0.03140203440219497,0.03340057188554917,0){1043504.847547499,1822017.581094539,1627711.140556734};
+ST(0.015,0.03100000000000718,0,0.01695330271573897,0.02978147206288096,0,0.01659942811445075,0.03140203440219502,0){314297.6201831806,700227.5226125802,452380.309287892};
+ST(0.03099999999999255,0.015,0,0.03140203440219678,0.01659942811444605,0,0.02980597193963426,0.01698781140710807,0){7710515.858036811,8734574.207694963,9129179.508373139};
+ST(0.015,0.01900000000000812,0,0.01659942811444598,0.0185979655978033,0,0.01698781140710817,0.02019402806036576,0){1043504.847547819,1627711.140554541,1838661.946960948};
+ST(0.035,0.03099999999999255,0,0.03340057188555399,0.03140203440219681,0,0.03301218859289196,0.02980597193963428,0){2322362.747301393,2427549.509934699,3491325.400350199};
+ST(0.01900000000000812,0.035,0,0.0185979655978033,0.03340057188555402,0,0.02019402806036574,0.03301218859289186,0){314297.6201842385,452380.3092871403,703170.2194260015};
+ST(0.01802267861832114,0.02781420225278695,0,0.01897469406756322,0.02955587648676389,0,0.01695330271573897,0.02978147206288096,0){1293569.239124938,1189893.288594959,700227.5226125802};
+ST(0.02781420225278688,0.03197732138167883,0,0.029555876486764,0.03102530593243676,0,0.02978147206288093,0.033046697284261,0){2271057.768509638,2959555.02311246,1822017.581094539};
+ST(0.0221857977472131,0.0180226786183216,0,0.02044412351323598,0.01897469406756427,0,0.02021852793711804,0.01695330271573968,0){5221942.849585891,3970630.174131208,3500792.923489964};
+ST(0.03197732138167837,0.02218579774721317,0,0.03102530593243569,0.0204441235132358,0,0.03304669728426029,0.02021852793711807,0){9167915.807342248,9875926.345166853,9109191.0106201};
+ST(0.01841005291005176,0.01841005291005398,0,0.01900766867126057,0.02025952709593554,0,0.01698781140710817,0.02019402806036576,0){2544761.773764699,2984994.992984212,1838661.946960948};
+ST(0.03158994708994604,0.01841005291005184,0,0.02974047290406433,0.01900766867126047,0,0.02980597193963426,0.01698781140710807,0){9506981.828722859,9895343.918024337,9129179.508373139};
+ST(0.01841005291005399,0.03158994708994826,0,0.02025952709593563,0.03099233132873949,0,0.02019402806036574,0.03301218859289186,0){681163.8648190071,1156621.769609863,703170.2194260015};
+ST(0.0315899470899482,0.03158994708994603,0,0.03099233132873955,0.02974047290406432,0,0.03301218859289196,0.02980597193963428,0){2544761.773765826,3834234.301988425,3491325.400350199};
+ST(0.035,0.01699999999999357,0,0.035,0.0189999999999931,0,0.0334005718855488,0.01859796559780545,0){7117703.227623346,7710515.858034595,8734574.207698569};
+ST(0.01699999999999357,0.015,0,0.0189999999999931,0.015,0,0.01859796559780539,0.01659942811445112,0){1437039.497774626,2322362.747294063,2427549.509936762};
+ST(0.0330000000000061,0.035,0,0.03100000000000718,0.035,0,0.03140203440219497,0.03340057188554917,0){963276.3823045766,1043504.847547499,1627711.140556734};
+ST(0.015,0.0330000000000061,0,0.015,0.03100000000000718,0,0.01659942811445075,0.03140203440219502,0){194482.1474535455,314297.6201831806,452380.309287892};
+ST(0.03099999999999255,0.015,0,0.03299999999999222,0.015,0,0.03140203440219678,0.01659942811444605,0){7710515.858036811,7117703.227629418,8734574.207694963};
+ST(0.035,0.03099999999999255,0,0.035,0.03299999999999222,0,0.03340057188555399,0.03140203440219681,0){2322362.747301393,1437039.497779934,2427549.509934699};
+ST(0.015,0.01900000000000812,0,0.015,0.01700000000000744,0,0.01659942811444598,0.0185979655978033,0){1043504.847547819,963276.3823053603,1627711.140554541};
+ST(0.01900000000000812,0.035,0,0.01700000000000744,0.035,0,0.0185979655978033,0.03340057188555402,0){314297.6201842385,194482.14745423,452380.3092871403};
+ST(0.01718800098674322,0.02423878267049852,0,0.01836860548693934,0.02621782439543604,0,0.01666349631902295,0.0263865043232969,0){1618401.802394528,1756069.912092069,1123094.041931162};
+ST(0.02423878267049847,0.03281199901325677,0,0.02621744120711934,0.0316307091801319,0,0.02638650432329687,0.03333650368097703,0){1389849.759568753,2240667.074053982,1481993.704088334};
+ST(0.02576121732950113,0.01718800098674297,0,0.02380323230261221,0.01834737080814979,0,0.02361349567670159,0.01666349631902287,0){7720189.599155798,6627806.876270836,5949981.737036595};
+ST(0.03281199901325701,0.02576121732950117,0,0.03165262919185018,0.02380323230261227,0,0.0333365036809771,0.02361349567670162,0){6629938.030430425,8420148.806340775,7851377.66252339};
+ST(0.035,0.01699999999999357,0,0.0334005718855488,0.01859796559780545,0,0.03300243786989122,0.01699756213010876,0){7117703.227623346,8734574.207698569,8350257.942833206};
+ST(0.01699999999999357,0.015,0,0.01859796559780539,0.01659942811445112,0,0.01699756213010875,0.01699756213010876,0){1437039.497774626,2427549.509936762,1685066.208253873};
+ST(0.0330000000000061,0.035,0,0.03140203440219497,0.03340057188554917,0,0.03300243786989124,0.03300243786989122,0){963276.3823045766,1627711.140556734,1685066.208253887};
+ST(0.015,0.0330000000000061,0,0.01659942811445075,0.03140203440219502,0,0.01699756213010876,0.03300243786989126,0){194482.1474535455,452380.309287892,340043.1634134258};
+ST(0.03299999999999222,0.015,0,0.03300243786989122,0.01699756213010876,0,0.03140203440219678,0.01659942811444605,0){7117703.227629418,8350257.942833206,8734574.207694963};
+ST(0.015,0.01700000000000744,0,0.01699756213010875,0.01699756213010876,0,0.01659942811444598,0.0185979655978033,0){963276.3823053603,1685066.208253873,1627711.140554541};
+ST(0.035,0.03299999999999222,0,0.03300243786989124,0.03300243786989122,0,0.03340057188555399,0.03140203440219681,0){1437039.497779934,1685066.208253887,2427549.509934699};
+ST(0.01700000000000744,0.035,0,0.01699756213010876,0.03300243786989126,0,0.0185979655978033,0.03340057188555402,0){194482.14745423,340043.1634134258,452380.3092871403};
+ST(0.01841005291005176,0.01841005291005398,0,0.02044412351323598,0.01897469406756427,0,0.01900766867126057,0.02025952709593554,0){2544761.773764699,3970630.174131208,2984994.992984212};
+ST(0.03158994708994604,0.01841005291005184,0,0.03102530593243569,0.0204441235132358,0,0.02974047290406433,0.01900766867126047,0){9506981.828722859,9875926.345166853,9895343.918024337};
+ST(0.01841005291005399,0.03158994708994826,0,0.01897469406756322,0.02955587648676389,0,0.02025952709593563,0.03099233132873949,0){681163.8648190071,1189893.288594959,1156621.769609863};
+ST(0.0315899470899482,0.03158994708994603,0,0.029555876486764,0.03102530593243676,0,0.03099233132873955,0.02974047290406432,0){2544761.773765826,2959555.02311246,3834234.301988425};
+ST(0.02219686411149813,0.02366027874564487,0,0.02102282424207957,0.02483178028678816,0,0.0200943405656333,0.02291758750152122,0){4757486.064843449,3536817.246116393,3442648.322572604};
+ST(0.02633972125435513,0.02219686411149813,0,0.0251688693519739,0.02080806878318098,0,0.02708241249847881,0.02009434056563312,0){8334025.300326461,7866845.121381707,9183174.245855754};
+ST(0.02355978257831163,0.02772127893890192,0,0.02481685402985528,0.02896618798686593,0,0.02291758750152112,0.02990565943436665,0){3638747.773110651,3421263.069148063,2269951.661359691};
+ST(0.02787134264717452,0.02636461712411311,0,0.02919889109015292,0.02517140974684725,0,0.02990565943436682,0.0270824124984789,0){6373761.262323579,7604442.001547791,6055036.612150463};
+ST(0.0192147217946296,0.02486333911250368,0,0.01718800098674322,0.02423878267049852,0,0.01865207083058131,0.02351995976535901,0){2466603.640555887,1618401.802394528,2437392.843189596};
+ST(0.0248633391125036,0.03078527820537037,0,0.02423878267049847,0.03281199901325677,0,0.02351995976535895,0.0313479291694187,0){2400098.990608844,1389849.759568753,1812887.565675998};
+ST(0.0251366608874963,0.01921472179462946,0,0.02576121732950113,0.01718800098674297,0,0.02648004023464091,0.0186520708305811,0){7845168.002306522,7720189.599155798,8675602.666294761};
+ST(0.03078527820537049,0.02513666088749638,0,0.03281199901325701,0.02576121732950117,0,0.03134792916941886,0.02648004023464098,0){7633646.319945103,6629938.030430425,6452752.268645341};
+ST(0.01802267861831516,0.02218579774721222,0,0.01865207083058131,0.02351995976535901,0,0.01718800098674322,0.02423878267049852,0){2271057.768506458,2437392.843189596,1618401.802394528};
+ST(0.02218579774721212,0.03197732138168483,0,0.02351995976535895,0.0313479291694187,0,0.02423878267049847,0.03281199901325677,0){1293569.239122899,1812887.565675998,1389849.759568753};
+ST(0.02781420225278781,0.01802267861831496,0,0.02648004023464091,0.0186520708305811,0,0.02576121732950113,0.01718800098674297,0){9167915.807340225,8675602.666294761,7720189.599155798};
+ST(0.03197732138168502,0.0278142022527879,0,0.03134792916941886,0.02648004023464098,0,0.03281199901325701,0.02576121732950117,0){5221942.849583705,6452752.268645341,6629938.030430425};
+ST(0.015,0.02900000000000975,0,0.01632789553061661,0.02799643572779862,0,0.01695330271573897,0.02978147206288096,0){468876.952199062,813744.2457889724,700227.5226125802};
+ST(0.02900000000000975,0.035,0,0.0279964357277986,0.03367210446938337,0,0.02978147206288093,0.033046697284261,0){1043504.847547853,1481682.088749314,1822017.581094539};
+ST(0.01802267861831516,0.02218579774721222,0,0.01643969820243331,0.02212372169562202,0,0.01698781140710817,0.02019402806036576,0){2271057.768506458,1519931.557947947,1838661.946960948};
+ST(0.02781420225278781,0.01802267861831496,0,0.02787627830437758,0.0164396982024332,0,0.02980597193963426,0.01698781140710807,0){9167915.807340225,8420975.407718396,9129179.508373139};
+ST(0.02218579774721212,0.03197732138168483,0,0.02212372169562198,0.03356030179756669,0,0.02019402806036574,0.03301218859289186,0){1293569.239122899,855054.256044248,703170.2194260015};
+ST(0.03197732138168502,0.0278142022527879,0,0.0335603017975668,0.02787627830437762,0,0.03301218859289196,0.02980597193963428,0){5221942.849583705,4737312.561715177,3491325.400350199};
+ST(0.01802267861832114,0.02781420225278695,0,0.01695330271573897,0.02978147206288096,0,0.01632789553061661,0.02799643572779862,0){1293569.239124938,700227.5226125802,813744.2457889724};
+ST(0.02781420225278688,0.03197732138167883,0,0.02978147206288093,0.033046697284261,0,0.0279964357277986,0.03367210446938337,0){2271057.768509638,1822017.581094539,1481682.088749314};
+ST(0.02099999999998355,0.015,0,0.022003564272199,0.01632789553061683,0,0.02021852793711804,0.01695330271573968,0){3464558.103290316,4610386.703738825,3500792.923489964};
+ST(0.035,0.02099999999998355,0,0.03367210446938315,0.02200356427219903,0,0.03304669728426029,0.02021852793711807,0){7710515.858038196,8394685.967355566,9109191.0106201};
+ST(0.015,0.02100000000001463,0,0.01698781140710817,0.02019402806036576,0,0.01643969820243331,0.02212372169562202,0){1043504.847547345,1838661.946960948,1519931.557947947};
+ST(0.02100000000001463,0.035,0,0.02019402806036574,0.03301218859289186,0,0.02212372169562198,0.03356030179756669,0){468876.9522011194,703170.2194260015,855054.256044248};
+ST(0.02899999999998503,0.015,0,0.02980597193963426,0.01698781140710807,0,0.02787627830437758,0.0164396982024332,0){7710515.858033354,9129179.508373139,8420975.407718396};
+ST(0.035,0.02899999999998503,0,0.03301218859289196,0.02980597193963428,0,0.0335603017975668,0.02787627830437762,0){3464558.103309908,3491325.400350199,4737312.561715177};
+ST(0.0221857977472131,0.0180226786183216,0,0.02021852793711804,0.01695330271573968,0,0.022003564272199,0.01632789553061683,0){5221942.849585891,3500792.923489964,4610386.703738825};
+ST(0.03197732138167837,0.02218579774721317,0,0.03304669728426029,0.02021852793711807,0,0.03367210446938315,0.02200356427219903,0){9167915.807342248,9109191.0106201,8394685.967355566};
+ST(0.0192147217946296,0.02486333911250368,0,0.0200943405656333,0.02291758750152122,0,0.02102282424207957,0.02483178028678816,0){2466603.640555887,3442648.322572604,3536817.246116393};
+ST(0.0248633391125036,0.03078527820537037,0,0.02291758750152112,0.02990565943436665,0,0.02481685402985528,0.02896618798686593,0){2400098.990608844,2269951.661359691,3421263.069148063};
+ST(0.0251366608874963,0.01921472179462946,0,0.02708241249847881,0.02009434056563312,0,0.0251688693519739,0.02080806878318098,0){7845168.002306522,9183174.245855754,7866845.121381707};
+ST(0.03078527820537049,0.02513666088749638,0,0.02990565943436682,0.0270824124984789,0,0.02919889109015292,0.02517140974684725,0){7633646.319945103,6055036.612150463,7604442.001547791};
+ST(0.02219686411149813,0.02366027874564487,0,0.02321297569192563,0.02525629976722429,0,0.02102282424207957,0.02483178028678816,0){4757486.064843449,4785838.434520846,3536817.246116393};
+ST(0.0218094427360229,0.0263682556397269,0,0.02102282424207957,0.02483178028678816,0,0.02321297569192563,0.02525629976722429,0){3408205.856754326,3536817.246116393,4785838.434520846};
+ST(0.02355978257831163,0.02772127893890192,0,0.0251772925710883,0.02680645346909019,0,0.02481685402985528,0.02896618798686593,0){3638747.773110651,4986438.53885695,3421263.069148063};
+ST(0.02630937339024055,0.02808524626481607,0,0.02481685402985528,0.02896618798686593,0,0.0251772925710883,0.02680645346909019,0){4538814.13250844,3421263.069148063,4986438.53885695};
+ST(0.03158994708994604,0.01841005291005184,0,0.0334005718855488,0.01859796559780545,0,0.03304669728426029,0.02021852793711807,0){9506981.828722859,8734574.207698569,9109191.0106201};
+ST(0.01841005291005176,0.01841005291005398,0,0.01859796559780539,0.01659942811445112,0,0.02021852793711804,0.01695330271573968,0){2544761.773764699,2427549.509936762,3500792.923489964};
+ST(0.0315899470899482,0.03158994708994603,0,0.03140203440219497,0.03340057188554917,0,0.02978147206288093,0.033046697284261,0){2544761.773765826,1627711.140556734,1822017.581094539};
+ST(0.01841005291005399,0.03158994708994826,0,0.01659942811445075,0.03140203440219502,0,0.01695330271573897,0.02978147206288096,0){681163.8648190071,452380.309287892,700227.5226125802};
+ST(0.03158994708994604,0.01841005291005184,0,0.02980597193963426,0.01698781140710807,0,0.03140203440219678,0.01659942811444605,0){9506981.828722859,9129179.508373139,8734574.207694963};
+ST(0.01841005291005176,0.01841005291005398,0,0.01698781140710817,0.02019402806036576,0,0.01659942811444598,0.0185979655978033,0){2544761.773764699,1838661.946960948,1627711.140554541};
+ST(0.0315899470899482,0.03158994708994603,0,0.03301218859289196,0.02980597193963428,0,0.03340057188555399,0.03140203440219681,0){2544761.773765826,3491325.400350199,2427549.509934699};
+ST(0.01841005291005399,0.03158994708994826,0,0.02019402806036574,0.03301218859289186,0,0.0185979655978033,0.03340057188555402,0){681163.8648190071,703170.2194260015,452380.3092871403};
+ST(0.0218094427360229,0.0263682556397269,0,0.02192345518181211,0.02793885538913906,0,0.02004288562623822,0.02710123186687632,0){3408205.856754326,2773264.32237239,2240880.53857153};
+ST(0.0209606099800624,0.02937278878959456,0,0.02004288562623822,0.02710123186687632,0,0.02192345518181211,0.02793885538913906,0){1834895.766186544,2240880.53857153,2773264.32237239};
+ST(0.02630937339024055,0.02808524626481607,0,0.02794496480441602,0.02806239939751506,0,0.02709387158569049,0.02994395049886663,0){4538814.13250844,5004465.210416885,3418862.023887427};
+ST(0.0238260250094999,0.02167653270546622,0,0.02215001361588054,0.02196095145324617,0,0.02292305321427705,0.02002627187241877,0){6641234.082622026,5196103.016175205,6060221.998352256};
+ST(0.02937278878959463,0.02903939001993757,0,0.02709387158569049,0.02994395049886663,0,0.02794496480441602,0.02806239939751506,0){4399739.940052692,3418862.023887427,5004465.210416885};
+ST(0.02062721121040539,0.02096060998006265,0,0.02292305321427705,0.02002627187241877,0,0.02215001361588054,0.02196095145324617,0){4115933.579576971,6060221.998352256,5196103.016175205};
+ST(0.02832346729453371,0.0238260250095,0,0.02803904854675385,0.0221500136158806,0,0.02997372812758119,0.02292305321427708,0){8398833.273041854,9188072.738583822,9180996.867303077};
+ST(0.0290393900199373,0.0206272112104053,0,0.02997372812758119,0.02292305321427708,0,0.02803904854675385,0.0221500136158806,0){9869245.814602781,9180996.867303077,9188072.738583822};
+ST(0.0209606099800624,0.02937278878959456,0,0.02025952709593563,0.03099233132873949,0,0.01897469406756322,0.02955587648676389,0){1834895.766186544,1156621.769609863,1189893.288594959};
+ST(0.02062721121040539,0.02096060998006265,0,0.01900766867126057,0.02025952709593554,0,0.02044412351323598,0.01897469406756427,0){4115933.579576971,2984994.992984212,3970630.174131208};
+ST(0.02937278878959463,0.02903939001993757,0,0.03099233132873955,0.02974047290406432,0,0.029555876486764,0.03102530593243676,0){4399739.940052692,3834234.301988425,2959555.02311246};
+ST(0.0290393900199373,0.0206272112104053,0,0.02974047290406433,0.01900766867126047,0,0.03102530593243569,0.0204441235132358,0){9869245.814602781,9895343.918024337,9875926.345166853};
+ST(0.0192147217946296,0.02486333911250368,0,0.01865207083058131,0.02351995976535901,0,0.0200943405656333,0.02291758750152122,0){2466603.640555887,2437392.843189596,3442648.322572604};
+ST(0.0248633391125036,0.03078527820537037,0,0.02351995976535895,0.0313479291694187,0,0.02291758750152112,0.02990565943436665,0){2400098.990608844,1812887.565675998,2269951.661359691};
+ST(0.0251366608874963,0.01921472179462946,0,0.02648004023464091,0.0186520708305811,0,0.02708241249847881,0.02009434056563312,0){7845168.002306522,8675602.666294761,9183174.245855754};
+ST(0.03078527820537049,0.02513666088749638,0,0.03134792916941886,0.02648004023464098,0,0.02990565943436682,0.0270824124984789,0){7633646.319945103,6452752.268645341,6055036.612150463};
+ST(0.02062721121040539,0.02096060998006265,0,0.0200943405656333,0.02291758750152122,0,0.01937302306416,0.02166131996144063,0){4115933.579576971,3442648.322572604,3144512.483378664};
+ST(0.0290393900199373,0.0206272112104053,0,0.02708241249847881,0.02009434056563312,0,0.02833868003855931,0.01937302306415976,0){9869245.814602781,9183174.245855754,9689610.723385418};
+ST(0.0209606099800624,0.02937278878959456,0,0.02291758750152112,0.02990565943436665,0,0.02166131996144063,0.03062697693584002,0){1834895.766186544,2269951.661359691,1612721.068750821};
+ST(0.02937278878959463,0.02903939001993757,0,0.02990565943436682,0.0270824124984789,0,0.03062697693584022,0.02833868003855934,0){4399739.940052692,6055036.612150463,4969495.094771407};
+ST(0.01802267861831516,0.02218579774721222,0,0.01900766867126057,0.02025952709593554,0,0.01937302306416,0.02166131996144063,0){2271057.768506458,2984994.992984212,3144512.483378664};
+ST(0.02781420225278781,0.01802267861831496,0,0.02974047290406433,0.01900766867126047,0,0.02833868003855931,0.01937302306415976,0){9167915.807340225,9895343.918024337,9689610.723385418};
+ST(0.02218579774721212,0.03197732138168483,0,0.02025952709593563,0.03099233132873949,0,0.02166131996144063,0.03062697693584002,0){1293569.239122899,1156621.769609863,1612721.068750821};
+ST(0.03197732138168502,0.0278142022527879,0,0.03099233132873955,0.02974047290406432,0,0.03062697693584022,0.02833868003855934,0){5221942.849583705,3834234.301988425,4969495.094771407};
+ST(0.01802267861832114,0.02781420225278695,0,0.01836860548693934,0.02621782439543604,0,0.02004288562623822,0.02710123186687632,0){1293569.239124938,1756069.912092069,2240880.53857153};
+ST(0.02781420225278688,0.03197732138167883,0,0.02621744120711934,0.0316307091801319,0,0.02709387158569049,0.02994395049886663,0){2271057.768509638,2240667.074053982,3418862.023887427};
+ST(0.0221857977472131,0.0180226786183216,0,0.02380323230261221,0.01834737080814979,0,0.02292305321427705,0.02002627187241877,0){5221942.849585891,6627806.876270836,6060221.998352256};
+ST(0.03197732138167837,0.02218579774721317,0,0.03165262919185018,0.02380323230261227,0,0.02997372812758119,0.02292305321427708,0){9167915.807342248,8420148.806340775,9180996.867303077};
+ST(0.0200943405656333,0.02291758750152122,0,0.01865207083058131,0.02351995976535901,0,0.01937302306416,0.02166131996144063,0){3442648.322572604,2437392.843189596,3144512.483378664};
+ST(0.02291758750152112,0.02990565943436665,0,0.02351995976535895,0.0313479291694187,0,0.02166131996144063,0.03062697693584002,0){2269951.661359691,1812887.565675998,1612721.068750821};
+ST(0.02708241249847881,0.02009434056563312,0,0.02648004023464091,0.0186520708305811,0,0.02833868003855931,0.01937302306415976,0){9183174.245855754,8675602.666294761,9689610.723385418};
+ST(0.02990565943436682,0.0270824124984789,0,0.03134792916941886,0.02648004023464098,0,0.03062697693584022,0.02833868003855934,0){6055036.612150463,6452752.268645341,4969495.094771407};
+ST(0.03158994708994604,0.01841005291005184,0,0.03300243786989122,0.01699756213010876,0,0.0334005718855488,0.01859796559780545,0){9506981.828722859,8350257.942833206,8734574.207698569};
+ST(0.01841005291005176,0.01841005291005398,0,0.01699756213010875,0.01699756213010876,0,0.01859796559780539,0.01659942811445112,0){2544761.773764699,1685066.208253873,2427549.509936762};
+ST(0.03158994708994604,0.01841005291005184,0,0.03140203440219678,0.01659942811444605,0,0.03300243786989122,0.01699756213010876,0){9506981.828722859,8734574.207694963,8350257.942833206};
+ST(0.01841005291005176,0.01841005291005398,0,0.01659942811444598,0.0185979655978033,0,0.01699756213010875,0.01699756213010876,0){2544761.773764699,1627711.140554541,1685066.208253873};
+ST(0.0315899470899482,0.03158994708994603,0,0.03300243786989124,0.03300243786989122,0,0.03140203440219497,0.03340057188554917,0){2544761.773765826,1685066.208253887,1627711.140556734};
+ST(0.01841005291005399,0.03158994708994826,0,0.01699756213010876,0.03300243786989126,0,0.01659942811445075,0.03140203440219502,0){681163.8648190071,340043.1634134258,452380.309287892};
+ST(0.0315899470899482,0.03158994708994603,0,0.03340057188555399,0.03140203440219681,0,0.03300243786989124,0.03300243786989122,0){2544761.773765826,2427549.509934699,1685066.208253887};
+ST(0.01841005291005399,0.03158994708994826,0,0.0185979655978033,0.03340057188555402,0,0.01699756213010876,0.03300243786989126,0){681163.8648190071,452380.3092871403,340043.1634134258};
+ST(0.01802267861831516,0.02218579774721222,0,0.01937302306416,0.02166131996144063,0,0.01865207083058131,0.02351995976535901,0){2271057.768506458,3144512.483378664,2437392.843189596};
+ST(0.02218579774721212,0.03197732138168483,0,0.02166131996144063,0.03062697693584002,0,0.02351995976535895,0.0313479291694187,0){1293569.239122899,1612721.068750821,1812887.565675998};
+ST(0.02781420225278781,0.01802267861831496,0,0.02833868003855931,0.01937302306415976,0,0.02648004023464091,0.0186520708305811,0){9167915.807340225,9689610.723385418,8675602.666294761};
+ST(0.03197732138168502,0.0278142022527879,0,0.03062697693584022,0.02833868003855934,0,0.03134792916941886,0.02648004023464098,0){5221942.849583705,4969495.094771407,6452752.268645341};
+ST(0.025,0.025,0,0.02321297569192563,0.02525629976722429,0,0.02417751621551469,0.023393463480524,0){6065306.597126335,4785838.434520846,6349753.383152546};
+ST(0.02219686411149813,0.02366027874564487,0,0.02417751621551469,0.023393463480524,0,0.02321297569192563,0.02525629976722429,0){4757486.064843449,6349753.383152546,4785838.434520846};
+ST(0.015,0.02700000000001855,0,0.01666349631902295,0.0263865043232969,0,0.01632789553061661,0.02799643572779862,0){645703.4689300079,1123094.041931162,813744.2457889724};
+ST(0.02700000000001855,0.035,0,0.02638650432329687,0.03333650368097703,0,0.0279964357277986,0.03367210446938337,0){963276.382306002,1481993.704088334,1481682.088749314};
+ST(0.02299999999997874,0.015,0,0.02361349567670159,0.01666349631902287,0,0.022003564272199,0.01632789553061683,0){4771139.155196142,5949981.737036595,4610386.703738825};
+ST(0.035,0.02299999999997874,0,0.0333365036809771,0.02361349567670162,0,0.03367210446938315,0.02200356427219903,0){7117703.227635174,7851377.66252339,8394685.967355566};
+ST(0.0221857977472131,0.0180226786183216,0,0.022003564272199,0.01632789553061683,0,0.02361349567670159,0.01666349631902287,0){5221942.849585891,4610386.703738825,5949981.737036595};
+ST(0.01802267861832114,0.02781420225278695,0,0.01632789553061661,0.02799643572779862,0,0.01666349631902295,0.0263865043232969,0){1293569.239124938,813744.2457889724,1123094.041931162};
+ST(0.02781420225278688,0.03197732138167883,0,0.0279964357277986,0.03367210446938337,0,0.02638650432329687,0.03333650368097703,0){2271057.768509638,1481682.088749314,1481993.704088334};
+ST(0.03197732138167837,0.02218579774721317,0,0.03367210446938315,0.02200356427219903,0,0.0333365036809771,0.02361349567670162,0){9167915.807342248,8394685.967355566,7851377.66252339};
+ST(0.02099999999998355,0.015,0,0.02299999999997874,0.015,0,0.022003564272199,0.01632789553061683,0){3464558.103290316,4771139.155196142,4610386.703738825};
+ST(0.035,0.02099999999998355,0,0.035,0.02299999999997874,0,0.03367210446938315,0.02200356427219903,0){7710515.858038196,7117703.227635174,8394685.967355566};
+ST(0.015,0.02900000000000975,0,0.015,0.02700000000001855,0,0.01632789553061661,0.02799643572779862,0){468876.952199062,645703.4689300079,813744.2457889724};
+ST(0.02900000000000975,0.035,0,0.02700000000001855,0.035,0,0.0279964357277986,0.03367210446938337,0){1043504.847547853,963276.382306002,1481682.088749314};
+ST(0.015,0.02300000000001896,0,0.015,0.02100000000001463,0,0.01643969820243331,0.02212372169562202,0){963276.3823038344,1043504.847547345,1519931.557947947};
+ST(0.02300000000001896,0.035,0,0.02100000000001463,0.035,0,0.02212372169562198,0.03356030179756669,0){645703.4689333984,468876.9522011194,855054.256044248};
+ST(0.02699999999997969,0.015,0,0.02899999999998503,0.015,0,0.02787627830437758,0.0164396982024332,0){7117703.227617423,7710515.858033354,8420975.407718396};
+ST(0.035,0.02699999999997969,0,0.035,0.02899999999998503,0,0.0335603017975668,0.02787627830437762,0){4771139.155223909,3464558.103309908,4737312.561715177};
+ST(0.01802267861832114,0.02781420225278695,0,0.01666349631902295,0.0263865043232969,0,0.01836860548693934,0.02621782439543604,0){1293569.239124938,1123094.041931162,1756069.912092069};
+ST(0.02781420225278688,0.03197732138167883,0,0.02638650432329687,0.03333650368097703,0,0.02621744120711934,0.0316307091801319,0){2271057.768509638,1481993.704088334,2240667.074053982};
+ST(0.0221857977472131,0.0180226786183216,0,0.02361349567670159,0.01666349631902287,0,0.02380323230261221,0.01834737080814979,0){5221942.849585891,5949981.737036595,6627806.876270836};
+ST(0.03197732138167837,0.02218579774721317,0,0.0333365036809771,0.02361349567670162,0,0.03165262919185018,0.02380323230261227,0){9167915.807342248,7851377.66252339,8420148.806340775};
+ST(0.02062721121040539,0.02096060998006265,0,0.01937302306416,0.02166131996144063,0,0.01900766867126057,0.02025952709593554,0){4115933.579576971,3144512.483378664,2984994.992984212};
+ST(0.0290393900199373,0.0206272112104053,0,0.02833868003855931,0.01937302306415976,0,0.02974047290406433,0.01900766867126047,0){9869245.814602781,9689610.723385418,9895343.918024337};
+ST(0.0209606099800624,0.02937278878959456,0,0.02166131996144063,0.03062697693584002,0,0.02025952709593563,0.03099233132873949,0){1834895.766186544,1612721.068750821,1156621.769609863};
+ST(0.02937278878959463,0.02903939001993757,0,0.03062697693584022,0.02833868003855934,0,0.03099233132873955,0.02974047290406432,0){4399739.940052692,4969495.094771407,3834234.301988425};
+ST(0.025,0.025,0,0.02646397875083406,0.02390163971325145,0,0.02634916017313363,0.02577951485829298,0){6065306.597126335,7578556.193563664,6266831.340785229};
+ST(0.025,0.025,0,0.02634916017313363,0.02577951485829298,0,0.0251772925710883,0.02680645346909019,0){6065306.597126335,6266831.340785229,4986438.53885695};
+};
diff --git a/Thermics/Qsource.opt b/Thermics/Qsource.opt
new file mode 100644
index 0000000000000000000000000000000000000000..5682ccd0dc9c3e5f89977f48bf70fb328350847a
--- /dev/null
+++ b/Thermics/Qsource.opt
@@ -0,0 +1,11 @@
+Plugin(NewView).View=-1;
+Plugin(NewView).Dimension=1;
+Plugin(NewView).Run;
+
+Plugin(MathEval).Expression0="1e7*Exp(-1e4*((x-0.03)^2+(y-0.02)^2))";
+Plugin(MathEval).PhysicalRegion=111;
+Plugin(MathEval).Run;
+
+View[1].NormalRaise = 4e-05;
+
+Save View[1] "Q.pos";
diff --git a/Thermics/brick.geo b/Thermics/brick.geo
new file mode 100644
index 0000000000000000000000000000000000000000..0dbaa0e01655ccbecefe9b4c3cdf1b951b2c9265
--- /dev/null
+++ b/Thermics/brick.geo
@@ -0,0 +1,111 @@
+/* -------------------------------------------------------------------
+ File "brick.geo"
+
+ This file is the geometrical description used by GMSH to produce
+ the file "brick.msh".
+ ------------------------------------------------------------------- */
+
+/* Gmsh options can be set directly in the .geo file.
+ Setting "Solver.AutoMesh" to 2 ensures that GMSH systematically
+ regenerates the mesh at each model execution, and does not reuse
+ the mesh on disk if it exists (which is the default option).
+ This option is needed in this model, because the interactive
+ model parameters "Flag_Regularization" has an effect
+ on the definition of the regions, which makes remeshing mandatory. */
+
+Solver.AutoMesh = 2;
+
+Include "brick_common.pro";
+
+/* The (very) simple geometry of this thermal model only contains
+ rectangles. It is the opportunity to illustrate
+ the function definition capabilities of Gmsh. */
+
+Function Def_Rectangle
+p1=newp; Point(newp)={xc_, yc_, 0, lc_};
+p2=newp; Point(newp)={xc_+dx_, yc_, 0, lc_};
+p3=newp; Point(newp)={xc_+dx_, yc_+dy_, 0, lc_};
+p4=newp; Point(newp)={xc_, yc_+dy_, 0, lc_};
+
+l1=newl; Line(newl)={p1,p2};
+l2=newl; Line(newl)={p2,p3};
+l3=newl; Line(newl)={p3,p4};
+l4=newl; Line(newl)={p4,p1};
+
+lines_[] = {l1, l2, l3, l4};
+
+ll1 = newll; Line Loop(ll1) = {lines_[]};
+s_ = news; Plane Surface(news) = {ll1, ll_Holes_[]};
+Return // end of Function Def_Rectangle
+
+/* Note that the code above is not parsed before it is called.
+ Good practice for readability (but by no means mandatory)
+ is to distinguish function variables from other variables.
+ This is done here with a trailing "_" in the variable name. */
+
+
+// Window 1 (with optional layer of thickness "e_layer")
+dx_ = dx_Win1;
+dy_ = dy_Win1;
+lc_ = lc_Win1;
+xc_ = xc_Win1-dx_/2; yc_=yc_Win1-dy_/2;
+ll_Holes_[]={};
+Call Def_Rectangle;
+s_Win1 = s_;
+l_Win1[] = lines_[];
+
+If( !Flag_Regularization )
+ ll_HolesForBrick[] += ll1;
+Else
+ ll_HolesForLayer[] += ll1;
+ dx_ = dx_Win1+2*e_layer;
+ dy_ = dy_Win1+2*e_layer;
+ lc_ = lc_Win1;
+ xc_ = xc_Win1-dx_/2; yc_=yc_Win1-dy_/2;
+
+ ll_Holes_[] = {ll_HolesForLayer[]};
+ Call Def_Rectangle;
+ s_Win1_Layer = s_;
+ l_Win1_Layer[] = lines_[];
+ ll_HolesForBrick[] += ll1;
+EndIf
+
+// Window 2
+dx_ = dx_Win2;
+dy_ = dy_Win2;
+lc_ = lc_Win2;
+xc_ = xc_Win2-dx_/2; yc_=yc_Win2-dy_/2;
+ll_Holes_[]={};
+Call Def_Rectangle;
+s_Win2 = s_;
+l_Win2[] = lines_[];
+ll_HolesForBrick[] += ll1;
+
+// Brick
+dx_ = dx_Brick; dy_ = dy_Brick;
+lc_ = lc_Brick;
+xc_ = 0; yc_ = 0;
+ll_Holes_[] = {ll_HolesForBrick[]};
+Call Def_Rectangle;
+s_Brick = s_;
+l_Brick[] = lines_[];
+
+//Printf("l_Brick", l_Brick[]);
+
+
+// Physical regions
+
+Physical Surface("Brick", 100) = {s_Brick};
+Physical Surface("Window1", 111) = {s_Win1};
+Physical Surface("Window2", 112) = {s_Win2};
+If( Flag_Regularization )
+ Physical Surface("LayerWindow1", 115) = {s_Win1_Layer};
+EndIf
+
+Physical Line("Surface1", 201) = { l_Brick[{3}] };
+Physical Line("Surface2", 202) = { l_Brick[{1}] };
+Physical Line("Surface3", 203) = { l_Brick[{0}], l_Brick[{2}] };
+
+Physical Line("SurfWindow1", 211) = {l_Win1[]};
+Physical Line("SurfWindow2", 212) = {l_Win2[]};
+
diff --git a/Thermics/brick.pro b/Thermics/brick.pro
new file mode 100644
index 0000000000000000000000000000000000000000..88694acd9155e4567fdd91a481162e2599eb0f03
--- /dev/null
+++ b/Thermics/brick.pro
@@ -0,0 +1,454 @@
+/* -------------------------------------------------------------------
+ Tutorial 5 : thermal problem with contact resistances
+
+ Features:
+ - Contact resistances: scalar FunctionSpace with a surface discontinuity
+ - Region with a uniform temperature (infinite thermal conductivity)
+ - Computation of the heat flux through surfaces
+ - Import a source field from a file
+
+ To compute the solution in a terminal:
+ getdp brick.pro -solve Thermal_T -pos Map_T
+
+ To compute the solution interactively from the Gmsh GUI:
+ File > Open > brick.pro
+ Run (button at the bottom of the left panel)
+ ------------------------------------------------------------------- */
+
+/* This model is a rectangular brick with two windows,
+ where various kinds of thermal constraints can be set.
+ Dirichlet, Neumann and convection boundary conditions are imposed
+ on different parts of the surface of the brick.
+ The model is rather academic but it
+ demonstrates some useful high-level GetDP features.
+
+ Governing eqations are
+
+ div ( -lambda[] grad T ) = Q in Vol_Lambda_The
+ -lambda[] grad T . n = qn = 0 on Sur_Neumann_The
+
+ Contact thermal resistance:
+ First, it is shown how to implement contact thermal resistances
+ with surface elements. The surface elements are associated with a
+ thickness and a thermal conductivity (typically much lower than that
+ of surrounding regions). The implementation takes advantage
+ of the powerful FunctionSpace definition in GetDP.
+
+ With the flag "Flag_Regularization", the contact surface
+ can be, for the sake of comparison, replaced
+ by a thin volume conducting region.
+
+ Thermal "electrode":
+ The floating potential idea (introduced in tutorial 4)
+ is reconsidered here in a thermal context to represent a region
+ with a very large thermal conductivity where, consequently,
+ the temperature field is uniform (exactly like the electric potential
+ is uniform on an electrode).
+ The dual quantity of this uniform temperature "T_electrode" [K]
+ (which is the "associated global quantity" in GetDP language)
+ is the heat flux "Q_electrode" [W] injected in the electrode
+ by the agent that maintains the temperature equal
+ to the prescribed value.
+
+ The value of Q_electrode is a by-product of the system resolution
+ provided the term
+
+ GlobalTerm { [-Dof{Q_electrode} , {T_electrode} ] ; In Tfloating_The ; }
+
+ is present in the "Resolution" section. This term triggers the writing
+ in the linear system of a supplementary equation associated
+ with the global basis function BF{T_electrode}.
+ All integrations are automatically done by getDP,
+ and the value of Q_electrode is obtained in postprocessing
+ with the PostOperation
+
+ Print[ Q_electrode, OnRegion Tfloating_The, ... ]
+
+ Heat flux through surfaces:
+ The purpose of a thermal simulation usually goes beyond
+ the mere calculation of a temperature distribution.
+ One is in general also interested in evaluating
+ the heat flux q(S) through some specific surface S:
+
+ q(S) = ( -lambda[] grad T . n )_S
+
+ This quantity cannot be computed from the temperature
+ distribution available on the surface S only.
+ As heat flux is related with the gradient of temperature
+ in the direction normal to the surface, its computation relies on
+ the temperature distribution in a neighborhood of the surface.
+ This means that volume elements in contact with the considered surface
+ need be involved in the computation.
+ To achieve this with getDP, a good method proceeds
+ by the definition a smooth auxiliary function g(S),
+ with g(S)=1 on S, and g(S)=0 outside a finite neighborhood of S.
+ Typically, g(S) is the sum of the shape functions of the nodes on S.
+ Let w(S) be the support of g(S),
+ and let dw(S) denote the boundary of w(S).
+ We then have, just adding and substracting dw(S)
+ to the surface of integration S
+
+ q(S) = ( -lambda[] grad T . n g(S) )_{ dw(S) - ( dw(S)-S ) }.
+
+ dw(S) being a boundary, Stokes theorem can be invoked and,
+ after an integration by part one ends up with
+
+ q(S) = ( -lambda[] grad T . grad g(S) )_w(S)
+ - ( -lambda[] grad T . n g(S) )_{dw(S)-S}.
+ + ( Q g(S) )_w(S)
+
+ Now, g(S) is zero on {dw(S)-S}, except maybe at some surface elements
+ adjacents to dS, but not in dS.
+ The second terme vanishes then if either S is closed,
+ or adjacent to a homogeneous Neumann boundary condition.
+
+ The third term also vanishes, except if a region
+ with a nonzero heatsource Q is in contact with the surface S.
+
+ So we have nearly always the following practical formula
+ to evaluate the heat flux across a surface S,
+ in terms of a well-chosen auxiliary scalar function g(S).
+
+ q(S) = ( -lambda[] grad T . grad g(S) )_{support of g(S)}
+
+
+ Particular cases:
+
+ For the heat flux through the boundary of a thermal electrode,
+ one uses g(S) = BF(T_electrode).
+ Note that this heat flux is equal to Q_electrode
+ in the stationary case.
+ This is the case for the heat flux through the boundary of Window2
+
+ Auxiliary functions g(S) are also generated
+ for the surfaces named "Surface_i", i=1,2,3 in the model.
+ Note that the flux computed through Surface_3 is incorrect
+ because this surface is not adjacent to surfaces
+ with homogeneous Neumann boundary conditions.
+*/
+
+Include "brick_common.pro";
+
+QWindow1 =
+ DefineNumber[1e3, Name "Window1/Heat source [W]"];
+
+/* The user is given the choice of setting either the global temperature
+ or the global heat flux in Window2.
+ Check how the variable "Flag_ConstraintWin2" is used at different
+ places to alter the model according to that choice
+ and to manage the visibility of related input and output data.*/
+QWindow2 =
+ DefineNumber[1e3, Name "Window2/Heat source [W]",
+ Visible Flag_ConstraintWin2];
+TWindow2 =
+ DefineNumber[50, Name "Window2/Temperature [degC]",
+ Visible !Flag_ConstraintWin2];
+outQWindow2 =
+ DefineNumber[0, Name "Output/Q window 2 [degC]",
+ Visible !Flag_ConstraintWin2, Highlight "Ivory"];
+outTWindow2 =
+ DefineNumber[0, Name "Output/T window 2 [degC]",
+ Visible Flag_ConstraintWin2, Highlight "Ivory"];
+
+
+ConvectionCoef =
+ DefineNumber[1000, Name"Surface2/hconv",
+ Label "Convection coefficient [W/(m^2K)]"];
+T_Ambiance =
+ DefineNumber[20, Name"Surface2/Ambiance temperature [degC]"];
+T_Dirichlet =
+ DefineNumber[20, Name"Surface1/Imposed temperature [degC]"];
+
+Group {
+ /* Geometrical regions: */
+ Brick = Region[100];
+ LayerWindow1 = Region[115];
+ Window1 = Region[111];
+ Window2 = Region[112];
+ SurfWindow1 = Region[211];
+ SurfWindow2 = Region[212];
+
+ Surface_1 = Region[201];
+ Surface_2 = Region[202];
+ Surface_3 = Region[203];
+ NbSurface = 3;
+
+ /* Abstract regions:
+
+ Vol_Lambda_Ele : volume regions with a thermal conductivity lambda[]
+ Sur_Dirichlet_The : Dirichlet bondary condition surface
+ Sur_Neu_Ele : Neumann bondary condition surface
+ Sur_Convection_The : convective surface q.n = h ( T-Tinf )
+ Vol_Qsource_The : volume heat source regions
+ Tfloating_The : thermal electrodes
+ */
+
+ Vol_Lambda_The = Region[ {Brick, Window1, Window2, LayerWindow1} ];
+ Sur_Dirichlet_The = Region[ Surface_1 ];
+ Sur_Neumann_The = Region[ Surface_3 ];
+ Sur_Convection_The = Region[ Surface_2 ];
+
+ Vol_Qsource_The = Region[ Window1 ];
+ Tfloating_The = Region[ Window2 ];
+
+ If( !Flag_Regularization )
+ Vol_OneSide_The = Region[ Brick ];
+ Sur_Tdisc_The = Region[ SurfWindow1 ];
+ Else
+ Vol_OneSide_The = Region[ {} ];
+ Sur_Tdisc_The = Region[ {} ];
+ EndIf
+}
+
+
+Function {
+ lambda_Brick = 50.; // steel
+ lambda[ Brick ] = lambda_Brick;
+ lambda[ Region[ {Window1, Window2} ] ] = lambda_Brick ; // comment
+
+ lambda_Layer = 1.0; // 50 time smaller than surrounding regions
+ If( Flag_Regularization )
+ lambda[ LayerWindow1 ] = lambda_Layer;
+ EndIf
+ lambda[ Sur_Tdisc_The ] = lambda_Layer;
+ thickness[ Sur_Tdisc_The ] = e_layer;
+
+ h[ Surface_2 ] = ConvectionCoef;
+ Tinf[ Surface_2 ] = T_Ambiance;
+
+ If( Flag_QFromFile )
+ Qsource[ Window1 ] = ScalarField[XYZ[],0,1]{1};
+ Else
+ Qsource[ Window1 ] = 0*QWindow1/SurfaceArea[];
+ EndIf
+}
+
+Constraint {
+ { Name Dirichlet_The ;
+ Case {
+ { Region Sur_Dirichlet_The ; Value T_Dirichlet ; }
+ }
+ }
+ { Name T_Discontinuity ;
+ Case {
+ { Region SurfWindow1 ; Value 20 ; }
+ }
+ }
+ { Name T_electrode;
+ Case {
+ If( !Flag_ConstraintWin2 )
+ { Region Window2 ; Value TWindow2 ; }
+ EndIf
+ }
+ }
+ { Name Q_electrode;
+ Case {
+ If( Flag_ConstraintWin2 )
+ { Region Window2 ; Value QWindow2 ; }
+ EndIf
+ }
+ }
+ For i In {1:NbSurface}
+ { Name FluxLayer~{i} ;
+ Case {
+ { Region Surface~{i} ; Value 1. ; }
+ }
+ }
+ EndFor
+}
+
+Integration {
+ { Name Int ;
+ Case {
+ {
+ Type Gauss ;
+ Case {
+ { GeoElement Triangle ; NumberOfPoints 4 ; }
+ { GeoElement Quadrangle ; NumberOfPoints 4 ; }
+ { GeoElement Line ; NumberOfPoints 4 ; }
+ }
+ }
+ }
+ }
+}
+
+Jacobian {
+ { Name Vol ;
+ Case {
+ { Region All ; Jacobian Vol ; }
+ }
+ }
+ { Name Sur ;
+ Case {
+ { Region All ; Jacobian Sur ; }
+ }
+ }
+}
+
+Group {
+ Dom_Hgrad_T = Region[ {Vol_Lambda_The,
+ Sur_Convection_The,
+ Sur_Tdisc_The} ];
+ DomainWithSurf_TL_The =
+ ElementsOf[ {Vol_OneSide_The, Sur_Tdisc_The},
+ OnOneSideOf Sur_Tdisc_The ];
+}
+
+FunctionSpace {
+ { Name Hgrad_T; Type Form0 ;
+ BasisFunction {
+ { Name sn ; NameOfCoef Tn ; Function BF_Node ;
+ Support Dom_Hgrad_T ;
+ Entity NodesOf[All, Not Tfloating_The] ; }
+ { Name sf ; NameOfCoef Tf ; Function BF_GroupOfNodes ;
+ Support Dom_Hgrad_T ; Entity GroupsOfNodesOf[ Tfloating_The ] ; }
+ { Name sdn ; NameOfCoef Tdn ; Function BF_Node ;
+ Support DomainWithSurf_TL_The ; Entity NodesOf[ Sur_Tdisc_The ] ; }
+ }
+ SubSpace {
+ { Name Tcont ; NameOfBasisFunction { sn, sf } ; }
+ { Name Tdisc ; NameOfBasisFunction { sdn } ; }
+ }
+ GlobalQuantity {
+ { Name T_electrode ; Type AliasOf ; NameOfCoef Tf ; }
+ { Name Q_electrode ; Type AssociatedWith ; NameOfCoef Tf ; }
+ }
+ Constraint {
+ { NameOfCoef Tn ; EntityType NodesOf ;
+ NameOfConstraint Dirichlet_The ; }
+ // { NameOfCoef Tdn ; EntityType NodesOf ; //
+ // NameOfConstraint T_Discontinuity ; } //
+ { NameOfCoef T_electrode ; EntityType GroupsOfNodesOf ;
+ NameOfConstraint T_electrode ; }
+ { NameOfCoef Q_electrode ; EntityType GroupsOfNodesOf ;
+ NameOfConstraint Q_electrode ; }
+ }
+ }
+ For i In {1:NbSurface}
+ { Name FluxLayer~{i} ; Type Form0 ;
+ BasisFunction {
+ { Name gn ; NameOfCoef un ; Function BF_GroupOfNodes;
+ Support Dom_Hgrad_T ; Entity GroupsOfNodesOf[ Surface~{i} ] ; }
+ }
+ Constraint {
+ { NameOfCoef un ; EntityType GroupsOfNodesOf ;
+ NameOfConstraint FluxLayer~{i} ; }
+ }
+ }
+ EndFor
+
+}
+
+Formulation {
+ { Name Thermal_T ; Type FemEquation ;
+ Quantity {
+ { Name T ; Type Local ; NameOfSpace Hgrad_T ; }
+ { Name Tcont ; Type Local ; NameOfSpace Hgrad_T[Tcont] ; }
+ { Name Tdisc ; Type Local ; NameOfSpace Hgrad_T[Tdisc] ; }
+ { Name Tglob ; Type Global ; NameOfSpace Hgrad_T[T_electrode] ; }
+ { Name Qglob ; Type Global ; NameOfSpace Hgrad_T[Q_electrode] ; }
+ For i In {1:NbSurface}
+ { Name un~{i} ; Type Local ; NameOfSpace FluxLayer~{i} ; }
+ EndFor
+
+ }
+ Equation {
+
+ Galerkin { [ lambda[] * Dof{d T} , {d T} ] ;
+ In Vol_Lambda_The; Integration Int ; Jacobian Vol ; }
+
+ Galerkin { [ ( lambda[]/thickness[] ) * Dof{Tdisc} , {Tdisc} ] ;
+ In Sur_Tdisc_The; Integration Int ; Jacobian Sur ; }
+
+ Galerkin { [ -Qsource[] , {T} ] ;
+ In Vol_Qsource_The ; Integration Int ; Jacobian Vol ; }
+
+ Galerkin { [ h[] * Dof{T} , {T} ] ;
+ In Sur_Convection_The ; Integration Int ; Jacobian Sur ; }
+
+ Galerkin { [ -h[] * Tinf[] , {T} ] ;
+ In Sur_Convection_The ; Integration Int ; Jacobian Sur ; }
+
+ GlobalTerm { [-Dof{Qglob} , {Tglob} ] ; In Tfloating_The ; }
+
+ For i In {1:NbSurface}
+ Galerkin { [ 0 * Dof{un~{i}} , {un~{i}} ] ;
+ In Vol_Lambda_The ; Integration Int ; Jacobian Vol ; }
+ EndFor
+ }
+ }
+}
+
+
+Resolution {
+ { Name Thermal_T ;
+ System {
+ { Name Sys_The ; NameOfFormulation Thermal_T ; }
+ }
+ Operation {
+ If( Flag_QFromFile )
+ GmshRead[ "Q.pos", 1];
+ EndIf
+ DeleteFile["output.txt"];
+ Generate Sys_The ; Solve Sys_The ; SaveSolution Sys_The ;
+ }
+ }
+}
+
+PostProcessing {
+ { Name Thermal_T ; NameOfFormulation Thermal_T ;
+ PostQuantity {
+ { Name T ; Value { Term { [ {T} ] ;
+ In Dom_Hgrad_T ; Jacobian Vol ; } } }
+ { Name q ; Value { Term { [ -lambda[] * {d T} ] ;
+ In Dom_Hgrad_T ; Jacobian Vol ; } } }
+ { Name Tcont ; Value { Term { [ {Tcont} ] ;
+ In Dom_Hgrad_T ; Jacobian Vol ; } } }
+ { Name Tdisc ; Value { Term { [ {Tdisc} ] ;
+ In Dom_Hgrad_T ; Jacobian Vol ; } } }
+ { Name Qglob; Value { Term { [ {Qglob} ]; In Tfloating_The; } } }
+ { Name Tglob; Value { Term { [ {Tglob} ]; In Tfloating_The; } } }
+
+ For i In {1:NbSurface}
+ { Name un~{i} ; Value { Local { [ {un~{i}} ] ;
+ In Vol_Lambda_The ; Jacobian Vol ; } } }
+ { Name IFlux~{i} ; Value { Integral { [ -lambda[]*{d T} * {d un~{i}} ];
+ In Vol_Lambda_The ; Jacobian Vol ; Integration Int ; } } }
+ EndFor
+
+ }
+ }
+}
+
+PostOperation Map_T UsingPost Thermal_T {
+ If( !Flag_Regularization )
+ Print[ Tcont, OnElementsOf Vol_Lambda_The, File "Tcont.pos"] ;
+ Print[ Tdisc, OnElementsOf Vol_Lambda_The, File "Tdisc.pos"] ;
+ EndIf
+
+ If(Flag_ConstraintWin2)
+ Print[ Tglob, OnRegion Tfloating_The, File > "output.txt", Color "Ivory",
+ Format Table, SendToServer "Output/T window 2 [degC]" ];
+ Else
+ Print[ Qglob, OnRegion Tfloating_The, File > "output.txt", Color "Ivory",
+ Format Table, SendToServer "Output/Q window 2 [degC]" ];
+ EndIf
+
+ For i In {1:NbSurface}
+ Print[un~{i}, OnElementsOf Vol_Lambda_The,
+ File Sprintf("FluxLayer_%g.pos",i)];
+ Print[ IFlux~{i}[Vol_Lambda_The], OnGlobal,
+ Format TimeTable, File > "Fluxes.dat", Color "Ivory",
+ SendToServer Sprintf("Output/Heat flux surface %g [W]", i)];
+ EndFor
+ Print[ T, OnElementsOf Vol_Lambda_The, File "T.pos" ] ;
+ Echo[ StrCat["l=PostProcessing.NbViews-1;",
+ "View[l].IntervalsType = 3;",
+ "View[l].NbIso = 30;",
+ "View[l].NormalRaise = 0.0005;"],
+ File "tmp.geo", LastTimeStepOnly] ;
+ Print [ q, OnLine {{0.02,0.0001,0}{0.02,0.05,0}} {200},
+ Format SimpleTable, File "Cut.txt" ];
+}
+
+
diff --git a/Thermics/brick_common.pro b/Thermics/brick_common.pro
new file mode 100644
index 0000000000000000000000000000000000000000..09a999c27200d6ec8a9f896d1f1eaac94c6f9440
--- /dev/null
+++ b/Thermics/brick_common.pro
@@ -0,0 +1,37 @@
+
+// Onelab parameters
+
+Flag_Regularization =
+ DefineNumber[0, Name "Options/Regularize field", Choices {0,1} ];
+Flag_QFromFile =
+ DefineNumber[0, Name "Options/Heat source from file", Choices {0,1} ];
+Flag_ConstraintWin2 =
+ DefineNumber[0, Name "Options/Constraint in Window2",
+ Choices {0="Fixed temperature", 1="Fixed heat flux"} ];
+MeshRefinement =
+ DefineNumber[1, Name "Options/0Mesh refinement",
+ Help "Choose 1 for a coarse mesh and 0.1 for a fine mesh."];
+
+// Geometrical dimensions
+
+mm=1e-3; // mm to m conversion factor
+
+dx_Brick=100*mm; dy_Brick= 50*mm;
+e_layer = 1*mm;
+dx_Win1 = 20*mm; dy_Win1 = 20*mm;
+If( Flag_Regularization )
+ dx_Win1 += e_layer;
+ dy_Win1 += e_layer;
+EndIf
+
+dx_Win2 = 20*mm; dy_Win2 = 20*mm;
+xc_Win1 = 25*mm; yc_Win1 = dy_Brick/2;
+xc_Win2 = 75*mm; yc_Win2 = dy_Brick/2;
+
+// Element sizes
+
+s=1;
+lc_Brick = dy_Brick/20 *s;
+lc_Win1 = dx_Win1/10 *s;
+lc_Win2 = dx_Win2/10 *s;
+
diff --git a/Thermics/screenshot1.png b/Thermics/screenshot1.png
new file mode 100644
index 0000000000000000000000000000000000000000..9e33f556f9188aad0bf142a03795b063f5cb7816
Binary files /dev/null and b/Thermics/screenshot1.png differ
diff --git a/Thermics/screenshot1_512.png b/Thermics/screenshot1_512.png
new file mode 100644
index 0000000000000000000000000000000000000000..a3ece9e6659eeb45a19eaa56f67e14f248a3f87c
Binary files /dev/null and b/Thermics/screenshot1_512.png differ
diff --git a/Thermics/screenshot2.png b/Thermics/screenshot2.png
new file mode 100644
index 0000000000000000000000000000000000000000..9b4757d415dcf8e0afcf3df89352d71263f3aaf1
Binary files /dev/null and b/Thermics/screenshot2.png differ
diff --git a/Thermics/screenshot2_512.png b/Thermics/screenshot2_512.png
new file mode 100644
index 0000000000000000000000000000000000000000..2fc975627f77f7ef548395a76b6c8e4561655e99
Binary files /dev/null and b/Thermics/screenshot2_512.png differ