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

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

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

revision 1.4 by dgoldberg, Tue Sep 18 17:06:48 2012 UTC revision 1.5 by dgoldberg, Fri May 10 18:38:07 2013 UTC
# Line 85  C        X DIRECTION AND 3 CELLS OUT IN Line 85  C        X DIRECTION AND 3 CELLS OUT IN
85       &            (STREAMICE_hmask(i-2,j,bi,bj).eq.1.0))       &            (STREAMICE_hmask(i-2,j,bi,bj).eq.1.0))
86       &             H0_valid(i,j,bi,bj)=.true.       &             H0_valid(i,j,bi,bj)=.true.
87                            
88                IF ((Gi.eq.1).and.(STREAMICE_hmask(i-1,j,bi,bj).eq.3.0))                IF ((STREAMICE_hmask(i-1,j,bi,bj).eq.3.0))
89       &         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
90                 hflux_x (i,j,bi,bj) = h(i-1,j,bi,bj) * uface                 hflux_x (i,j,bi,bj) = h(i-1,j,bi,bj) * uface
91    !               PRINT *, "BOUNDARY FLUX RIGHT", hflux_x (i,j,bi,bj),
92    !     &          h(i-1,j,bi,bj),uface,i,j
93                ELSEIF (H0_valid(i,j,bi,bj)) THEN                ELSEIF (H0_valid(i,j,bi,bj)) THEN
94                 phi = SLOPE_LIMITER (                 phi = SLOPE_LIMITER (
95       &          stencil(0)-stencil(-1),       &          stencil(0)-stencil(-1),
# Line 107  C        X DIRECTION AND 3 CELLS OUT IN Line 109  C        X DIRECTION AND 3 CELLS OUT IN
109       &            (STREAMICE_hmask(i+1,j,bi,bj).eq.1.0))       &            (STREAMICE_hmask(i+1,j,bi,bj).eq.1.0))
110       &             H0_valid(i,j,bi,bj)=.true.       &             H0_valid(i,j,bi,bj)=.true.
111    
112                IF ((Gi.eq.Nx).and.(STREAMICE_hmask(i+1,j,bi,bj).eq.3.0))                IF ((STREAMICE_hmask(i,j,bi,bj).eq.3.0))
113       &         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
114                 hflux_x (i,j,bi,bj) = h(i+1,j,bi,bj) * uface                 hflux_x (i,j,bi,bj) = h(i,j,bi,bj) * uface
115    !               PRINT *, "BOUNDARY FLUX LEFT", hflux_x(i,j,bi,bj),
116    !     &          h(i,j,bi,bj),uface
117                ELSEIF (H0_valid(i,j,bi,bj)) THEN                ELSEIF (H0_valid(i,j,bi,bj)) THEN
118                 phi = SLOPE_LIMITER (                 phi = SLOPE_LIMITER (
119       &          stencil(0)-stencil(-1),       &          stencil(0)-stencil(-1),
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22