calculate ip-gradients with respect to old positions if MUDFLAT

src/3d/internal_pressure.F90

@@ -67,7 +67,12 @@
 #else
    use domain, only: dx,dy
 #endif
-   use variables_3d, only: kmin,hn,hun,hvn,idpdx,idpdy,buoy,ssun,ssvn,ssen
+   use variables_3d, only: kmin,hun,hvn,idpdx,idpdy,buoy,ssun,ssvn,ssen
+#ifdef MUDFLAT
+   use variables_3d, only: hn=>ho
+#else
+   use variables_3d, only: hn
+#endif
    IMPLICIT NONE
 !
 ! !PUBLIC DATA MEMBERS:

src/3d/ip_shchepetkin_mcwilliams.F90

@@ -38,7 +38,7 @@
 !
 ! !USES:
    use internal_pressure
-   use variables_3d, only: hn,buoy,sseo
+   use variables_3d, only: buoy,sseo
    use domain, only: H,az,au,av
 !$ use omp_lib
    IMPLICIT NONE