link ../../libf_gcm/phymars/lwdiff.F