added +debug variant to ISSM build as pre-release

README.md

@@ -54,3 +54,4 @@ module avail access-issm
 For users of ACCESS-ISSM model configurations released by ACCESS-NRI, the exact location of the model executables is not required. Model configurations will be updated with new model components when necessary.
 
 For information on contributing your own fixes to the ACCESS-ISSM `spack.yaml`, see the [Hive Docs article](https://docs.access-hive.org.au/models/run-a-model/create-a-prerelease/) on creating a prerelease.
+

config/versions.json

@@ -1,5 +1,5 @@
 {
     "$schema": "https://raw.githubusercontent.com/ACCESS-NRI/schema/main/au.org.access-nri/model/deployment/config/versions/3-0-0.json",
     "spack": "0.22",
-    "spack-packages": "2025.04.001"
+    "spack-packages": "justin/debug"
 }

spack.yaml

@@ -1,11 +1,14 @@
 spack:
   specs:
-    - access-issm@git.2025.03.0
+    - access-issm@git.2025.06.0
+
   packages:
     issm:
       require:
-        - '@git.2025.04.11'
+        - '@git.2025.06.19'
         - +wrappers
+        - +debug
+
     # Mark Python 3.9.2 as external so Spack reuses the system module
     python:
       externals:
@@ -27,3 +30,14 @@ spack:
   # Make the solver unify dependencies
   concretizer:
     unify: true
+
+  # Enable generating module files (TCL by default).
+  modules:
+    default:
+      tcl:
+        include:
+          - access-issm
+          - issm
+          - access-triangle
+        projections:
+          access-issm: '{name}/2025.06.0-{hash:7}'

test/Mismip.py

@@ -0,0 +1,76 @@
+import numpy as np
+from SetIceShelfBC import SetIceShelfBC
+from mismipbasalforcings import mismipbasalforcings
+
+# Creating thickness
+print('      creating thickness')
+bx = -150 - 728.8*(md.mesh.x/300000)**2 + 343.91*(md.mesh.x/300000)**4 - 50.57*(md.mesh.x/300000)**6
+by = (500./(1 + np.exp(-2./4000.*(md.mesh.y - 80000./2 - 24000))) +
+      500./(1 + np.exp( 2./4000.*(md.mesh.y - 80000./2 + 24000))))
+by0 = (500./(1 + np.exp(-2./4000.*(0     - 80000./2 - 24000))) +
+       500./(1 + np.exp( 2./4000.*(0     - 80000./2 + 24000))))
+md.geometry.bed     = np.maximum(bx + by, -720)
+md.geometry.surface = np.maximum(bx + by0 + 100, 10)
+md.geometry.base    = np.maximum(md.geometry.bed, -90)
+md.geometry.thickness = md.geometry.surface - md.geometry.base
+
+# Creating drag
+print('      creating drag')
+md.friction.coefficient = np.sqrt(3.160e6) * np.ones(md.mesh.numberofvertices)
+md.friction.p           =            3 * np.ones(md.mesh.numberofelements)
+md.friction.q           =            0 * np.ones(md.mesh.numberofelements)
+
+# Creating flow law parameter
+print('      creating flow law parameter')
+md.materials.rheology_B   = 1/((6.338e-25)**(1/3)) * np.ones(md.mesh.numberofvertices)
+md.materials.rheology_n   =            3 * np.ones(md.mesh.numberofelements)
+md.materials.rheology_law = 'None'
+
+# Boundary conditions for diagnostic model
+print('      boundary conditions for diagnostic model')
+md = SetIceShelfBC(md, './Front.exp')
+md.mask.ice_levelset[:]   = -1
+md.mask.ocean_levelset[:] = -1
+pos = np.where((md.mesh.x < 640000.1) & (md.mesh.x > 639999.9))[0]
+md.mask.ice_levelset[pos] = 0
+md.stressbalance.spcvx[:]  = np.nan
+md.stressbalance.spcvy[:]  = np.nan
+pos = np.where(((md.mesh.y <  80000.1) & (md.mesh.y >  79999.9)) |
+               ((md.mesh.y <      0.1) & (md.mesh.y >     -0.1)))[0]
+md.stressbalance.spcvy[pos] = 0
+pos2 = np.where((md.mesh.x <      0.1) & (md.mesh.x >     -0.1))[0]
+md.stressbalance.spcvx[pos2] = 0
+md.stressbalance.spcvy[pos2] = 0
+
+# Forcing conditions
+print('      forcing conditions')
+# --- basal melt and buttressing conditions -----------------------
+md.basalforcings = mismipbasalforcings()      #  ← remove md
+
+
+md.basalforcings.meltrate_factor        =   0
+md.basalforcings.threshold_thickness    =  75
+md.basalforcings.upperdepth_melt        = -100
+md.smb.mass_balance                    = 0.3 * np.ones(md.mesh.numberofvertices)
+md.basalforcings.geothermalflux        = 0.5 * np.ones(md.mesh.numberofvertices)
+md.basalforcings.groundedice_melting_rate = 0. * np.ones(md.mesh.numberofvertices)
+
+# Other model parameters
+md.thermal.spctemperature        = np.nan * np.ones(md.mesh.numberofvertices)
+md.groundingline.migration       = 'SubelementMigration'
+md.materials.rho_ice             = 918
+md.materials.rho_water           = 1028
+md.constants.g                   = 9.8
+md.constants.yts                 = 31556926
+md.transient.isthermal           = 0
+md.transient.isgroundingline     = 1
+md.stressbalance.isnewton        = 0
+
+# Initialization
+md.initialization.vx          = np.ones(md.mesh.numberofvertices)
+md.initialization.vy          = np.ones(md.mesh.numberofvertices)
+md.initialization.vz          = np.ones(md.mesh.numberofvertices)
+md.initialization.vel         = np.sqrt(2) * np.ones(md.mesh.numberofvertices)
+md.initialization.pressure    = md.constants.g * md.materials.rho_ice * md.geometry.thickness
+md.initialization.temperature = 273 * np.ones(md.mesh.numberofvertices)
+