pkg/salt_plume/salt_plume_init_varia.F

e4775240e5 Dimi*0001 #include "SALT_PLUME_OPTIONS.h"
                0002
                0003 CBOP 0
                0004 C     !ROUTINE: SALT_PLUME_INIT_VARIA
                0005 C     !INTERFACE:
                0006       SUBROUTINE SALT_PLUME_INIT_VARIA( myThid )
                0007
                0008 C     !DESCRIPTION:
                0009 C     Initialize SALT_PLUME variables
                0010
                0011 C     !USES:
                0012       IMPLICIT NONE
                0013 #include "SIZE.h"
                0014 #include "EEPARAMS.h"
                0015 #include "SALT_PLUME.h"
                0016
                0017 C     !INPUT/OUTPUT PARAMETERS:
                0018 C     myThid ::  my Thread Id number
                0019       INTEGER myThid
                0020 CEOP
                0021
762bde3a3e Dimi*0022 #ifdef ALLOW_SALT_PLUME
                0023
e4775240e5 Dimi*0024 C     !LOCAL VARIABLES:
ed6012c5a0 Jean*0025       INTEGER i,j, bi,bj,Nrp1
                0026 #ifdef SALT_PLUME_VOLUME
                0027       INTEGER k
                0028 #endif
e4775240e5 Dimi*0029
1f89baba18 Patr*0030       Nrp1=Nr+1
e4775240e5 Dimi*0031       DO bj = myByLo(myThid), myByHi(myThid)
                0032        DO bi = myBxLo(myThid), myBxHi(myThid)
                0033
                0034 C     Initialize arrays in common blocks :
                0035 C     2-dim.
ed6012c5a0 Jean*0036         DO j=1-OLy,sNy+OLy
                0037          DO i=1-OLx,sNx+OLx
e4775240e5 Dimi*0038           SaltPlumeDepth (i,j,bi,bj) = 0. _d 0
                0039           saltPlumeFlux  (i,j,bi,bj) = 0. _d 0
1f89baba18 Patr*0040 #ifdef SALT_PLUME_VOLUME
                0041           SPbrineVolFlux (i,j,bi,bj) = 0. _d 0
                0042           SPforcS1       (i,j,bi,bj) = 0. _d 0
                0043           DO k=1,Nr
                0044            dSPvolSurf2kLev (i,j,k,bi,bj) = 0. _d 0
                0045            dSPvolkLev2Above(i,j,k,bi,bj) = 0. _d 0
                0046            SPforcingS      (i,j,k,bi,bj) = 0. _d 0
                0047            SPforcingT      (i,j,k,bi,bj) = 0. _d 0
                0048            SPplumek        (i,j,k,bi,bj) = 1. _d 0
                0049           ENDDO
                0050           dSPvolkLev2Above (i,j,Nrp1,bi,bj) = 0. _d 0
                0051           SPplumek         (i,j,Nrp1,bi,bj) = 1. _d 0
                0052 #endif
e4775240e5 Dimi*0053          ENDDO
                0054         ENDDO
                0055
                0056 C- end bi,bj loops
                0057        ENDDO
                0058       ENDDO
                0059
762bde3a3e Dimi*0060 #endif /* ALLOW_SALT_PLUME */
                0061
e4775240e5 Dimi*0062       RETURN
                0063       END