Daniel Shapero

I could also use this; it looks like the main concern is what assumptions gmsh makes internally. The API really seems like it wants you to work on only one...

I'd also welcome this change in icepack.

Hm that is odd! Can you try changing `opts` to this: ```python opts = { "dirichlet_ids": boundary_ids, "diagnostic_solver_type": "petsc", "diagnostic_solver_parameters": { "ksp_type": "preonly", "pc_type": "lu", "pc_factor_mat_solver_type": "mumps", }, } ```...

Thanks for sharing! I'm able to run your code and working through this but I'm a little behind on writing a grant proposal so it might be a few more...

Hey David, thanks for sending this in and sorry for the slow reply! Been out of work some with a broken leg (funny story!) but back at it now. This...

Thanks for catching this! I suspect the problem is happening because the temperature crosses the slope break at -10C. That could certainly cause the interpolated or projected fluidity field to...

Hi Meghana, thanks for reporting this! Are you able to attach the complete runnable script? I'll have an easier time getting you a fix if I can try it myself.

I'm pretty sure this is because scipy messed up in the 1.10.0 release. I've just pinned icepack to use 1.9.3. Can you rebase from master? Your changes didn't touch interpolation...

Hi Jack, thanks for taking a look! While it would be very useful to include more general data types, C++ will always beat Fortran on this point via template meta-programming....

Hi Victor, thanks for your interest! SiGMA is admittedly lacking in any documentation. I can make a UML diagram of the various classes with some notes if that will be...