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

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

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

revision 1.25 by jmc, Thu Aug 30 15:02:51 2007 UTC revision 1.26 by jmc, Mon Oct 1 13:38:08 2007 UTC
# Line 268  CADJ STORE  icemask(i,j,bi,bj) = comlev1 Line 268  CADJ STORE  icemask(i,j,bi,bj) = comlev1
268            icFlxAtm(i,j,bi,bj) = icFrac*icFlxAtm(i,j,bi,bj)            icFlxAtm(i,j,bi,bj) = icFrac*icFlxAtm(i,j,bi,bj)
269       &                        - opFrac*Qnet(i,j,bi,bj)       &                        - opFrac*Qnet(i,j,bi,bj)
270            icFrwAtm(i,j,bi,bj) = icFrac*icFrwAtm(i,j,bi,bj)            icFrwAtm(i,j,bi,bj) = icFrac*icFrwAtm(i,j,bi,bj)
271       &                        + opFrac*rhofw*EmPmR(i,j,bi,bj)       &                        + opFrac*EmPmR(i,j,bi,bj)
272            Qnet(i,j,bi,bj) = -icFrac*flx2oc(i,j) + opFrac*Qnet(i,j,bi,bj)            Qnet(i,j,bi,bj) = -icFrac*flx2oc(i,j) + opFrac*Qnet(i,j,bi,bj)
273            EmPmR(i,j,bi,bj)= -icFrac*frw2oc(i,j)/rhofw            EmPmR(i,j,bi,bj)= -icFrac*frw2oc(i,j)
274       &                    +  opFrac*EmPmR(i,j,bi,bj)       &                    +  opFrac*EmPmR(i,j,bi,bj)
275            saltFlux(i,j,bi,bj) = -icFrac*fsalt(i,j)            saltFlux(i,j,bi,bj) = -icFrac*fsalt(i,j)
276  #endif  #endif
# Line 282  CADJ STORE  icemask(i,j,bi,bj) = comlev1 Line 282  CADJ STORE  icemask(i,j,bi,bj) = comlev1
282       &           snowHeight(i,j,bi,bj), Qnet(i,j,bi,bj)       &           snowHeight(i,j,bi,bj), Qnet(i,j,bi,bj)
283  #endif  #endif
284    
285          ELSEIF (hOceMxL(i,j,bi,bj).gt.0. _d 0) THEN          ELSEIF (hOceMxL(i,j,bi,bj).GT.0. _d 0) THEN
286            icFlxAtm(i,j,bi,bj) = -Qnet(i,j,bi,bj)            icFlxAtm(i,j,bi,bj) = -Qnet(i,j,bi,bj)
287            icFrwAtm(i,j,bi,bj) = rhofw*EmPmR(i,j,bi,bj)            icFrwAtm(i,j,bi,bj) = EmPmR(i,j,bi,bj)
288          ELSE          ELSE
289            icFlxAtm(i,j,bi,bj) = 0. _d 0            icFlxAtm(i,j,bi,bj) = 0. _d 0
290            icFrwAtm(i,j,bi,bj) = 0. _d 0            icFrwAtm(i,j,bi,bj) = 0. _d 0
# Line 322  C--    Net fluxes : Line 322  C--    Net fluxes :
322            sFluxFromIce(i,j)= sFluxFromIce(i,j) - fsalt(i,j)            sFluxFromIce(i,j)= sFluxFromIce(i,j) - fsalt(i,j)
323  #else  #else
324            Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) - flx2oc(i,j)            Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) - flx2oc(i,j)
325            EmPmR(i,j,bi,bj)= EmPmR(i,j,bi,bj)- frw2oc(i,j)/rhofw            EmPmR(i,j,bi,bj)= EmPmR(i,j,bi,bj)- frw2oc(i,j)
326            saltFlux(i,j,bi,bj)=saltFlux(i,j,bi,bj) - fsalt(i,j)            saltFlux(i,j,bi,bj)=saltFlux(i,j,bi,bj) - fsalt(i,j)
327  #endif  #endif
328    
# Line 374  C--   note: those fluxes should to be ad Line 374  C--   note: those fluxes should to be ad
374          DO i = iMin, iMax          DO i = iMin, iMax
375           IF ( hOceMxL(i,j,bi,bj).GT.0. _d 0 ) THEN           IF ( hOceMxL(i,j,bi,bj).GT.0. _d 0 ) THEN
376            Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) - oceQnet(i,j,bi,bj)            Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) - oceQnet(i,j,bi,bj)
377            EmPmR(i,j,bi,bj)= EmPmR(i,j,bi,bj)- oceFWfx(i,j,bi,bj)/rhofw            EmPmR(i,j,bi,bj)= EmPmR(i,j,bi,bj)- oceFWfx(i,j,bi,bj)
378            saltFlux(i,j,bi,bj)=saltFlux(i,j,bi,bj) - oceSflx(i,j,bi,bj)            saltFlux(i,j,bi,bj)=saltFlux(i,j,bi,bj) - oceSflx(i,j,bi,bj)
379           ENDIF           ENDIF
380          ENDDO          ENDDO
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22