Hi, all!
I am doing sensitivity tests in the CESM2.1.3 MARBL. I am trying to add an additional time-varying term in the equation of computing no3 and nh4 interior tendency in the marbl_interior_tendency_mod.F90 to simulate nitrate and ammonium fertilization in the ocean. I wonder how to add a time-varying term in the marbl.
Take no3 for example, the default interior tendency for no3 in the marbl is like in Line 3459:
interior_tendencies(no3_ind,k) = nitrif(k) - denitrif(k) - sed_denitrif(k) - sum(NO3_V(:,k))
While my goal is:
interior_tendencies(no3_ind,k) = nitrif(k) - denitrif(k) - sed_denitrif(k) - sum(NO3_V(:,k)) + nitrate(t)
where nitrate(t) is a time-varying additional numerical term (e.g. nitrate(t) = ln(t+1)), t is the model run time.
I wonder if it is feasible in the marbl and how to add this new variable in the output .nc files.
Thank you for any suggestions!
I am doing sensitivity tests in the CESM2.1.3 MARBL. I am trying to add an additional time-varying term in the equation of computing no3 and nh4 interior tendency in the marbl_interior_tendency_mod.F90 to simulate nitrate and ammonium fertilization in the ocean. I wonder how to add a time-varying term in the marbl.
Take no3 for example, the default interior tendency for no3 in the marbl is like in Line 3459:
interior_tendencies(no3_ind,k) = nitrif(k) - denitrif(k) - sed_denitrif(k) - sum(NO3_V(:,k))
While my goal is:
interior_tendencies(no3_ind,k) = nitrif(k) - denitrif(k) - sed_denitrif(k) - sum(NO3_V(:,k)) + nitrate(t)
where nitrate(t) is a time-varying additional numerical term (e.g. nitrate(t) = ln(t+1)), t is the model run time.
I wonder if it is feasible in the marbl and how to add this new variable in the output .nc files.
Thank you for any suggestions!