c prof to test s/r profgen_jones real*8 x,pxpr,r,w0,w character*8 modx,id common r(6)/ww/id(10),w0,w(400) common / xx/ modx(2),x,pxpr open (unit=11,file='data.out',status='old') w(6) = 1.d0 w(2) = 6371.2 do 100 ih = 80,400,5 hh = ih r(1) = dble (hh) + w(2) call electx c plasma freq (recall f=1) x = dsqrt(x) write(6,200) ih,x,pxpr write(11,200) ih,x,pxpr 200 format(i10,2f10.3) 100 continue end