There was a bug in the heat flux (1433) calculation because tmp1d was allocated, but never initialized to zero. Often it would be zero and so the code would run fine, but occasionally NaNs were output in 1433. The initialization should fix this issue.
There was a bug in the heat flux (1433) calculation because tmp1d was allocated, but never initialized to zero. Often it would be zero and so the code would run fine, but occasionally NaNs were output in 1433. The initialization should fix this issue.