@@ -134,24 +134,6 @@ def _unrolled_cone_curvilinear_grid():
134134 )
135135
136136
137- def simple_UV_dataset (dims = (360 , 2 , 30 , 4 ), maxdepth = 1 , mesh_type = "spherical" ):
138- max_lon = 180.0 if mesh_type == "spherical" else 1e6
139-
140- return xr .Dataset (
141- {"U" : (["time" , "depth" , "YG" , "XG" ], np .zeros (dims )), "V" : (["time" , "depth" , "YG" , "XG" ], np .zeros (dims ))},
142- coords = {
143- "time" : (["time" ], xr .date_range ("2000" , "2001" , dims [0 ]), {"axis" : "T" }),
144- "depth" : (["depth" ], np .linspace (0 , maxdepth , dims [1 ]), {"axis" : "Z" }),
145- "YC" : (["YC" ], np .arange (dims [2 ]) + 0.5 , {"axis" : "Y" }),
146- "YG" : (["YG" ], np .arange (dims [2 ]), {"axis" : "Y" , "c_grid_axis_shift" : - 0.5 }),
147- "XC" : (["XC" ], np .arange (dims [3 ]) + 0.5 , {"axis" : "X" }),
148- "XG" : (["XG" ], np .arange (dims [3 ]), {"axis" : "X" , "c_grid_axis_shift" : - 0.5 }),
149- "lat" : (["YG" ], np .linspace (- 90 , 90 , dims [2 ]), {"axis" : "Y" , "c_grid_axis_shift" : 0.5 }),
150- "lon" : (["XG" ], np .linspace (- max_lon , max_lon , dims [3 ]), {"axis" : "X" , "c_grid_axis_shift" : - 0.5 }),
151- },
152- )
153-
154-
155137datasets = {
156138 "2d_left_rotated" : _rotated_curvilinear_grid (),
157139 "ds_2d_left" : xr .Dataset ( # MITgcm indexing style
0 commit comments