dwinf Subroutine

   pure subroutine dwinf &
      (n, m, np, nq, ldwe, ld2we, isodr, &
       deltai, epsi, xplusi, fni, sdi, vcvi, &
       rvari, wssi, wssdei, wssepi, rcondi, etai, &
       olmavi, taui, alphai, actrsi, pnormi, rnorsi, prersi, &
       partli, sstoli, taufci, epsmai, &
       beta0i, betaci, betasi, betani, si, ssi, ssfi, qrauxi, ui, &
       fsi, fjacbi, we1i, diffi, &
       deltsi, deltni, ti, tti, omegai, fjacdi, &
       wrk1i, wrk2i, wrk3i, wrk4i, wrk5i, wrk6i, wrk7i, &
       loweri, upperi, &
       lwkmn)
   !! Get storage locations within REAL (wp) work space.
      ! Routines Called  (NONE)
      ! Date Written   860529   (YYMMDD)
      ! Revision Date  920619   (YYMMDD)
      integer, intent(in) :: n
         !! The number of observations.
      integer, intent(in) :: m
         !! The number of columns of data in the explanatory variable.
      integer, intent(in) :: np
         !! The number of function parameters.
      integer, intent(in) :: nq
         !! The number of responses per observation.
      integer, intent(in) :: ldwe
         !! The leading dimension of array `we`.
      integer, intent(in) :: ld2we
         !! The second dimension of array `we`.
      logical, intent(in) :: isodr
         !! The variable designating whether the solution is by ODR (`isodr`=.true.) or by OLS (`isodr`=.false.).
      integer, intent(out) :: deltai
         !! The starting location in array `work` of array `delta`.
      integer, intent(out) :: epsi
         !! The starting location in array `work` of array `eps`.
      integer, intent(out) :: xplusi
         !! The starting location in array `work` of array `xplusd`.
      integer, intent(out) :: fni
         !! The starting location in array `work` of array `fn`.
      integer, intent(out) :: sdi
         !! The starting location in array `work` of array `sd`.
      integer, intent(out) :: vcvi
         !! The starting location in array `work` of array `vcv`.
      integer, intent(out) :: rvari
         !! The location in array `work` of variable `rvar`.
      integer, intent(out) :: wssi
         !! The location in array `work` of variable `wss`.
      integer, intent(out) :: wssdei
         !! The location in array `work` of variable `wssdel`.
      integer, intent(out) :: wssepi
         !! The location in array `work` of variable `wsseps`.
      integer, intent(out) :: rcondi
         !! The location in array `work` of variable `rcondi`.
      integer, intent(out) :: etai
         !! The location in array `work` of variable `eta`.
      integer, intent(out) :: olmavi
         !! The location in array `work` of variable `olmavg`.
      integer, intent(out) :: taui
         !! The location in array `work` of variable `tau`.
      integer, intent(out) :: alphai
         !! The location in array `work` of variable `alpha`.
      integer, intent(out) :: actrsi
         !! The location in array `work` of variable `actrs`.
      integer, intent(out) :: pnormi
         !! The location in array `work` of variable `pnorm`.
      integer, intent(out) :: rnorsi
         !! The location in array `work` of variable `rnorms`.
      integer, intent(out) :: prersi
         !! The location in array `work` of variable `prers`.
      integer, intent(out) :: partli
         !! The location in array `work` of variable `partol`.
      integer, intent(out) :: sstoli
         !! The location in array `work` of variable `sstol`.
      integer, intent(out) :: taufci
         !! The location in array `work` of variable `taufac`.
      integer, intent(out) :: epsmai
         !! The location in array `work` of variable `epsmac`.
      integer, intent(out) :: beta0i
         !! The starting location in array `work` of array `beta0`.
      integer, intent(out) :: betaci
         !! The starting location in array `work` of array `betac`.
      integer, intent(out) :: betasi
         !! The starting location in array `work` of array `betas`.
      integer, intent(out) :: betani
         !! The starting location in array `work` of array `betan`.
      integer, intent(out) :: si
         !! The starting location in array `work` of array `s`.
      integer, intent(out) :: ssi
         !! The starting location in array `work` of array `ss`.
      integer, intent(out) :: ssfi
         !! The starting location in array `work` of array `ssf`.
      integer, intent(out) :: qrauxi
         !! The starting location in array `work` of array `qraux`.
      integer, intent(out) :: ui
         !! The starting location in array `work` of array `u`.
      integer, intent(out) :: fsi
         !! The starting location in array `work` of array `fs`.
      integer, intent(out) :: fjacbi
         !! The starting location in array `work` of array `fjacb`.
      integer, intent(out) :: we1i
         !! The starting location in array `work` of array `we1`.
      integer, intent(out) :: diffi
         !! The starting location in array `work` of array `diff`.
      integer, intent(out) :: deltsi
         !! The starting location in array `work` of array `deltas`.
      integer, intent(out) :: deltni
         !! The starting location in array `work` of array `deltan`.
      integer, intent(out) :: ti
         !! The starting location in array `work` of array `t`.
      integer, intent(out) :: tti
         !! The starting location in array `work` of array `tt`.
      integer, intent(out) :: omegai
         !! The starting location in array `work` of array `omega`.
      integer, intent(out) :: fjacdi
         !! The starting location in array `work` of array `fjacd`.
      integer, intent(out) :: wrk1i
         !! The starting location in array `work` of array `wrk1`.
      integer, intent(out) :: wrk2i
         !! The starting location in array `work` of array `wrk2`.
      integer, intent(out) :: wrk3i
         !! The starting location in array `work` of array `wrk3`.
      integer, intent(out) :: wrk4i
         !! The starting location in array `work` of array `wrk4`.
      integer, intent(out) :: wrk5i
         !! The starting location in array `work` of array `wrk5`.
      integer, intent(out) :: wrk6i
         !! The starting location in array `work` of array `wrk6`.
      integer, intent(out) :: wrk7i
         !! The starting location in array `work` of array `wrk7`.
      integer, intent(out) :: loweri
         !! The starting location in array `work` of array `lower`.
      integer, intent(out) :: upperi
         !! The starting location in array `work` of array `upper`.
      integer, intent(out) :: lwkmn
         !! The minimum acceptable length of vector `work`.

      ! Local scalars
      integer :: next

      ! Variable Definitions (alphabetically)
      !  ACTRSI:  The location in array WORK of variable ACTRS.
      !  ALPHAI:  The location in array WORK of variable ALPHA.
      !  BETACI:  The starting location in array WORK of array BETAC.
      !  BETANI:  The starting location in array WORK of array BETAN.
      !  BETASI:  The starting location in array WORK of array BETAS.
      !  BETA0I:  The starting location in array WORK of array BETA0.
      !  DELTAI:  The starting location in array WORK of array DELTA.
      !  DELTNI:  The starting location in array WORK of array DELTAN.
      !  DELTSI:  The starting location in array WORK of array DELTAS.
      !  DIFFI:   The starting location in array WORK of array DIFF.
      !  EPSI:    The starting location in array WORK of array EPS.
      !  EPSMAI:  The location in array WORK of variable EPSMAC.
      !  ETAI:    The location in array WORK of variable ETA.
      !  FJACBI:  The starting location in array WORK of array FJACB.
      !  FJACDI:  The starting location in array WORK of array FJACD.
      !  FNI:     The starting location in array WORK of array FN.
      !  FSI:     The starting location in array WORK of array FS.
      !  ISODR:   The variable designating whether the solution is by ODR (ISODR=TRUE) or
      !           by OLS (ISODR=FALSE).
      !  LDWE:    The leading dimension of array WE.
      !  LD2WE:   The second dimension of array WE.
      !  LWKMN:   The minimum acceptable length of vector work.
      !  M:       The number of columns of data in the explanatory variable.
      !  N:       The number of observations.
      !  NEXT:    The next available location with WORK.
      !  NP:      The number of function parameters.
      !  NQ:      The number of responses per observation.
      !  OLMAVI:  The location in array WORK of variable OLMAVG.
      !  OMEGAI:  The starting location in array WORK of array OMEGA.
      !  PARTLI:  The location in array WORK of variable PARTOL.
      !  PNORMI:  The location in array WORK of variable PNORM.
      !  PRERSI:  The location in array WORK of variable PRERS.
      !  QRAUXI:  The starting location in array WORK of array QRAUX.
      !  RCONDI:  The location in array WORK of variable RCONDI.
      !  RNORSI:  The location in array WORK of variable RNORMS.
      !  RVARI:   The location in array WORK of variable RVAR.
      !  SDI:     The starting location in array WORK of array SD.
      !  SI:      The starting location in array WORK of array S.
      !  SSFI:    The starting location in array WORK of array SSF.
      !  SSI:     The starting location in array WORK of array SS.
      !  SSTOLI:  The location in array WORK of variable SSTOL.
      !  TAUFCI:  The location in array WORK of variable TAUFAC.
      !  TAUI:    The location in array WORK of variable TAU.
      !  TI:      The starting location in array WORK of array T.
      !  TTI:     The starting location in array WORK of array TT.
      !  UI:      The starting location in array WORK of array U.
      !  VCVI:    The starting location in array WORK of array VCV.
      !  WE1I:    The starting location in array WORK of array WE1.
      !  WRK1I:   The starting location in array WORK of array WRK1.
      !  WRK2I:   The starting location in array WORK of array WRK2.
      !  WRK3I:   The starting location in array WORK of array WRK3.
      !  WRK4I:   The starting location in array WORK of array WRK4.
      !  WRK5I:   The starting location in array WORK of array WRK5.
      !  WRK6I:   The starting location in array WORK of array WRK6.
      !  WRK7I:   The starting location in array WORK of array WRK7.
      !  WSSI:    The location in array WORK of variable WSS.
      !  WSSDEI:  The location in array WORK of variable WSSDEL.
      !  WSSEPI:  The location in array work of variable WSSEPS.
      !  XPLUSI:  The starting location in array WORK of array XPLUSD.

      if (n >= 1 .and. m >= 1 .and. np >= 1 .and. nq >= 1 .and. &
          ldwe >= 1 .and. ld2we >= 1) then

         deltai = 1
         epsi = deltai + n*m
         xplusi = epsi + n*nq
         fni = xplusi + n*m
         sdi = fni + n*nq
         vcvi = sdi + np
         rvari = vcvi + np*np

         wssi = rvari + 1
         wssdei = wssi + 1
         wssepi = wssdei + 1
         rcondi = wssepi + 1
         etai = rcondi + 1
         olmavi = etai + 1

         taui = olmavi + 1
         alphai = taui + 1
         actrsi = alphai + 1
         pnormi = actrsi + 1
         rnorsi = pnormi + 1
         prersi = rnorsi + 1
         partli = prersi + 1
         sstoli = partli + 1
         taufci = sstoli + 1
         epsmai = taufci + 1
         beta0i = epsmai + 1

         betaci = beta0i + np
         betasi = betaci + np
         betani = betasi + np
         si = betani + np
         ssi = si + np
         ssfi = ssi + np
         qrauxi = ssfi + np
         ui = qrauxi + np
         fsi = ui + np

         fjacbi = fsi + n*nq

         we1i = fjacbi + n*np*nq

         diffi = we1i + ldwe*ld2we*nq

         next = diffi + nq*(np + m)

         if (isodr) then
            deltsi = next
            deltni = deltsi + n*m
            ti = deltni + n*m
            tti = ti + n*m
            omegai = tti + n*m
            fjacdi = omegai + nq*nq
            wrk1i = fjacdi + n*m*nq
            next = wrk1i + n*m*nq
         else
            deltsi = deltai
            deltni = deltai
            ti = deltai
            tti = deltai
            omegai = deltai
            fjacdi = deltai
            wrk1i = deltai
         end if

         wrk2i = next
         wrk3i = wrk2i + n*nq
         wrk4i = wrk3i + np
         wrk5i = wrk4i + m*m
         wrk6i = wrk5i + m
         wrk7i = wrk6i + n*nq*np
         loweri = wrk7i + 5*nq
         upperi = loweri + np
         next = upperi + np

         lwkmn = next
      else
         deltai = 1
         epsi = 1
         xplusi = 1
         fni = 1
         sdi = 1
         vcvi = 1
         rvari = 1
         wssi = 1
         wssdei = 1
         wssepi = 1
         rcondi = 1
         etai = 1
         olmavi = 1
         taui = 1
         alphai = 1
         actrsi = 1
         pnormi = 1
         rnorsi = 1
         prersi = 1
         partli = 1
         sstoli = 1
         taufci = 1
         epsmai = 1
         beta0i = 1
         betaci = 1
         betasi = 1
         betani = 1
         si = 1
         ssi = 1
         ssfi = 1
         qrauxi = 1
         fsi = 1
         ui = 1
         fjacbi = 1
         we1i = 1
         diffi = 1
         deltsi = 1
         deltni = 1
         ti = 1
         tti = 1
         fjacdi = 1
         omegai = 1
         wrk1i = 1
         wrk2i = 1
         wrk3i = 1
         wrk4i = 1
         wrk5i = 1
         wrk6i = 1
         wrk7i = 1
         loweri = 1
         upperi = 1
         lwkmn = 1
      end if

   end subroutine dwinf