573 |
IF ( saltWtrIce(I,J) .GE. 0.0 ) THEN |
IF ( saltWtrIce(I,J) .GE. 0.0 ) THEN |
574 |
saltFlux(I,J,bi,bj) = HEFFM(I,J,bi,bj)*saltWtrIce(I,J)* |
saltFlux(I,J,bi,bj) = HEFFM(I,J,bi,bj)*saltWtrIce(I,J)* |
575 |
& ICE2WATR*rhoConstFresh*SEAICE_salinity* |
& ICE2WATR*rhoConstFresh*SEAICE_salinity* |
576 |
& salt(I,J,1,bi,bj)/SEAICE_deltaTtherm |
& salt(I,j,kSurface,bi,bj)/SEAICE_deltaTtherm |
577 |
|
#ifdef ALLOW_SALT_PLUME |
578 |
|
C saltPlumeFlux is defined only during freezing: |
579 |
|
saltPlumeFlux(I,J,bi,bj)=HEFFM(I,J,bi,bj)*saltWtrIce(I,J)* |
580 |
|
& ICE2WATR*rhoConstFresh*(1-SEAICE_salinity)* |
581 |
|
& salt(I,j,kSurface,bi,bj)/SEAICE_deltaTtherm |
582 |
|
#endif /* ALLOW_SALT_PLUME */ |
583 |
C saltWtrIce < 0 : m of sea ice that is melted |
C saltWtrIce < 0 : m of sea ice that is melted |
584 |
ELSE |
ELSE |
585 |
saltFlux(I,J,bi,bj) = HEFFM(I,J,bi,bj)*saltWtrIce(I,J)* |
saltFlux(I,J,bi,bj) = HEFFM(I,J,bi,bj)*saltWtrIce(I,J)* |