polar plots using ncl

From: Phil Austin (paustin AT XXXXXX)
Date: Fri Feb 22 2002 - 12:20:37 MST

I would like to convert a conpack routine that does a simple polar plot of
a cloud-top radiance field to ncl or hlu, and am looking for an example
to get me started. I see the discussion of base plots at
http://ngwww.ucar.edu/ngdoc/ng/ref/hlu/obj/ContourPlot.obj.html
as well as the discussion of cnConpackParams, but I'm
uncertain about the preferred way to do this.

Here are the pertinent conpack calls we currently use:
(rad is the array holding the radiance field)

      
      call cpseti('MAP',2)
      call cpsetr('XC1',0.)
      call cpsetr('XCM',90.)
      call cpsetr('YC1',0.)
      call cpsetr('YCN',180.)
      call cpseti('SET',0)
      call set(0.,1.,0.,1.,-90.,90.,-90.,90.,1)
      call getset(fl,fr,fb,ft,ul,ur,ub,ut)
      
      
      call cpmpxy(map,xinp,yinp,xotp,yotp)
      
      call cpsetr('CWM',1.)
      
      call cpseti('CLS',0)
      call cpseti('NCL',40)
      
      do i = 1, 40
         call cpseti('PAI',i)
         call cpsetr('CLV',float(i)*dp)
         call cpseti('CLU',3)
         call cpsetr('CLL',2)
         call cpseti('AIA',i+1)
         call cpseti('AIB',i)
         call cpsetr('LLS',.5)
      end do
      
      call cprect(rad,iterx,dimsiz(3),dimsiz(2),
     $ rwrk,5000,iwrk,1000)
      call cpback(rad,rwrk,iwrk)
      call cpcldr(rad,rwrk,iwrk)
      

Thanks -- Phil Austin

This archive was generated by hypermail 2b29 : Fri Feb 22 2002 - 17:22:51 MST