icebalance module
Performs ice balance for the non-RGH case.
- param jmx:
Number of grid cells (integer)
- param ice:
1D numpy array of grid cell indices (0-indexed) where ice exists
- param notice:
1D numpy array of grid cell indices (0-indexed) where no ice is present (ocean)
- param conduct:
Conduction constant (scalar)
- param h:
Sea ice thickness array (length jmx)
- param Tfrz:
Freezing temperature (scalar, e.g., -2)
- param rprimew:
Radiative forcing for water (array of length jmx)
- param Cw_delt:
Ocean scaling factor (scalar; typically Cw/delt)
- param M:
Matrix of size (2*jmx, 2*jmx)
- param r:
Right-hand side vector of length (2*jmx,)
- param Diff_Op:
Diffusion operator matrix of size (jmx, jmx)
- param B:
Parameter (scalar)
- param nu_fw:
Array of water nu values (length jmx)
- param fw:
Ocean fraction array (length jmx)
- param delx:
Spatial grid spacing (scalar)
- param Cw:
Parameter for ocean (scalar)
- param W:
Water temperature array (length jmx)
- returns:
Full temperature vector (length 2*jmx; interleaved land and ocean) L : Land temperature array (length jmx) W_new : Updated water temperature array (length jmx) Fnet : Net flux array (flattened to 1D, length jmx)
- rtype:
T