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.6 | C $Header: $ |
2 | C $Name: $ | ||
3 | heimbach | 1.1 | |
4 | edhill | 1.6 | #include "EXF_OPTIONS.h" |
5 | heimbach | 1.1 | |
6 | subroutine exf_filter_rl( | ||
7 | I arr, | ||
8 | I ckind, | ||
9 | I mythid | ||
10 | & ) | ||
11 | |||
12 | c ================================================================== | ||
13 | c SUBROUTINE exf_filter_rl | ||
14 | c ================================================================== | ||
15 | c | ||
16 | c o Read a flux record for external forcing. | ||
17 | c | ||
18 | c started: Ralf.Giering@FastOpt.de 24-Mai-2000 | ||
19 | dimitri | 1.4 | c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002 |
20 | heimbach | 1.1 | c |
21 | c ================================================================== | ||
22 | c SUBROUTINE exf_filter_rl | ||
23 | c ================================================================== | ||
24 | |||
25 | implicit none | ||
26 | |||
27 | c == global variables == | ||
28 | |||
29 | #include "EEPARAMS.h" | ||
30 | #include "SIZE.h" | ||
31 | #include "GRID.h" | ||
32 | dimitri | 1.3 | #include "PARAMS.h" |
33 | heimbach | 1.1 | #include "exf_constants.h" |
34 | heimbach | 1.2 | #include "exf.h" |
35 | #include "exf_param.h" | ||
36 | heimbach | 1.1 | |
37 | c == routine arguments == | ||
38 | |||
39 | _RL arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) | ||
40 | character*1 ckind | ||
41 | integer mythid | ||
42 | |||
43 | c == local variables == | ||
44 | |||
45 | integer bi,bj | ||
46 | integer i,j | ||
47 | integer itlo,ithi | ||
48 | integer jtlo,jthi | ||
49 | |||
50 | c == end of interface == | ||
51 | |||
52 | jtlo = mybylo(mythid) | ||
53 | jthi = mybyhi(mythid) | ||
54 | itlo = mybxlo(mythid) | ||
55 | ithi = mybxhi(mythid) | ||
56 | |||
57 | edhill | 1.5 | c Do not filter with pkg/seaice because of B/C-grid interpolation |
58 | dimitri | 1.3 | IF ( .NOT. useSEAICE ) THEN |
59 | |||
60 | heimbach | 1.1 | c filter forcing field array |
61 | do bj = jtlo,jthi | ||
62 | do bi = itlo,ithi | ||
63 | |||
64 | c Set undefined values to zero. | ||
65 | crg not necessary and | ||
66 | crg would require additional intermediate results in adjoint | ||
67 | crg do j = 1,sny | ||
68 | crg do i = 1,snx | ||
69 | crg if (arr(i,j,bi,bj) .le. exf_undef) then | ||
70 | crg arr(i,j,bi,bj) = 0. _d 0 | ||
71 | crg endif | ||
72 | crg enddo | ||
73 | crg enddo | ||
74 | |||
75 | c Set land points to zero | ||
76 | if (ckind .eq. 's') then | ||
77 | |||
78 | do j = 1,sny | ||
79 | do i = 1,snx | ||
80 | if ( _hFacC(i,j,1,bi,bj) .eq. 0. ) then | ||
81 | arr(i,j,bi,bj) = 0. _d 0 | ||
82 | endif | ||
83 | enddo | ||
84 | enddo | ||
85 | |||
86 | else if (ckind .eq. 'u') then | ||
87 | |||
88 | do j = 1,sny | ||
89 | do i = 1,snx | ||
90 | if ( maskw(i,j,1,bi,bj) .eq. 0. ) then | ||
91 | arr(i,j,bi,bj) = 0. _d 0 | ||
92 | endif | ||
93 | enddo | ||
94 | enddo | ||
95 | |||
96 | else if (ckind .eq. 'v') then | ||
97 | |||
98 | do j = 1,sny | ||
99 | do i = 1,snx | ||
100 | if ( masks(i,j,1,bi,bj) .eq. 0. ) then | ||
101 | arr(i,j,bi,bj) = 0. _d 0 | ||
102 | endif | ||
103 | enddo | ||
104 | enddo | ||
105 | |||
106 | end if | ||
107 | |||
108 | enddo | ||
109 | enddo | ||
110 | dimitri | 1.3 | |
111 | ENDIF | ||
112 | dimitri | 1.4 | c END IF ( .NOT. useSEAICE ) |
113 | heimbach | 1.1 | |
114 | end |
ViewVC Help | |
Powered by ViewVC 1.1.22 |