Hello
I'm coupling TELEMAC3D with aed2 and I want to output the CO2 flux at the air/water interface (in regards with some water input rich in dissolved organic carbon (DIC) in my study field).
I succeeded in output variables car_DIC, car_CH4 and car_pH explained in the aed manual (the "car" term refer to carbon module of aed2) and make them appear in bluekenue. But I didn't succeed in output the variables called "Diagnostics" in the manual, which are notably CAR_pco2 (CO2 partial pressure in water) and CAR_atm_co2_flux (CO2 flux at the air/water interface).
I put here an example of aed.nml file that I found in telemac, in which we should be able to configure the output of co2 flux at the air/water interface...
! aed2_carbon - AED carbon model
!
! dic_initial [float, unit = mmol/m^3]
! initial DIC concentration (can be overwritten by GLM or TFV IC)
! Fsed_dic [float, unit = mmol/m^2/day]
! sediment CO2 flux
! Ksed_dic [float, unit = mmol/m^3]
! half-saturation oxygen concentration controlling CO2 flux
! theta_sed_dic [float, unit = - ]
! Arrhenius temperature multiplier for sediment CO2 flux
! pH_initial [float, unit = - ]
! initial water column pH
! atmco2 [float, unit = ppm]
! atmospheric CO2 concentration
! ionic [float, unit = meq]
! average ionic strength of the water column
!
! ch4_initial [float, unit = mmol/m^3]
! initial CH4 concentration (can be overwritten by GLM or TFV IC)
! Rch4ox [float, unit = - ]
! maximum reaction rate of CH4 oxidation @ 20C
! Kch4ox [float, unit = ppm]
! half-saturation oxygen concentration for CH4 oxidation
! vTch4ox [float, unit = meq]
! Arrhenius temperature multiplier for CH4 oxidation
! Fsed_ch4 [float, unit = mmol/m^2/day]
! sediment CH4 flux
! Ksed_ch4 [float, unit = mmol/m^3]
! half-saturation oxygen concentration controlling CH4 flux
! theta_sed_ch4 [float, unit = - ]
! Arrhenius temperature multiplier for sediment CH4 flux
! methane_reactant_variable
! state variable to be consumed during CH4 oxidation
!
! NOTE: Processes such as photosynthesis and bacterial respiration from other
! AED modules should be configured to use the following names:
! CAR_dic
! CAR_ch4
! CAR_pH
!
! DIAGNOSTICS:
! -
!
&aed2_carbon
!-- DIC and pH --
dic_initial = 1600.5
Fsed_dic = 3.0
Ksed_dic = 20.0
theta_sed_dic = 1.08
Fsed_dic_variable='SDF_Fsed_dic'
pH_initial = 7.5
atmco2 = 370e-6
ionic = 0.1
!-- CH4 --
ch4_initial = 27.6
Rch4ox = 0.01
Kch4ox = 0.5
vTch4ox = 1.08
Fsed_ch4 = 0.5
Ksed_ch4 = 100.0
theta_sed_ch4 = 1.08
methane_reactant_variable='OXY_oxy'
! Fsed_ch4_variable='SDF_Fsed_ch4'
/
I tried to add the line below in the aed.nml file but it doesn't work (error message in link with the use of new terms in aed2_carbon) :
! DIAGNOSTICS:
CAR_atm_co2_flux
!
&aed2_carbon
!-- DIC and pH --
co2_model = 1
alk_code = 1
ionic = 0.1
co2_piston_model = 1
Thanks a lot for your help
Raphaël
