/[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.6 by dgoldberg, Wed Aug 27 19:29:13 2014 UTC
# Line 39  C     LOCAL VARIABLES Line 39  C     LOCAL VARIABLES
39       &        Gi, Gj       &        Gi, Gj
40        LOGICAL at_west_bdry, at_east_bdry,        LOGICAL at_west_bdry, at_east_bdry,
41       &        at_north_bdry, at_south_bdry       &        at_north_bdry, at_south_bdry
42        _RL sx, sy, diffx, diffy, neu_val        _RL sx, sy, diffx, diffy, geom_fac
43                
44        IF (myXGlobalLo.eq.1) at_west_bdry = .true.        IF (myXGlobalLo.eq.1) at_west_bdry = .true.
45        IF (myYGlobalLo.eq.1) at_south_bdry = .true.        IF (myYGlobalLo.eq.1) at_south_bdry = .true.
# Line 50  C     LOCAL VARIABLES Line 50  C     LOCAL VARIABLES
50    
51        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
52         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
53          DO j=1-OLy,sNy+OLy          DO j=1-OLy+1,sNy+OLy-1
54           DO i=1-OLx,sNx+OLx           DO i=1-OLx+1,sNx+OLx-1
55            taudx_SI(i,j,bi,bj) = 0. _d 0            taudx_SI(i,j,bi,bj) = 0. _d 0
56            taudy_SI(i,j,bi,bj) = 0. _d 0            taudy_SI(i,j,bi,bj) = 0. _d 0
57           ENDDO           ENDDO
# Line 62  C     LOCAL VARIABLES Line 62  C     LOCAL VARIABLES
62        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
63         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
64                    
65          DO i=0,sNx+1          DO i=1,sNx
66           DO j=0,sNy+1           DO j=1,sNy
67    
68            diffx = 0. _d 0            diffx = 0. _d 0
69            diffy = 0. _d 0            diffy = 0. _d 0
# Line 73  C     LOCAL VARIABLES Line 73  C     LOCAL VARIABLES
73            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i
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_umask(i,j,bi,bj).eq.1.0) THEN
77             IF (Gi .eq. 1) 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.AND.
79               sx = (surf_el_streamice(i+1,j,bi,bj)-       &      streamice_hmask(i,j,bi,bj).eq.1.0) THEN
80       &             surf_el_streamice(i,j,bi,bj))/dxC(i+1,j,bi,bj)          
81              ELSE  !             geom_fac = sqrt(rA(i-1,j,bi,bj)*recip_rA(i,j,bi,bj)*
82               sx = 0. _d 0  !      &        dxF(i,j,bi,bj)*recip_dxF(i-1,j,bi,bj))
83              ENDIF  
84             ELSEIF (Gi .eq. Nx) THEN              taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
85              IF (STREAMICE_hmask(i-1,j,bi,bj).eq.1.0) THEN       &       0.25 * dyG(i,j,bi,bj) *
86               sx = (surf_el_streamice(i,j,bi,bj)-       &       streamice_density * gravity *
87       &             surf_el_streamice(i-1,j,bi,bj))/dxC(i,j,bi,bj)       &       (H_streamice(i,j,bi,bj)+
88              ELSE       &        H_streamice(i-1,j,bi,bj)) *
89               sx = 0. _d 0       &       (surf_el_streamice(i,j,bi,bj)-
90              ENDIF       &              surf_el_streamice(i-1,j,bi,bj))
91             ELSE  
92              IF (STREAMICE_hmask(i+1,j,bi,bj).eq.1.0) THEN              taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
93               diffx = diffx + dxC(i+1,j,bi,bj)       &       streamice_density * gravity *
94               sx = surf_el_streamice(i+1,j,bi,bj)       &       streamice_bg_surf_slope_x * .25 * rA(i,j,bi,bj) *
95         &       (H_streamice(i-1,j,bi,bj) + H_streamice(i,j,bi,bj))
96    
97               ELSE IF (streamice_hmask(i-1,j,bi,bj).eq.1.0) THEN
98    
99                IF (float_frac_streamice(i-1,j,bi,bj) .eq. 1.0) THEN
100    
101                 taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) +
102         &        0.25 * dyG(i,j,bi,bj) *
103         &        gravity *
104         &        (streamice_density * H_streamice(i-1,j,bi,bj)**2 -
105    #ifdef USE_ALT_RLOW
106         &         streamice_density_ocean_avg * R_low_si(i-1,j,bi,bj)**2)
107    #else          
108         &         streamice_density_ocean_avg * R_low(i-1,j,bi,bj)**2)
109    #endif
110    
111              ELSE              ELSE
112               sx = surf_el_streamice(i,j,bi,bj)  
113                 taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) +
114         &        0.25 * dyG(i,j,bi,bj) *
115         &        streamice_density * gravity *
116         &        (1-streamice_density/streamice_density_ocean_avg) *
117         &         H_streamice(i-1,j,bi,bj)**2
118    
119              ENDIF              ENDIF
120              IF (STREAMICE_hmask(i-1,j,bi,bj).eq.1.0) THEN  
121               diffx = diffx + dxC(i,j,bi,bj)             ELSE IF (streamice_hmask(i,j,bi,bj).eq.1.0) THEN
122               sx = sx - surf_el_streamice(i-1,j,bi,bj)  
123                IF (float_frac_streamice(i,j,bi,bj) .eq. 1.0) THEN
124    
125                 taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
126         &        0.25 * dyG(i,j,bi,bj) *
127         &        gravity *
128         &        (streamice_density * H_streamice(i,j,bi,bj)**2 -
129    #ifdef USE_ALT_RLOW
130         &         streamice_density_ocean_avg * R_low_si(i,j,bi,bj)**2)
131    #else          
132         &         streamice_density_ocean_avg * R_low(i,j,bi,bj)**2)
133    #endif
134    
135              ELSE              ELSE
136               sx = sx - surf_el_streamice(i,j,bi,bj)  
137              ENDIF               taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
138              IF (diffx .gt. 0. _d 0) THEN       &        0.25 * dyG(i,j,bi,bj) *
139               sx = sx / diffx       &        streamice_density * gravity *
140         &        (1-streamice_density/streamice_density_ocean_avg) *
141         &         H_streamice(i,j,bi,bj)**2
142    
143                END IF
144               END IF
145    
146    ! cells below
147    
148               IF (streamice_hmask(i-1,j-1,bi,bj).eq.1.0.AND.
149         &      streamice_hmask(i,j-1,bi,bj).eq.1.0) THEN
150            
151    !             geom_fac = sqrt(rA(i-1,j-1,bi,bj)*recip_rA(i,j-1,bi,bj)*
152    !      &             dxF(i,j-1,bi,bj)*recip_dxF(i-1,j-1,bi,bj))
153    
154                taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
155         &       0.25 * dyg(i,j-1,bi,bj) *
156         &       streamice_density * gravity *
157         &       (H_streamice(i,j-1,bi,bj)+
158         &        H_streamice(i-1,j-1,bi,bj)) *
159         &       (surf_el_streamice(i,j-1,bi,bj)-
160         &              surf_el_streamice(i-1,j-1,bi,bj))
161    
162                taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
163         &       streamice_density * gravity *
164         &       streamice_bg_surf_slope_x * .25 * rA(i,j-1,bi,bj) *
165         &       (H_streamice(i-1,j-1,bi,bj) + H_streamice(i,j-1,bi,bj))
166    
167               ELSE IF (streamice_hmask(i-1,j-1,bi,bj).eq.1.0) THEN
168    
169                IF (float_frac_streamice(i-1,j-1,bi,bj) .eq. 1.0) THEN
170    
171                 taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) +
172         &        0.25 * dyg(i,j-1,bi,bj) *
173         &        gravity *
174         &        (streamice_density * H_streamice(i-1,j-1,bi,bj)**2 -
175    #ifdef USE_ALT_RLOW
176         &         streamice_density_ocean_avg*R_low_si(i-1,j-1,bi,bj)**2)
177    #else          
178         &         streamice_density_ocean_avg * R_low(i-1,j-1,bi,bj)**2)
179    #endif
180    
181              ELSE              ELSE
182               sx = 0. _d 0  
183                 taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) +
184         &        0.25 * dyg(i,j-1,bi,bj) *
185         &        streamice_density * gravity *
186         &        (1-streamice_density/streamice_density_ocean_avg) *
187         &         H_streamice(i-1,j-1,bi,bj)**2
188    
189              ENDIF              ENDIF
            ENDIF  
190    
191                         ELSE IF (streamice_hmask(i,j-1,bi,bj).eq.1.0) THEN
192    
193                IF (float_frac_streamice(i,j-1,bi,bj) .eq. 1.0) THEN
194    
195                 taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
196         &        0.25 * dyg(i,j-1,bi,bj) *
197         &        gravity *
198         &        (streamice_density * H_streamice(i,j-1,bi,bj)**2 -
199    #ifdef USE_ALT_RLOW
200         &         streamice_density_ocean_avg * R_low_si(i,j-1,bi,bj)**2)
201    #else          
202         &         streamice_density_ocean_avg * R_low(i,j-1,bi,bj)**2)
203    #endif
204    
            IF (Gj .eq. 1) THEN  
             IF (STREAMICE_hmask(i,j+1,bi,bj).eq.1.0) THEN  
              sy = (surf_el_streamice(i,j+1,bi,bj)-  
      &             surf_el_streamice(i,j,bi,bj))/dyC(i,j+1,bi,bj)  
205              ELSE              ELSE
206               sy = 0. _d 0  
207              ENDIF               taudx_si(i,j,bi,bj) = taudx_si(i,j,bi,bj) -
208             ELSEIF (Gj .eq. Ny) THEN       &        0.25 * dyg(i,j-1,bi,bj) *
209              IF (STREAMICE_hmask(i,j-1,bi,bj).eq.1.0) THEN       &        streamice_density * gravity *
210               sy = (surf_el_streamice(i,j,bi,bj)-       &        (1-streamice_density/streamice_density_ocean_avg) *
211       &             surf_el_streamice(i,j-1,bi,bj))/dyC(i,j,bi,bj)       &         H_streamice(i,j-1,bi,bj)**2
212    
213                END IF
214               END IF
215              END IF       ! if umask==1
216    
217    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
218    
219              IF (streamice_vmask(i,j,bi,bj).eq.1.0) THEN
220    
221               IF (streamice_hmask(i,j-1,bi,bj).eq.1.0.AND.
222         &      streamice_hmask(i,j,bi,bj).eq.1.0) THEN
223            
224    !             geom_fac = sqrt(rA(i,j-1,bi,bj)*recip_rA(i,j,bi,bj)*
225    !      &        dxF(i,j,bi,bj)*recip_dyF(i,j-1,bi,bj))
226    
227                taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
228         &       0.25 * dxG(i,j,bi,bj) *
229         &       streamice_density * gravity *
230         &       (H_streamice(i,j,bi,bj)+
231         &        H_streamice(i,j-1,bi,bj)) *
232         &       (surf_el_streamice(i,j,bi,bj)-
233         &              surf_el_streamice(i,j-1,bi,bj))
234    
235                taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
236         &       streamice_density * gravity *
237         &       streamice_bg_surf_slope_y * .25 * rA(i,j,bi,bj) *
238         &       (H_streamice(i,j-1,bi,bj) + H_streamice(i,j,bi,bj))
239    
240               ELSE IF (streamice_hmask(i,j-1,bi,bj).eq.1.0) THEN
241    
242                IF (float_frac_streamice(i,j-1,bi,bj) .eq. 1.0) THEN
243    
244                 taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) +
245         &        0.25 * dxG(i,j,bi,bj) *
246         &        gravity *
247         &        (streamice_density * H_streamice(i,j-1,bi,bj)**2 -
248    #ifdef USE_ALT_RLOW
249         &         streamice_density_ocean_avg * R_low_si(i,j-1,bi,bj)**2)
250    #else          
251         &         streamice_density_ocean_avg * R_low(i,j-1,bi,bj)**2)
252    #endif
253    
254              ELSE              ELSE
255               sy = 0. _d 0  
256                 taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) +
257         &        0.25 * dxG(i,j,bi,bj) *
258         &        streamice_density * gravity *
259         &        (1-streamice_density/streamice_density_ocean_avg) *
260         &         H_streamice(i,j-1,bi,bj)**2
261    
262              ENDIF              ENDIF
263             ELSE  
264              IF (STREAMICE_hmask(i,j+1,bi,bj).eq.1.0) THEN             ELSE IF (streamice_hmask(i,j,bi,bj).eq.1.0) THEN
265                
266               diffy = diffy + dyC(i,j+1,bi,bj)              IF (float_frac_streamice(i,j,bi,bj) .eq. 1.0) THEN
267               sy = surf_el_streamice(i,j+1,bi,bj)  
268                 taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
269         &        0.25 * dxG(i,j,bi,bj) *
270         &        gravity *
271         &        (streamice_density * H_streamice(i,j,bi,bj)**2 -
272    #ifdef USE_ALT_RLOW
273         &         streamice_density_ocean_avg * R_low_si(i,j,bi,bj)**2)
274    #else          
275         &         streamice_density_ocean_avg * R_low(i,j,bi,bj)**2)
276    #endif
277    
278              ELSE              ELSE
279               sy = surf_el_streamice(i,j,bi,bj)  
280              ENDIF               taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
281              IF (STREAMICE_hmask(i,j-1,bi,bj).eq.1.0) THEN       &        0.25 * dxG(i,j,bi,bj) *
282               diffy = diffy + dyC(i,j,bi,bj)       &        streamice_density * gravity *
283               sy = sy - surf_el_streamice(i,j-1,bi,bj)       &        (1-streamice_density/streamice_density_ocean_avg) *
284         &         H_streamice(i,j,bi,bj)**2
285    
286                END IF
287               END IF
288    
289    ! cells to left
290    
291               IF (streamice_hmask(i-1,j-1,bi,bj).eq.1.0.AND.
292         &      streamice_hmask(i-1,j,bi,bj).eq.1.0) THEN
293            
294    !             geom_fac = sqrt(rA(i-1,j-1,bi,bj)*recip_rA(i-1,j,bi,bj)*
295    !      &        dxF(i-1,j,bi,bj)*recip_dxF(i-1,j-1,bi,bj))
296    
297                taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
298         &       0.25 * dxG(i-1,j,bi,bj) *
299         &       streamice_density * gravity *
300         &       (H_streamice(i-1,j,bi,bj)+
301         &        H_streamice(i-1,j-1,bi,bj)) *
302         &       (surf_el_streamice(i-1,j,bi,bj)-
303         &              surf_el_streamice(i-1,j-1,bi,bj))
304    
305                taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
306         &       streamice_density * gravity *
307         &       streamice_bg_surf_slope_y * .25 * rA(i-1,j,bi,bj) *
308         &       (H_streamice(i-1,j-1,bi,bj) + H_streamice(i-1,j,bi,bj))
309    
310    
311               ELSE IF (streamice_hmask(i-1,j-1,bi,bj).eq.1.0) THEN
312    
313                IF (float_frac_streamice(i-1,j-1,bi,bj) .eq. 1.0) THEN
314    
315                 taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) +
316         &        0.25 * dxG(i-1,j,bi,bj) *
317         &        gravity *
318         &        (streamice_density * H_streamice(i-1,j-1,bi,bj)**2 -
319    #ifdef USE_ALT_RLOW
320         &         streamice_density_ocean_avg*R_low_si(i-1,j-1,bi,bj)**2)
321    #else          
322         &         streamice_density_ocean_avg * R_low(i-1,j-1,bi,bj)**2)
323    #endif
324    
325              ELSE              ELSE
326               sy = sy - surf_el_streamice(i,j,bi,bj)  
327                 taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) +
328         &        0.25 * dxG(i-1,j,bi,bj) *
329         &        streamice_density * gravity *
330         &        (1-streamice_density/streamice_density_ocean_avg) *
331         &         H_streamice(i-1,j-1,bi,bj)**2
332    
333              ENDIF              ENDIF
334              IF (diffy .gt. 0. _d 0) THEN  
335               sy = sy / diffy             ELSE IF (streamice_hmask(i-1,j,bi,bj).eq.1.0) THEN
336    
337                IF (float_frac_streamice(i-1,j,bi,bj) .eq. 1.0) THEN
338    
339                 taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
340         &        0.25 * dxG(i-1,j,bi,bj) *
341         &        gravity *
342         &        (streamice_density * H_streamice(i-1,j,bi,bj)**2 -
343    #ifdef USE_ALT_RLOW
344         &         streamice_density_ocean_avg * R_low_si(i-1,j,bi,bj)**2)
345    #else          
346         &         streamice_density_ocean_avg * R_low(i-1,j,bi,bj)**2)
347    #endif
348    
349              ELSE              ELSE
              sy = 0. _d 0  
             ENDIF  
            ENDIF  
350    
351             DO k=0,1               taudy_si(i,j,bi,bj) = taudy_si(i,j,bi,bj) -
352              DO l=0,1       &        0.25 * dxG(i-1,j,bi,bj) *
353               IF (STREAMICE_umask(i+k,j+l,bi,bj).eq.1.0) THEN       &        streamice_density * gravity *
354                taudx_SI(i+k,j+l,bi,bj) = taudx_SI(i+k,j+l,bi,bj) -       &        (1-streamice_density/streamice_density_ocean_avg) *
355       &         0.25 * streamice_density * gravity * sx *       &         H_streamice(i-1,j,bi,bj)**2
356       &         H_streamice(i,j,bi,bj) * rA(i,j,bi,bj)  
357                taudy_SI(i+k,j+l,bi,bj) = taudy_SI(i+k,j+l,bi,bj) -              END IF
358       &         0.25 * streamice_density * gravity * sy *             END IF
359       &         H_streamice(i,j,bi,bj) * rA(i,j,bi,bj)            END IF      ! if vmask ==1
                 
              ENDIF  
             ENDDO  
            ENDDO  
   
            IF (float_frac_streamice(i,j,bi,bj) .eq. 1.0) then  
             neu_val = .5 * gravity *  
      &       (streamice_density * H_streamice (i,j,bi,bj) ** 2 -  
      &        streamice_density_ocean_avg * R_low(i,j,bi,bj) ** 2)  
            ELSE  
             neu_val = .5 * gravity *  
      &       (1-streamice_density/streamice_density_ocean_avg) *  
      &        streamice_density * H_streamice(i,j,bi,bj) ** 2  
            ENDIF  
   
            IF ((STREAMICE_ufacemask(i,j,bi,bj) .eq. 2)  
      &      .OR. (STREAMICE_hmask(i-1,j,bi,bj) .eq. 0)  
      &      .OR. (STREAMICE_hmask(i-1,j,bi,bj) .eq. 2) ) THEN ! left face of the cell is at a stress boundary  
           ! the depth-integrated longitudinal stress is equal to the difference of depth-integrated pressure on either side of the face  
           ! on the ice side, it is rho g h^2 / 2  
           ! on the ocean side, it is rhow g (delta OD)^2 / 2  
           ! OD can be zero under the ice; but it is ASSUMED on the ice-free side of the face, topography elevation is not above the base of the  
           !     ice in the current cell  
               
              taudx_SI(i,j,bi,bj) = taudx_SI(i,j,bi,bj) -  
      &        .5 * dyG(i,j,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector  
              taudx_SI(i,j+1,bi,bj) = taudx_SI(i,j+1,bi,bj) -  
      &        .5 * dyG(i,j,bi,bj) * neu_val  
            ENDIF  
   
            IF ((STREAMICE_ufacemask(i+1,j,bi,bj) .eq. 2)  
      &      .OR. (STREAMICE_hmask(i+1,j,bi,bj) .eq. 0)  
      &      .OR. (STREAMICE_hmask(i+1,j,bi,bj) .eq. 2) ) THEN  
               
              taudx_SI(i+1,j,bi,bj) = taudx_SI(i+1,j,bi,bj) +  
      &        .5 * dyG(i+1,j,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector  
              taudx_SI(i+1,j+1,bi,bj) = taudx_SI(i+1,j+1,bi,bj) +  
      &        .5 * dyG(i+1,j,bi,bj) * neu_val  
            ENDIF  
   
            IF ((STREAMICE_vfacemask(i,j,bi,bj) .eq. 2)  
      &      .OR. (STREAMICE_hmask(i,j-1,bi,bj) .eq. 0)  
      &      .OR. (STREAMICE_hmask(i,j-1,bi,bj) .eq. 2) ) THEN  
               
              taudy_SI(i,j,bi,bj) = taudy_SI(i,j,bi,bj) -  
      &        .5 * dxG(i,j,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector  
              taudy_SI(i+1,j,bi,bj) = taudy_SI(i+1,j,bi,bj) -  
      &        .5 * dxG(i,j,bi,bj) * neu_val  
            ENDIF  
   
            IF ((STREAMICE_vfacemask(i,j+1,bi,bj) .eq. 2)  
      &      .OR. (STREAMICE_hmask(i,j+1,bi,bj) .eq. 0)  
      &      .OR. (STREAMICE_hmask(i,j+1,bi,bj) .eq. 2) ) THEN  
               
              taudy_SI(i,j+1,bi,bj) = taudy_SI(i,j+1,bi,bj) +  
      &        .5 * dxG(i,j+1,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector  
              taudy_SI(i+1,j+1,bi,bj) = taudy_SI(i+1,j+1,bi,bj) +  
      &        .5 * dxG(i,j+1,bi,bj) * neu_val  
            ENDIF  
360    
           ENDIF  
361           ENDDO           ENDDO
362          ENDDO          ENDDO
363         ENDDO         ENDDO
364        ENDDO        ENDDO
365    
366    !      taudx_SI (1,1,1,1) = taudx_SI (1,1,1,1) +
367    !     & streamice_v_normal_pert (1,1,1,1)
368    
369    
370    
371  #endif  #endif
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22