I have an array containing an area map (map of the southermost state of
brazil : rio grande do sul), as given in world coordinates attached as
rgs.dat. I modified the NCAR example: cpex08.f, v 1.5 1996/07/24 20:16:56
kennison Exp, to plot the data in MAPROJ('AE',-30.,-55.,0.).
The question is how to present the plotting only over the above mentioned
state ( rio grande do sul ) ?
- I included the following commands to BE A MASK
c OPEN AREA MAP OF RIO GRANDE DO SUL STATE : READ- LAT AND LON DATA
OPEN(97,FILE='rgs.dat',STATUS='unknown')
DO 97 I=1,120
READ(97,970, END= 98) LIXO,ALARS(I),ALORS(I)
97 CONTINUE
CLOSE(97)
98 CONTINUE
C LAST POINT = FIRST POINT IN THE CONTOUR
alars(121) = alars(1)
alors(121) = alors(1)
970 FORMAT(I1,2F8.2)
C INCLUDE 121 POINTS OF THE RIO GRANDE DO SUL STATE AREA MAP
C INTO GROUP 1 --> LAND INSIDE(2) & TAKE OCEN OUTSIDE (1)
NPT = 121
CALL AREDAM(IAMA,ALARS,ALORS,NPT,1,2,1)
and other minor changes in the maim programme ( attached as cpexy.f)
but I am unable to get the plotting only over RIO GANDE DO SUL STATE.
Actually I am trying to solve the above mentioned problem to better
understand how plotting works with mask routine. Actually I am trying to
plot the data collected by raingauges in the rio grande do sul state
(distributed as 20 scattered points of observations); use of IDSFFT
routine to get a regular grid distribution ; and plot the data over the
entire defined regular grid area. But I would like to present the results
only over the rio grande do sul state, excluding the plotting out of the
state.
best regards,
yyamazaki
This archive was generated by hypermail 2b29 : Wed Aug 01 2001 - 14:09:45 MDT