$NetBSD: patch-aa,v 1.2 2008/10/14 22:41:44 markd Exp $

--- src/lmm.f.orig	2003-12-04 09:19:01.000000000 +1300
+++ src/lmm.f
@@ -549,7 +549,8 @@ C     3 = inadequate information to obta
 C     4 = value of xi or inv(Ui) became non-pos.def. during iterations
 C     5 = t(X)%*%W%*%X became non-pos.def. during iterations
 C     6 = supplied xihat is non-pos.def.
-      integer gmax/50/
+      integer gmax
+      data gmax/50/
       double precision estarhat(0:50)
       integer ntot,subj(ntot),m,ist(m),ifin(m),occ(ntot),nmax,
      /     pcol,q,zcol(q),iflag,err,msg,iter,p,xcol(p),i,j,