MODEL: Pump_Tank_2EQUATIONS FOR BERKELEY MADONNAMETHOD RK4STARTTIME = 0STOPTIME = 10DT = 0.02 {Reservoirs} d/dt (V) = + IV_in - IV_out INIT V = A*h_init {Flows} IV_in = IF (TIME < 213) THEN IF (delta_pR > p_crit) THEN k_total*SQRT(delta_pR) ELSE delta_pR/RV_total ELSE 0 IV_out = IF (TIME > 213) THEN IF (delta_pC > p_crit_2) THEN k_pipe*SQRT(delta_pC) ELSE delta_pC/RV_pipe ELSE 0 {Functions} h_init = 0.07 k_total = 1.6e-7 delta_p_Pump = 15800 delta_pR = delta_p_Pump - delta_pC delta_pC = 1000*9.81*h h = V/A k_pipe = 3e-7 dpC = delta_pC/1000+dp dp = 100.8 A = PI*(0.02^2 - 0.015^2) RV_total = 5.5e8 p_crit = 5000 RV_pipe = 1.6e8 p_crit_2 = 2000