/[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.1 by dgoldberg, Thu Sep 27 23:37:39 2012 UTC revision 1.4 by dgoldberg, Wed Aug 27 19:29:13 2014 UTC
# Line 350  C     LOCAL VARIABLES Line 350  C     LOCAL VARIABLES
350  #else  #else
351              neu_val = .5 * gravity *              neu_val = .5 * gravity *
352       &       (streamice_density * H_streamice (i,j,bi,bj) ** 2 -       &       (streamice_density * H_streamice (i,j,bi,bj) ** 2 -
353       &        streamice_density_ocean_avg * R_low(i,j,bi,bj) ** 2       &        streamice_density_ocean_avg * R_low(i,j,bi,bj) ** 2)
354  #endif  #endif
355             ELSE             ELSE
356              neu_val = .5 * gravity *              neu_val = .5 * gravity *
# 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.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22