polar plots using ncl

From: Phil Austin (paustin AT unknown)
Date: Thu Feb 21 2002 - 18:34:25 MST

  • Next message: Erika Harnett: "semi-successful compile of NCARG with Intel Fortran compiler for Linux"

    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
    http://ngwww.ucar.edu/ngdoc/ng/ref/hlu/obj/ContourPlot.obj.html
    at, 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:20:35 MST