Problem with pdsolve command

I have a PDE system that relates four functions: sht1, svt1, Lt1, Jirt1.

I'm trying to solve this system numerically, but the pdsolve command returns an error (this error does not make sense to me).

Where am I going wrong?

Thanks!

Test.mw