Parent Directory | Revision Log | Revision Graph
o first check-in for the "branch-genmake2" merge o verification suite as run on shelley (gcc 3.2.2): Wed Oct 8 23:42:29 EDT 2003 T S U V G D M c m s m s m s m s E p a R g m m e . m m e . m m e . m m e . N n k u 2 i a a d i a a d i a a d i a a d 2 d e n d n x n . n x n . n x n . n x n . OPTFILE=NONE Y Y Y Y 13 16 16 16 0 16 16 16 16 16 16 16 16 13 12 0 0 pass adjustment.128x64x1 Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 0 0 16 16 0 0 pass adjustment.cs-32x32x1 Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 22 0 16 16 22 0 pass adjust_nlfs.cs-32x32x1 Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O advect_cs Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O advect_xy Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O advect_xz Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass aim.5l_cs Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass aim.5l_Equatorial_Channel Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass aim.5l_LatLon Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass exp0 Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass exp1 Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass exp2 Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass exp4 Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass exp5 Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass front_relax Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass global_ocean.90x40x15 Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL global_ocean.cs32x15 Y Y Y Y 6 11 12 13 13 12 13 16 13 9 9 9 9 10 9 9 11 FAIL global_ocean_pressure Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass global_with_exf Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16 9 16 pass hs94.128x64x5 Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass hs94.1x64x5 Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass hs94.cs-32x32x5 Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL ideal_2D_oce Y Y Y Y 8 16 16 16 16 16 16 16 16 13 13 8 16 16 16 16 16 FAIL internal_wave Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass inverted_barometer Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL lab_sea Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL natl_box Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass plume_on_slope Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass solid-body.cs-32x32x1
1 | C $Header: /u/u3/gcmpack/MITgcm/pkg/timeave/timeave_cumul_dif_1t.F,v 1.1.4.1 2003/10/02 18:30:08 adcroft Exp $ |
2 | C $Name: $ |
3 | #include "TIMEAVE_OPTIONS.h" |
4 | |
5 | CStartofinterface |
6 | SUBROUTINE TIMEAVE_CUMUL_DIF_1T( |
7 | O fldtave, |
8 | I fld1, fld2, Ksize, dir, deltaT, |
9 | I bi, bj, myThid ) |
10 | C /==========================================================* |
11 | C | SUBROUTINE TIMEAVE_CUMUL_DIF_1T |
12 | C | o Sum over time the product: |
13 | C | minus gradient of 1rst field x 2nd field (1 tile) |
14 | C \==========================================================* |
15 | IMPLICIT NONE |
16 | |
17 | C == Global variables === |
18 | #include "SIZE.h" |
19 | #include "EEPARAMS.h" |
20 | #include "GRID.h" |
21 | |
22 | C == Routine arguments == |
23 | C myThid - Thread number for this instance of the routine. |
24 | C fldtave - time averaged Field |
25 | C fld1,fld2 - Input Field |
26 | C dir - define the gradient direction: 1: X dir ; 2: Y dir ; 3: R dir |
27 | C Ksize - 3rd dimension of local arrays (Input and Output fields) |
28 | INTEGER Ksize, dir |
29 | _RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy) |
30 | _RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize) |
31 | _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy) |
32 | _RL deltaT |
33 | INTEGER bi, bj, myThid |
34 | |
35 | CEndofinterface |
36 | |
37 | #ifdef ALLOW_TIMEAVE |
38 | C == Local variables == |
39 | C i,j,k,bi,bj - Loop counters |
40 | INTEGER i, j, k |
41 | |
42 | IF ( dir.eq.1 ) THEN |
43 | |
44 | DO k=1,Ksize |
45 | DO j=1,sNy |
46 | DO i=1,sNx |
47 | fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) |
48 | & - ( fld1(i,j,k,bi,bj) - fld1(i-1,j,k,bi,bj) ) |
49 | & * recip_dxC(i,j,bi,bj) |
50 | & * fld2(i,j,k) * deltaT |
51 | ENDDO |
52 | ENDDO |
53 | ENDDO |
54 | |
55 | ELSEIF ( dir.eq.2 ) THEN |
56 | |
57 | DO k=1,Ksize |
58 | DO j=1,sNy |
59 | DO i=1,sNx |
60 | fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) |
61 | & - ( fld1(i,j,k,bi,bj) - fld1(i,j-1,k,bi,bj) ) |
62 | & * recip_dyC(i,j,bi,bj) |
63 | & * fld2(i,j,k) * deltaT |
64 | ENDDO |
65 | ENDDO |
66 | ENDDO |
67 | |
68 | ELSEIF ( dir.eq.3 ) THEN |
69 | |
70 | DO k=2,Ksize |
71 | DO j=1,sNy |
72 | DO i=1,sNx |
73 | fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) |
74 | & + ( fld1(i,j,k,bi,bj) - fld1(i,j,k-1,bi,bj) ) |
75 | & * recip_drC(k) |
76 | & * fld2(i,j,k) * deltaT |
77 | ENDDO |
78 | ENDDO |
79 | ENDDO |
80 | |
81 | ENDIF |
82 | |
83 | #endif /* ALLOW_TIMEAVE */ |
84 | |
85 | RETURN |
86 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |