sorry for that
the subroutine is kind of big .so i just give part of them ,i check the main program . it's correct . it should be something wrong with subroutine calrespknt .
the subroutine like :
subroutine calrespknt (modela,ipors,dt,t,pr,complt,modcnv)
.......
common /coff/...