/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_cg_wrapper.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_cg_wrapper.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by dgoldberg, Thu Jul 19 18:55:26 2012 UTC revision 1.4 by dgoldberg, Wed Aug 27 19:29:13 2014 UTC
# Line 13  CBOP Line 13  CBOP
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     /============================================================\
# Line 37  C     cg_Bu, cg_Bv - driving stress Line 38  C     cg_Bu, cg_Bv - driving stress
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)
# Line 73  C     DIRICHLET BOUNDARY VALUES ADDED TO Line 75  C     DIRICHLET BOUNDARY VALUES ADDED TO
75       O    ubd_SI,       O    ubd_SI,
76       O    vbd_SI)       O    vbd_SI)
77    
78          _EXCH_XY_RL(ubd_SI, myThid)
79          _EXCH_XY_RL(vbd_SI, myThid)
80    
81    !      CALL WRITE_FLD_XY_RL ( "ubd_SI", "",
82    !     &   ubd_SI, 0, myThid )
83    
84    !      CALL WRITE_FLD_XY_RL ( "vbd_SI", "",
85    !     &   STREAMICE_vmask, 0, myThid )
86    
87        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
88         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
89          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
# Line 100  C     GET DIAGONAL OF MATRIX Line 111  C     GET DIAGONAL OF MATRIX
111    
112  C     ccccc  C     ccccc
113    
114          DO bj = myByLo(myThid), myByHi(myThid)
115           DO bi = myBxLo(myThid), myBxHi(myThid)
116            DO j=1-OLy,sNy+OLy
117             DO i=1-OLy,sNx+OLy
118              IF (STREAMICE_umask(i,j,bi,bj).ne.1.0)
119         &     cg_Uin(i,j,bi,bj)=0.0
120              IF (STREAMICE_vmask(i,j,bi,bj).ne.1.0)
121         &     cg_Vin(i,j,bi,bj)=0.0
122    
123    !           print *, "rhs", i,j,RHSu_SI(i,j,bi,bj)
124    
125             ENDDO
126            ENDDO
127           ENDDO
128          ENDDO
129    
130    
131  #ifdef STREAMICE_CONSTRUCT_MATRIX  #ifdef STREAMICE_CONSTRUCT_MATRIX
# Line 107  C     ccccc Line 133  C     ccccc
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,
# Line 134  C     ccccc Line 237  C     ccccc
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 )
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22