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_2vfc.F,v 1.1.12.1 2003/10/02 18:30:08 adcroft Exp $ |
| 2 | C $Name: $ |
| 3 | #include "TIMEAVE_OPTIONS.h" |
| 4 | |
| 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 |