...=pack(q(:,i), q(:,i) miss)) sd(i) = getsd (x=pack(q(:,i), q(:,i) miss),...= sum(x)/nr end function getmean ! ====================================================== pure function getsd (x, m) IMPLICIT NONE REAL(rk), INTENT(IN) :: x(:), m REAL(rk) :: getsd INTEGER :: i real(rk) :: nr1, resid(size(x))...(i=1:size(x)) resid(i) = x(i) - m getsd = dot_product(resid,resid)/nr1 getsd = sqrt(...

... a bug report? gfortran -mcmodel=medium -shared -O -c geek.f90 gfortran -L/home/rcrump/lib -o geek geek.o -lran -llapack -lblas /home/rcrump/lib/libran.a(getsd.o)(.text+0x6b): In function `getsd_': : undefined reference to `s_wsle' These are routines in libg2c.a You need to do something like gfortran -mcmodel=medium -shared -O -c -ff2c geek.f90 gfortran -L/...

... -lblas /home/rcrump/lib/libran.a(getsd.o)(.text+0x6b): In function `getsd_': ...`s_wsle' /home/rcrump/lib/libran.a(getsd.o)(.text+0x95): In function `getsd_': ...`do_lio' /home/rcrump/lib/libran.a(getsd.o)(.text+0xbf): In function `getsd_': ...`do_lio' /home/rcrump/lib/libran.a(getsd.o)(.text+0xc6): In function `getsd_':...e_wsle' /home/rcrump/lib/libran.a(getsd.o)(.text+0xdc): In function `getsd_':...