_RL aim_abs_pCO2, aim_ref_pCO2
& aim_dragStrato, aim_abs_pCO2, aim_ref_pCO2,
aim_ref_pCO2 = 320. _d -6
aim_fixed_pCO2 = aim_ref_pCO2
atmpCO2init = aim_ref_pCO2
& + aim_abs_pCO2*LOG( aim_pCO2/aim_ref_pCO2 )