/[MITgcm]/MITgcm/pkg/thsice/thsice_advdiff.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_advdiff.F

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

revision 1.13 by jmc, Mon Jan 21 22:54:54 2013 UTC revision 1.14 by heimbach, Tue Jan 22 23:31:09 2013 UTC
# Line 89  C     msgBuf    :: Informational/error m Line 89  C     msgBuf    :: Informational/error m
89        _RL iceFrc    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL iceFrc    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
90        _RL iceVol    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL iceVol    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
91        _RL oldVol    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL oldVol    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
92          _RL iceTmp    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
93        _RL r_minArea        _RL r_minArea
94        _RL meanCellArea, areaEpsil, vol_Epsil        _RL meanCellArea, areaEpsil, vol_Epsil
95  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
# Line 167  C-     Ice-Volume weighted quantities: Line 168  C-     Ice-Volume weighted quantities:
168       &       DIAGNOSTICS_IS_ON('SI_AdvQ2',myThid) ) THEN       &       DIAGNOSTICS_IS_ON('SI_AdvQ2',myThid) ) THEN
169           DO j=1,sNy           DO j=1,sNy
170            DO i=1,sNx            DO i=1,sNx
171             iceVol(i,j) = iceMask(i,j,bi,bj)*iceHeight(i,j,bi,bj)             iceTmp(i,j) = iceMask(i,j,bi,bj)*iceHeight(i,j,bi,bj)
172            ENDDO            ENDDO
173           ENDDO           ENDDO
174           CALL DIAGNOSTICS_FRACT_FILL(           CALL DIAGNOSTICS_FRACT_FILL(
175       I                   Qice1(1-OLx,1-OLy,bi,bj),       I                   Qice1(1-OLx,1-OLy,bi,bj),
176       I                   iceVol,tmpFac,1,'SI_AdvQ1',       I                   iceTmp,tmpFac,1,'SI_AdvQ1',
177       I                   0,1,2,bi,bj,myThid)       I                   0,1,2,bi,bj,myThid)
178           CALL DIAGNOSTICS_FRACT_FILL(           CALL DIAGNOSTICS_FRACT_FILL(
179       I                   Qice2(1-OLx,1-OLy,bi,bj),       I                   Qice2(1-OLx,1-OLy,bi,bj),
180       I                   iceVol,tmpFac,1,'SI_AdvQ2',       I                   iceTmp,tmpFac,1,'SI_AdvQ2',
181       I                   0,1,2,bi,bj,myThid)       I                   0,1,2,bi,bj,myThid)
182          ENDIF          ENDIF
183        ENDIF        ENDIF

Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14

  ViewVC Help
Powered by ViewVC 1.1.22