C
C There are a total of 13 entries in the algebraic variable array.
C There are a total of 2 entries in each of the rate and state variable arrays.
C There are a total of 8 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is PPC in component pulmonary_fluid_dynamics (mmHg).
C CONSTS(2) is PPA in component pulmonary_fluid_dynamics (mmHg).
C CONSTS(3) is PLA in component pulmonary_fluid_dynamics (mmHg).
C CONSTS(4) is CPP in component pulmonary_fluid_dynamics (gram_per_L).
C CONSTS(5) is RPV in component pulmonary_fluid_dynamics (mmHg_minute_per_L).
C CONSTS(6) is RPA in component pulmonary_fluid_dynamics (mmHg_minute_per_L).
C CONSTS(8) is PCP in component pulmonary_capillary_pressure (mmHg).
C ALGBRC(5) is POS in component colloid_osmotic_pressure_of_pulmonary_interstitium (mmHg).
C ALGBRC(2) is PPI in component pulmonary_interstitial_fluid_pressure (mmHg).
C ALGBRC(6) is PFI in component fluid_filtration_into_pulmonary_interstitium (L_per_minute).
C CONSTS(7) is CPF in component parameter_values (L_per_minute_per_mmHg).
C ALGBRC(8) is PLF in component lung_lymphatic_protein_flow (L_per_minute).
C ALGBRC(11) is DFP in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
C ALGBRC(1) is VPF in component pulmonary_interstitial_free_fluid_volume (litre).
C ALGBRC(9) is DFZ in component pulmonary_interstitial_free_fluid_volume (L_per_minute).
C STATES(1) is VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
C ALGBRC(10) is PPO in component lung_lymphatic_protein_flow (gram_per_minute).
C ALGBRC(7) is PPN in component protein_leakage_into_pulmonary_interstitium (gram_per_minute).
C ALGBRC(13) is PPD in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
C ALGBRC(4) is CPN in component concentration_of_protein_in_pulmonary_interstitium (gram_per_L).
C ALGBRC(12) is PPZ in component concentration_of_protein_in_pulmonary_interstitium (gram_per_minute).
C STATES(2) is PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
C ALGBRC(3) is PPR in component concentration_of_protein_in_pulmonary_interstitium (gram).
C RATES(1) is d/dt VPF1 in component pulmonary_interstitial_free_fluid_volume (litre).
C RATES(2) is d/dt PPR1 in component concentration_of_protein_in_pulmonary_interstitium (gram).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 29.9941
      CONSTS(2) = 15.6376
      CONSTS(3) = 2
      CONSTS(4) = 71.9719
      CONSTS(5) = 1.55719
      CONSTS(6) = 1.5683
      CONSTS(7) = 0.0003
      STATES(1) = 0.0123238
      STATES(2) = 0.419998
      CONSTS(8) = ( (CONSTS(2) - CONSTS(3))*CONSTS(5))/(CONSTS(5)+CONSTS(6))+CONSTS(3)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(STATES(1).LT.0.00100000, 0.00100000, STATES(1))
      ALGBRC(3) = TERNRY(STATES(2).LT.0.0250000, 0.0250000, STATES(2))
      ALGBRC(4) = ALGBRC(3)/ALGBRC(1)
      ALGBRC(5) =  ALGBRC(4)*0.400000
      ALGBRC(2) = 2.00000 - 0.150000/ALGBRC(1)
      ALGBRC(6) =  (((CONSTS(8) - ALGBRC(2))+ALGBRC(5)) - CONSTS(1))*CONSTS(7)
      ALGBRC(8) =  (ALGBRC(2)+11.0000)*0.000300000
      ALGBRC(9) = ALGBRC(6) - ALGBRC(8)
      ALGBRC(11) = ALGBRC(9)
      RATES(1) = ALGBRC(11)
      ALGBRC(10) =  ALGBRC(8)*ALGBRC(4)
      ALGBRC(7) =  (CONSTS(4) - ALGBRC(4))*0.000225000
      ALGBRC(12) = ALGBRC(7) - ALGBRC(10)
      ALGBRC(13) = ALGBRC(12)
      RATES(2) = ALGBRC(13)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(STATES(1).LT.0.00100000, 0.00100000, STATES(1))
      ALGBRC(3) = TERNRY(STATES(2).LT.0.0250000, 0.0250000, STATES(2))
      ALGBRC(4) = ALGBRC(3)/ALGBRC(1)
      ALGBRC(5) =  ALGBRC(4)*0.400000
      ALGBRC(2) = 2.00000 - 0.150000/ALGBRC(1)
      ALGBRC(6) =  (((CONSTS(8) - ALGBRC(2))+ALGBRC(5)) - CONSTS(1))*CONSTS(7)
      ALGBRC(8) =  (ALGBRC(2)+11.0000)*0.000300000
      ALGBRC(9) = ALGBRC(6) - ALGBRC(8)
      ALGBRC(11) = ALGBRC(9)
      ALGBRC(10) =  ALGBRC(8)*ALGBRC(4)
      ALGBRC(7) =  (CONSTS(4) - ALGBRC(4))*0.000225000
      ALGBRC(12) = ALGBRC(7) - ALGBRC(10)
      ALGBRC(13) = ALGBRC(12)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END