use "../../data/processed/co2Up_1980_2024.nc" co2Up_ref
  use "../../data/processed/co2Up_1980_2024.nc" co2Up_linw

se wi/quality=high/thicken=2/textprom=2 1 
 se viewport ur
  plot/nolab/vlimit=-0.05:0.45 co2Up_ref
  plot/nolab/ov co2Up_linw

se wi/quality=high/thicken=2/textprom=2 2 
 let p=t[gt=co2Up_ref]
 let q=co2Up_ref
 set grid q
 go regresst 
 PLOT/SYM q
 PLOT/OVER qhat
 !list qhat => (0.2932-0.2474)/(1980-2024) ==> decrease of 0.001041 (Pg/yr^2) 

se wi/quality=high/thicken=2/textprom=2 3 
 let p=t[gt=co2Up_linw]
 let q=co2Up_linw
 set grid q
 go regresst 
 PLOT/SYM q
 PLOT/OVER qhat
 !list qhat => (0.2932-0.2474)/(1980-2024) ==> decrease of 0.005011 (Pg/yr^2) = 5.0PgC/1000yr^2


==> wind induces reduction of carbon uptake = 0.005011-0.001041=0.003970 (PgC/yr^2) = 4.0PgC/1000yr^2
==> responses from the Dietze et al