link ../../libf_gcm/phymars/nlthermeq.F