/[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.3 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC
# Line 73  C     DIRICHLET BOUNDARY VALUES ADDED TO Line 73  C     DIRICHLET BOUNDARY VALUES ADDED TO
73       O    ubd_SI,       O    ubd_SI,
74       O    vbd_SI)       O    vbd_SI)
75    
76          _EXCH_XY_RL (ubd_SI, myThid)
77          _EXCH_XY_RL (vbd_SI, myThid)
78    
79    !      CALL WRITE_FLD_XY_RL ( "ubd_SI", "",
80    !     &   ubd_SI, 0, myThid )
81    
82    !      CALL WRITE_FLD_XY_RL ( "vbd_SI", "",
83    !     &   STREAMICE_vmask, 0, myThid )
84    
85        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
86         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
87          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
# Line 100  C     GET DIAGONAL OF MATRIX Line 109  C     GET DIAGONAL OF MATRIX
109    
110  C     ccccc  C     ccccc
111    
112          DO bj = myByLo(myThid), myByHi(myThid)
113           DO bi = myBxLo(myThid), myBxHi(myThid)
114            DO j=1-OLy,sNy+OLy
115             DO i=1-OLy,sNx+OLy
116              IF (STREAMICE_umask(i,j,bi,bj).ne.1.0)
117         &     cg_Uin(i,j,bi,bj)=0.0
118              IF (STREAMICE_vmask(i,j,bi,bj).ne.1.0)
119         &     cg_Vin(i,j,bi,bj)=0.0
120    
121    !           print *, "rhs", i,j,RHSu_SI(i,j,bi,bj)
122    
123             ENDDO
124            ENDDO
125           ENDDO
126          ENDDO
127    
128    
129  #ifdef STREAMICE_CONSTRUCT_MATRIX  #ifdef STREAMICE_CONSTRUCT_MATRIX
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22