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

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

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

revision 1.2 by heimbach, Thu May 17 21:30:17 2012 UTC revision 1.3 by dgoldberg, Fri May 10 18:38:07 2013 UTC
# Line 81  C         Y DIRECTION Line 81  C         Y DIRECTION
81       &            (STREAMICE_hmask(i,j-2,bi,bj).eq.1.0))       &            (STREAMICE_hmask(i,j-2,bi,bj).eq.1.0))
82       &             H0_valid(i,j,bi,bj)=.true.       &             H0_valid(i,j,bi,bj)=.true.
83                            
84                IF ((Gj.eq.1).and.(STREAMICE_hmask(i,j-1,bi,bj).eq.3.0))                IF ((STREAMICE_hmask(i,j-1,bi,bj).eq.3.0))
85       &         THEN  ! we are at western bdry and there is a thick. bdry cond       &         THEN  ! we are at western bdry and there is a thick. bdry cond
86                 hflux_y (i,j,bi,bj) = h(i,j-1,bi,bj) * vface                 hflux_y (i,j,bi,bj) = h(i,j-1,bi,bj) * vface
87    !               PRINT *, "BOUNDARY FLUX UP", hflux_y (i,j,bi,bj),
88    !     &          h(i,j-1,bi,bj),vface
89                ELSEIF (H0_valid(i,j,bi,bj)) THEN                ELSEIF (H0_valid(i,j,bi,bj)) THEN
90                 phi = SLOPE_LIMITER (                 phi = SLOPE_LIMITER (
91       &          stencil(0)-stencil(-1),       &          stencil(0)-stencil(-1),
# Line 101  C         Y DIRECTION Line 103  C         Y DIRECTION
103       &            (STREAMICE_hmask(i,j+1,bi,bj).eq.1.0))       &            (STREAMICE_hmask(i,j+1,bi,bj).eq.1.0))
104       &             H0_valid(i,j,bi,bj)=.true.       &             H0_valid(i,j,bi,bj)=.true.
105    
106                IF ((Gj.eq.Ny).and.(STREAMICE_hmask(i,j+1,bi,bj).eq.3.0))                IF ((STREAMICE_hmask(i,j,bi,bj).eq.3.0))
107       &         THEN  ! we are at western bdry and there is a thick. bdry cond       &         THEN  ! we are at western bdry and there is a thick. bdry cond
108                 hflux_y (i,j,bi,bj) = h(i,j+1,bi,bj) * vface                 hflux_y (i,j,bi,bj) = h(i,j,bi,bj) * vface
109    !               PRINT *, "BOUNDARY FLUX DOWN", hflux_y (i,j,bi,bj),
110    !     &          h(i,j,bi,bj),vface
111                ELSEIF (H0_valid(i,j,bi,bj)) THEN                ELSEIF (H0_valid(i,j,bi,bj)) THEN
112                 phi = SLOPE_LIMITER (                 phi = SLOPE_LIMITER (
113       &          stencil(0)-stencil(-1),       &          stencil(0)-stencil(-1),
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22