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

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

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

revision 1.4 by cnh, Wed Jun 10 01:44:03 1998 UTC revision 1.11 by adcroft, Mon Mar 27 22:25:44 2000 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  CStartOfInterface  CStartOfInterface
6        SUBROUTINE INI_SALT ( myThid )        SUBROUTINE INI_SALT ( myThid )
# Line 18  C     |     integration. Line 18  C     |     integration.
18  C     | In addition to setting the salinity field we also        |  C     | In addition to setting the salinity field we also        |
19  C     | set the initial salinity tendency term here.             |  C     | set the initial salinity tendency term here.             |
20  C     \==========================================================/  C     \==========================================================/
21          IMPLICIT NONE
22    
23  C     === Global variables ===  C     === Global variables ===
24  #include "SIZE.h"  #include "SIZE.h"
# Line 32  C     myThid -  Number of this instance Line 33  C     myThid -  Number of this instance
33  CEndOfInterface  CEndOfInterface
34    
35  C     == Local variables ==  C     == Local variables ==
 C     iC, jC - Center of domain  
 C     iD, jD - Disitance from domain center.  
 C     rad    - Radius of initial patch  
 C     rD     - Radial displacement of point I,J  
 C     iG, jG - Global coordinate index  
36  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
37  C     I,J,K  C     I,J,K
       INTEGER iC, jC, iD, jD  
       INTEGER iG, jG  
38        INTEGER bi, bj        INTEGER bi, bj
39        INTEGER  I,  J, K        INTEGER  I,  J, K
       REAL rad, rD  
40    
41        _BARRIER        _BARRIER
42    
# Line 51  C     I,J,K Line 44  C     I,J,K
44  C--    Initialise salinity field to the vertical reference profile  C--    Initialise salinity field to the vertical reference profile
45         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
46          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
47           DO K=1,Nz           DO K=1,Nr
48            DO J=1,sNy            DO J=1,sNy
49             DO I=1,sNx             DO I=1,sNx
50              salt(I,J,K,bi,bj) = sRef(K)              salt(I,J,K,bi,bj) = sRef(K)
# Line 68  C--    Initialise salinity field to the Line 61  C--    Initialise salinity field to the
61  C     Set initial tendency terms  C     Set initial tendency terms
62        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
63         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
64          DO K=1,Nz          DO K=1,Nr
65           DO J=1,sNy           DO J=1,sNy
66            DO I=1,sNx            DO I=1,sNx
67             gs   (I,J,K,bi,bj) = 0. _d 0             gs   (I,J,K,bi,bj) = 0. _d 0
68             gsNM1(I,J,K,bi,bj) = 0. _d 0             gsNM1(I,J,K,bi,bj) = 0. _d 0
69               IF (hFacC(I,J,K,bi,bj).EQ.0) salt(I,J,K,bi,bj) = 0.
70            ENDDO            ENDDO
71           ENDDO           ENDDO
72          ENDDO          ENDDO
# Line 82  C Line 76  C
76        _EXCH_XYZ_R8(salt  , myThid )        _EXCH_XYZ_R8(salt  , myThid )
77        _EXCH_XYZ_R8(gs , myThid )        _EXCH_XYZ_R8(gs , myThid )
78        _EXCH_XYZ_R8(gsNM1 , myThid )        _EXCH_XYZ_R8(gsNM1 , myThid )
79    
80          CALL PLOT_FIELD_XYZRL( salt, 'Initial Salinity' , Nr, 1, myThid )
81    
82        RETURN        RETURN
83        END        END
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22