Hey @swensosc @erik ,
After the tile, excess ice and lateral snow redistribution implementations, we are now testing lateral heat flux between tiles in CTSM.
We placed the lateral heat flux function after the melting/freezing calculation but before the phase change calculation.
In our function, the temperature of each soil layer in Tile1 and Tile2 is updated based on the heat fluxes between connected layers.
Then we update errsoi_patch(p) in SoilFluxesMod.F90 --> errsoi_patch(p) = errsoi_patch(p)+eflx_lateral_col(c) (CTSM/SoilFluxesMod.F90 at be956e4139213c44ef739091eebd03134400b0a2 · mariuslam/CTSM)
The code seems to work when printing variables, but I don't manage to understand how to fix this: ERROR: ERROR in BalanceCheckMod.F90 at line 1021
The model works fine without the lateral heat fluxes. Do you think it is the updated temperatures that cause problems ?
Regards,
Marius Lambert
After the tile, excess ice and lateral snow redistribution implementations, we are now testing lateral heat flux between tiles in CTSM.
We placed the lateral heat flux function after the melting/freezing calculation but before the phase change calculation.
CTSM/SoilTemperatureMod.F90 at be956e4139213c44ef739091eebd03134400b0a2 · mariuslam/CTSM
Community Terrestrial Systems Model (includes the Community Land Model of CESM) - CTSM/SoilTemperatureMod.F90 at be956e4139213c44ef739091eebd03134400b0a2 · mariuslam/CTSM
github.com
In our function, the temperature of each soil layer in Tile1 and Tile2 is updated based on the heat fluxes between connected layers.
Then we update errsoi_patch(p) in SoilFluxesMod.F90 --> errsoi_patch(p) = errsoi_patch(p)+eflx_lateral_col(c) (CTSM/SoilFluxesMod.F90 at be956e4139213c44ef739091eebd03134400b0a2 · mariuslam/CTSM)
The code seems to work when printing variables, but I don't manage to understand how to fix this: ERROR: ERROR in BalanceCheckMod.F90 at line 1021
The model works fine without the lateral heat fluxes. Do you think it is the updated temperatures that cause problems ?
Regards,
Marius Lambert