link ../../libf_gcm/phymars/mkstat.F90