13 |
I cg_tauV, |
I cg_tauV, |
14 |
I tolerance, |
I tolerance, |
15 |
O iters, |
O iters, |
16 |
|
I maxIter, |
17 |
I myThid ) |
I myThid ) |
18 |
|
|
19 |
C /============================================================\ |
C /============================================================\ |
38 |
|
|
39 |
INTEGER myThid |
INTEGER myThid |
40 |
INTEGER iters |
INTEGER iters |
41 |
|
INTEGER maxIter |
42 |
_RL tolerance |
_RL tolerance |
43 |
|
|
44 |
_RL cg_Uin (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL cg_Uin (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
75 |
O ubd_SI, |
O ubd_SI, |
76 |
O vbd_SI) |
O vbd_SI) |
77 |
|
|
78 |
_EXCH_XY_RL (ubd_SI, myThid) |
_EXCH_XY_RL(ubd_SI, myThid) |
79 |
_EXCH_XY_RL (vbd_SI, myThid) |
_EXCH_XY_RL(vbd_SI, myThid) |
80 |
|
|
81 |
! CALL WRITE_FLD_XY_RL ( "ubd_SI", "", |
! CALL WRITE_FLD_XY_RL ( "ubd_SI", "", |
82 |
! & ubd_SI, 0, myThid ) |
! & ubd_SI, 0, myThid ) |
133 |
|
|
134 |
CALL STREAMICE_CG_MAKE_A(myThid) |
CALL STREAMICE_CG_MAKE_A(myThid) |
135 |
|
|
136 |
|
! call write_fld_xy_rl ("streamicb_cg_A1_m1_m1","", |
137 |
|
! & streamice_cg_A1(:,:,1,1,-1,-1),0,myThid) |
138 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_m1_0","", |
139 |
|
c & streamice_cg_A1(:,:,1,1,-1,0),0,myThid) |
140 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_m1_p1","", |
141 |
|
c & streamice_cg_A1(:,:,1,1,-1,1),0,myThid) |
142 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_0_m1","", |
143 |
|
c & streamice_cg_A1(:,:,1,1,0,-1),0,myThid) |
144 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_0_0","", |
145 |
|
c & streamice_cg_A1(:,:,1,1,0,0),0,myThid) |
146 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_0_p1","", |
147 |
|
c & streamice_cg_A1(:,:,1,1,0,1),0,myThid) |
148 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_p1_m1","", |
149 |
|
c & streamice_cg_A1(:,:,1,1,1,-1),0,myThid) |
150 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_p1_0","", |
151 |
|
c & streamice_cg_A1(:,:,1,1,1,0),0,myThid) |
152 |
|
c call write_fld_xy_rl ("streamicb_cg_A1_p1_p1","", |
153 |
|
c & streamice_cg_A1(:,:,1,1,1,1),0,myThid) |
154 |
|
c |
155 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_m1_m1","", |
156 |
|
c & streamice_cg_A2(:,:,1,1,-1,-1),0,myThid) |
157 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_m1_0","", |
158 |
|
c & streamice_cg_A2(:,:,1,1,-1,0),0,myThid) |
159 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_m1_p1","", |
160 |
|
c & streamice_cg_A2(:,:,1,1,-1,1),0,myThid) |
161 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_0_m1","", |
162 |
|
c & streamice_cg_A2(:,:,1,1,0,-1),0,myThid) |
163 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_0_0","", |
164 |
|
c & streamice_cg_A2(:,:,1,1,0,0),0,myThid) |
165 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_0_p1","", |
166 |
|
c & streamice_cg_A2(:,:,1,1,0,1),0,myThid) |
167 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_p1_m1","", |
168 |
|
c & streamice_cg_A2(:,:,1,1,1,-1),0,myThid) |
169 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_p1_0","", |
170 |
|
c & streamice_cg_A2(:,:,1,1,1,0),0,myThid) |
171 |
|
c call write_fld_xy_rl ("streamicb_cg_A2_p1_p1","", |
172 |
|
c & streamice_cg_A2(:,:,1,1,1,1),0,myThid) |
173 |
|
c |
174 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_m1_m1","", |
175 |
|
c & streamice_cg_A3(:,:,1,1,-1,-1),0,myThid) |
176 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_m1_0","", |
177 |
|
c & streamice_cg_A3(:,:,1,1,-1,0),0,myThid) |
178 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_m1_p1","", |
179 |
|
c & streamice_cg_A3(:,:,1,1,-1,1),0,myThid) |
180 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_0_m1","", |
181 |
|
c & streamice_cg_A3(:,:,1,1,0,-1),0,myThid) |
182 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_0_0","", |
183 |
|
c & streamice_cg_A3(:,:,1,1,0,0),0,myThid) |
184 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_0_p1","", |
185 |
|
c & streamice_cg_A3(:,:,1,1,0,1),0,myThid) |
186 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_p1_m1","", |
187 |
|
c & streamice_cg_A3(:,:,1,1,1,-1),0,myThid) |
188 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_p1_0","", |
189 |
|
c & streamice_cg_A3(:,:,1,1,1,0),0,myThid) |
190 |
|
c call write_fld_xy_rl ("streamicb_cg_A3_p1_p1","", |
191 |
|
c & streamice_cg_A3(:,:,1,1,1,1),0,myThid) |
192 |
|
c |
193 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_m1_m1","", |
194 |
|
c & streamice_cg_A4(:,:,1,1,-1,-1),0,myThid) |
195 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_m1_0","", |
196 |
|
c & streamice_cg_A4(:,:,1,1,-1,0),0,myThid) |
197 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_m1_p1","", |
198 |
|
c & streamice_cg_A4(:,:,1,1,-1,1),0,myThid) |
199 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_0_m1","", |
200 |
|
c & streamice_cg_A4(:,:,1,1,0,-1),0,myThid) |
201 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_0_0","", |
202 |
|
c & streamice_cg_A4(:,:,1,1,0,0),0,myThid) |
203 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_0_p1","", |
204 |
|
c & streamice_cg_A4(:,:,1,1,0,1),0,myThid) |
205 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_p1_m1","", |
206 |
|
c & streamice_cg_A4(:,:,1,1,1,-1),0,myThid) |
207 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_p1_0","", |
208 |
|
c & streamice_cg_A4(:,:,1,1,1,0),0,myThid) |
209 |
|
c call write_fld_xy_rl ("streamicb_cg_A4_p1_p1","", |
210 |
|
c & streamice_cg_A4(:,:,1,1,1,1),0,myThid) |
211 |
|
c |
212 |
|
|
213 |
! print *, "MATRIX 1" |
! print *, "MATRIX 1" |
214 |
! do i=1,sNx |
! do i=1,sNx |
215 |
! print *, i, |
! print *, i, |
237 |
& streamice_cg_A4, |
& streamice_cg_A4, |
238 |
& tolerance, |
& tolerance, |
239 |
& iters, |
& iters, |
240 |
|
& maxIter, |
241 |
& myThid ) |
& myThid ) |
242 |
|
|
243 |
_EXCH_XY_RL( RHSu_SI, myThid ) |
_EXCH_XY_RL( RHSu_SI, myThid ) |