subroutine IDIGT (input,iunit,itens,ihund,ithou) c c returns the values of the digits making up the supplied integer, c currently up to a maximum size for input of 9999 c ithou=0 ihund=0 itens=0 iunit=0 ipart=iabs(input) ndig=alog10(float(ipart))+1 if (ndig.gt.4) go to 50 go to (40,30,20,10),ndig c 10 ithou=ipart/1000 ipart=ipart-ithou*1000 c 20 ihund=ipart/100 ipart=ipart-ihund*100 c 30 itens=ipart/10 ipart=ipart-itens*10 c 40 iunit=ipart return c 50 type 90 90 format(' ?IDIGT-F-Error: input too large') call exit end