pkg/salt_plume/salt_plume_forcing_surf.F

C $Header: /u/gcmpack/MITgcm/pkg/salt_plume/salt_plume_forcing_surf.F,v 1.1 2007/11/28 00:18:17 dimitri Exp $
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
CEOP

#ifdef ALLOW_SALT_PLUME

C     !LOCAL VARIABLES:
C     i,j                  :: loop indices
C     ks                   :: surface level index
      INTEGER i, j, ks

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