1 SUBROUTINE OUTZ
2 1 (ASTAT, JMX, KMX , AKWD, ATTL, AUNT, IPWR, RLAT, AREA)
3 C----------------------------------------------------------------------
4 C ASTAT(JMX,KMX) ZONAL MEAN VALUE
5 C AKWD 要素についての KEYWORD
6 C ATTL 要素についてのコメント
7 C IPWR べき乗
8 C RLAT(JMX) GAUSSIAN LATITUDE
9 C AREA 'HEAD' HEADER だけ出力
10 C 'LAT ' LAT だけ出力
11 C 'KMAX' 鉛直レベル毎の値を出力
12 C 'LAND' 陸の値を出力
13 C 'SEA ' 海の値を出力
14 C 'BOTH' 海陸別・全体の値を出力
15 C 'ALL ' 全体の値を出力
16 C----------------------------------------------------------------------
17 C
18 DIMENSION ASTAT(JMX,KMX)
19 DIMENSION RLAT(JMX)
20 CHARACTER*4 AKWD
21 CHARACTER*32 ATTL
22 CHARACTER*16 AUNT
23 CHARACTER*4 AREA
24 C
25 CHARACTER*3 LEVL(3)
26 DATA LEVL/'LND','SEA','TOT'/
27 C
28 C ------------------------------------------------------
29 C >>> AREA='LAT ' なら RLAT だけを出力しておわり <<<
30 C ------------------------------------------------------
31 IF(AREA.EQ.'LAT ') THEN
32 WRITE(6,801) ( RLAT(J), J = 1, JMX )
33 RETURN
34 END IF
35 801 FORMAT(1H ,'LAT ',16F8.1)
36 C
37 C ------------------
38 C >>> HEADER <<<
39 C ------------------
40 IF(AREA.EQ.'HEAD') THEN
41 WRITE(6,802) AKWD, ATTL, AUNT, IPWR
42 RETURN
43 END IF
44 802 FORMAT(1X,A4,4X,A32,'(',A16,') **',I2)
45 C
46 C ==================================================================
47 C >>> 出力 <<<
48 C ==================================================================
49 C
50 C -------------------------------------
51 C >>> AREA='KMAX' なら LAT を出力する <<<
52 C -------------------------------------
53 IF(AREA.EQ.'KMAX') THEN
54 WRITE(6,802) AKWD, ATTL, AUNT, IPWR
55 WRITE(6,801) ( RLAT(J), J = 1, JMX )
56 END IF
57 C
58 C -----------------------------------------
59 C >>> 海陸別値を出力する <<<
60 C -----------------------------------------
61 IF( AREA.EQ.'LAND'.OR.AREA.EQ.'BOTH' ) THEN
62 WRITE(6,803) LEVL(1), ( ASTAT( J, 1 )*10**IPWR, J = 1, JMX )
63 ENDIF
64 IF( AREA.EQ.'SEA '.OR.AREA.EQ.'BOTH' ) THEN
65 WRITE(6,803) LEVL(2), ( ASTAT( J, 2 )*10**IPWR, J = 1, JMX )
66 ENDIF
67 IF( AREA.EQ.'ALL '.OR.AREA.EQ.'BOTH' ) THEN
68 WRITE(6,803) LEVL(3), ( ASTAT( J, 3 )*10**IPWR, J = 1, JMX )
69 ENDIF
70 803 FORMAT(1H ,A3,1X,16F8.1)
71 IF( AREA.EQ.'KMAX' ) THEN
72 C ------------------------------------------------------
73 C >>> AREA='KMAX' なら鉛直レベル毎の値を出力する <<<
74 C ------------------------------------------------------
75 DO 120 K=1,KMX
76 WRITE(6,804) K, ( ASTAT( J, K )*10**IPWR, J = 1, JMX )
77 120 CONTINUE
78 804 FORMAT(1H ,I2,2X,16F8.1)
79 ENDIF
80 C
81 C ===================================================================
82 C
83 RETURN
84 END SUBROUTINE OUTZ