1*59599516SKenneth E. Jansen subroutine scatnu (ien, strl, xmudmi, xnut, shp) 2*59599516SKenneth E. Jansen 3*59599516SKenneth E. Jansen include "common.h" 4*59599516SKenneth E. Jansen 5*59599516SKenneth E. Jansen dimension ien(npro,nshl), strl(npro,ngauss), 6*59599516SKenneth E. Jansen & xmudmi(npro,ngauss), shp(nshl,ngauss) 7*59599516SKenneth E. Jansen dimension xnut(numnp) 8*59599516SKenneth E. Jansen 9*59599516SKenneth E. Jansen xmudmi=zero 10*59599516SKenneth E. Jansen 11*59599516SKenneth E. Jansen if(iLES.eq.5) return ! Debugging with zero-ed model 12*59599516SKenneth E. Jansen 13*59599516SKenneth E. Jansen do in = 1,nshl 14*59599516SKenneth E. Jansen do int = 1, ngauss 15*59599516SKenneth E. Jansen xmudmi(:,int) = xmudmi(:,int) + xnut(ien(:,in)) * strl(:,int) 16*59599516SKenneth E. Jansen & *shp(in,int) 17*59599516SKenneth E. Jansen enddo 18*59599516SKenneth E. Jansen enddo 19*59599516SKenneth E. Jansenc 20*59599516SKenneth E. Jansenc local clipping 21*59599516SKenneth E. Jansenc 22*59599516SKenneth E. Jansen rmu=datmat(1,2,1) 23*59599516SKenneth E. Jansen xmudmi=min(xmudmi,1000.0*rmu) !don't let it get larger than 1000 mu 24*59599516SKenneth E. Jansen xmudmi=max(xmudmi, zero) ! don't let (xmudmi) < 0 25*59599516SKenneth E. Jansenc xmudmi=max(xmudmi, -rmu) ! don't let (xmudmi + mu) < 0 26*59599516SKenneth E. Jansen return 27*59599516SKenneth E. Jansen end 28