diff --git a/run_model_on_snellius/exe/STEMMUS_SCOPE b/run_model_on_snellius/exe/STEMMUS_SCOPE
index 4ecb2fd0..4b8d877b 100755
Binary files a/run_model_on_snellius/exe/STEMMUS_SCOPE and b/run_model_on_snellius/exe/STEMMUS_SCOPE differ
diff --git a/src/ebal.m b/src/ebal.m
index 4796caf5..13d3b80e 100644
--- a/src/ebal.m
+++ b/src/ebal.m
@@ -511,7 +511,7 @@
     fluxes.aPAR_Wm2 = Rntot_PAR; % [W m-2]      absorbed PAR
     fluxes.aPAR_Cab_eta = aPAR_Cab_eta;
     fluxes.GPP    = Actot * 12 / 1000000000;  % [kg C m-2 s-1]      gross primary production
-    fluxes.NEE      = (Actot - Resp) * 12 / 1000000000; % [kg C m-2 s-1]      net primary production
+    fluxes.NEE      = (Resp - Actot) * 12 / 1000000000; % [kg C m-2 s-1]      net primary production
 
     thermal.Ta    = Ta;       % [oC]                air temperature (as in input)
     thermal.Ts    = Ts;       % [oC]                soil temperature, sunlit and shaded [2x1]