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/update_surf_dr.F,v 1.3.6.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: UPDATE_SURF_DR |
9 | C !INTERFACE: |
10 | SUBROUTINE UPDATE_SURF_DR( myTime, myIter, myThid ) |
11 | C !DESCRIPTION: \bv |
12 | C *==========================================================* |
13 | C | SUBROUTINE UPDATE_SURF_DR |
14 | C | o Update the surface-level thickness fraction (hFacC,W,S) |
15 | C | according to the surface r-position = Non-Linear FrSurf |
16 | C *==========================================================* |
17 | C \ev |
18 | |
19 | C !USES: |
20 | IMPLICIT NONE |
21 | C == Global variables |
22 | #include "SIZE.h" |
23 | #include "EEPARAMS.h" |
24 | #include "PARAMS.h" |
25 | c #include "DYNVARS.h" |
26 | #include "GRID.h" |
27 | #include "SURFACE.h" |
28 | |
29 | C !INPUT/OUTPUT PARAMETERS: |
30 | C == Routine arguments == |
31 | C myTime - Current time in simulation |
32 | C myIter - Current iteration number in simulation |
33 | C myThid - Thread number for this instance of the routine. |
34 | _RL myTime |
35 | INTEGER myIter |
36 | INTEGER myThid |
37 | |
38 | C !LOCAL VARIABLES: |
39 | #ifdef NONLIN_FRSURF |
40 | C Local variables |
41 | C i,j,k,bi,bj - loop counter |
42 | INTEGER i,j,k,bi,bj |
43 | INTEGER ks |
44 | CEOP |
45 | |
46 | DO bj=myByLo(myThid), myByHi(myThid) |
47 | DO bi=myBxLo(myThid), myBxHi(myThid) |
48 | |
49 | C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
50 | |
51 | #ifdef ALLOW_OBCS |
52 | C-- Apply OBC to hFac_surfW,S before updating hFacW,S |
53 | IF (useOBCS) CALL OBCS_APPLY_SURF_DR( |
54 | I bi, bj, ksurfC, ksurfW, ksurfS, |
55 | U hFac_surfC, hFac_surfW, hFac_surfS, |
56 | I myThid ) |
57 | #endif /* ALLOW_OBCS */ |
58 | |
59 | C-- Update the fractional thickness "hFacC" of the surface level ksurfC : |
60 | DO j=1-Oly,sNy+Oly |
61 | DO i=1-Olx,sNx+Olx |
62 | ks = ksurfC(i,j,bi,bj) |
63 | IF (ks.LE.Nr) THEN |
64 | hFacC(i,j,ks,bi,bj) = hFac_surfC(i,j,bi,bj) |
65 | recip_hFacC(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfC(i,j,bi,bj) |
66 | ENDIF |
67 | ENDDO |
68 | ENDDO |
69 | |
70 | C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
71 | C-- Update fractional thickness "hFacW" & "hFacS" (at U and V points) |
72 | |
73 | DO j=1-Oly,sNy+Oly |
74 | DO i=2-Olx,sNx+Olx |
75 | ks = ksurfW(i,j,bi,bj) |
76 | IF (ks.LE.Nr) THEN |
77 | hFacW(i,j,ks,bi,bj) = hFac_surfW(i,j,bi,bj) |
78 | recip_hFacW(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfW(i,j,bi,bj) |
79 | ENDIF |
80 | ENDDO |
81 | ENDDO |
82 | DO j=2-Oly,sNy+Oly |
83 | DO i=1-Olx,sNx+Olx |
84 | ks = ksurfS(i,j,bi,bj) |
85 | IF (ks.LE.Nr) THEN |
86 | hFacS(i,j,ks,bi,bj) = hFac_surfS(i,j,bi,bj) |
87 | recip_hFacS(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfS(i,j,bi,bj) |
88 | ENDIF |
89 | ENDDO |
90 | ENDDO |
91 | |
92 | C- end bi,bj loop |
93 | ENDDO |
94 | ENDDO |
95 | |
96 | C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
97 | #endif /* NONLIN_FRSURF */ |
98 | |
99 | RETURN |
100 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |