I have a subroutine with a complex part and need guidance
on how I could use it with wrapit.
Below is the section which have a complex part
do 10 i=1,n
omg=2.*pi/float(i)
omg=omg*dt
z=exp(cmplx(0,-omg))
w=a*(1-z**2)/(1.+b1*z+b2*z**2)
h(i)=abs(w)**2
10 continue
return
end
Thanks
CS
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Aug 6 11:55:17 2011
This archive was generated by hypermail 2.1.8 : Fri Aug 12 2011 - 11:48:30 MDT