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), |
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), |