← Back to branch summary

~asc/meshing/bastion

« back to all changes in this revision

Viewing changes to examples/valid/filchner-ronne/shorelines.geo

  • Committer: Adam Candy
  • Date: 2012-05-21 17:55:22 UTC
  • Revision ID: adam.candy@imperial.ac.uk-20120521175522-kmoa2r62pdg3gckf
It is now possible to generate the Southern Ocean simulation with this version, which should (awaiting confirmation) run without the failures we've seen recently.  Note the changes to the boundary IDs - the open and contour boundaries are now set as prescribed by the boundary class object - and are fixed, where they were previously variable depending on what the Gmsh command 'newll' returned.

Show diffs side-by-side

added added

removed removed

Lines of Context:
examples/valid/filchner-ronne/shorelines.geo
1
 
// Arguments: -f filchner-ronne/shorelines.geo -n ../RTopo105b_50S.nc -no -b -85.0:-20.0,-89.0:-75.0 -64.0:-30.0,-89.0:-70.0 -30.0:-20.0,-89.0:-75.0 -lat -65.0
 
1
// Arguments:  -f filchner-ronne/shorelines.geo -n ../RTopo105b_50S.nc -no -b -85.0:-20.0,-89.0:-75.0 -64.0:-30.0,-89.0:-70.0 -30.0:-20.0,-89.0:-75.0 -lat -65.0
2
2
// Source netCDF located at ../RTopo105b_50S.nc
3
3
// Output to filchner-ronne/shorelines.geo
4
4
// Region defined by ((longitude >= -85.0) and (longitude <= -20.0) and (latitude >= -89.0) and (latitude <= -75.0)) or ((longitude >= -64.0) and (longitude <= -30.0) and (latitude >= -89.0) and (latitude <= -70.0)) or ((longitude >= -30.0) and (longitude <= -20.0) and (latitude >= -89.0) and (latitude <= -75.0))
12345
12345
LoopStart0 = IP + 2;
12346
12346
LoopEnd0 = IP + 12324;
12347
12347
BSpline ( IL + 0 ) = { IP + 2 : IP + 12324 };
12348
 
Physical Line( IL + 0 ) = { IL + 0 };
 
12348
Physical Line( 3 ) = { IL + 0 };
12349
12349
 
12350
12350
// Closing path with parallels and merdians, from (-62.15833282, -70.00000000) to  (-20.00000000, -75.52499771)
12351
12351
// Drawing meridian to max latitude index 12324 at -62.158333.2, -70.000000.2 (to match -65.000000.2)
12929
12929
LoopStart1 = IP + 12324;
12930
12930
LoopEnd1 = IP + 12899;
12931
12931
BSpline ( IL + 1 ) = { IP + 12324 : IP + 12899, IP + 2 };
12932
 
Physical Line( IL + 1 ) = { IL + 1 };
 
12932
Physical Line( 4 ) = { IL + 1 };
12933
12933
 
12934
12934
Line Loop( ILL + 0 ) = { IL + 0, IL + 1 };
12935
12935
// **** END end of loop False12323/12323True
14148
14148
LoopStart2 = IP + 12900;
14149
14149
LoopEnd2 = IP + 14103;
14150
14150
BSpline ( IL + 2 ) = { IP + 12900 : IP + 14103, IP + 12900 };
14151
 
Physical Line( IL + 2 ) = { IL + 2 };
 
14151
Physical Line( 3 ) = { IL + 2 };
14152
14152
 
14153
14153
Line Loop( ILL + 1 ) = { IL + 2 };
14154
14154
// Ice-Land mass number 5
16348
16348
LoopStart3 = IP + 14104;
16349
16349
LoopEnd3 = IP + 16295;
16350
16350
BSpline ( IL + 3 ) = { IP + 14104 : IP + 16295, IP + 14104 };
16351
 
Physical Line( IL + 3 ) = { IL + 3 };
 
16351
Physical Line( 3 ) = { IL + 3 };
16352
16352
 
16353
16353
Line Loop( ILL + 2 ) = { IL + 3 };
16354
16354
// Ice-Land mass number 6
16920
16920
LoopStart4 = IP + 16296;
16921
16921
LoopEnd4 = IP + 16859;
16922
16922
BSpline ( IL + 4 ) = { IP + 16296 : IP + 16859, IP + 16296 };
16923
 
Physical Line( IL + 4 ) = { IL + 4 };
 
16923
Physical Line( 3 ) = { IL + 4 };
16924
16924
 
16925
16925
Line Loop( ILL + 3 ) = { IL + 4 };
16926
16926
// Ice-Land mass number 7
17765
17765
LoopStart5 = IP + 16860;
17766
17766
LoopEnd5 = IP + 17693;
17767
17767
BSpline ( IL + 5 ) = { IP + 16860 : IP + 17693, IP + 16860 };
17768
 
Physical Line( IL + 5 ) = { IL + 5 };
 
17768
Physical Line( 3 ) = { IL + 5 };
17769
17769
 
17770
17770
Line Loop( ILL + 4 ) = { IL + 5 };
17771
17771
// Ice-Land mass number 9
18123
18123
LoopStart6 = IP + 17694;
18124
18124
LoopEnd6 = IP + 18043;
18125
18125
BSpline ( IL + 6 ) = { IP + 17694 : IP + 18043, IP + 17694 };
18126
 
Physical Line( IL + 6 ) = { IL + 6 };
 
18126
Physical Line( 3 ) = { IL + 6 };
18127
18127
 
18128
18128
Line Loop( ILL + 5 ) = { IL + 6 };
18129
18129
// Ice-Land mass number 10
18279
18279
LoopStart7 = IP + 18044;
18280
18280
LoopEnd7 = IP + 18191;
18281
18281
BSpline ( IL + 7 ) = { IP + 18044 : IP + 18191, IP + 18044 };
18282
 
Physical Line( IL + 7 ) = { IL + 7 };
 
18282
Physical Line( 3 ) = { IL + 7 };
18283
18283
 
18284
18284
Line Loop( ILL + 6 ) = { IL + 7 };
18285
18285
// Ice-Land mass number 11
18484
18484
LoopStart8 = IP + 18192;
18485
18485
LoopEnd8 = IP + 18259;
18486
18486
BSpline ( IL + 8 ) = { IP + 18192 : IP + 18259, IP + 18192 };
18487
 
Physical Line( IL + 8 ) = { IL + 8 };
 
18487
Physical Line( 3 ) = { IL + 8 };
18488
18488
 
18489
18489
Line Loop( ILL + 7 ) = { IL + 8 };
18490
18490
// Ice-Land mass number 55
18587
18587
LoopStart9 = IP + 18260;
18588
18588
LoopEnd9 = IP + 18351;
18589
18589
BSpline ( IL + 9 ) = { IP + 18260 : IP + 18351, IP + 18260 };
18590
 
Physical Line( IL + 9 ) = { IL + 9 };
 
18590
Physical Line( 3 ) = { IL + 9 };
18591
18591
 
18592
18592
Line Loop( ILL + 8 ) = { IL + 9 };
18593
18593
// Ice-Land mass number 57
19370
19370
 
19371
19371
Delete { Point{1}; }
19372
19372
Delete { Point{2}; }
19373
 
// Open boundaries   (id 10010): 2
19374
 
// Closed boundaries (id 10009): 1 3 4 5 6 7 8 9 10
 
19373
// Open boundaries   (id 4): 2
 
19374
// Closed boundaries (id 3): 1 3 4 5 6 7 8 9 10
19375
19375
 
19376
 
Plane Surface( 10012 ) = { ILL + 0, ILL + 1, ILL + 2, ILL + 3, ILL + 4, ILL + 5, ILL + 6, ILL + 7, ILL + 8 };
19377
 
Physical Surface( 10012 ) = { 10012 };
 
19376
Plane Surface( 10 ) = { ILL + 0, ILL + 1, ILL + 2, ILL + 3, ILL + 4, ILL + 5, ILL + 6, ILL + 7, ILL + 8 };
 
19377
Physical Surface( 10 ) = { 10 };
19378
19378
 
19379
19379
Printf("Assigning characteristic mesh sizes...");
19380
19380
// Field[ IFI + 1] = Attractor;