/[MITgcm]/MITgcm/model/src/external_forcing.F
ViewVC logotype

Diff of /MITgcm/model/src/external_forcing.F

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

revision 1.43 by dimitri, Sun Jul 22 23:51:16 2007 UTC revision 1.44 by dimitri, Mon Jul 23 00:17:44 2007 UTC
# Line 398  C     == Global data == Line 398  C     == Global data ==
398  #include "SURFACE.h"  #include "SURFACE.h"
399  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
400  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
401  #include "SEAICE.h"  #include "SEAICE_PARAMS.h"
402  #endif /* ALLOW_SEAICE */  #endif /* ALLOW_SEAICE */
403  #endif /* ALLOW_SALT_PLUME */  #endif /* ALLOW_SALT_PLUME */
404    
# Line 422  C     kSurface  :: index of surface laye Line 422  C     kSurface  :: index of surface laye
422        INTEGER kSurface        INTEGER kSurface
423  CEOP  CEOP
424  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
425       _RL saltPlume        _RL saltPlume
426  #endif /* ALLOW_SALT_PLUME */  #endif /* ALLOW_SALT_PLUME */
427    
428        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
# Line 492  C Duffy et al. (GRL 1999) Line 492  C Duffy et al. (GRL 1999)
492            saltPlume = 0.            saltPlume = 0.
493  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
494            IF ( saltFlux(i,j,bi,bj) .GT. 0. .AND.            IF ( saltFlux(i,j,bi,bj) .GT. 0. .AND.
495       &         salt(i,j,ks,bi,bj)  .GT. SEAICE_salinity ) THEN       &         salt(i,j,kSurface,bi,bj)  .GT. SEAICE_salinity ) THEN
496             saltPlume = (salt(i,j,ks,bi,bj)-SEAICE_salinity) *             saltPlume = (salt(i,j,kSurface,bi,bj)-SEAICE_salinity) *
497       &          saltFlux(i,j,bi,bj) / salt(i,j,ks,bi,bj)       &          saltFlux(i,j,bi,bj) / salt(i,j,kSurface,bi,bj)
498            ENDIF            ENDIF
499  #endif /* ALLOW_SEAICE */  #endif /* ALLOW_SEAICE */
500            IF ( SaltPlumeDepth(i,j,bi,bj) .GT. -rF(kLev) ) THEN            IF ( SaltPlumeDepth(i,j,bi,bj) .GT. -rF(kLev) ) THEN
# Line 503  C Duffy et al. (GRL 1999) Line 503  C Duffy et al. (GRL 1999)
503       &          *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)       &          *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)
504       &          *max(drF(kLev),SaltPlumeDepth(i,j,bi,bj)+rF(kLev))       &          *max(drF(kLev),SaltPlumeDepth(i,j,bi,bj)+rF(kLev))
505       &          /SaltPlumeDepth(i,j,bi,bj)       &          /SaltPlumeDepth(i,j,bi,bj)
506              ENDIF
507          ENDDO          ENDDO
508         ENDDO         ENDDO
509  #endif /* ALLOW_SALT_PLUME */  #endif /* ALLOW_SALT_PLUME */
Legend:
Removed from v.1.43  
changed lines
  Added in v.1.44
  ViewVC Help
Powered by ViewVC 1.1.22