/[MITgcm]/MITgcm_contrib/verification_other/shelfice_remeshing/code/shelfice_step_icemass.F
ViewVC logotype

Diff of /MITgcm_contrib/verification_other/shelfice_remeshing/code/shelfice_step_icemass.F

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

revision 1.2 by dgoldberg, Thu May 5 18:16:04 2016 UTC revision 1.3 by ksnow, Mon Nov 14 11:17:53 2016 UTC
# Line 53  C  i,j                  :: loop indices Line 53  C  i,j                  :: loop indices
53           DO bi = myBxLo(myThid), myBxHi(myThid)           DO bi = myBxLo(myThid), myBxHi(myThid)
54            DO j=1-OLy,sNy+OLy-1            DO j=1-OLy,sNy+OLy-1
55             DO i=1-OLx+1,sNx+OLx-1             DO i=1-OLx+1,sNx+OLx-1
56              if (streamice_hmask(i,j,bi,bj).eq.1 .or.              if (streamice_hmask(i,j,bi,bj).eq.1) then
      &          streamice_hmask(i,j,bi,bj).eq.2) then  
   
57               shelficeMass(i,j,bi,bj) =               shelficeMass(i,j,bi,bj) =
58       &        H_streamice(I,J,bi,bj) * streamice_density       &        H_streamice(I,J,bi,bj) * streamice_density
59                elseif (streamice_hmask(i,j,bi,bj).eq.2) then
60                 shelficeMass(i,j,bi,bj) =
61         &        H_streamice(I,J,bi,bj) * streamice_density
62         &        * area_shelf_streamice(i,j,bi,bj)*recip_rA(i,j,bi,bj)
63              endif              endif
64             ENDDO             ENDDO
65            ENDDO            ENDDO
# Line 80  C  i,j                  :: loop indices Line 81  C  i,j                  :: loop indices
81    
82              shelficeMass(i,j,bi,bj) = shelficeMass(i,j,bi,bj)              shelficeMass(i,j,bi,bj) = shelficeMass(i,j,bi,bj)
83       &      + shelfIceMassDynTendency(I,J,bi,bj)*deltaT       &      + shelfIceMassDynTendency(I,J,bi,bj)*deltaT
84                
85             ENDDO             ENDDO
86            ENDDO            ENDDO
87           ENDDO           ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22