1*-----------------------------------------------------------------------
2*     Copyright (C) 2000-2004 GFD Dennou Club. All rights reserved.
3*-----------------------------------------------------------------------
4      PROGRAM U2DF03
5
6      PARAMETER ( NX=19, NY=19 )
7      PARAMETER ( XMIN=0, XMAX=360, DX1=10, DX2=60 )
8      PARAMETER ( YMIN=-90, YMAX=90, MY=7, NC=3 )
9      PARAMETER ( PI=3.141592, DRAD=PI/180, DZ=0.05 )
10
11      PARAMETER ( XMNZ=60, XMXZ=300, NXZ=13, IXZ1=4 )
12      PARAMETER ( IYZ1=4, NYZ=13 )
13
14      REAL      P(NX, NY), UY1(NY), UY2(MY)
15      CHARACTER CH(MY)*(NC)
16
17      DATA      CH/ 'SP ', '60S', '30S', 'EQ ', '30N', '60N', 'NP ' /
18
19
20      DO 20 J = 1, NY
21        DO 10 I = 1, NX
22          ALON = ( XMIN + (XMAX-XMIN) * (I-1) / (NX-1) ) * DRAD
23          ALAT = ( YMIN + (YMAX-YMIN) * (J-1) / (NY-1) ) * DRAD
24          SLAT = SIN(ALAT)
25          P(I,J) = COS(ALON) * (1-SLAT**2) * SIN(2*PI*SLAT) + DZ
26   10   CONTINUE
27   20 CONTINUE
28
29      DO 30 J = 1, NY
30        UY1(J) = SIN( ( YMIN + (YMAX-YMIN) * (J-1) / (NY-1) ) * DRAD )
31   30 CONTINUE
32      DO 40 J = 1, MY
33        UY2(J) = SIN( ( YMIN + (YMAX-YMIN) * (J-1) / (MY-1) ) * DRAD )
34   40 CONTINUE
35
36      WRITE(*,*) ' WORKSTATION ID (I)  ? ;'
37      CALL SGPWSN
38      READ (*,*) IWS
39
40      CALL GROPN( IWS )
41      CALL GRFRM
42
43      CALL GRSWND( XMIN, XMAX, UY1(1), UY1(NY) )
44      CALL GRSVPT( 0.2, 0.8, 0.2, 0.8 )
45      CALL GRSTRN( 1 )
46      CALL GRSTRF
47
48      CALL UXAXDV( 'B', DX1, DX2 )
49      CALL UXAXDV( 'T', DX1, DX2 )
50      CALL UXSTTL( 'B', 'LONGITUDE', 0.0 )
51
52      CALL UYAXLB( 'L', UY1, NY, UY2, CH, NC, MY )
53      CALL UYAXLB( 'R', UY1, NY, UY2, CH, NC, MY )
54      CALL UYSTTL( 'L', 'LATITUDE', 0.0 )
55
56      CALL UWSGXB( XMNZ, XMXZ, NXZ )
57      CALL UWSGYA( UY1(IYZ1), NYZ )
58
59      CALL UDCNTR( P(IXZ1,IYZ1), NX, NXZ, NYZ )
60
61      CALL GRCLS
62
63      END
64