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