jimmyc wrote:Well, you have a semi big task ahead of you. You need to start in solve_em.F in dyn_em/

Check the registry to see if any of the variables are specified before adding new ones.

Otherwise, it will be a matter of tracking down the subroutine where the terms are calculated. They will probably be coupled to mass, and applied as tendency terms. You will have to decide if you want coupled or uncoupled terms, and if you want instantaneous values or accumulated values.

Good luck.

Dear Sir!

I am a student, I also need to do this work, The answer that you provided has given me a big help，But i also don't know how to solve this problem. Could you describe in detail? could you give a modified solve_em.F to me? I know it is easy for you to do this work, However it is a huge task for me to do.

I am looking forward to your reply.

Thank you very much!!

cheers

tkgdewxf1983