Slab ocean Q-flux is read in from file, but not propagated to `Sfcprop%qfluxadj`
Closed this issue · 0 comments
Describe the bug
Within sfcsub.F
, infrastructure is in place to read in an ocean Q-flux from a file into a locally defined qflux
variable in the clima
subroutine:
SHiELD_physics/gsmphys/sfcsub.F
Lines 7521 to 7533 in d44f50f
However, data from this local qflux
variable is not propagated to the output qfluxadj
variable, and therefore the Q-flux assigned to Sfcprop%qfluxadj
(for eventual use elsewhere in the model) is always zero.
To Reproduce
Run SHiELD with a prescribed Q-flux file, i.e. with a file specified for the namsfc.fnqfluxc
namelist parameter. It is useful to define a diagnostic pointing to Sfcprop%qfluxadj
to assess this (I have done this in a fork).
Expected behavior
qfluxadj
should be populated in a similar manner to mldclm
, i.e. interpolated between monthly means loaded in from a file:
SHiELD_physics/gsmphys/sfcsub.F
Lines 8008 to 8012 in d44f50f
I have implemented this in a fork, and it seems to address this issue. I will make a PR after a little more testing.