Martin Losch

I use gfortran on MacOS with `darwin_arm64_gfortran`. The problem is related to the division by `pLoc = ABS(R_Shelfice)` in the icefront part, which can be zero. I can add an...

I don't understand either, why this is happening. BTW in #772, this does not happen, and I certainly did not change `stic_thermodynamics.F` or `stic_solve4fluxes.F` in the PR.

I did this both on my Mac and on my Linux box: ```` (base) bkli04l026::review (branch_controls_snap_tc)> git remote add owang01 [email protected]:owang01/MITgcm.git (base) bkli04l026::review (branch_controls_snap_tc)> git fetch owang01 X11 forwarding request...

I find that the problem is related to `ALLOW_shiTransCoeff_3d`. When this is undefined, the model runs also on my Apple MacBook: ``` Y Y Y Y>12

@owang01 thanks, but I am afraid it still does not work for me ... bummer.

I see the problem with `iceConductionDistance`. I print in stic_thermodynamics before line 505 ``` print *, 'ml-stic:', i,j,bi,bj,k,ploc,tloc,sloc, & gammaTloc, gammaSloc, & thermalConductionDistance, thermalConductionTemp ml-stic: -1 -1 1 1 7...

With the latest change, I can compile and run the verification experiment isomip.stic, but the AD-experiment fails with a compile error with addummy_in_stepping.f I used this idea to update #773

@owang01 Now the AD-test passes and I get this: ``` Y Y Y Y 16>15< 9 16 16 13 16 16 16 13 16 12 13 13 13 14 13...

I am mentioning (draft) PR #574 where I tried to implement a quadratic bottom drag according to the law of the wall. I guess there's a lot of overlap with...

@averdy, I didn't know if you still used this code for AD-simulations, that's why I didn't just go ahead and suggest something. In our verification experiments, `pkg/bling` is only used...