[Gmsh] constant mesh despite background field

Tue Apr 5 17:34:46 CEST 2016

Hi everyone,
my first post here, I hope it is not too silly a question.

I am trying to generate a gis-based 2D mesh for some large-scale 
geophysical fluid dynamics model. My .geo file is attached, and has been 
generated by this QGIS plugin:
  https://github.com/ccorail/qgis-gmsh/wiki

The essential part is this:
NF = newf;
Field[NF] = Structured;
Field[NF].TextFormat = 0;
Field[NF].FileName = "/home/cimatori/Data/Mesh/mesh.dat";
Background Field = NF;

As far as I understand, this should mean that the mesh follows the 
scales given in the mesh.dat file. However, the mesh I get is almost 
entirely fixed at the smallest scale (apart from some small region close 
to complex geometry).
I did check the content of this file, and it contains the numbers I 
expect to find (range 30-300m). The data in mesh.dat seems to be 
completely ignored. Could it be related to the fact that on all the 
boundaries the scale is expected to be the same (30m)? But how to fix it?

Thanks for your help.

-- 

Andrea Cimatoribus
post-doctoral researcher
EPFL ENAC IIE ECOL
people.epfl.ch/andrea.cimatoribus
-------------- next part --------------
IP = newp;
IL = newl;
IS = news;
ILL = newll;
Point(IP+0) = {500044.1749674345, 117961.3523919472, 0};
Point(IP+1) = {556366.0496117522, 137244.5617184398, 0};
Point(IP+2) = {500399, 118091, 0};
Point(IP+3) = {500795, 118539, 0};
Point(IP+4) = {500698, 119270, 0};
Point(IP+5) = {500870, 119672, 0};
Point(IP+6) = {500528, 120288, 0};
Point(IP+7) = {500685, 121027, 0};
Point(IP+8) = {500660, 121394, 0};
Point(IP+9) = {500849, 122076, 0};
Point(IP+10) = {500812, 122227, 0};
Point(IP+11) = {500917, 122408, 0};
Point(IP+12) = {500739, 122642, 0};
Point(IP+13) = {500746, 122985, 0};
Point(IP+14) = {500924, 123252, 0};
Point(IP+15) = {501257, 123864, 0};
Point(IP+16) = {501570, 124289, 0};
Point(IP+17) = {501783, 124362, 0};
Point(IP+18) = {502071, 124947, 0};
Point(IP+19) = {502269, 125377, 0};
Point(IP+20) = {502234, 125453, 0};
Point(IP+21) = {502314, 125586, 0};
Point(IP+22) = {502009, 126354, 0};
Point(IP+23) = {502249, 127165, 0};
Point(IP+24) = {502237, 127538, 0};
Point(IP+25) = {502304, 127923, 0};
Point(IP+26) = {502657, 128183, 0};
Point(IP+27) = {502880, 128429, 0};
Point(IP+28) = {503156, 128817, 0};
Point(IP+29) = {503218, 129107, 0};
Point(IP+30) = {503684, 129784, 0};
Point(IP+31) = {504135, 130144, 0};
Point(IP+32) = {504403, 130532, 0};
Point(IP+33) = {504841, 131696, 0};
Point(IP+34) = {504839, 132215, 0};
Point(IP+35) = {505197, 132721, 0};
Point(IP+36) = {505225, 133668, 0};
Point(IP+37) = {505561, 134068, 0};
Point(IP+38) = {505711, 134472, 0};
Point(IP+39) = {505871, 134554, 0};
Point(IP+40) = {505953, 134898, 0};
Point(IP+41) = {506183, 135061, 0};
Point(IP+42) = {506322, 135415, 0};
Point(IP+43) = {506649, 135625, 0};
Point(IP+44) = {507253, 136285, 0};
Point(IP+45) = {507383, 136789, 0};
Point(IP+46) = {507914, 137211, 0};
Point(IP+47) = {508139, 137574, 0};
Point(IP+48) = {508234, 137610, 0};
Point(IP+49) = {508351, 137966, 0};
Point(IP+50) = {509214, 138638, 0};
Point(IP+51) = {509686, 138600, 0};
Point(IP+52) = {509895, 138410, 0};
Point(IP+53) = {510399, 138314, 0};
Point(IP+54) = {510724, 138931, 0};
Point(IP+55) = {511028, 139044, 0};
Point(IP+56) = {511113, 139403, 0};
Point(IP+57) = {511028, 140038, 0};
Point(IP+58) = {511672, 140931, 0};
Point(IP+59) = {511694, 141618, 0};
Point(IP+60) = {511819, 141992, 0};
Point(IP+61) = {512168, 142047, 0};
Point(IP+62) = {512318, 142320, 0};
Point(IP+63) = {512972, 142674, 0};
Point(IP+64) = {513845, 143750, 0};
Point(IP+65) = {513950, 144161, 0};
Point(IP+66) = {514202, 144249, 0};
Point(IP+67) = {514681, 144723, 0};
Point(IP+68) = {514996, 144859, 0};
Point(IP+69) = {515315, 145262, 0};
Point(IP+70) = {515465, 145589, 0};
Point(IP+71) = {515883, 145980, 0};
Point(IP+72) = {516171, 146146, 0};
Point(IP+73) = {516341, 146078, 0};
Point(IP+74) = {516818, 146222, 0};
Point(IP+75) = {517371, 146257, 0};
Point(IP+76) = {518269, 146567, 0};
Point(IP+77) = {518964, 146459, 0};
Point(IP+78) = {519367, 146121, 0};
Point(IP+79) = {519555, 145940, 0};
Point(IP+80) = {519862, 145831, 0};
Point(IP+81) = {520387, 145803, 0};
Point(IP+82) = {520825, 146005, 0};
Point(IP+83) = {520930, 146247, 0};
Point(IP+84) = {521828, 146844, 0};
Point(IP+85) = {523154, 147104, 0};
Point(IP+86) = {523934, 147396, 0};
Point(IP+87) = {524039, 147469, 0};
Point(IP+88) = {524064, 147623, 0};
Point(IP+89) = {524454, 147963, 0};
Point(IP+90) = {524952, 147908, 0};
Point(IP+91) = {524919, 148268, 0};
Point(IP+92) = {524994, 148674, 0};
Point(IP+93) = {525570, 149284, 0};
Point(IP+94) = {525792, 149339, 0};
Point(IP+95) = {526087, 149243, 0};
Point(IP+96) = {526533, 149392, 0};
Point(IP+97) = {526610, 150007, 0};
Point(IP+98) = {526890, 150428, 0};
Point(IP+99) = {527038, 150428, 0};
Point(IP+100) = {527600, 150750, 0};
Point(IP+101) = {528111, 151446, 0};
Point(IP+102) = {528063, 151645, 0};
Point(IP+103) = {528563, 152045, 0};
Point(IP+104) = {528693, 152071, 0};
Point(IP+105) = {528959, 151940, 0};
Point(IP+106) = {529206, 152068, 0};
Point(IP+107) = {529709, 151904, 0};
Point(IP+108) = {530722, 151123, 0};
Point(IP+109) = {530954, 151050, 0};
Point(IP+110) = {531300, 151133, 0};
Point(IP+111) = {532125, 151058, 0};
Point(IP+112) = {532608, 151165, 0};
Point(IP+113) = {532970, 151614, 0};
Point(IP+114) = {533563, 151745, 0};
Point(IP+115) = {534038, 152227, 0};
Point(IP+116) = {534348, 152295, 0};
Point(IP+117) = {534686, 152162, 0};
Point(IP+118) = {534898, 152313, 0};
Point(IP+119) = {535298, 152114, 0};
Point(IP+120) = {535784, 151556, 0};
Point(IP+121) = {536001, 151483, 0};
Point(IP+122) = {536964, 151122, 0};
Point(IP+123) = {537264, 150998, 0};
Point(IP+124) = {537602, 150743, 0};
Point(IP+125) = {537719, 150705, 0};
Point(IP+126) = {537859, 150771, 0};
Point(IP+127) = {538122, 150995, 0};
Point(IP+128) = {538974, 150695, 0};
Point(IP+129) = {539205, 150824, 0};
Point(IP+130) = {539900, 150758, 0};
Point(IP+131) = {540240, 150773, 0};
Point(IP+132) = {540430, 150748, 0};
Point(IP+133) = {540645, 150672, 0};
Point(IP+134) = {540822, 150602, 0};
Point(IP+135) = {541185, 150561, 0};
Point(IP+136) = {541528, 150622, 0};
Point(IP+137) = {541748, 150576, 0};
Point(IP+138) = {542095, 150327, 0};
Point(IP+139) = {542380, 150314, 0};
Point(IP+140) = {542686, 150094, 0};
Point(IP+141) = {542848, 150107, 0};
Point(IP+142) = {543343, 149913, 0};
Point(IP+143) = {544976, 148585, 0};
Point(IP+144) = {545646, 148615, 0};
Point(IP+145) = {545944, 148878, 0};
Point(IP+146) = {546204, 148840, 0};
Point(IP+147) = {546312, 148812, 0};
Point(IP+148) = {546664, 148890, 0};
Point(IP+149) = {547242, 148522, 0};
Point(IP+150) = {547967, 147812, 0};
Point(IP+151) = {549262, 147138, 0};
Point(IP+152) = {549647, 147070, 0};
Point(IP+153) = {550132, 147226, 0};
Point(IP+154) = {551013, 146871, 0};
Point(IP+155) = {551646, 146779, 0};
Point(IP+156) = {551758, 146749, 0};
Point(IP+157) = {552466, 146678, 0};
Point(IP+158) = {552706, 146454, 0};
Point(IP+159) = {553324, 146480, 0};
Point(IP+160) = {553650, 145997, 0};
Point(IP+161) = {553657, 145664, 0};
Point(IP+162) = {554090, 145436, 0};
Point(IP+163) = {554185, 145471, 0};
Point(IP+164) = {554928, 145257, 0};
Point(IP+165) = {555143, 144720, 0};
Point(IP+166) = {555308, 144658, 0};
Point(IP+167) = {555440, 144475, 0};
Point(IP+168) = {555568, 144126, 0};
Point(IP+169) = {555921, 144291, 0};
Point(IP+170) = {556273, 144194, 0};
Point(IP+171) = {556363, 144249, 0};
Point(IP+172) = {556596, 144209, 0};
Point(IP+173) = {556858, 144043, 0};
Point(IP+174) = {556981, 143897, 0};
Point(IP+175) = {556993, 143672, 0};
Point(IP+176) = {557316, 143383, 0};
Point(IP+177) = {557856, 143210, 0};
Point(IP+178) = {558642, 143388, 0};
Point(IP+179) = {558814, 143310, 0};
Point(IP+180) = {559287, 142831, 0};
Point(IP+181) = {559347, 142608, 0};
Point(IP+182) = {559284, 142292, 0};
Point(IP+183) = {559442, 142124, 0};
Point(IP+184) = {559572, 142061, 0};
Point(IP+185) = {559825, 142186, 0};
Point(IP+186) = {560082, 142103, 0};
Point(IP+187) = {560410, 141549, 0};
Point(IP+188) = {560390, 141185, 0};
Point(IP+189) = {560685, 140879, 0};
Point(IP+190) = {560665, 140420, 0};
Point(IP+191) = {560863, 140222, 0};
Point(IP+192) = {560938, 139725, 0};
Point(IP+193) = {560765, 139341, 0};
Point(IP+194) = {560497, 139100, 0};
Point(IP+195) = {560530, 138902, 0};
Point(IP+196) = {560372, 138541, 0};
Point(IP+197) = {560092, 138413, 0};
Point(IP+198) = {560055, 138355, 0};
Point(IP+199) = {559775, 138307, 0};
Point(IP+200) = {559617, 138352, 0};
Point(IP+201) = {559164, 138468, 0};
Point(IP+202) = {558592, 138470, 0};
Point(IP+203) = {558231, 138250, 0};
Point(IP+204) = {557714, 138367, 0};
Point(IP+205) = {557619, 138468, 0};
Point(IP+206) = {557376, 138373, 0};
Point(IP+207) = {556886, 138388, 0};
Point(IP+208) = {556626, 138247, 0};
Point(IP+209) = {556491, 138307, 0};
Point(IP+210) = {555840.2028041825, 138385.7866801331, 0};
Point(IP+211) = {555558.0201996198, 138289.8398288973, 0};
Line(IL+0) = {IP+0, IP+2, IP+3, IP+4, IP+5, IP+6, IP+7, IP+8, IP+9, IP+10, IP+11, IP+12, IP+13, IP+14, IP+15, IP+16, IP+17, IP+18, IP+19, IP+20, IP+21, IP+22, IP+23, IP+24, IP+25, IP+26, IP+27, IP+28, IP+29, IP+30, IP+31, IP+32, IP+33, IP+34, IP+35, IP+36, IP+37, IP+38, IP+39, IP+40, IP+41, IP+42, IP+43, IP+44, IP+45, IP+46, IP+47, IP+48, IP+49, IP+50, IP+51, IP+52, IP+53, IP+54, IP+55, IP+56, IP+57, IP+58, IP+59, IP+60, IP+61, IP+62, IP+63, IP+64, IP+65, IP+66, IP+67, IP+68, IP+69, IP+70, IP+71, IP+72, IP+73, IP+74, IP+75, IP+76, IP+77, IP+78, IP+79, IP+80, IP+81, IP+82, IP+83, IP+84, IP+85, IP+86, IP+87, IP+88, IP+89, IP+90, IP+91, IP+92, IP+93, IP+94, IP+95, IP+96, IP+97, IP+98, IP+99, IP+100, IP+101, IP+102, IP+103, IP+104, IP+105, IP+106, IP+107, IP+108, IP+109, IP+110, IP+111, IP+112, IP+113, IP+114, IP+115, IP+116, IP+117, IP+118, IP+119, IP+120, IP+121, IP+122, IP+123, IP+124, IP+125, IP+126, IP+127, IP+128, IP+129, IP+130, IP+131, IP+132, IP+133, IP+134, IP+135, IP+136, IP+137, IP+138, IP+139, IP+140, IP+141, IP+142, IP+143, IP+144, IP+145, IP+146, IP+147, IP+148, IP+149, IP+150, IP+151, IP+152, IP+153, IP+154, IP+155, IP+156, IP+157, IP+158, IP+159, IP+160, IP+161, IP+162, IP+163, IP+164, IP+165, IP+166, IP+167, IP+168, IP+169, IP+170, IP+171, IP+172, IP+173, IP+174, IP+175, IP+176, IP+177, IP+178, IP+179, IP+180, IP+181, IP+182, IP+183, IP+184, IP+185, IP+186, IP+187, IP+188, IP+189, IP+190, IP+191, IP+192, IP+193, IP+194, IP+195, IP+196, IP+197, IP+198, IP+199, IP+200, IP+201, IP+202, IP+203, IP+204, IP+205, IP+206, IP+207, IP+208, IP+209, IP+210, IP+211, IP+1};
Point(IP+212) = {556297.8924934099, 137187.4198919509, 0};
Point(IP+213) = {500076.8093026768, 117840.9429481225, 0};
Point(IP+214) = {555403.740256654, 138188.3997148289, 0};
Point(IP+215) = {555265.0017526142, 138122.8059588878, 0};
Point(IP+216) = {554395, 137607, 0};
Point(IP+217) = {553204, 137321, 0};
Point(IP+218) = {552342, 137635, 0};
Point(IP+219) = {552099, 137841, 0};
Point(IP+220) = {551686, 137901, 0};
Point(IP+221) = {551411, 138094, 0};
Point(IP+222) = {551281, 138287, 0};
Point(IP+223) = {550168, 138149, 0};
Point(IP+224) = {548918, 138694, 0};
Point(IP+225) = {548753, 138679, 0};
Point(IP+226) = {547805, 139014, 0};
Point(IP+227) = {547587, 139260, 0};
Point(IP+228) = {546012, 139521, 0};
Point(IP+229) = {544994, 139882, 0};
Point(IP+230) = {544123, 139841, 0};
Point(IP+231) = {542528, 139918, 0};
Point(IP+232) = {542047, 139788, 0};
Point(IP+233) = {541484, 139836, 0};
Point(IP+234) = {540756, 139562, 0};
Point(IP+235) = {539637, 139613, 0};
Point(IP+236) = {539234, 139713, 0};
Point(IP+237) = {538523, 139647, 0};
Point(IP+238) = {538318, 139726, 0};
Point(IP+239) = {537694, 139554, 0};
Point(IP+240) = {536615, 139424, 0};
Point(IP+241) = {536423, 139465, 0};
Point(IP+242) = {535567, 139296, 0};
Point(IP+243) = {535053, 139265, 0};
Point(IP+244) = {534588, 139370, 0};
Point(IP+245) = {531282, 138519, 0};
Point(IP+246) = {530525, 138860, 0};
Point(IP+247) = {529429, 139604, 0};
Point(IP+248) = {529212, 139667, 0};
Point(IP+249) = {528739, 139687, 0};
Point(IP+250) = {528291, 139511, 0};
Point(IP+251) = {528036, 139294, 0};
Point(IP+252) = {527191, 138981, 0};
Point(IP+253) = {526538, 138368, 0};
Point(IP+254) = {526463, 138108, 0};
Point(IP+255) = {526473, 137275, 0};
Point(IP+256) = {526303, 136801, 0};
Point(IP+257) = {526068, 136513, 0};
Point(IP+258) = {525743, 136384, 0};
Point(IP+259) = {525368, 136137, 0};
Point(IP+260) = {524230, 135670, 0};
Point(IP+261) = {523265, 135113, 0};
Point(IP+262) = {522515, 134805, 0};
Point(IP+263) = {521107, 134777, 0};
Point(IP+264) = {520447, 133854, 0};
Point(IP+265) = {520257, 133710, 0};
Point(IP+266) = {519921, 133619, 0};
Point(IP+267) = {519466, 132769, 0};
Point(IP+268) = {519266, 132602, 0};
Point(IP+269) = {518574, 132746, 0};
Point(IP+270) = {518093, 133031, 0};
Point(IP+271) = {517431, 133238, 0};
Point(IP+272) = {517168, 133468, 0};
Point(IP+273) = {517008, 133490, 0};
Point(IP+274) = {516736, 134008, 0};
Point(IP+275) = {516758, 134237, 0};
Point(IP+276) = {516368, 134775, 0};
Point(IP+277) = {516383, 134946, 0};
Point(IP+278) = {516243, 135153, 0};
Point(IP+279) = {516293, 135691, 0};
Point(IP+280) = {515733, 136026, 0};
Point(IP+281) = {515420, 135961, 0};
Point(IP+282) = {515015, 136049, 0};
Point(IP+283) = {514708, 135937, 0};
Point(IP+284) = {514366, 136167, 0};
Point(IP+285) = {514187, 136167, 0};
Point(IP+286) = {513966.3413133933, 135984.0811300198, 0};
Point(IP+287) = {513563, 135876, 0};
Point(IP+288) = {512989, 135610, 0};
Point(IP+289) = {512722, 135633, 0};
Point(IP+290) = {512580, 135528, 0};
Point(IP+291) = {512552, 135407, 0};
Point(IP+292) = {511814, 134966, 0};
Point(IP+293) = {511255, 134477, 0};
Point(IP+294) = {510863, 133985, 0};
Point(IP+295) = {510577, 133811, 0};
Point(IP+296) = {510577, 133292, 0};
Point(IP+297) = {510206, 132570, 0};
Point(IP+298) = {509630, 131731, 0};
Point(IP+299) = {509437, 131565, 0};
Point(IP+300) = {509375, 131351, 0};
Point(IP+301) = {508966, 131001, 0};
Point(IP+302) = {509089, 130530, 0};
Point(IP+303) = {508713, 129593, 0};
Point(IP+304) = {508318, 129004, 0};
Point(IP+305) = {507734, 128769, 0};
Point(IP+306) = {507864, 128160, 0};
Point(IP+307) = {507273, 127075, 0};
Point(IP+308) = {506707, 126707, 0};
Point(IP+309) = {506206, 125931, 0};
Point(IP+310) = {506001, 125465, 0};
Point(IP+311) = {505280, 124523, 0};
Point(IP+312) = {504551, 124428, 0};
Point(IP+313) = {504285, 124148, 0};
Point(IP+314) = {504052, 123642, 0};
Point(IP+315) = {504072, 123423, 0};
Point(IP+316) = {503952, 123126, 0};
Point(IP+317) = {503987, 122771, 0};
Point(IP+318) = {503835, 122348, 0};
Point(IP+319) = {503867, 122101, 0};
Point(IP+320) = {503999, 122003, 0};
Point(IP+321) = {503874, 121329, 0};
Point(IP+322) = {503707, 120907, 0};
Point(IP+323) = {503092, 120087, 0};
Point(IP+324) = {502501, 118928, 0};
Point(IP+325) = {502200, 118835, 0};
Point(IP+326) = {501622, 118259, 0};
Point(IP+327) = {501525, 118239, 0};
Point(IP+328) = {501136, 117985, 0};
Point(IP+329) = {500852, 117789, 0};
Point(IP+330) = {500591, 117882, 0};
Line(IL+1) = {IP+212, IP+214, IP+215, IP+216, IP+217, IP+218, IP+219, IP+220, IP+221, IP+222, IP+223, IP+224, IP+225, IP+226, IP+227, IP+228, IP+229, IP+230, IP+231, IP+232, IP+233, IP+234, IP+235, IP+236, IP+237, IP+238, IP+239, IP+240, IP+241, IP+242, IP+243, IP+244, IP+245, IP+246, IP+247, IP+248, IP+249, IP+250, IP+251, IP+252, IP+253, IP+254, IP+255, IP+256, IP+257, IP+258, IP+259, IP+260, IP+261, IP+262, IP+263, IP+264, IP+265, IP+266, IP+267, IP+268, IP+269, IP+270, IP+271, IP+272, IP+273, IP+274, IP+275, IP+276, IP+277, IP+278, IP+279, IP+280, IP+281, IP+282, IP+283, IP+284, IP+285, IP+286, IP+287, IP+288, IP+289, IP+290, IP+291, IP+292, IP+293, IP+294, IP+295, IP+296, IP+297, IP+298, IP+299, IP+300, IP+301, IP+302, IP+303, IP+304, IP+305, IP+306, IP+307, IP+308, IP+309, IP+310, IP+311, IP+312, IP+313, IP+314, IP+315, IP+316, IP+317, IP+318, IP+319, IP+320, IP+321, IP+322, IP+323, IP+324, IP+325, IP+326, IP+327, IP+328, IP+329, IP+330, IP+213};
Line(IL+2) = {IP+212, IP+1};
Line(IL+3) = {IP+0, IP+213};
Line Loop(ILL+0) = {-IL-3, IL+0, -IL-2, IL+1};
NF = newf;
Field[NF] = Structured;
Field[NF].TextFormat = 0;
Field[NF].FileName = "/home/cimatori/Data/Mesh/mesh.dat";
Background Field = NF;
Plane Surface(IS) = {ILL:ILL+0};
Physical Surface("Domain") = {IS};
Physical Line("coast_n") = {IL + 0};
Physical Line("geneva_out") = {IL + 3};
Physical Line("rhone_in") = {IL + 2};
Physical Line("coast_s") = {IL + 1};