/[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.2 by dgoldberg, Wed Jan 9 21:56:18 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    
77    !      CALL WRITE_FLD_XY_RL ( "ubd_SI", "",
78    !     &   ubd_SI, 0, myThid )
79    
80    !      CALL WRITE_FLD_XY_RL ( "vbd_SI", "",
81    !     &   STREAMICE_vmask, 0, myThid )
82    
83        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
84         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
85          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
# Line 100  C     GET DIAGONAL OF MATRIX Line 107  C     GET DIAGONAL OF MATRIX
107    
108  C     ccccc  C     ccccc
109    
110          DO bj = myByLo(myThid), myByHi(myThid)
111           DO bi = myBxLo(myThid), myBxHi(myThid)
112            DO j=1-OLy,sNy+OLy
113             DO i=1-OLy,sNx+OLy
114              IF (STREAMICE_umask(i,j,bi,bj).ne.1.0)
115         &     cg_Uin(i,j,bi,bj)=0.0
116              IF (STREAMICE_vmask(i,j,bi,bj).ne.1.0)
117         &     cg_Vin(i,j,bi,bj)=0.0
118    
119    !           print *, "rhs", i,j,RHSu_SI(i,j,bi,bj)
120    
121             ENDDO
122            ENDDO
123           ENDDO
124          ENDDO
125    
126    
127  #ifdef STREAMICE_CONSTRUCT_MATRIX  #ifdef STREAMICE_CONSTRUCT_MATRIX
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22