When you say iTEBD, do you mean to find the ground state or to do real-time evolution? Since you referrred to MBL, I'm assuming you mean real time evolution but I'll comment on both.
Here is a "code formula" page that shows how to make a code to do time evolution of MPS using Trotter gates, which is the best approach if you have a short ranged Hamiltonian:
To include disorder, just define the individual Hamiltonian terms ("hterm" in the code) appropriately, with parameters depending on the bond number "b". Be careful that you use the same parameters for the reverse set of gates defined in the second loop going from N-1 down to 1.
If you pass the Gate::tReal flag as shown in the sample code, it will produce real-time evolving gates of the form exp(-i h tau). If you pass the Gate::tImag flag it will do imaginary time evolution and you can use that to project into the ground state.
But for finding the ground state, a more powerful approach is to use DMRG. Be careful, though, as DMRG can get stuck in an eigenstate which isn't the ground state if disorder is extremely strong (and depending on the initial state you give).