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 | edhill | 1.2 | C $Header: /u/u3/gcmpack/MITgcm/pkg/timeave/timeave_cumul_2vfc.F,v 1.1.12.1 2003/10/02 18:30:08 adcroft Exp $ |
2 | jmc | 1.1 | C $Name: $ |
3 | edhill | 1.2 | #include "TIMEAVE_OPTIONS.h" |
4 | jmc | 1.1 | |
5 | CStartofinterface | ||
6 | SUBROUTINE TIMEAVE_CUMUL_2VFC( | ||
7 | O fldtave, | ||
8 | I fld1, fld2, fld3, Ksize, dir, deltaT, | ||
9 | I bi, bj, myThid ) | ||
10 | C /==========================================================\ | ||
11 | C | SUBROUTINE TIMEAVE_CUMUL_2VFC | | ||
12 | C | o Sum over time a product of 2 fields and 1 factor array | | ||
13 | C | depending on the relative position on the grid of the | | ||
14 | C | 2nd field & factor (tracer_grid, u_grid, v_grid). | | ||
15 | C | The first array is on tracer_grid or at same position | | ||
16 | C | as the 2 others ; 2nd and 3rd always at same position | | ||
17 | C \==========================================================/ | ||
18 | IMPLICIT NONE | ||
19 | |||
20 | C == Global variables === | ||
21 | #include "SIZE.h" | ||
22 | #include "EEPARAMS.h" | ||
23 | |||
24 | C == Routine arguments == | ||
25 | C myThid - Thread number for this instance of the routine. | ||
26 | C fldtave - time averaged Field | ||
27 | C fld1,fld2 - Input fields | ||
28 | C fld3 - Input factor | ||
29 | C dir - type of grid for 2nd field and factor, relatively to the 1rst : | ||
30 | C 0: same location ; 1: zonal,West ; 2: merid, South | ||
31 | C Ksize - 3rd dimension of local arrays (Input and Output fields) | ||
32 | INTEGER Ksize, dir | ||
33 | _RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy) | ||
34 | _RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy) | ||
35 | _RS fld3(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy) | ||
36 | _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy) | ||
37 | _RL deltaT | ||
38 | INTEGER bi, bj, myThid | ||
39 | |||
40 | CEndofinterface | ||
41 | |||
42 | C == Local variables == | ||
43 | C i,j,k,bi,bj - Loop counters | ||
44 | INTEGER i, j, k | ||
45 | |||
46 | IF ( dir.eq.0 ) THEN | ||
47 | |||
48 | C DO bj = myByLo(myThid), myByHi(myThid) | ||
49 | C DO bi = myBxLo(myThid), myBxHi(myThid) | ||
50 | DO k=1,Ksize | ||
51 | DO j=1,sNy | ||
52 | DO i=1,sNx | ||
53 | fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaT* | ||
54 | & fld1(i,j,k,bi,bj)*fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj) | ||
55 | ENDDO | ||
56 | ENDDO | ||
57 | ENDDO | ||
58 | C ENDDO | ||
59 | C ENDDO | ||
60 | |||
61 | ELSEIF ( dir.eq.1 ) THEN | ||
62 | |||
63 | C DO bj = myByLo(myThid), myByHi(myThid) | ||
64 | C DO bi = myBxLo(myThid), myBxHi(myThid) | ||
65 | DO k=1,Ksize | ||
66 | DO j=1,sNy | ||
67 | DO i=1,sNx | ||
68 | fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaT* | ||
69 | & 0.5*( fld1(i-1,j,k,bi,bj) + fld1(i,j,k,bi,bj) ) | ||
70 | & *fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj) | ||
71 | ENDDO | ||
72 | ENDDO | ||
73 | ENDDO | ||
74 | C ENDDO | ||
75 | C ENDDO | ||
76 | |||
77 | ELSEIF ( dir.eq.2 ) THEN | ||
78 | |||
79 | C DO bj = myByLo(myThid), myByHi(myThid) | ||
80 | C DO bi = myBxLo(myThid), myBxHi(myThid) | ||
81 | DO k=1,Ksize | ||
82 | DO j=1,sNy | ||
83 | DO i=1,sNx | ||
84 | fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaT* | ||
85 | & 0.5*( fld1(i,j-1,k,bi,bj) + fld1(i,j,k,bi,bj) ) | ||
86 | & *fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj) | ||
87 | ENDDO | ||
88 | ENDDO | ||
89 | ENDDO | ||
90 | C ENDDO | ||
91 | C ENDDO | ||
92 | |||
93 | ENDIF | ||
94 | |||
95 | RETURN | ||
96 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |