rotation of modelgrid

[hjia1005]

I have been struggling with a model grid rotation of an mf6 disv modelgrid. The model grid need to be rotated by -45 degree(clockwise) to be along with the principal structure direction. I have to use gridgen to set the active domain and the refinement. Here are the procedures I followed:

1. I rotated all the shapefiles that will be used to define the model extent and refinement by 45 degree. Sorry I really don't know how to get the markdown work here.
   

2. create the gridgen object, using the above xll and yll as the xoff and yoff, do the add_active_domain and add_refinment_features, all shapefiles are from the rotated ones. The gridgen object was created as expected.

3. create a mf6-disv model using the gridprops from the gridgen object. And apply the rotation of -45 degree.
   
   Here the proplem starts, I used the same xll and yll, but the rotated modelgrid run away to another continent. If I want to get it to the right position, apparently something wrong with the xoff and yoff. How should I apply xoff and yoff, why I can't use the original xoff and yoff from the basemodel? Even in each and every step I applied the right epsg it's still the same.

Hope someone could help.

[hjia1005]

It's finally solved by removing the real-world coordinates from the base model and the gridgen object, and making the model origin as (0,0).

To ensure gridgen can still work with real-world shapefiles, they have to be rotated AND shifted beforehand.

Then the code in the second screenshot in my previous post works perfectly. Hope it could help someone with the same needs.