c...in order to include this file, one must have NPT,NZ among arguments dimension csy(1), csyc(1), dxp(1), dyp(1) dimension slx(npt,1), sly(npt,1), gtrz(npt,1), gtr(1), * trx(npt,1), try(npt,1), trz(npt,1), * sigx(npt,1),sigy(npt,1),sigz(npt,1) real Kmap(1), Kaspect(1) pointer (p_csy, csy),(p_csyc, csyc), * (p_dxp,dxp),(p_dyp,dyp),(p_gtr, gtr), * (p_slx,slx), (p_sly, sly), (p_gtrz,gtrz), * (p_trx,trx), (p_try, try), (p_trz, trz), * (p_sigx,sigx),(p_sigy,sigy),(p_sigz,sigz), * (p_Kmap,Kmap), (p_Kaspect,Kaspect) common /diff_data/ p_gtr, p_gtrz, p_slx, p_sly, * p_trx, p_try, p_trz, * p_sigx, p_sigy, p_sigz, p_Kmap, p_Kaspect common /metrics_other/ p_csy, p_csyc, p_dxp, p_dyp