/[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.3 by dgoldberg, Thu Jul 26 16:13:18 2012 UTC revision 1.5 by dgoldberg, Fri May 10 18:38:07 2013 UTC
# Line 39  C     LOCAL VARIABLES Line 39  C     LOCAL VARIABLES
39                          ! there are valid cells to calculate a                          ! there are valid cells to calculate a
40                          ! slope-limited 2nd order flux                          ! slope-limited 2nd order flux
41        _RL SLOPE_LIMITER        _RL SLOPE_LIMITER
42        _RL total_vol_out  !       _RL total_vol_out
43        external SLOPE_LIMITER        external SLOPE_LIMITER
44    
45        total_vol_out = 0.0  !       total_vol_out = 0.0
46    
47        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
48         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# 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),
# Line 124  C        X DIRECTION AND 3 CELLS OUT IN Line 128  C        X DIRECTION AND 3 CELLS OUT IN
128    
129              ENDIF              ENDIF
130    
131              if (streamice_ufacemask(i,j,bi,bj).eq.2.0) THEN  !             if (streamice_ufacemask(i,j,bi,bj).eq.2.0) THEN
132               total_vol_out = total_vol_out +  !              total_vol_out = total_vol_out +
133       &         hflux_x (i,j,bi,bj) * time_step  !      &         hflux_x (i,j,bi,bj) * time_step
134              ENDIF  !             ENDIF
135    
136             ENDIF             ENDIF
137            ENDDO            ENDDO
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22