@@ -15,39 +15,39 @@ println("Loading Chmy-GMG tools")
1515
1616Creates a GMG `CartGrid` data structure from a `Chmy` grid object
1717"""
18- function create_CartGrid (grid:: StructuredGrid ; ylevel= 0.0 )
18+ function create_CartGrid (grid:: StructuredGrid ; ylevel = 0.0 )
1919
20- coord1D = Vector .(coords (grid, Vertex ()))
20+ coord1D = Vector .(coords (grid, Vertex ()))
2121 coord1D_cen = Vector .(coords (grid, Center ()))
22- N = length .(coord1D)
23- L = extent (grid, Vertex ())
24- X₁ = origin (grid, Vertex ())
25- Δ = spacing (grid)
26- ConstantΔ = false ;
22+ N = length .(coord1D)
23+ L = extent (grid, Vertex ())
24+ X₁ = origin (grid, Vertex ())
25+ Δ = spacing (grid)
26+ ConstantΔ = false
2727 if isa (grid, UniformGrid)
2828 ConstantΔ = true
2929 end
30- if ndims (grid)== 2
30+ if ndims (grid) == 2
3131 # we need a special treatment of this, as all GMG routines work with 3D coordinates
32- X₁ = (X₁[1 ], ylevel, X₁[2 ])
33- L = (L[1 ], 0.0 , L[2 ])
34- Δ = (Δ[1 ], 0.0 , Δ[2 ])
35- N = (N[1 ],1 , N[2 ])
32+ X₁ = (X₁[1 ], ylevel, X₁[2 ])
33+ L = (L[1 ], 0.0 , L[2 ])
34+ Δ = (Δ[1 ], 0.0 , Δ[2 ])
35+ N = (N[1 ], 1 , N[2 ])
3636 coord1D = (coord1D[1 ], [0.0 ], coord1D[2 ])
3737 coord1D_cen = (coord1D_cen[1 ], [0.0 ], coord1D_cen[2 ])
3838 end
39- Xₙ = X₁ .+ L
40-
39+ Xₙ = X₁ .+ L
4140
42- return CartGrid (ConstantΔ,N,Δ,L,X₁,Xₙ,coord1D, coord1D_cen)
41+
42+ return CartGrid (ConstantΔ, N, Δ, L, X₁, Xₙ, coord1D, coord1D_cen)
4343end
4444
4545# all functions to be ported
4646function_names = (:add_box! , :add_sphere! , :add_ellipsoid! , :add_cylinder! , :add_layer! , :add_polygon! , :add_slab! , :add_stripes! , :add_volcano! )
4747
4848for fn in function_names
4949
50- @eval begin
50+ @eval begin
5151 """
5252 $($ fn) ( Phase::Field,
5353 Temp::Field,
@@ -56,22 +56,24 @@ for fn in function_names
5656
5757 Sets `$($ fn) ` function for `Chmy` fields and grids.
5858 """
59- function $fn ( Phase:: Field ,
60- Temp:: Field ,
61- Grid:: StructuredGrid ; # required input
62- kwargs... )
59+ function $fn (
60+ Phase:: Field ,
61+ Temp:: Field ,
62+ Grid:: StructuredGrid ; # required input
63+ kwargs...
64+ )
6365
6466 CartGrid = create_CartGrid (Grid)
6567
6668 cell = false
67- if all (location (Phase).== Center ())
69+ if all (location (Phase) .== Center ())
6870 cell = true
6971 end
70-
71- return ($ fn)(Phase, Temp, CartGrid; cell= cell, kwargs... )
72+
73+ return ($ fn)(Phase, Temp, CartGrid; cell = cell, kwargs... )
7274 end
7375 end
74-
76+
7577end
7678
7779
@@ -80,28 +82,30 @@ function_names = (:above_surface, :below_surface)
8082
8183for fn in function_names
8284
83- @eval begin
85+ @eval begin
8486 """
8587 $($ fn) ( Grid::StructuredGrid, field::Field, DataSurface_Cart::CartData; kwargs...)
8688
8789 Sets `$($ fn) ` function for `Chmy` grids and the field `field` which can be either on vertices or centers
8890 """
89- function $fn ( Grid:: StructuredGrid ,
90- field:: Field ,
91- DataSurface_Cart:: CartData ;
92- kwargs... )
91+ function $fn (
92+ Grid:: StructuredGrid ,
93+ field:: Field ,
94+ DataSurface_Cart:: CartData ;
95+ kwargs...
96+ )
9397
9498 CartGrid = create_CartGrid (Grid)
95-
99+
96100 cell = false
97- if all (location (field).== Center ())
101+ if all (location (field) .== Center ())
98102 cell = true
99103 end
100104
101- return ($ fn)(CartGrid, DataSurface_Cart; cell= cell, kwargs... )
105+ return ($ fn)(CartGrid, DataSurface_Cart; cell = cell, kwargs... )
102106 end
103107 end
104-
108+
105109end
106110
107111
0 commit comments