344 |
gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj) |
gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj) |
345 |
& -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj) |
& -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj) |
346 |
& -swfracb(2)*maskC(i,j,kp1, bi,bj)) |
& -swfracb(2)*maskC(i,j,kp1, bi,bj)) |
347 |
& *recip_Cp*recip_rhoConst |
& *recip_Cp*mass2rUnit |
348 |
& *recip_drF(klev)*_recip_hFacC(i,j,kLev,bi,bj) |
& *recip_drF(klev)*_recip_hFacC(i,j,kLev,bi,bj) |
349 |
ENDDO |
ENDDO |
350 |
ENDDO |
ENDDO |
352 |
#endif |
#endif |
353 |
|
|
354 |
#ifdef ALLOW_RBCS |
#ifdef ALLOW_RBCS |
355 |
if (useRBCS) then |
IF (useRBCS) THEN |
356 |
call RBCS_ADD_TENDENCY(bi,bj,klev, 1, |
CALL RBCS_ADD_TENDENCY(bi,bj,klev, 1, |
357 |
& myTime, myThid ) |
& myTime, myThid ) |
358 |
endif |
ENDIF |
359 |
#endif |
#endif |
360 |
|
|
361 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
501 |
#endif /* ALLOW_SEAICE */ |
#endif /* ALLOW_SEAICE */ |
502 |
IF ( SaltPlumeDepth(i,j,bi,bj) .GT. -rF(kLev) ) THEN |
IF ( SaltPlumeDepth(i,j,bi,bj) .GT. -rF(kLev) ) THEN |
503 |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
504 |
& +saltPlume*horiVertRatio*recip_rhoConst |
& +saltPlume*mass2rUnit |
505 |
& *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj) |
& *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj) |
506 |
& *min(drF(kLev),SaltPlumeDepth(i,j,bi,bj)+rF(kLev)) |
& *min(drF(kLev),SaltPlumeDepth(i,j,bi,bj)+rF(kLev)) |
507 |
& /SaltPlumeDepth(i,j,bi,bj) |
& /SaltPlumeDepth(i,j,bi,bj) |
511 |
#endif /* ALLOW_SALT_PLUME */ |
#endif /* ALLOW_SALT_PLUME */ |
512 |
|
|
513 |
#ifdef ALLOW_RBCS |
#ifdef ALLOW_RBCS |
514 |
if (useRBCS) then |
IF (useRBCS) THEN |
515 |
call RBCS_ADD_TENDENCY(bi,bj,klev, 2, |
CALL RBCS_ADD_TENDENCY(bi,bj,klev, 2, |
516 |
& myTime, myThid ) |
& myTime, myThid ) |
517 |
endif |
ENDIF |
518 |
#endif /* ALLOW_RBCS */ |
#endif /* ALLOW_RBCS */ |
519 |
|
|
520 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |