1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
|
C Copyright 1981-2007 ECMWF
C
C Licensed under the GNU Lesser General Public License which
C incorporates the terms and conditions of version 3 of the GNU
C General Public License.
C See LICENSE and gpl-3.0.txt for details.
C
SUBROUTINE JNUMGG( KNUM, HTYPE, KPTS, KRET)
C
C---->
C**** JNUMGG
C
C PURPOSE
C _______
C
C This routine returns the number of points in each latitude of
C a gaussian grid.
C
C
C INTERFACE
C _________
C
C CALL JNUMGG( KNUM, HTYPE, KPTS, KRET)
C
C
C Input parameters
C ________________
C
C KNUM - Gaussian grid number
C HTYPE - Gaussian grid type
C = 'R' for reduced, = 'F' for full
C
C
C Output parameters
C ________________
C
C KPTS - Array giving number of points at each gaussian latitude
C (both hemispheres)
C KRET - Return status code
C 0 = OK
C
C Common block usage
C __________________
C
C None.
C
C Method
C ______
C
C Simplified interface to JGETGG.
C
C
C Externals
C _________
C
C JGETGG - Get the gaussian grid definition.
C INTLOG - Output log message
C
C
C Reference
C _________
C
C None
C
C
C Comments
C ________
C
C None
C
C
C AUTHOR
C ______
C
C J.D.Chambers *ECMWF* Feb 1994
C
C
C MODIFICATIONS
C _____________
C
C None.
C
C----<
C
IMPLICIT NONE
#include "jparams.h"
#include "parim.h"
C
C Subroutine arguments
C
INTEGER KNUM, KPTS, KRET
CHARACTER*1 HTYPE
DIMENSION KPTS(*)
C
C Parameters
INTEGER JPROUTINE
PARAMETER ( JPROUTINE = 30700 )
C
C Local variables
C
REAL ZLAT
DIMENSION ZLAT(JPMAX)
C _______________________________________________________
C
C* Section 1. Get the information.
C _______________________________________________________
C
100 CONTINUE
C
KRET = 0
C
C Make sure local array is big enough for return values.
C
IF ( KNUM .GT. (JPMAX/2) ) THEN
CALL INTLOG(JP_ERROR,
X 'JNUMGG - Local array too small for return values.',JPMAX/2)
KRET = JPROUTINE + 1
GOTO 900
ENDIF
C
CALL JGETGG( KNUM, HTYPE, ZLAT, KPTS, KRET)
C
C _______________________________________________________
C
C* Section 9. Return to calling routine. Format statements
C _______________________________________________________
C
900 CONTINUE
C
RETURN
END
|
