/[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.5 by heimbach, Mon Apr 16 22:38:24 2007 UTC revision 1.6 by jmc, Sun Apr 29 23:48:44 2007 UTC
# Line 124  C     if ice area (=ice fraction * grid- Line 124  C     if ice area (=ice fraction * grid-
124  C     thickness * grid-cell area) are too small (i.e.: < areaEpsil,vol_Epsil)  C     thickness * grid-cell area) are too small (i.e.: < areaEpsil,vol_Epsil)
125  C     will assume that ice is gone, and will loose mass or energy.  C     will assume that ice is gone, and will loose mass or energy.
126  C     However, if areaEpsil,vol_Epsil are much smaller than minimun ice area  C     However, if areaEpsil,vol_Epsil are much smaller than minimun ice area
127  C     (iceMaskMin*rAc) and minimum ice volume (iceMaskMin*himin*rAc),  C     (iceMaskMin*rAc) and minimum ice volume (iceMaskMin*hIceMin*rAc),
128  C     good chance that this will never happen within 1 time step.  C     good chance that this will never happen within 1 time step.
129    
130        dBugFlag = debugLevel.GE.debLevB        dBugFlag = debugLevel.GE.debLevB
# Line 433  C      and adjust Ice thickness and snow Line 433  C      and adjust Ice thickness and snow
433  C-     adjust sea-ice state if not enough ice.  C-     adjust sea-ice state if not enough ice.
434          DO j=1,sNy          DO j=1,sNy
435           DO i=1,sNx           DO i=1,sNx
436            IF ( iceHeight(i,j,bi,bj).LT.himin ) THEN            IF ( iceHeight(i,j,bi,bj).LT.hIceMin ) THEN
437  C-    Not enough ice, melt the tiny amount of snow & ice:  C-    Not enough ice, melt the tiny amount of snow & ice:
438  C     and return fresh-water, salt & energy to the ocean (flx > 0 = into ocean)  C     and return fresh-water, salt & energy to the ocean (flx > 0 = into ocean)
439  C- -  Note: using 1rst.Order Upwind, I can get the same results as when  C- -  Note: using 1rst.Order Upwind, I can get the same results as when
# Line 443  C- - Line 443  C- -
443              oceFWfx(i,j,bi,bj) =  ( rhos*snowHeight(i,j,bi,bj)              oceFWfx(i,j,bi,bj) =  ( rhos*snowHeight(i,j,bi,bj)
444       &                             +rhoi*iceHeight(i,j,bi,bj)       &                             +rhoi*iceHeight(i,j,bi,bj)
445       &                            )*iceMask(i,j,bi,bj)/thSIce_deltaT       &                            )*iceMask(i,j,bi,bj)/thSIce_deltaT
446              oceSflx(i,j,bi,bj) =    rhoi*iceHeight(i,j,bi,bj)*saltice              oceSflx(i,j,bi,bj) =    rhoi*iceHeight(i,j,bi,bj)*saltIce
447       &                             *iceMask(i,j,bi,bj)/thSIce_deltaT       &                             *iceMask(i,j,bi,bj)/thSIce_deltaT
448              oceQnet(i,j,bi,bj) = -( rhos*snowHeight(i,j,bi,bj)*qsnow              oceQnet(i,j,bi,bj) = -( rhos*snowHeight(i,j,bi,bj)*qsnow
449       &                             +rhoi*iceHeight(i,j,bi,bj)       &                             +rhoi*iceHeight(i,j,bi,bj)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22