<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Good day, colleagues,</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">I have hit a roadblock. I'm attempting to plot zonal and depth-integrated v4r4 variables as a function of latitude. I'm using the xgcm python package and jupyter notebook. I did find LLCMapper and have attempted to use that information to create something that will work for my needs. My problems are</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">
1) 'faces' are called 'tile' in v4r4 variables</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">2) I'm getting the error that I must have 'i' coordinate (I have 'i_g', perhaps shift to 'i'?)</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">3) I'm not entirely sure LLCMapper is going to allow me to do what I eventually want to do.</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Would someone mind taking a look at my code below and giving feedback/suggestions? Thanks in advance!</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Warm regards,</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><div dir="ltr"><span class="gmail_default"></span>Deb Gardner (she/her/ella)</div><div dir="ltr">NSF Graduate Research Fellow</div><div dir="ltr">Best Waffle Topping, 2019 | Scripps Institution of Oceanography<br></div><div dir="ltr">SIO Graduate Student Council, Junior Co-Chair<span class="gmail_default"> 2020-2021</span></div><div>--</div><div>I respectfully acknowledge that UC San Diego and its surrounding communities occupy the stolen territory of the Kumeyaay/Kumiai Nation and to whom I owe honor and gratitude. / <a href="https://native-land.ca/" target="_blank">Whose land are you on?<br></a></div></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"># calculate pressure gradient in x-dir<br>dP = gcm_grid.derivative(P,'X',boundary='extend')<br><br># zonal and depth integral of pressure gradient in x<br>dP_int = gcm_grid.integrate(dP,['Z','X'])<br># dP coords - k, tile, j, i_g<br><br># zonal integration of windstress<br>taue_int = gcm_grid.cumsum(oceTAUE,'X',boundary='extend')</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"># note: grid.integrate does not work with oceTAUE<br>#taue_int has coordinates tile, j, i_g<br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Pmapper = LLCMapper(P_coords)<br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">tauemapper = LLCMapper(taue_coords)</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">#calling either of the mappers with the integrated values produces the same error (below)</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">---------------------------------------------------------------------------<br>AssertionError Traceback (most recent call last)<br><ipython-input-34-1562da53e232> in <module><br>----> 1 Pmapper(dP_int)<br><br><ipython-input-8-efce9e1da44c> in __call__(self, da, ax, projection, lon_0, **plt_kwargs)<br> 14 self.new_grid = pyresample.geometry.GridDefinition(lons=self.new_grid_lon,lats=self.new_grid_lat)<br> 15 def __call__(self, da, ax=None, projection=cart.crs.Robinson(), lon_0=-60, **plt_kwargs):<br>---> 16 assert set(da.dims) == set(['face', 'j', 'i']), "da must have dimensions ['face', 'j', 'i']"<br> 17 if ax is None:<br> 18 fig, ax = plt.subplots(figsize=(12, 6))<br><br>AssertionError: da must have dimensions ['face', 'j', 'i']<br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"></span></div><div><a href="https://native-land.ca/" target="_blank"></a><br></div></div></div></div></div></div></div></div>