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

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

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

revision 1.2 by heimbach, Fri Sep 28 02:55:40 2012 UTC revision 1.3 by dgoldberg, Mon Dec 10 02:34:45 2012 UTC
# Line 368  C     LOCAL VARIABLES Line 368  C     LOCAL VARIABLES
368            !     ice in the current cell            !     ice in the current cell
369                            
370               taudx_SI(i,j,bi,bj) = taudx_SI(i,j,bi,bj) -               taudx_SI(i,j,bi,bj) = taudx_SI(i,j,bi,bj) -
371       &        .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)
372               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) -
373       &        .5 * dyG(i,j,bi,bj) * neu_val       &        .5 * dyG(i,j,bi,bj)*(neu_val+streamice_addl_backstress)
374             ENDIF             ENDIF
375    
376             IF ((STREAMICE_ufacemask(i+1,j,bi,bj) .eq. 2)             IF ((STREAMICE_ufacemask(i+1,j,bi,bj) .eq. 2)
# Line 378  C     LOCAL VARIABLES Line 378  C     LOCAL VARIABLES
378       &      .OR. (STREAMICE_hmask(i+1,j,bi,bj) .eq. 2) ) THEN       &      .OR. (STREAMICE_hmask(i+1,j,bi,bj) .eq. 2) ) THEN
379                            
380               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) +
381       &        .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
382               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) +
383       &        .5 * dyG(i+1,j,bi,bj) * neu_val       &        .5 * dyG(i+1,j,bi,bj)*(neu_val+streamice_addl_backstress)
384             ENDIF             ENDIF
385    
386             IF ((STREAMICE_vfacemask(i,j,bi,bj) .eq. 2)             IF ((STREAMICE_vfacemask(i,j,bi,bj) .eq. 2)
# Line 388  C     LOCAL VARIABLES Line 388  C     LOCAL VARIABLES
388       &      .OR. (STREAMICE_hmask(i,j-1,bi,bj) .eq. 2) ) THEN       &      .OR. (STREAMICE_hmask(i,j-1,bi,bj) .eq. 2) ) THEN
389                            
390               taudy_SI(i,j,bi,bj) = taudy_SI(i,j,bi,bj) -               taudy_SI(i,j,bi,bj) = taudy_SI(i,j,bi,bj) -
391       &        .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)
392               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) -
393       &        .5 * dxG(i,j,bi,bj) * neu_val       &        .5 * dxG(i,j,bi,bj)*(neu_val+streamice_addl_backstress)
394             ENDIF             ENDIF
395    
396             IF ((STREAMICE_vfacemask(i,j+1,bi,bj) .eq. 2)             IF ((STREAMICE_vfacemask(i,j+1,bi,bj) .eq. 2)
# Line 398  C     LOCAL VARIABLES Line 398  C     LOCAL VARIABLES
398       &      .OR. (STREAMICE_hmask(i,j+1,bi,bj) .eq. 2) ) THEN       &      .OR. (STREAMICE_hmask(i,j+1,bi,bj) .eq. 2) ) THEN
399                            
400               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) +
401       &        .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)
402               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) +
403       &        .5 * dxG(i,j+1,bi,bj) * neu_val       &        .5 * dxG(i,j+1,bi,bj)*(neu_val+streamice_addl_backstress)
404             ENDIF             ENDIF
405    
406            ENDIF            ENDIF
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22