pkg/salt_plume/salt_plume_forcing_surf.F

C $Header:  $
C $Name:  $

#include "SALT_PLUME_OPTIONS.h"

CBOP
C     !ROUTINE: SALT_PLUME_FORCING_SURF
C     !INTERFACE:
      SUBROUTINE SALT_PLUME_FORCING_SURF(
     I                            bi, bj, iMin, iMax, jMin, jMax,
     I                            myTime,myIter,myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R SALT_PLUME_FORCING_SURF
C     | o saltPlume is the amount of salt rejected by ice while freezing;
C     |   it is here subtracted from surfaceForcingS and will be redistributed
C     |   to multiple vertical levels later on as per Duffy et al. (GRL 1999)
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "FFIELDS.h"
#include "DYNVARS.h"
#include "SALT_PLUME.h"

C !INPUT PARAMETERS: ===================================================
C  bi,bj                :: tile indices
C  myTime               :: model time
C  myIter               :: time-step number
C  myThid               :: thread number
      INTEGER bi, bj, iMin, iMax, jMin, jMax
      _RL myTime
      INTEGER myIter
      INTEGER myThid

C     !LOCAL VARIABLES:
C     == Local variables ==
C     i,j                  :: loop indices
C     ks                   :: surface level index
      INTEGER i, j, ks
CEOP

      IF ( usingPCoords ) THEN
         ks = Nr
      ELSE
         ks = 1
      ENDIF

      DO j = jMin, jMax
         DO i = iMin, iMax
            surfaceForcingS(i,j,bi,bj) = surfaceForcingS(i,j,bi,bj)
     &           - saltPlumeFlux(i,j,bi,bj) * mass2rUnit
         ENDDO
      ENDDO
     
      RETURN
      END
1	dimitri	1.1	C $Header: $
2			C $Name: $
3
4			#include "SALT_PLUME_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: SALT_PLUME_FORCING_SURF
8			C !INTERFACE:
9			SUBROUTINE SALT_PLUME_FORCING_SURF(
10			I bi, bj, iMin, iMax, jMin, jMax,
11			I myTime,myIter,myThid )
12			C !DESCRIPTION: \bv
13			C ==========================================================
14			C \| S/R SALT_PLUME_FORCING_SURF
15			C \| o saltPlume is the amount of salt rejected by ice while freezing;
16			C \| it is here subtracted from surfaceForcingS and will be redistributed
17			C \| to multiple vertical levels later on as per Duffy et al. (GRL 1999)
18			C ==========================================================
19			C \ev
20
21			C !USES:
22			IMPLICIT NONE
23			C == Global variables ==
24			#include "SIZE.h"
25			#include "EEPARAMS.h"
26			#include "PARAMS.h"
27			#include "FFIELDS.h"
28			#include "DYNVARS.h"
29			#include "SALT_PLUME.h"
30
31			C !INPUT PARAMETERS: ===================================================
32			C bi,bj :: tile indices
33			C myTime :: model time
34			C myIter :: time-step number
35			C myThid :: thread number
36			INTEGER bi, bj, iMin, iMax, jMin, jMax
37			_RL myTime
38			INTEGER myIter
39			INTEGER myThid
40
41			C !LOCAL VARIABLES:
42			C == Local variables ==
43			C i,j :: loop indices
44			C ks :: surface level index
45			INTEGER i, j, ks
46			CEOP
47
48			IF ( usingPCoords ) THEN
49			ks = Nr
50			ELSE
51			ks = 1
52			ENDIF
53
54			DO j = jMin, jMax
55			DO i = iMin, iMax
56			surfaceForcingS(i,j,bi,bj) = surfaceForcingS(i,j,bi,bj)
57			& - saltPlumeFlux(i,j,bi,bj) * mass2rUnit
58			ENDDO
59			ENDDO
60
61			RETURN
62			END