probably a fairly simple question but I lost my head a little bit and I need a fresh look. Using Python and Xarray to load and manipulate my data. I have gridded data of atmospheric field. I've calculated its gradient by using "data.differentiate('longitude')" and "data.differentiate('latitude')". Then I rotated it by using "tx = -ddata_dlat" and ty = ddata_dlon. Now I have component of vectors I guess, that are tangent to and point along the isolines of my data at every grid cell. The next step is to calculate a 10 km step in the direction of the vectors at each point, so that eventually I will have new set of pairs of longitude and latitude that are moved by 10 km in the direction of my isolines. What would be the correct way to do that? I assume I need to calculate the unit vectors for the x and y direction, by I just can't wrap my head around what code should be written to accomplish this 10 km step, broken into some movement in the x and y directions. AI chatbots didn't helped so I seek some human insight here. Thank you those who read this far. Continue reading...
