1*-----------------------------------------------------------------------
2*     Copyright (C) 2000-2004 GFD Dennou Club. All rights reserved.
3*-----------------------------------------------------------------------
4      PROGRAM MAP3D3
5      PARAMETER( NX=16, NY=16 )
6      PARAMETER( XMIN=120, XMAX=150, YMIN=20, YMAX=50 )
7      PARAMETER( FACT=0.2 )
8      REAL U(NX,NY), V(NX,NY), ALON(NX), ALAT(NY)
9
10      DO 10 I=1,NX
11        ALON(I) = XMIN + (XMAX-XMIN)*(I-1)/(NX-1)
12   10 CONTINUE
13      DO 20 J=1,NY
14        ALAT(J) = YMIN + (YMAX-YMIN)*(J-1)/(NY-1)
15   20 CONTINUE
16
17      DO 30 J=1,NY
18      DO 30 I=1,NX
19        U(I,J) =   (J-1-(NY-1)/2.) * FACT
20        V(I,J) = - (I-1-(NX-1)/2.) * FACT
21   30 CONTINUE
22
23      WRITE(*,*) ' WORKSTATION ID (I)  ? ;'
24      CALL SGPWSN
25      READ (*,*) IWS
26
27      CALL GROPN( IWS )
28      CALL GRFRM
29      CALL SGLSET( 'LCLIP', .TRUE. )
30
31      CALL GRSWND( XMIN, XMAX, YMIN, YMAX )
32      CALL GRSVPT( 0.1, 0.9, 0.1, 0.9 )
33      CALL GRSTRN( 10 )
34      CALL UMPFIT
35      CALL GRSTRF
36
37      CALL SLPVPR( 3 )
38      CALL UMRSET( 'DGRIDMJ', 10. )
39      CALL UMRSET( 'DGRIDMN',  2. )
40      CALL UMPMAP( 'coast_world' )
41      CALL UMPGLB
42
43      CALL SGSLAI( 3 )
44      DO 40 J=1,NY
45      DO 40 I=1,NX
46        CALL SGLAU( ALON(I), ALAT(J), ALON(I)+U(I,J), ALAT(J)+V(I,J) )
47   40   CONTINUE
48
49      CALL GRCLS
50
51      END
52