/[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.3 by dgoldberg, Thu Sep 27 20:29:00 2012 UTC revision 1.5 by dgoldberg, Mon Dec 10 02:34:45 2012 UTC
# Line 204  C     LOCAL VARIABLES Line 204  C     LOCAL VARIABLES
204  #else  #else
205              neu_val = .5 * gravity *              neu_val = .5 * gravity *
206       &       (streamice_density * H_streamice (i,j,bi,bj) ** 2 -       &       (streamice_density * H_streamice (i,j,bi,bj) ** 2 -
207       &        streamice_density_ocean_avg * R_low(i,j,bi,bj) ** 2       &        streamice_density_ocean_avg * R_low(i,j,bi,bj) ** 2)
208  #endif  #endif
209             ELSE             ELSE
210              neu_val = .5 * gravity *              neu_val = .5 * gravity *
# Line 222  C     LOCAL VARIABLES Line 222  C     LOCAL VARIABLES
222            !     ice in the current cell            !     ice in the current cell
223                            
224               taudx_SI(i,j,bi,bj) = taudx_SI(i,j,bi,bj) -               taudx_SI(i,j,bi,bj) = taudx_SI(i,j,bi,bj) -
225       &        .5 * dyG(i,j,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector       &        .5 * dyG(i,j,bi,bj)*(neu_val+streamice_addl_backstress)  ! note negative sign is due to direction of normal vector
226               taudx_SI(i,j+1,bi,bj) = taudx_SI(i,j+1,bi,bj) -               taudx_SI(i,j+1,bi,bj) = taudx_SI(i,j+1,bi,bj) -
227       &        .5 * dyG(i,j,bi,bj) * neu_val       &        .5 * dyG(i,j,bi,bj)*(neu_val+streamice_addl_backstress)
228             ENDIF             ENDIF
229    
230             IF ((STREAMICE_ufacemask(i+1,j,bi,bj) .eq. 2)             IF ((STREAMICE_ufacemask(i+1,j,bi,bj) .eq. 2)
# Line 232  C     LOCAL VARIABLES Line 232  C     LOCAL VARIABLES
232       &      .OR. (STREAMICE_hmask(i+1,j,bi,bj) .eq. 2) ) THEN       &      .OR. (STREAMICE_hmask(i+1,j,bi,bj) .eq. 2) ) THEN
233                            
234               taudx_SI(i+1,j,bi,bj) = taudx_SI(i+1,j,bi,bj) +               taudx_SI(i+1,j,bi,bj) = taudx_SI(i+1,j,bi,bj) +
235       &        .5 * dyG(i+1,j,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector       &        .5 * dyG(i+1,j,bi,bj)*(neu_val+streamice_addl_backstress)  ! note negative sign is due to direction of normal vector
236               taudx_SI(i+1,j+1,bi,bj) = taudx_SI(i+1,j+1,bi,bj) +               taudx_SI(i+1,j+1,bi,bj) = taudx_SI(i+1,j+1,bi,bj) +
237       &        .5 * dyG(i+1,j,bi,bj) * neu_val       &        .5 * dyG(i+1,j,bi,bj)*(neu_val+streamice_addl_backstress)
238             ENDIF             ENDIF
239    
240             IF ((STREAMICE_vfacemask(i,j,bi,bj) .eq. 2)             IF ((STREAMICE_vfacemask(i,j,bi,bj) .eq. 2)
# Line 242  C     LOCAL VARIABLES Line 242  C     LOCAL VARIABLES
242       &      .OR. (STREAMICE_hmask(i,j-1,bi,bj) .eq. 2) ) THEN       &      .OR. (STREAMICE_hmask(i,j-1,bi,bj) .eq. 2) ) THEN
243                            
244               taudy_SI(i,j,bi,bj) = taudy_SI(i,j,bi,bj) -               taudy_SI(i,j,bi,bj) = taudy_SI(i,j,bi,bj) -
245       &        .5 * dxG(i,j,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector       &        .5 * dxG(i,j,bi,bj)*(neu_val+streamice_addl_backstress)  ! note negative sign is due to direction of normal vector
246               taudy_SI(i+1,j,bi,bj) = taudy_SI(i+1,j,bi,bj) -               taudy_SI(i+1,j,bi,bj) = taudy_SI(i+1,j,bi,bj) -
247       &        .5 * dxG(i,j,bi,bj) * neu_val       &        .5 * dxG(i,j,bi,bj)*(neu_val+streamice_addl_backstress)
248             ENDIF             ENDIF
249    
250             IF ((STREAMICE_vfacemask(i,j+1,bi,bj) .eq. 2)             IF ((STREAMICE_vfacemask(i,j+1,bi,bj) .eq. 2)
# Line 252  C     LOCAL VARIABLES Line 252  C     LOCAL VARIABLES
252       &      .OR. (STREAMICE_hmask(i,j+1,bi,bj) .eq. 2) ) THEN       &      .OR. (STREAMICE_hmask(i,j+1,bi,bj) .eq. 2) ) THEN
253                            
254               taudy_SI(i,j+1,bi,bj) = taudy_SI(i,j+1,bi,bj) +               taudy_SI(i,j+1,bi,bj) = taudy_SI(i,j+1,bi,bj) +
255       &        .5 * dxG(i,j+1,bi,bj) * neu_val  ! note negative sign is due to direction of normal vector       &        .5 * dxG(i,j+1,bi,bj)*(neu_val+streamice_addl_backstress)  ! note negative sign is due to direction of normal vector
256               taudy_SI(i+1,j+1,bi,bj) = taudy_SI(i+1,j+1,bi,bj) +               taudy_SI(i+1,j+1,bi,bj) = taudy_SI(i+1,j+1,bi,bj) +
257       &        .5 * dxG(i,j+1,bi,bj) * neu_val       &        .5 * dxG(i,j+1,bi,bj)*(neu_val+streamice_addl_backstress)
258             ENDIF             ENDIF
259    
260            ENDIF            ENDIF
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22