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/model/src/correction_step.F,v 1.20.2.1 2003/10/02 18:10:45 edhill Exp $ |
2 | C $Name: $ |
3 | |
4 | #include "PACKAGES_CONFIG.h" |
5 | #include "CPP_OPTIONS.h" |
6 | |
7 | CBOP |
8 | C !ROUTINE: CORRECTION_STEP |
9 | C !INTERFACE: |
10 | SUBROUTINE CORRECTION_STEP( bi, bj, iMin, iMax, jMin, jMax, |
11 | I K, phiSurfX, phiSurfY, |
12 | I myCurrentTime, myThid ) |
13 | C !DESCRIPTION: \bv |
14 | C *==========================================================* |
15 | C | S/R CORRECTION_STEP |
16 | C | o Corrects the horizontal flow fields with the surface |
17 | C | slope. |
18 | C *==========================================================* |
19 | C \ev |
20 | |
21 | C !USES: |
22 | IMPLICIT NONE |
23 | C == Global variables == |
24 | #include "SIZE.h" |
25 | #include "DYNVARS.h" |
26 | #include "EEPARAMS.h" |
27 | #include "PARAMS.h" |
28 | #include "GRID.h" |
29 | #ifdef ALLOW_NONHYDROSTATIC |
30 | #include "SOLVE_FOR_PRESSURE3D.h" |
31 | #endif |
32 | |
33 | C !INPUT/OUTPUT PARAMETERS: |
34 | C == Routine Arguments == |
35 | C phiSurfX, phiSurfY - Surface Potential gradient |
36 | C bi,bj,iMin,iMax,jMin,jMax, K - Loop counters |
37 | C myThid - Instance number for |
38 | C this call to S/R CORRECTION_STEP |
39 | C myCurrentTime - Current simulation time for this instance. |
40 | _RL phiSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
41 | _RL phiSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
42 | INTEGER bi,bj,iMin,iMax,jMin,jMax |
43 | INTEGER K |
44 | INTEGER myThid |
45 | _RL myCurrentTime |
46 | |
47 | C !LOCAL VARIABLES: |
48 | C == Local variables == |
49 | C i,j :: Loop counters |
50 | C hxFac, hyFac :: Tracer parameters for supressing gradients |
51 | C hx3dFac,hy3dFac |
52 | INTEGER i,j |
53 | _RL hxFac,hyFac |
54 | _RL hx3dFac,hy3dFac |
55 | CEOP |
56 | |
57 | C On/off scaling paramters |
58 | hxFac = pfFacMom |
59 | hyFac = pfFacMom |
60 | IF ( nonHydrostatic ) THEN |
61 | hx3dFac = pfFacMom |
62 | hy3dFac = pfFacMom |
63 | ELSE |
64 | hx3dFac = 0. |
65 | hy3dFac = 0. |
66 | ENDIF |
67 | |
68 | C Step forward zonal velocity |
69 | DO j=jMin,jMax |
70 | DO i=iMin,iMax |
71 | uVel(i,j,k,bi,bj)=( gU(i,j,k,bi,bj) |
72 | & -deltaTmom*hxFac*implicSurfPress*phiSurfX(i,j) |
73 | #ifdef ALLOW_NONHYDROSTATIC |
74 | ceh3 needs an IF ( useNONHYDROSTATIC ) THEN |
75 | & -deltaTmom*hx3dFac*_recip_dxC(i,j,bi,bj)* |
76 | & (phi_nh(i,j,k,bi,bj)-phi_nh(i-1,j,k,bi,bj)) |
77 | #endif |
78 | & )*_maskW(i,j,k,bi,bj) |
79 | ENDDO |
80 | ENDDO |
81 | |
82 | C Step forward meridional velocity |
83 | DO j=jMin,jMax |
84 | DO i=iMin,iMax |
85 | vVel(i,j,k,bi,bj)=( gV(i,j,k,bi,bj) |
86 | & -deltaTmom*hyFac*implicSurfPress*phiSurfY(i,j) |
87 | #ifdef ALLOW_NONHYDROSTATIC |
88 | ceh3 needs an IF ( useNONHYDROSTATIC ) THEN |
89 | & -deltaTmom*hy3dFac*_recip_dyC(i,j,bi,bj)* |
90 | & (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj)) |
91 | #endif |
92 | & )*_maskS(i,j,k,bi,bj) |
93 | ENDDO |
94 | ENDDO |
95 | |
96 | RETURN |
97 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |