There are 3 subroutine used in a lagrangian two phase transient simulation-----drmast.f, drheat.f, sorent.f. The last one is used to add source term to continuou phase, and in this subroutine, a parameter difined in 'drmast.f'-----'EM' was used to difine the parameter S1P:

S1P = (-8240+7157+29326)*1000*EM

S2P = 0

In order to use "EM" in 'sorent.f', I difined this in subroutine:

CALL DRMAST(EM,DEMU)

However, it seems that the 'sorent.f' subroutine didn't work.IS the "call" wrong? How can I use 'EM' in this subroutine?

Sorent is called per cell, but drmast is called per droplet. EM in drmast is for a specific droplet, the one for which drmast is called. So it's not well posed to try to access EM in sorent - which droplet do you mean? There could be any number of droplets in the cell, all with a different EM.