I am trying to do the coupling NEMO4.2.1-FABM-ERGOM2.1 together with an external bio-optical module OASIM. The module is compiled through FABM but it needs additional meteorological forcing parameters, such as cloud liquid water content and water vapour.
Physical ocean parameters such as SST and pressure can be passed via trcsms_fabm.F90 to be used by the biogeochemical model. Also surface air pressure, apr, is being called:
CALL model%link_horizontal_data(fabm_standard_variables%surface_air_pressure, apr(:,:))
Is there a possibility to also pass surface air temperature, surface specific humidity and additionally load additional meteorological forcing parameters to be used by the biogeochemical model and the bio-optical module?
1 Like
Potentially the some of the fields : theta_zt,q_zt or theta_zu, q_zu or ssq which are similar to the fields you are looking for should be available for integration into FABM with the code you mentioned.
These variables are only available if the computation of them is activated. I assume that for these ln_blk=true is needed.
An alternative might be to use either the variables (for 3D different call for FABM needed) t_abl,q_abl (or tz1_abl, qz1_abl, tz1_dta,qz1_dta)
In my namelist ln_blk=TRUE is activated. I tried both:
CALL model%link_horizontal_data(fabm_standard_variables%surface_temperature, theta_zt(:,:))
alternatively
CALL model%link_horizontal_data(fabm_standard_variables%surface_temperature, theta_zt)
and each time I receive following error message:
link_horizontal_data_by_variable: surface_temperature: extents of provided array do not match domain extents.
What is the correct way to link surface temperature with FABM?