@@ -108,7 +108,6 @@ contains
108
108
case default
109
109
call error_stop("ERROR (moment): wrong dimension")
110
110
end select
111
-
112
111
res = res / n
113
112
114
113
end function ${RName}$
@@ -143,7 +142,8 @@ contains
143
142
case(${fi}$)
144
143
if (present(center)) then
145
144
do i = 1, size(x, ${fi}$)
146
- res = res + (real(x${select_subarray(rank, [(fi, 'i')])}$, dp) - center)**order
145
+ res = res + (real(x${select_subarray(rank, [(fi, 'i')])}$, dp) -&
146
+ center)**order
147
147
end do
148
148
else
149
149
do i = 1, size(x, ${fi}$)
@@ -154,7 +154,6 @@ contains
154
154
case default
155
155
call error_stop("ERROR (moment): wrong dimension")
156
156
end select
157
-
158
157
res = res / n
159
158
160
159
end function ${RName}$
@@ -234,7 +233,8 @@ contains
234
233
case(${fi}$)
235
234
if (present(center)) then
236
235
do i = 1, size(x, ${fi}$)
237
- res = res + merge( (x${select_subarray(rank, [(fi, 'i')])}$ - center)**order,&
236
+ res = res + merge( (x${select_subarray(rank, [(fi, 'i')])}$ -&
237
+ center)**order,&
238
238
#:if t1[0] == 'r'
239
239
0._${k1}$,&
240
240
#:else
@@ -257,7 +257,6 @@ contains
257
257
case default
258
258
call error_stop("ERROR (moment): wrong dimension")
259
259
end select
260
-
261
260
res = res / n
262
261
263
262
end function ${RName}$
@@ -281,26 +280,27 @@ contains
281
280
282
281
n = count(mask, dim)
283
282
284
- res = 0
285
- select case(dim)
286
- #:for fi in range(1, rank+1)
287
- case(${fi}$)
288
- if (present(center)) then
289
- do i = 1, size(x, ${fi}$)
290
- res = res + merge((real(x${select_subarray(rank, [(fi, 'i')])}$, dp) - center)**order,&
291
- 0._dp, mask${select_subarray(rank, [(fi, 'i')])}$)
292
- end do
293
- else
294
- do i = 1, size(x, ${fi}$)
295
- res = res + merge((real(x${select_subarray(rank, [(fi, 'i')])}$, dp))**order,&
296
- 0._dp, mask${select_subarray(rank, [(fi, 'i')])}$)
297
- end do
298
- end if
299
- #:endfor
300
- case default
301
- call error_stop("ERROR (moment): wrong dimension")
302
- end select
303
-
283
+ res = 0
284
+ select case(dim)
285
+ #:for fi in range(1, rank+1)
286
+ case(${fi}$)
287
+ if (present(center)) then
288
+ do i = 1, size(x, ${fi}$)
289
+ res = res + merge((real(x${select_subarray(rank, [(fi, 'i')])}$, dp) -&
290
+ center)**order,&
291
+ 0._dp, mask${select_subarray(rank, [(fi, 'i')])}$)
292
+ end do
293
+ else
294
+ do i = 1, size(x, ${fi}$)
295
+ res = res + merge((real(x${select_subarray(rank, [(fi, 'i')])}$, dp))&
296
+ **order,&
297
+ 0._dp, mask${select_subarray(rank, [(fi, 'i')])}$)
298
+ end do
299
+ end if
300
+ #:endfor
301
+ case default
302
+ call error_stop("ERROR (moment): wrong dimension")
303
+ end select
304
304
res = res / n
305
305
306
306
end function ${RName}$
0 commit comments