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/seaice/diffus.F,v 1.7.2.1 2003/10/02 18:18:34 adcroft Exp $ |
2 | C $Name: $ |
3 | |
4 | #include "SEAICE_OPTIONS.h" |
5 | |
6 | CStartOfInterface |
7 | SUBROUTINE diffus( HEFF, DIFFA, HEFFM, DELTT, myThid ) |
8 | C /==========================================================\ |
9 | C | SUBROUTINE diffus | |
10 | C | o Add diffusion terms to ice mass conservation equations | |
11 | C |==========================================================| |
12 | C \==========================================================/ |
13 | IMPLICIT NONE |
14 | |
15 | C === Global variables === |
16 | #include "SIZE.h" |
17 | #include "EEPARAMS.h" |
18 | #include "SEAICE_GRID.h" |
19 | |
20 | C === Routine arguments === |
21 | C myThid - Thread no. that called this routine. |
22 | _RL HEFF (1-OLx:sNx+OLx,1-OLy:sNy+OLy,3,nSx,nSy) |
23 | _RL HEFFM (1-OLx:sNx+OLx,1-OLy:sNy+OLy, nSx,nSy) |
24 | _RL DIFFA (1-OLx:sNx+OLx,1-OLy:sNy+OLy, nSx,nSy) |
25 | _RL DELTT |
26 | INTEGER myThid |
27 | CEndOfInterface |
28 | |
29 | #ifdef ALLOW_SEAICE |
30 | |
31 | C === Local variables === |
32 | C i,j,bi,bj - Loop counters |
33 | |
34 | INTEGER i, j, bi, bj |
35 | _RL DELTXX1, DELTYY1, DELTXX, DELTYY |
36 | _RL HEFF1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
37 | |
38 | DO bj=myByLo(myThid),myByHi(myThid) |
39 | DO bi=myBxLo(myThid),myBxHi(myThid) |
40 | |
41 | DO j=1-OLy,sNy+OLy |
42 | DO i=1-OLx,sNx+OLx |
43 | HEFF1(I,J,bi,bj)=0.0 _d 0 |
44 | ENDDO |
45 | ENDDO |
46 | |
47 | C NOW DO DIFFUSION WITH NUIT CONVERSION |
48 | DO j=1,sNy |
49 | DO i=1,sNx |
50 | DELTXX1=DELTT*DIFFA(I,J,bi,bj) |
51 | DELTYY1=DELTT*DIFFA(I,J,bi,bj) |
52 | DELTXX=DELTXX1/(DXTICE(I,J,bi,bj)*DXTICE(I,J,bi,bj) |
53 | & *CSTICE(I,J,bi,bj)*CSTICE(I,J,bi,bj)) |
54 | DELTYY=DELTYY1/(DYTICE(I,J,bi,bj)*DYTICE(I,J,bi,bj) |
55 | & *CSTICE(I,J,bi,bj)) |
56 | HEFF1(I,J,bi,bj)=DELTXX*((HEFF(I+1,J,3,bi,bj) |
57 | & -HEFF(I,J,3,bi,bj))*HEFFM(I+1,J,bi,bj) |
58 | & -(HEFF(I,J,3,bi,bj)-HEFF(I-1,J,3,bi,bj))*HEFFM(I-1,J,bi,bj)) |
59 | & +DELTYY*((HEFF(I,J+1,3,bi,bj)-HEFF(I,J,3,bi,bj)) |
60 | & *CSUICE(I+1,J+1,bi,bj)*HEFFM(I,J+1,bi,bj) |
61 | & -(HEFF(I,J,3,bi,bj)-HEFF(I,J-1,3,bi,bj)) |
62 | & *CSUICE(I+1,J,bi,bj)*HEFFM(I,J-1,bi,bj)) |
63 | ENDDO |
64 | ENDDO |
65 | |
66 | ENDDO |
67 | ENDDO |
68 | |
69 | _EXCH_XY_R8(HEFF1, myThid) |
70 | |
71 | DO bj=myByLo(myThid),myByHi(myThid) |
72 | DO bi=myBxLo(myThid),myBxHi(myThid) |
73 | DO j=1-OLy,sNy+OLy |
74 | DO i=1-OLx,sNx+OLx |
75 | HEFF(I,J,3,bi,bj)=HEFF1(I,J,bi,bj) |
76 | ENDDO |
77 | ENDDO |
78 | ENDDO |
79 | ENDDO |
80 | |
81 | #endif /* ALLOW_SEAICE */ |
82 | |
83 | RETURN |
84 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |