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

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

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

revision 1.1 by heimbach, Thu Mar 29 15:59:21 2012 UTC revision 1.2 by dgoldberg, Tue Sep 18 17:06:48 2012 UTC
# Line 74  C     LOCAL VARIABLES Line 74  C     LOCAL VARIABLES
74            Gj = (myYGlobalLo-1)+(bj-1)*sNy+j            Gj = (myYGlobalLo-1)+(bj-1)*sNy+j
75    
76            IF (STREAMICE_hmask(i,j,bi,bj).eq.1.0) THEN            IF (STREAMICE_hmask(i,j,bi,bj).eq.1.0) THEN
77             IF (Gi .eq. 1) THEN             IF (Gi.eq.1.AND..NOT.STREAMICE_EW_periodic) THEN
78              IF (STREAMICE_hmask(i+1,j,bi,bj).eq.1.0) THEN              IF (STREAMICE_hmask(i+1,j,bi,bj).eq.1.0) THEN
79               sx = (surf_el_streamice(i+1,j,bi,bj)-               sx = (surf_el_streamice(i+1,j,bi,bj)-
80       &             surf_el_streamice(i,j,bi,bj))/dxC(i+1,j,bi,bj)       &             surf_el_streamice(i,j,bi,bj))/dxC(i+1,j,bi,bj)
81              ELSE              ELSE
82               sx = 0. _d 0               sx = 0. _d 0
83              ENDIF              ENDIF
84             ELSEIF (Gi .eq. Nx) THEN             ELSEIF (Gi.eq.Nx.AND..NOT.STREAMICE_EW_periodic) THEN
85              IF (STREAMICE_hmask(i-1,j,bi,bj).eq.1.0) THEN              IF (STREAMICE_hmask(i-1,j,bi,bj).eq.1.0) THEN
86               sx = (surf_el_streamice(i,j,bi,bj)-               sx = (surf_el_streamice(i,j,bi,bj)-
87       &             surf_el_streamice(i-1,j,bi,bj))/dxC(i,j,bi,bj)       &             surf_el_streamice(i-1,j,bi,bj))/dxC(i,j,bi,bj)
# Line 110  C     LOCAL VARIABLES Line 110  C     LOCAL VARIABLES
110    
111                        
112    
113             IF (Gj .eq. 1) THEN             IF (Gj.eq.1.AND..NOT.STREAMICE_NS_periodic) THEN
114              IF (STREAMICE_hmask(i,j+1,bi,bj).eq.1.0) THEN              IF (STREAMICE_hmask(i,j+1,bi,bj).eq.1.0) THEN
115               sy = (surf_el_streamice(i,j+1,bi,bj)-               sy = (surf_el_streamice(i,j+1,bi,bj)-
116       &             surf_el_streamice(i,j,bi,bj))/dyC(i,j+1,bi,bj)       &             surf_el_streamice(i,j,bi,bj))/dyC(i,j+1,bi,bj)
117              ELSE              ELSE
118               sy = 0. _d 0               sy = 0. _d 0
119              ENDIF              ENDIF
120             ELSEIF (Gj .eq. Ny) THEN             ELSEIF (Gj.eq.Ny.AND..NOT.STREAMICE_NS_periodic) THEN
121              IF (STREAMICE_hmask(i,j-1,bi,bj).eq.1.0) THEN              IF (STREAMICE_hmask(i,j-1,bi,bj).eq.1.0) THEN
122               sy = (surf_el_streamice(i,j,bi,bj)-               sy = (surf_el_streamice(i,j,bi,bj)-
123       &             surf_el_streamice(i,j-1,bi,bj))/dyC(i,j,bi,bj)       &             surf_el_streamice(i,j-1,bi,bj))/dyC(i,j,bi,bj)
# Line 149  C     LOCAL VARIABLES Line 149  C     LOCAL VARIABLES
149              DO l=0,1              DO l=0,1
150               IF (STREAMICE_umask(i+k,j+l,bi,bj).eq.1.0) THEN               IF (STREAMICE_umask(i+k,j+l,bi,bj).eq.1.0) THEN
151                taudx_SI(i+k,j+l,bi,bj) = taudx_SI(i+k,j+l,bi,bj) -                taudx_SI(i+k,j+l,bi,bj) = taudx_SI(i+k,j+l,bi,bj) -
152       &         0.25 * streamice_density * gravity * sx *       &         0.25 * streamice_density * gravity *
153         &         (streamice_bg_surf_slope_x+sx) *
154       &         H_streamice(i,j,bi,bj) * rA(i,j,bi,bj)       &         H_streamice(i,j,bi,bj) * rA(i,j,bi,bj)
155    !     &          (streamice_bg_surf_slope_x) *
156    !     &         1000. * rA(i,j,bi,bj)
157                taudy_SI(i+k,j+l,bi,bj) = taudy_SI(i+k,j+l,bi,bj) -                taudy_SI(i+k,j+l,bi,bj) = taudy_SI(i+k,j+l,bi,bj) -
158       &         0.25 * streamice_density * gravity * sy *       &         0.25 * streamice_density * gravity *
159         &         (streamice_bg_surf_slope_y+sy) *
160       &         H_streamice(i,j,bi,bj) * rA(i,j,bi,bj)       &         H_streamice(i,j,bi,bj) * rA(i,j,bi,bj)
161                                
162               ENDIF               ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22