Coupled PDEs in Coefficient form PDEs

Hello to all,

I am working on ultrashort pulsed laser heating modeling.
Two Temperature Model is used in this model. which coupled with two Partial differential equations.
Ce*(d/dt)*Te= d/dz(k*(d/dz)*Te)- G(Te-Tl)+S(z,t)
Cl*(d/dt)*Ti=G(Te-Ti)
I want to solve these equation by using coefficient form PDEs in COMSOL.
How can i coupled these two equations in my model ?

please guide me to solve this problem.

Thanks
Kiran