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
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
|
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
INTEGER FUNCTION IGGMEM (KPR, KERR)
C
C---->
C**** *IGGMEM*
C
C PURPOSE
C _______
C
C This routine acquires heap space.
C
C INTERFACE
C _________
C
C IERR = IGGMEM (KPR, KERR)
C
C Input parameters
C ________________
C
C KPR - The debug print switch.
C 0 , No debugging output.
C 1 , Produce debugging output.
C
C KERR - The error control flag.
C -ve, No error message. Return error code.
C 0 , Hard failure with error message.
C +ve, Print error message. Return error code.
C
C Output parameters
C ________________
C
C None
C
C Return value
C ____________
C
C The error indicator (INTEGER).
C
C Error and Warning Return Values
C _______________________________
C
C None
C
C Common block usage
C __________________
C
C nifld.common - This file contains all the input field
C definition variables.
C
C NINS - Number of grid points in NS direction for input
C field (used in grspace.h).
C NIWE - Number of grid points in WE direction for input
C field (used in grspace.h).
C
C nofld.common - This file contains all the output field
C definition variables.
C
C NONS - Number of grid points in NS direction for output
C field (used in grspace.h).
C NOWE - Number of grid points in WE direction for output
C field (used in grspace.h).
C
C grspace.h - This file contains all the work space array
C definitions for grid point to grid point
C interpolation.
C
C GRID_POINT
C
C NILATGP - POINTER to array MILATG.
C NILONGP - POINTER to array MILONG.
C NILSMP - POINTER to array MILSM.
C NINPNTP - POINTER to array RINPNT.
C NISTRTP - POINTER to array MISTRT.
C NNSDISTP - POINTER to array MNSDIST.
C NNSINDP - POINTER to array MNSIND.
C NMAXP - POINTER to array RMAX.
C NOLATGP - POINTER to array MOLATG.
C NOLONGP - POINTER to array MOLONG.
C NOLSMP - POINTER to array MOLSM.
C NWEDISTP - POINTER to array MWEDIST.
C NWEINDP - POINTER to array MWEIND.
C NWFACTP - POINTER to array WFACT.
C
C MILATG - Dynamic array of length (NINS) (see nifld.common)
C which holds the values of the input field
C latitudes.
C MILONG - Dynamic array of length (NIWE + 1) (see
C nifld.common) which holds the values of the
C input field longitudes.
C MILSM - Dynamic array of length (NIWE, 2) (see
C nifld.common) used in calculating the effects of
C the land sea mask on interpolation.
C MISTRT - Dynamic array of length (NINS) (see nifld.common)
C which holds the array offsets of the start of
C each latitude line for a quasi regular Gaussian
C input field. No space is assigned to this array
C for a regular input field.
C MNSDIST - Dynamic array of length (2, NONS) (see
C nofld.common) which holds the distances to
C neighbouring latitude lines of the input field
C from the associated line of latitude in the
C output field.
C MNSIND - Dynamic array of length (2, NONS) (see
C nofld.common) which holds the latitude line
C numbers (array offset) of the input field
C associated with each line of latitude in the
C output field.
C MOLATG - Dynamic array of length (NONS) (see nofld.common)
C which holds the values of the output field
C latitudes.
C MOLONG - Dynamic array of length (NOWE) (see nofld.common)
C which holds the values of the output field
C longitudes.
C MOLSM - Dynamic array of length (NOWE) (see nofld.common)
C used in calculating the effects of the land sea
C mask on interpolation.
C MWEDIST - Dynamic array of length (2, NONS) for a regular
C input field and of length (2, 2 * NONS * NOWE)
C (see nofld.common) for a quasi regular Gaussian
C input field. This array holds the distances to
C neighbouring longitude points of the input field
C from the associated longitude points in the
C output field.
C MWEIND - Dynamic array of length (2, NONS) for a regular
C input field and of length (2, 2 * NONS * NOWE)
C (see nofld.common) for a quasi regular Gaussian
C input field. This array holds the longitude
C points (array offset) from the input field
C associated with each longitude point in the
C output field.
C RINPNT - Dynamic array of length (NOWE) (see nofld.common)
C used to aid vectorisation in processing
C precipitation fields.
C RMAX - Dynamic array of length (NOWE) (see nofld.common)
C used to aid vectorisation in processing
C precipitation fields.
C WFACT - Dynamic array of length (4, NOWE * NONS) (see
C nofld.common) which holds the interpolation
C weights for each point in the output field.
C
C EXTERNALS
C _________
C
C IGALLOC - Get a heap allocation for array space.
C INTLOG - Logs messages.
C
C METHOD
C ______
C
C This routine calculates how much heap space is required. It
C then calls IGALLOC to acquire the required space. Finally the
C POINTERs to the arrays in grspace.h are initialised to allow
C the use of the dynamic arrays.
C
C REFERENCE
C _________
C
C None
C
C COMMENTS
C ________
C
C Program contains sections 0 to 3 and 9
C
C AUTHOR
C ______
C
C K. Fielding *ECMWF* Jan 1994
C
C MODIFICATIONS
C _____________
C
C None
C
C----<
C -----------------------------------------------------------------|
C* Section 0. Definition of variables.
C -----------------------------------------------------------------|
C
IMPLICIT NONE
C
#include "parim.h"
#include "nifld.common"
#include "nofld.common"
#include "grspace.h"
C
C Dummy arguments
C
INTEGER KPR, KERR
C
C Local variables
C
#ifdef POINTER_64
INTEGER*8 IBASE
#else
INTEGER IBASE
#endif
INTEGER IREQUEST, IERR
INTEGER JPROUTINE
PARAMETER (JPROUTINE = 22100)
C
C External functions
C
INTEGER IGALLOC
C
C -----------------------------------------------------------------|
C* Section 1. Initialisation
C -----------------------------------------------------------------|
C
100 CONTINUE
C
IF( KPR.GE.1 ) CALL INTLOG(JP_DEBUG,'IGGMEM: Section 1.',JPQUIET)
C
IGGMEM = 0
C
C -----------------------------------------------------------------|
C* Section 2. Calculate memory required and get heap
C -----------------------------------------------------------------|
C
200 CONTINUE
C
IF( KPR.GE.1 ) CALL INTLOG(JP_DEBUG,'IGGMEM: Section 2.',JPQUIET)
C
C REAL arrays (This may change for double precision version)
C
IREQUEST = 4 * NOWE * NONS * JPRLEN
C
C REAL arrays for precipitation calculations
C
IREQUEST = IREQUEST + 2 * NOWE * JPRLEN
C
C Fixed integer fields
C
IREQUEST = IREQUEST + (NIWE + 1 + NINS + NOWE + 5 * NONS) *
1 JPILEN
C
C Regular input fields (These change for reduced Gaussian fields)
C
IREQUEST = IREQUEST + 4 * NOWE * JPILEN
C
C Fields for LSM calculations
C
IREQUEST = IREQUEST + (2 * NIWE + NOWE) * JPILEN
C
C Get the required memory
C
IERR = IGALLOC (JPINNER, IREQUEST, IBASE, KPR, KERR)
C
IF( IERR .GT. 0) THEN
IGGMEM = IERR
GOTO 900
ENDIF
C
C -----------------------------------------------------------------|
C* Section 3. Initialise base addresses for POINTER arrays
C -----------------------------------------------------------------|
C
300 CONTINUE
C
IF( KPR.GE.1 ) CALL INTLOG(JP_DEBUG,'IGGMEM: Section 3.',JPQUIET)
C
C REAL arrays first to ease conversion to DOUBLE PRECISION
C
NWFACTP = IBASE
C
C Precipitation REAL arrays
C
NMAXP = NWFACTP + 4 * NOWE * NONS * JPRLEN
NINPNTP = NMAXP + NOWE * JPRLEN
C
C This offset may change if DOUBLE PRECISION
C
NILONGP = NINPNTP + NOWE * JPRLEN
C
NILATGP = NILONGP + (NIWE + 1) * JPILEN
NOLONGP = NILATGP + NINS * JPILEN
NOLATGP = NOLONGP + NOWE * JPILEN
NNSINDP = NOLATGP + NONS * JPILEN
NNSDISTP = NNSINDP + 2 * NONS * JPILEN
C
C These fields change for quasi regular Gaussian input
C
NWEINDP = NNSDISTP + 2 * NONS * JPILEN
NWEDISTP = NWEINDP + 2 * NOWE * JPILEN
NISTRTP = NWEDISTP + 2 * NOWE * JPILEN
C
C MISTRT is not use for regular input grids
C LSM arrays
C
NILSMP = NISTRTP
C
NOLSMP = NILSMP + 2 * NIWE * JPILEN
C
IF( KPR.GE.1 ) THEN
CALL INTLOG(JP_DEBUG,'IGGMEM: Output base addresses.',JPQUIET)
CALL INTLOG(JP_DEBUG,'IGGMEM: WFACT address = ',NWFACTP)
CALL INTLOG(JP_DEBUG,'IGGMEM: RMAX address = ',NMAXP )
CALL INTLOG(JP_DEBUG,'IGGMEM: RINPNT address = ',NINPNTP)
CALL INTLOG(JP_DEBUG,'IGGMEM: Input long. address = ',NILONGP)
CALL INTLOG(JP_DEBUG,'IGGMEM: Input lat. address = ',NILATGP)
CALL INTLOG(JP_DEBUG,'IGGMEM: Output long. address = ',NOLONGP)
CALL INTLOG(JP_DEBUG,'IGGMEM: Output lat. address = ',NOLATGP)
CALL INTLOG(JP_DEBUG,'IGGMEM: NS offset = ',NNSINDP)
CALL INTLOG(JP_DEBUG,'IGGMEM: NS address = ',NNSDISTP)
CALL INTLOG(JP_DEBUG,'IGGMEM: WE offset = ',NWEINDP)
CALL INTLOG(JP_DEBUG,'IGGMEM: WE address = ',NWEDISTP)
CALL INTLOG(JP_DEBUG,
X 'IGGMEM: Line start array address = ',NISTRTP)
CALL INTLOG(JP_DEBUG,'IGGMEM: Input LSM address = ',NILSMP)
CALL INTLOG(JP_DEBUG,'IGGMEM: Output LSM address = ',NOLSMP)
ENDIF
C
C -----------------------------------------------------------------|
C* Section 9. Return to calling routine. Format statements
C -----------------------------------------------------------------|
C
900 CONTINUE
C
IF( KPR.GE.1 ) CALL INTLOG(JP_DEBUG,'IGGMEM: Section 9.',JPQUIET)
C
RETURN
END
|
