/[MITgcm]/MITgcm/pkg/seaice/seaice_growth.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_growth.F

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

revision 1.8 by mlosch, Wed Dec 20 20:49:12 2006 UTC revision 1.9 by mlosch, Thu Dec 21 10:50:46 2006 UTC
# Line 477  C     freezing, add precip as snow Line 477  C     freezing, add precip as snow
477            ELSE            ELSE
478  C     add precip as rain, water converted into equivalent m of  C     add precip as rain, water converted into equivalent m of
479  C     ice by 1/ICE_DENS.  C     ice by 1/ICE_DENS.
480  C     Sign issue? frWtrIce (>0 for more ice)  C     Do not get confused by the sign:
481  C     has the opposite sign of precip (>0 for downward flux of water)  C     precip   > 0 for downward flux of fresh water
482    C     frWtrIce > 0 for more ice (corresponds to an upward "fresh water flux"),
483    C     so that here the rain is added *as if* it is melted ice (which is not
484    C     true, but just a trick; physically the rain just runs as water
485    C     through the ice into the ocean)
486             frWtrIce(I,J) = frWtrIce(I,J)             frWtrIce(I,J) = frWtrIce(I,J)
487       &            -PRECIP(I,J,bi,bj)*AREA(I,J,2,bi,bj)*       &            -PRECIP(I,J,bi,bj)*AREA(I,J,2,bi,bj)*
488       &            SEAICE_deltaTtherm/ICE_DENS       &            SEAICE_deltaTtherm/ICE_DENS
489            ENDIF            ENDIF
490    
 C     now add in precip over open water directly into ocean as negative salt  
 C     Sign issue? frWtrIce (>0 for more ice)  
 C     has the opposite sign of precip (>0 for downward flux of water)  
           frWtrIce(I,J) = frWtrIce(I,J)  
      &         -PRECIP(I,J,bi,bj)*(ONE-AREA(I,J,2,bi,bj))  
      &         *SEAICE_deltaTtherm/ICE_DENS  
   
491  C     Now melt snow if there is residual heat left in surface level  C     Now melt snow if there is residual heat left in surface level
492  C     Note that units of YNEG and frWtrIce are m of ice  C     Note that units of YNEG and frWtrIce are m of ice
493  cph( very sensitive bit here by JZ  cph( very sensitive bit here by JZ
# Line 506  cph) Line 503  cph)
503    
504  C NOW GET FRESH WATER FLUX  C NOW GET FRESH WATER FLUX
505            EmPmR(I,J,bi,bj)  = maskC(I,J,kSurface,bi,bj)*(            EmPmR(I,J,bi,bj)  = maskC(I,J,kSurface,bi,bj)*(
506       &         EVAP(I,J,bi,bj)*(ONE-AREA(I,J,2,bi,bj))       &         ( EVAP(I,J,bi,bj)-PRECIP(I,J,bi,bj) )
507       &         -RUNOFF(I,J,bi,bj)       &         * ( ONE - AREA(I,J,2,bi,bj) )
508       &         +frWtrIce(I,J)*ICE_DENS/SEAICE_deltaTtherm       &         - RUNOFF(I,J,bi,bj)
509         &         + frWtrIce(I,J)*ICE_DENS/SEAICE_deltaTtherm
510       &         )       &         )
511    
512  C NOW GET TOTAL QNET AND QSW  C NOW GET TOTAL QNET AND QSW
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22