subroutine collis(ws3,wc3,vs3,vc3,tpars3,tparc3,ts3,tc3, * qs3,qc3,qpars3,qparc3,dt) c c Advance collisional part of equations in time. c implicit none include 'itg.par' include 'itg.cmn' real wtmp1(lz,mz,nz,nspecz),wtmp2(lz,mz,nz,nspecz), * vc3(lz,mz,nz,nspecz),vs3(lz,mz,nz,nspecz) real ws3(lz,mz,nz,nspecz),wc3(lz,mz,nz,nspecz), * tpars3(lz,mz,nz,nspecz),tparc3(lz,mz,nz,nspecz), * ts3(lz,mz,nz,nspecz),tc3(lz,mz,nz,nspecz), * qs3(lz,mz,nz,nspecz),qc3(lz,mz,nz,nspecz),dt, * qparc3(lz,mz,nz,nspecz),qpars3(lz,mz,nz,nspecz) if(rmu1.lt.1.0e-33) goto 23 c call poismu(wtmp1,ws3,-1,rmu1,dt) call poismu(wtmp2,wc3,-1,rmu1,dt) call copy(wtmp1,ws3) call copy(wtmp2,wc3) 23 continue c if(tkap.lt.1.0e-33) goto 241 call poismu(wtmp1,ts3,-1,tkap,dt) call poismu(wtmp2,tc3,-1,tkap,dt) call copy(wtmp1,ts3) call copy(wtmp2,tc3) 241 continue c if(qkap.lt.1.0e-33) goto 242 call poismu(wtmp1,qs3,-1,qkap,dt) call poismu(wtmp2,qc3,-1,qkap,dt) call copy(wtmp1,qs3) call copy(wtmp2,qc3) 242 continue c if(qkappar.lt.1.0e-33) goto 243 call poismu(wtmp1,qpars3,-1,qkappar,dt) call poismu(wtmp2,qparc3,-1,qkappar,dt) call copy(wtmp1,qpars3) call copy(wtmp2,qparc3) 243 continue c if(rkap.lt.1.0e-33) goto 24 call poismu(wtmp1,tpars3,-1,rkap,dt) call poismu(wtmp2,tparc3,-1,rkap,dt) call copy(wtmp1,tpars3) call copy(wtmp2,tparc3) 24 continue c if(rmu2.lt.1.0e-33) goto 25 call poismu(wtmp1,vs3,-1,rmu2,dt) call poismu(wtmp2,vc3,-1,rmu2,dt) call copy(wtmp1,vs3) call copy(wtmp2,vc3) 25 continue c return end