diff --git a/src/frun_CEMANEIGE.f90 b/src/frun_CEMANEIGE.f90
index 2d659015e9d55bf3563b12a058123f21a3e15947..a83b5dbd8e4e2395a1dd8db3c790881c1a6e26b3 100644
--- a/src/frun_CEMANEIGE.f90
+++ b/src/frun_CEMANEIGE.f90
@@ -201,17 +201,17 @@
! Storage of outputs
DO I=1,NOutputs
- IF(IndOutputs(I).EQ.1) Outputs(k,I)=Pliq ! Pliq ! observed liquid precipitation [mm/time step]
- IF(IndOutputs(I).EQ.2) Outputs(k,I)=Psol ! Psol ! observed solid precipitation [mm/time step]
- IF(IndOutputs(I).EQ.3) Outputs(k,I)=G ! SnowPack ! snow pack [mm]
- IF(IndOutputs(I).EQ.4) Outputs(k,I)=eTG ! ThermalState ! thermal state [°C]
- IF(IndOutputs(I).EQ.5) Outputs(k,I)=Gratio ! Gratio ! Gratio [-]
- IF(IndOutputs(I).EQ.6) Outputs(k,I)=PotMelt ! PotMelt ! potential snow melt [mm/time step]
- IF(IndOutputs(I).EQ.7) Outputs(k,I)=Melt ! Melt ! melt [mm/time step]
- IF(IndOutputs(I).EQ.8) Outputs(k,I)=PliqAndMelt ! PliqAndMelt ! liquid precipitation + melt [mm/time step]
- IF(IndOutputs(I).EQ.9) Outputs(k,I)=InputsTemp(k) ! Temp ! air temperature [°C]
- IF(IndOutputs(I).EQ.10) Outputs(k,I)=Gthreshold ! Gthreshold ! melt threshold [mm]
- IF(IndOutputs(I).EQ.11) Outputs(k,I)=Glocalmax ! Glocalmax ! local melt threshold for hysteresis [mm]
+ IF(IndOutputs(I).EQ.1) Outputs(k,I)=Pliq ! Pliq ! [numeric] observed liquid precipitation [mm/time step]
+ IF(IndOutputs(I).EQ.2) Outputs(k,I)=Psol ! Psol ! [numeric] observed solid precipitation [mm/time step]
+ IF(IndOutputs(I).EQ.3) Outputs(k,I)=G ! SnowPack ! [numeric] snow pack [mm]
+ IF(IndOutputs(I).EQ.4) Outputs(k,I)=eTG ! ThermalState ! [numeric] thermal state [°C]
+ IF(IndOutputs(I).EQ.5) Outputs(k,I)=Gratio ! Gratio ! [numeric] Gratio [-]
+ IF(IndOutputs(I).EQ.6) Outputs(k,I)=PotMelt ! PotMelt ! [numeric] potential snow melt [mm/time step]
+ IF(IndOutputs(I).EQ.7) Outputs(k,I)=Melt ! Melt ! [numeric] melt [mm/time step]
+ IF(IndOutputs(I).EQ.8) Outputs(k,I)=PliqAndMelt ! PliqAndMelt ! [numeric] liquid precipitation + melt [mm/time step]
+ IF(IndOutputs(I).EQ.9) Outputs(k,I)=InputsTemp(k) ! Temp ! [numeric] air temperature [°C]
+ IF(IndOutputs(I).EQ.10) Outputs(k,I)=Gthreshold ! Gthreshold ! [numeric] melt threshold [mm]
+ IF(IndOutputs(I).EQ.11) Outputs(k,I)=Glocalmax ! Glocalmax ! [numeric] local melt threshold for hysteresis [mm]
ENDDO
ENDDO
diff --git a/src/frun_GR4J.f90 b/src/frun_GR4J.f90
index bfba7078519d593158c70235006461b65d3f8644..c169991ac7d19fde92fb13dbda388c95b2325ccc 100644
--- a/src/frun_GR4J.f90
+++ b/src/frun_GR4J.f90
@@ -292,24 +292,24 @@
IF(Q.LT.0.) Q=0.
! Variables storage
- MISC( 1)=E ! PE ! observed potential evapotranspiration [mm/day]
- MISC( 2)=P1 ! Precip ! observed total precipitation [mm/day]
- MISC( 3)=St(1) ! Prod ! production store level (St(1)) [mm]
- MISC( 4)=PN ! Pn ! net rainfall [mm/day]
- MISC( 5)=PS ! Ps ! part of Ps filling the production store [mm/day]
- MISC( 6)=AE ! AE ! actual evapotranspiration [mm/day]
- MISC( 7)=PERC ! Perc ! percolation (PERC) [mm/day]
- MISC( 8)=PR ! PR ! PR=PN-PS+PERC [mm/day]
- MISC( 9)=StUH1(1) ! Q9 ! outflow from UH1 (Q9) [mm/day]
- MISC(10)=StUH2(1) ! Q1 ! outflow from UH2 (Q1) [mm/day]
- MISC(11)=St(2) ! Rout ! routing store level (St(2)) [mm]
- MISC(12)=EXCH ! Exch ! potential semi-exchange between catchments (EXCH) [mm/day]
- MISC(13)=AEXCH1 ! AExch1 ! actual exchange between catchments from branch 1 (AEXCH1) [mm/day]
- MISC(14)=AEXCH2 ! AExch2 ! actual exchange between catchments from branch 2 (AEXCH2) [mm/day]
- MISC(15)=AEXCH1+AEXCH2 ! AExch ! actual total exchange between catchments (AEXCH1+AEXCH2) [mm/day]
- MISC(16)=QR ! QR ! outflow from routing store (QR) [mm/day]
- MISC(17)=QD ! QD ! outflow from UH2 branch after exchange (QD) [mm/day]
- MISC(18)=Q ! Qsim ! simulated outflow at catchment outlet [mm/day]
+ MISC( 1)=E ! PE ! [numeric] observed potential evapotranspiration [mm/day]
+ MISC( 2)=P1 ! Precip ! [numeric] observed total precipitation [mm/day]
+ MISC( 3)=St(1) ! Prod ! [numeric] production store level (St(1)) [mm]
+ MISC( 4)=PN ! Pn ! [numeric] net rainfall [mm/day]
+ MISC( 5)=PS ! Ps ! [numeric] part of Ps filling the production store [mm/day]
+ MISC( 6)=AE ! AE ! [numeric] actual evapotranspiration [mm/day]
+ MISC( 7)=PERC ! Perc ! [numeric] percolation (PERC) [mm/day]
+ MISC( 8)=PR ! PR ! [numeric] PR=PN-PS+PERC [mm/day]
+ MISC( 9)=StUH1(1) ! Q9 ! [numeric] outflow from UH1 (Q9) [mm/day]
+ MISC(10)=StUH2(1) ! Q1 ! [numeric] outflow from UH2 (Q1) [mm/day]
+ MISC(11)=St(2) ! Rout ! [numeric] routing store level (St(2)) [mm]
+ MISC(12)=EXCH ! Exch ! [numeric] potential semi-exchange between catchments (EXCH) [mm/day]
+ MISC(13)=AEXCH1 ! AExch1 ! [numeric] actual exchange between catchments from branch 1 (AEXCH1) [mm/day]
+ MISC(14)=AEXCH2 ! AExch2 ! [numeric] actual exchange between catchments from branch 2 (AEXCH2) [mm/day]
+ MISC(15)=AEXCH1+AEXCH2 ! AExch ! [numeric] actual total exchange between catchments (AEXCH1+AEXCH2) [mm/day]
+ MISC(16)=QR ! QR ! [numeric] outflow from routing store (QR) [mm/day]
+ MISC(17)=QD ! QD ! [numeric] outflow from UH2 branch after exchange (QD) [mm/day]
+ MISC(18)=Q ! Qsim ! [numeric] simulated outflow at catchment outlet [mm/day]
diff --git a/src/frun_GR5H.f90 b/src/frun_GR5H.f90
index e5ec1a1cbd991fe81b03c7cbb1fea59e60de3daf..b6f04fb0c1993ca2f0172182c8d76f94284f764b 100644
--- a/src/frun_GR5H.f90
+++ b/src/frun_GR5H.f90
@@ -340,27 +340,27 @@
IF(Q.LT.0.) Q=0.
! Variables storage
- MISC( 1)=E ! PE ! observed potential evapotranspiration [mm/h]
- MISC( 2)=P1 ! Precip ! observed total precipitation [mm/h]
- MISC( 3)=St(3) ! Interc ! interception store level (St(3)) [mm]
- MISC( 4)=St(1) ! Prod ! production store level (St(1)) [mm]
- MISC( 5)=PN ! Pn ! net rainfall [mm/h]
- MISC( 6)=PS ! Ps ! part of Ps filling the production store [mm/h]
- MISC( 7)=AE ! AE ! actual evapotranspiration [mm/h]
- MISC( 8)=EI ! EI ! evapotranspiration from rainfall neutralisation or interception store [mm/h]
- MISC( 9)=ES ! ES ! evapotranspiration from production store [mm/h]
- MISC(10)=PERC ! Perc ! percolation (PERC) [mm/h]
- MISC(11)=PR ! PR ! PR=PN-PS+PERC [mm/h]
- MISC(12)=Q9 ! Q9 ! outflow from UH1 (Q9) [mm/h]
- MISC(13)=Q1 ! Q1 ! outflow from UH2 (Q1) [mm/h]
- MISC(14)=St(2) ! Rout ! routing store level (St(2)) [mm]
- MISC(15)=EXCH ! Exch ! potential semi-exchange between catchments (EXCH) [mm/h]
- MISC(16)=AEXCH1 ! AExch1 ! actual exchange between catchments from branch 1 (AEXCH1) [mm/h]
- MISC(17)=AEXCH2 ! AExch2 ! actual exchange between catchments from branch 2 (AEXCH2) [mm/h]
- MISC(18)=AEXCH1+AEXCH2 ! AExch ! actual total exchange between catchments (AEXCH1+AEXCH2) [mm/h]
- MISC(19)=QR ! QR ! outflow from routing store (QR) [mm/h]
- MISC(20)=QD ! QD ! outflow from UH2 branch after exchange (QD) [mm/h]
- MISC(21)=Q ! Qsim ! simulated outflow at catchment outlet [mm/h]
+ MISC( 1)=E ! PE ! [numeric] observed potential evapotranspiration [mm/h]
+ MISC( 2)=P1 ! Precip ! [numeric] observed total precipitation [mm/h]
+ MISC( 3)=St(3) ! Interc ! [numeric] interception store level (St(3)) [mm]
+ MISC( 4)=St(1) ! Prod ! [numeric] production store level (St(1)) [mm]
+ MISC( 5)=PN ! Pn ! [numeric] net rainfall [mm/h]
+ MISC( 6)=PS ! Ps ! [numeric] part of Ps filling the production store [mm/h]
+ MISC( 7)=AE ! AE ! [numeric] actual evapotranspiration [mm/h]
+ MISC( 8)=EI ! EI ! [numeric] evapotranspiration from rainfall neutralisation or interception store [mm/h]
+ MISC( 9)=ES ! ES ! [numeric] evapotranspiration from production store [mm/h]
+ MISC(10)=PERC ! Perc ! [numeric] percolation (PERC) [mm/h]
+ MISC(11)=PR ! PR ! [numeric] PR=PN-PS+PERC [mm/h]
+ MISC(12)=Q9 ! Q9 ! [numeric] outflow from UH1 (Q9) [mm/h]
+ MISC(13)=Q1 ! Q1 ! [numeric] outflow from UH2 (Q1) [mm/h]
+ MISC(14)=St(2) ! Rout ! [numeric] routing store level (St(2)) [mm]
+ MISC(15)=EXCH ! Exch ! [numeric] potential semi-exchange between catchments (EXCH) [mm/h]
+ MISC(16)=AEXCH1 ! AExch1 ! [numeric] actual exchange between catchments from branch 1 (AEXCH1) [mm/h]
+ MISC(17)=AEXCH2 ! AExch2 ! [numeric] actual exchange between catchments from branch 2 (AEXCH2) [mm/h]
+ MISC(18)=AEXCH1+AEXCH2 ! AExch ! [numeric] actual total exchange between catchments (AEXCH1+AEXCH2) [mm/h]
+ MISC(19)=QR ! QR ! [numeric] outflow from routing store (QR) [mm/h]
+ MISC(20)=QD ! QD ! [numeric] outflow from UH2 branch after exchange (QD) [mm/h]
+ MISC(21)=Q ! Qsim ! [numeric] simulated outflow at catchment outlet [mm/h]
ENDSUBROUTINE
diff --git a/src/frun_GR5J.f90 b/src/frun_GR5J.f90
index 47d3c79657a5a5f674902b9b0acb3c514c1fa573..9f12016f086b934880ed5a86f1af1f5871fdb856 100644
--- a/src/frun_GR5J.f90
+++ b/src/frun_GR5J.f90
@@ -279,24 +279,24 @@
IF(Q.LT.0.) Q=0.
! Variables storage
- MISC( 1)=E ! PE ! observed potential evapotranspiration [mm/day]
- MISC( 2)=P1 ! Precip ! observed total precipitation [mm/day]
- MISC( 3)=St(1) ! Prod ! production store level (St(1)) [mm]
- MISC( 4)=PN ! Pn ! net rainfall [mm/day]
- MISC( 5)=PS ! Ps ! part of Ps filling the production store [mm/day]
- MISC( 6)=AE ! AE ! actual evapotranspiration [mm/day]
- MISC( 7)=PERC ! Perc ! percolation (PERC) [mm/day]
- MISC( 8)=PR ! PR ! PR=PN-PS+PERC [mm/day]
- MISC( 9)=Q9 ! Q9 ! outflow from UH1 (Q9) [mm/day]
- MISC(10)=Q1 ! Q1 ! outflow from UH2 (Q1) [mm/day]
- MISC(11)=St(2) ! Rout ! routing store level (St(2)) [mm]
- MISC(12)=EXCH ! Exch ! potential semi-exchange between catchments (EXCH) [mm/day]
- MISC(13)=AEXCH1 ! AExch1 ! actual exchange between catchments from branch 1 (AEXCH1) [mm/day]
- MISC(14)=AEXCH2 ! AExch2 ! actual exchange between catchments from branch 2 (AEXCH2) [mm/day]
- MISC(15)=AEXCH1+AEXCH2 ! AExch ! actual total exchange between catchments (AEXCH1+AEXCH2) [mm/day]
- MISC(16)=QR ! QR ! outflow from routing store (QR) [mm/day]
- MISC(17)=QD ! QD ! outflow from UH2 branch after exchange (QD) [mm/day]
- MISC(18)=Q ! Qsim ! simulated outflow at catchment outlet [mm/day]
+ MISC( 1)=E ! PE ! [numeric] observed potential evapotranspiration [mm/day]
+ MISC( 2)=P1 ! Precip ! [numeric] observed total precipitation [mm/day]
+ MISC( 3)=St(1) ! Prod ! [numeric] production store level (St(1)) [mm]
+ MISC( 4)=PN ! Pn ! [numeric] net rainfall [mm/day]
+ MISC( 5)=PS ! Ps ! [numeric] part of Ps filling the production store [mm/day]
+ MISC( 6)=AE ! AE ! [numeric] actual evapotranspiration [mm/day]
+ MISC( 7)=PERC ! Perc ! [numeric] percolation (PERC) [mm/day]
+ MISC( 8)=PR ! PR ! [numeric] PR=PN-PS+PERC [mm/day]
+ MISC( 9)=Q9 ! Q9 ! [numeric] outflow from UH1 (Q9) [mm/day]
+ MISC(10)=Q1 ! Q1 ! [numeric] outflow from UH2 (Q1) [mm/day]
+ MISC(11)=St(2) ! Rout ! [numeric] routing store level (St(2)) [mm]
+ MISC(12)=EXCH ! Exch ! [numeric] potential semi-exchange between catchments (EXCH) [mm/day]
+ MISC(13)=AEXCH1 ! AExch1 ! [numeric] actual exchange between catchments from branch 1 (AEXCH1) [mm/day]
+ MISC(14)=AEXCH2 ! AExch2 ! [numeric] actual exchange between catchments from branch 2 (AEXCH2) [mm/day]
+ MISC(15)=AEXCH1+AEXCH2 ! AExch ! [numeric] actual total exchange between catchments (AEXCH1+AEXCH2) [mm/day]
+ MISC(16)=QR ! QR ! [numeric] outflow from routing store (QR) [mm/day]
+ MISC(17)=QD ! QD ! [numeric] outflow from UH2 branch after exchange (QD) [mm/day]
+ MISC(18)=Q ! Qsim ! [numeric] simulated outflow at catchment outlet [mm/day]
END SUBROUTINE
diff --git a/src/frun_GR6J.f90 b/src/frun_GR6J.f90
index f5a6e8a966bd0eb936d1448243972247b5ef2479..4616075333f8d26ba6f21ec204cec34af494d27f 100644
--- a/src/frun_GR6J.f90
+++ b/src/frun_GR6J.f90
@@ -316,26 +316,26 @@
IF(Q.LT.0.) Q=0.
! Variables storage
- MISC( 1)=E ! PE ! observed potential evapotranspiration [mm/day]
- MISC( 2)=P1 ! Precip ! observed total precipitation [mm/day]
- MISC( 3)=St(1) ! Prod ! production store level (St(1)) [mm]
- MISC( 4)=PN ! Pn ! net rainfall [mm/day]
- MISC( 5)=PS ! Ps ! part of Ps filling the production store [mm/day]
- MISC( 6)=AE ! AE ! actual evapotranspiration [mm/day]
- MISC( 7)=PERC ! Perc ! percolation (PERC) [mm/day]
- MISC( 8)=PR ! PR ! PR=PN-PS+PERC [mm/day]
- MISC( 9)=StUH1(1) ! Q9 ! outflow from UH1 (Q9) [mm/day]
- MISC(10)=StUH2(1) ! Q1 ! outflow from UH2 (Q1) [mm/day]
- MISC(11)=St(2) ! Rout ! routing store level (St(2)) [mm]
- MISC(12)=EXCH ! Exch ! potential third-exchange between catchments (EXCH) [mm/day]
- MISC(13)=AEXCH1 ! AExch1 ! actual exchange between catchments from routing store (AEXCH1) [mm/day]
- MISC(14)=AEXCH2 ! AExch2 ! actual exchange between catchments from direct branch (after UH2) (AEXCH2) [mm/day]
- MISC(15)=AEXCH1+AEXCH2+EXCH ! AExch ! actual total exchange between catchments (AEXCH1+AEXCH2+EXCH) [mm/day]
- MISC(16)=QR ! QR ! outflow from routing store (QR) [mm/day]
- MISC(17)=QRExp ! QRExp ! outflow from exponential store (QRExp) [mm/day]
- MISC(18)=St(3) ! Exp ! exponential store level (St(3)) (negative) [mm]
- MISC(19)=QD ! QD ! outflow from UH2 branch after exchange (QD) [mm/day]
- MISC(20)=Q ! Qsim ! simulated outflow at catchment outlet [mm/day]
+ MISC( 1)=E ! PE ! [numeric] observed potential evapotranspiration [mm/day]
+ MISC( 2)=P1 ! Precip ! [numeric] observed total precipitation [mm/day]
+ MISC( 3)=St(1) ! Prod ! [numeric] production store level (St(1)) [mm]
+ MISC( 4)=PN ! Pn ! [numeric] net rainfall [mm/day]
+ MISC( 5)=PS ! Ps ! [numeric] part of Ps filling the production store [mm/day]
+ MISC( 6)=AE ! AE ! [numeric] actual evapotranspiration [mm/day]
+ MISC( 7)=PERC ! Perc ! [numeric] percolation (PERC) [mm/day]
+ MISC( 8)=PR ! PR ! [numeric] PR=PN-PS+PERC [mm/day]
+ MISC( 9)=StUH1(1) ! Q9 ! [numeric] outflow from UH1 (Q9) [mm/day]
+ MISC(10)=StUH2(1) ! Q1 ! [numeric] outflow from UH2 (Q1) [mm/day]
+ MISC(11)=St(2) ! Rout ! [numeric] routing store level (St(2)) [mm]
+ MISC(12)=EXCH ! Exch ! [numeric] potential third-exchange between catchments (EXCH) [mm/day]
+ MISC(13)=AEXCH1 ! AExch1 ! [numeric] actual exchange between catchments from routing store (AEXCH1) [mm/day]
+ MISC(14)=AEXCH2 ! AExch2 ! [numeric] actual exchange between catchments from direct branch (after UH2) (AEXCH2) [mm/day]
+ MISC(15)=AEXCH1+AEXCH2+EXCH ! AExch ! [numeric] actual total exchange between catchments (AEXCH1+AEXCH2+EXCH) [mm/day]
+ MISC(16)=QR ! QR ! [numeric] outflow from routing store (QR) [mm/day]
+ MISC(17)=QRExp ! QRExp ! [numeric] outflow from exponential store (QRExp) [mm/day]
+ MISC(18)=St(3) ! Exp ! [numeric] exponential store level (St(3)) (negative) [mm]
+ MISC(19)=QD ! QD ! [numeric] outflow from UH2 branch after exchange (QD) [mm/day]
+ MISC(20)=Q ! Qsim ! [numeric] simulated outflow at catchment outlet [mm/day]
END SUBROUTINE