subroutine nonline(e_density1,e_p1,e_r1,phi_ba1) implicit none include 'itg.par' include 'itg.cmn' include 'itg.mom' complex phi_ba1(kz,mz,nz) integer i c c We must calculate {a,b} for a,b=: c uv,w uv,tpar uv,v uv,T c uv,q uv,qpar uflr,w uflr,T c uflr,v uflr2,T uflr2,q c c In each case we want to pass the two variables' real and c imaginary parts to a subroutine which will return the real c and imaginary parts of {a,b} in (x,ky,kz) space. c c call nlpsce(0,phi_ba,e_density1,nl_edensity) call nlpsce(1,phi_ba,e_p1,nl_ep) call nlpsce(1,phi_ba,e_r1,nl_er) return end