wrapit and complex number

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