I have experienced problems using CONPACK to display data on a world map
whenever I supply data that goes outside the area which actually appears
on the map. Below is a section from my program showing the routines I am
calling.
I find that unless the lat, lon limits (ylat1,ylatn, xlon1,xlonm) are
totally contained in the area of the globe on the map, the colouring
can be completely wrong although the contours are plotted correctly.
This means that I am not able to plot a view of the world that is not
bounded by lat-lon lines.
Does anyone know if there is some variable, flag or whatever that I
should be setting that would solve this problem?
If it matters, we have NCARG 4.0.
Thanks for any advice you can offer!!
Rosemarie Drummond
Section of program:
CALL MAPROJ (PROJ,plat,plon,rota)
CALL MAPSET ('AN - ANGLES',angl,angr,angb,angt)
CALL MAPSTI ('EL - ELLIPTICAL BOUNDARY',1)
CALL MAPSTC ('OU - OUTLINE DATASET','CO')
CALL MAPSTI ('VS - VERTICAL STRIPPING',20)
CALL MAPINT
CALL MAPBLA (IAMA)
C
C Tell CONPACK not to do the SET call (since it's already been done),
C to use mapping function 1 (EZMAP background), and what range of X and
C Y coordinates to send into the mapping function. The X coordinates
C will be treated as longitudes, and the Y coordinates will
C be treated as latitudes.
C
CALL CPSETI ('SET - DO-SET-CALL FLAG',0)
CALL CPSETI ('MAP - MAPPING FLAG',1)
CALL CPSETR ('XC1 - X COORDINATE AT I=1',xlon1)
CALL CPSETR ('XCM - X COORDINATE AT I=M',xlonm)
C
C Initialize the drawing of the contour plot.
C
CALL CPSETR ('YC1 - Y COORDINATE AT J=1',ylat1)
CALL CPSETR ('YCN - Y COORDINATE AT J=N',ylatn)
CALL CPRECT (zdat(jl1,kl1),ILON+1,nlon+1,nlat,
. RWRK,iwdim,IWRK,iwdim)
C
C Add contour lines to the area map.
C
CALL CPCLAM (zdat(jl1,kl1),RWRK,IWRK,IAMA)
C
C Color the map.
C
CALL ARSCAM (IAMA,XCRA,YCRA,ixdim,IAIA,IGIA,10,COLRAM3)
C
C put black contour lines over the color map
C
if (icline.ne.0) then
call gslwsc(1.)
CALL GSPLCI (icline)
CALL CPCLDM (zdat(jl1,kl1),RWRK,IWRK,IAMA,COLRCL)
endif
Rosemarie Drummond.
Dept of Physics, University of Toronto,
60 St. George St., Toronto, Ontario, CANADA, M5S 1A7
Tel (416)-978-8218 Fax (416)-978-8905
email: root@atmosp.physics.utoronto.ca
This archive was generated by hypermail 2b29 : Wed Jun 28 2000 - 09:45:38 MDT