Creating a common scaling factor for all instances

[Jorrit Mesman]

Dear all,

In a FABM model I’m working on, I am trying to create a single scaling factor for all instances of a submodel. For example, I have a growth rate “r0” in the phytoplankton-submodel, and I want to scale “r0” with a parameter “sc_r0”, but “sc_r0” must be the same for all instances. The ultimate goal is that I can calibrate “sc_r0” using parsac, and that the different phytoplankton groups retain their functional differences compared to each other.

So far, I didn’t manage to do this in FABM. I tried to make “sc_r0” a parameter in another submodule (which sets nutrient concentrations, and is initialised only once) and adding it to an aggregate_variable, but I don’t know how to retrieve this value in phytoplankton.F90 (making it a diagnostic_variable would have the same problem). I also tried to make it a dependency, but I didn’t see how to set the value.

Does anyone have an idea how to do this? A parsac-only solution would work for me as well.

Cheers,

Jorrit