maxdepth = k2
if(maxdepth.eq.k) goto 1000
maxdepth=k2
Pda(maxdepth,ktr) = Pd(maxdepth,ktr)
kmx = maxdepth-1
do k2=k,maxdepth-1
tda(maxdepth) = td(maxdepth)
sda(maxdepth) = sd(maxdepth)
do k2=1,maxdepth-1
do k2=k,maxdepth
parameter (maxcoord=3000, maxdepth=100)
real delZ(maxdepth)