link ../../libf_gcm/phymars/inistats.F