Skip to content

Commit 3c348f8

Browse files
author
Fredrik Jansson
committed
replacing sign(1,x)*x by abs(x) in advection routines
1 parent e3cd8b4 commit 3c348f8

2 files changed

Lines changed: 128 additions & 128 deletions

File tree

‎src/advec_52.f90‎

Lines changed: 36 additions & 36 deletions
Original file line numberDiff line numberDiff line change
@@ -222,21 +222,21 @@ subroutine advecu_52(putin,putout)
222222
(&
223223
(u0(i+1,j,k)+u0(i,j,k))/60.&
224224
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
225-
-sign(1.,(u0(i+1,j,k)+u0(i,j,k)))*(u0(i+1,j,k)+u0(i,j,k))/60.&
225+
-abs(u0(i+1,j,k)+u0(i,j,k))/60.&
226226
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
227227
-(u0(i,j,k)+u0(i-1,j,k))/60.&
228228
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
229-
+sign(1.,(u0(i,j,k)+u0(i-1,j,k)))*(u0(i,j,k)+u0(i-1,j,k))/60.&
229+
+abs(u0(i,j,k)+u0(i-1,j,k))/60.&
230230
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
231231
)*dxi5 &
232232
+(&
233233
(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
234234
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
235-
-sign(1.,(v0(i,j+1,k)+v0(i-1,j+1,k)))*(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
235+
-abs(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
236236
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
237237
-(v0(i,j,k)+v0(i-1,j,k))/60.&
238238
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
239-
+sign(1.,(v0(i,j,k)+v0(i-1,j,k)))*(v0(i,j,k)+v0(i-1,j,k))/60.&
239+
+abs(v0(i,j,k)+v0(i-1,j,k))/60.&
240240
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
241241
)* dyi5 &
242242
+(1./rhobf(k))*( &
@@ -257,21 +257,21 @@ subroutine advecu_52(putin,putout)
257257
( &
258258
(u0(i+1,j,k)+u0(i,j,k))/60.&
259259
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
260-
-sign(1.,(u0(i+1,j,k)+u0(i,j,k)))*(u0(i+1,j,k)+u0(i,j,k))/60.&
260+
-abs(u0(i+1,j,k)+u0(i,j,k))/60.&
261261
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
262262
-(u0(i,j,k)+u0(i-1,j,k))/60.&
263263
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
264-
+sign(1.,(u0(i,j,k)+u0(i-1,j,k)))*(u0(i,j,k)+u0(i-1,j,k))/60.&
264+
+abs(u0(i,j,k)+u0(i-1,j,k))/60.&
265265
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
266266
)*dxi5&
267267
+(&
268268
(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
269269
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
270-
-sign(1.,(v0(i,j+1,k)+v0(i-1,j+1,k)))*(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
270+
-abs(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
271271
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
272272
-(v0(i,j,k)+v0(i-1,j,k))/60.&
273273
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
274-
+sign(1.,(v0(i,j,k)+v0(i-1,j,k)))*(v0(i,j,k)+v0(i-1,j,k))/60.&
274+
+abs(v0(i,j,k)+v0(i-1,j,k))/60.&
275275
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
276276
)* dyi5 &
277277
+ ( &
@@ -293,21 +293,21 @@ subroutine advecu_52(putin,putout)
293293
(&
294294
(u0(i+1,j,k)+u0(i,j,k))/60.&
295295
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
296-
-sign(1.,(u0(i+1,j,k)+u0(i,j,k)))*(u0(i+1,j,k)+u0(i,j,k))/60.&
296+
-abs(u0(i+1,j,k)+u0(i,j,k))/60.&
297297
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
298298
-(u0(i,j,k)+u0(i-1,j,k))/60.&
299299
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
300-
+sign(1.,(u0(i,j,k)+u0(i-1,j,k)))*(u0(i,j,k)+u0(i-1,j,k))/60.&
300+
+abs(u0(i,j,k)+u0(i-1,j,k))/60.&
301301
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
302302
)*dxi5 &
303303
+(&
304304
(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
305305
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
306-
-sign(1.,(v0(i,j+1,k)+v0(i-1,j+1,k)))*(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
306+
-abs(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
307307
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
308308
-(v0(i,j,k)+v0(i-1,j,k))/60.&
309309
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
310-
+sign(1.,(v0(i,j,k)+v0(i-1,j,k)))*(v0(i,j,k)+v0(i-1,j,k))/60.&
310+
+abs(v0(i,j,k)+v0(i-1,j,k))/60.&
311311
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
312312
)* dyi5 &
313313
+(1./rhobf(k))*( &
@@ -328,21 +328,21 @@ subroutine advecu_52(putin,putout)
328328
( &
329329
(u0(i+1,j,k)+u0(i,j,k))/60.&
330330
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
331-
-sign(1.,(u0(i+1,j,k)+u0(i,j,k)))*(u0(i+1,j,k)+u0(i,j,k))/60.&
331+
-abs(u0(i+1,j,k)+u0(i,j,k))/60.&
332332
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
333333
-(u0(i,j,k)+u0(i-1,j,k))/60.&
334334
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
335-
+sign(1.,(u0(i,j,k)+u0(i-1,j,k)))*(u0(i,j,k)+u0(i-1,j,k))/60.&
335+
+abs(u0(i,j,k)+u0(i-1,j,k))/60.&
336336
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
337337
)*dxi5&
338338
+(&
339339
(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
340340
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
341-
-sign(1.,(v0(i,j+1,k)+v0(i-1,j+1,k)))*(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
341+
-abs(v0(i,j+1,k)+v0(i-1,j+1,k))/60.&
342342
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
343343
-(v0(i,j,k)+v0(i-1,j,k))/60.&
344344
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
345-
+sign(1.,(v0(i,j,k)+v0(i-1,j,k)))*(v0(i,j,k)+v0(i-1,j,k))/60.&
345+
+abs(v0(i,j,k)+v0(i-1,j,k))/60.&
346346
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
347347
)* dyi5 &
348348
+ ( &
@@ -379,21 +379,21 @@ subroutine advecv_52(putin, putout)
379379
( &
380380
(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
381381
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
382-
-sign(1.,(u0(i+1,j,k)+u0(i+1,j-1,k)))*(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
382+
-abs(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
383383
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
384384
-(u0(i,j,k)+u0(i,j-1,k))/60.&
385385
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
386-
+sign(1.,(u0(i,j,k)+u0(i,j-1,k)))*(u0(i,j,k)+u0(i,j-1,k))/60.&
386+
+abs(u0(i,j,k)+u0(i,j-1,k))/60.&
387387
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
388388
)*dxi5&
389389
+(&
390390
(v0(i,j+1,k)+v0(i,j,k))/60.&
391391
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
392-
-sign(1.,(v0(i,j+1,k)+v0(i,j,k)))*(v0(i,j+1,k)+v0(i,j,k))/60.&
392+
-abs(v0(i,j+1,k)+v0(i,j,k))/60.&
393393
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
394394
-(v0(i,j,k)+v0(i,j-1,k))/60.&
395395
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
396-
+sign(1.,(v0(i,j,k)+v0(i,j-1,k)))*(v0(i,j,k)+v0(i,j-1,k))/60.&
396+
+abs(v0(i,j,k)+v0(i,j-1,k))/60.&
397397
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
398398
)* dyi5 &
399399
+(1./rhobf(k))*( &
@@ -414,21 +414,21 @@ subroutine advecv_52(putin, putout)
414414
( &
415415
(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
416416
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
417-
-sign(1.,(u0(i+1,j,k)+u0(i+1,j-1,k)))*(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
417+
-abs(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
418418
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
419419
-(u0(i,j,k)+u0(i,j-1,k))/60.&
420420
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
421-
+sign(1.,(u0(i,j,k)+u0(i,j-1,k)))*(u0(i,j,k)+u0(i,j-1,k))/60.&
421+
+abs(u0(i,j,k)+u0(i,j-1,k))/60.&
422422
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
423423
)*dxi5&
424424
+(&
425425
(v0(i,j+1,k)+v0(i,j,k))/60.&
426426
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
427-
-sign(1.,(v0(i,j+1,k)+v0(i,j,k)))*(v0(i,j+1,k)+v0(i,j,k))/60.&
427+
-abs(v0(i,j+1,k)+v0(i,j,k))/60.&
428428
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
429429
-(v0(i,j,k)+v0(i,j-1,k))/60.&
430430
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
431-
+sign(1.,(v0(i,j,k)+v0(i,j-1,k)))*(v0(i,j,k)+v0(i,j-1,k))/60.&
431+
+abs(v0(i,j,k)+v0(i,j-1,k))/60.&
432432
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
433433
)* dyi5 &
434434
+ ( &
@@ -447,21 +447,21 @@ subroutine advecv_52(putin, putout)
447447
( &
448448
(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
449449
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
450-
-sign(1.,(u0(i+1,j,k)+u0(i+1,j-1,k)))*(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
450+
-abs(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
451451
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
452452
-(u0(i,j,k)+u0(i,j-1,k))/60.&
453453
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
454-
+sign(1.,(u0(i,j,k)+u0(i,j-1,k)))*(u0(i,j,k)+u0(i,j-1,k))/60.&
454+
+abs(u0(i,j,k)+u0(i,j-1,k))/60.&
455455
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
456456
)*dxi5&
457457
+(&
458458
(v0(i,j+1,k)+v0(i,j,k))/60.&
459459
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
460-
-sign(1.,(v0(i,j+1,k)+v0(i,j,k)))*(v0(i,j+1,k)+v0(i,j,k))/60.&
460+
-abs(v0(i,j+1,k)+v0(i,j,k))/60.&
461461
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
462462
-(v0(i,j,k)+v0(i,j-1,k))/60.&
463463
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
464-
+sign(1.,(v0(i,j,k)+v0(i,j-1,k)))*(v0(i,j,k)+v0(i,j-1,k))/60.&
464+
+abs(v0(i,j,k)+v0(i,j-1,k))/60.&
465465
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
466466
)* dyi5 &
467467
+(1./rhobf(k))*( &
@@ -482,21 +482,21 @@ subroutine advecv_52(putin, putout)
482482
( &
483483
(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
484484
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
485-
-sign(1.,(u0(i+1,j,k)+u0(i+1,j-1,k)))*(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
485+
-abs(u0(i+1,j,k)+u0(i+1,j-1,k))/60.&
486486
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
487487
-(u0(i,j,k)+u0(i,j-1,k))/60.&
488488
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
489-
+sign(1.,(u0(i,j,k)+u0(i,j-1,k)))*(u0(i,j,k)+u0(i,j-1,k))/60.&
489+
+abs(u0(i,j,k)+u0(i,j-1,k))/60.&
490490
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
491491
)*dxi5&
492492
+(&
493493
(v0(i,j+1,k)+v0(i,j,k))/60.&
494494
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
495-
-sign(1.,(v0(i,j+1,k)+v0(i,j,k)))*(v0(i,j+1,k)+v0(i,j,k))/60.&
495+
-abs(v0(i,j+1,k)+v0(i,j,k))/60.&
496496
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
497497
-(v0(i,j,k)+v0(i,j-1,k))/60.&
498498
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
499-
+sign(1.,(v0(i,j,k)+v0(i,j-1,k)))*(v0(i,j,k)+v0(i,j-1,k))/60.&
499+
+abs(v0(i,j,k)+v0(i,j-1,k))/60.&
500500
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
501501
)* dyi5 &
502502
+(1./rhobf(k))*( &
@@ -535,21 +535,21 @@ subroutine advecw_52(putin, putout)
535535
(&
536536
(u0(i+1,j,k)+u0(i+1,j,k-1))/60.&
537537
*(37.*(putin(i+1,j,k)+putin(i,j,k))-8.*(putin(i+2,j,k)+putin(i-1,j,k))+(putin(i+3,j,k)+putin(i-2,j,k)))&
538-
-sign(1.,(u0(i+1,j,k)+u0(i+1,j,k-1)))*(u0(i+1,j,k)+u0(i+1,j,k-1))/60.&
538+
-abs(u0(i+1,j,k)+u0(i+1,j,k-1))/60.&
539539
*(10.*(putin(i+1,j,k)-putin(i,j,k))-5.*(putin(i+2,j,k)-putin(i-1,j,k))+(putin(i+3,j,k)-putin(i-2,j,k)))&
540540
-(u0(i,j,k)+u0(i,j,k-1))/60.&
541541
*(37.*(putin(i,j,k)+putin(i-1,j,k))-8.*(putin(i+1,j,k)+putin(i-2,j,k))+(putin(i+2,j,k)+putin(i-3,j,k)))&
542-
+sign(1.,(u0(i,j,k)+u0(i,j,k-1)))*(u0(i,j,k)+u0(i,j,k-1))/60.&
542+
+abs(u0(i,j,k)+u0(i,j,k-1))/60.&
543543
*(10.*(putin(i,j,k)-putin(i-1,j,k))-5.*(putin(i+1,j,k)-putin(i-2,j,k))+(putin(i+2,j,k)-putin(i-3,j,k)))&
544544
)*dxi5&
545545
+ (&
546546
(v0(i,j+1,k)+v0(i,j+1,k-1))/60.&
547547
*(37.*(putin(i,j+1,k)+putin(i,j,k))-8.*(putin(i,j+2,k)+putin(i,j-1,k))+(putin(i,j+3,k)+putin(i,j-2,k)))&
548-
-sign(1.,(v0(i,j+1,k)+v0(i,j+1,k-1)))*(v0(i,j+1,k)+v0(i,j+1,k-1))/60.&
548+
-abs(v0(i,j+1,k)+v0(i,j+1,k-1))/60.&
549549
*(10.*(putin(i,j+1,k)-putin(i,j,k))-5.*(putin(i,j+2,k)-putin(i,j-1,k))+(putin(i,j+3,k)-putin(i,j-2,k)))&
550550
-(v0(i,j,k)+v0(i,j,k-1))/60.&
551551
*(37.*(putin(i,j,k)+putin(i,j-1,k))-8.*(putin(i,j+1,k)+putin(i,j-2,k))+(putin(i,j+2,k)+putin(i,j-3,k)))&
552-
+sign(1.,(v0(i,j,k)+v0(i,j,k-1)))*(v0(i,j,k)+v0(i,j,k-1))/60.&
552+
+abs(v0(i,j,k)+v0(i,j,k-1))/60.&
553553
*(10.*(putin(i,j,k)-putin(i,j-1,k))-5.*(putin(i,j+1,k)-putin(i,j-2,k))+(putin(i,j+2,k)-putin(i,j-3,k)))&
554554
)* dyi5 &
555555
+ (1./rhobh(k))*( &

0 commit comments

Comments
 (0)