/[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.1 by cnh, Wed Apr 22 19:15:30 1998 UTC revision 1.9 by adcroft, Wed Dec 9 16:11:52 1998 UTC
# Line 1  Line 1 
1  C $Id$  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 45  C     I,J,K Line 46  C     I,J,K
46        INTEGER  I,  J, K        INTEGER  I,  J, K
47        REAL rad, rD        REAL rad, rD
48    
49  C--   Initialise salinity field to the vertical reference profile        _BARRIER
50        DO bj = myByLo(myThid), myByHi(myThid)  
51         DO bi = myBxLo(myThid), myBxHi(myThid)        IF ( hydrogSaltFile .EQ. ' ' ) THEN
52          DO K=1,Nz  C--    Initialise salinity field to the vertical reference profile
53           DO J=1,sNy         DO bj = myByLo(myThid), myByHi(myThid)
54            DO I=1,sNx          DO bi = myBxLo(myThid), myBxHi(myThid)
55             salt(I,J,K,bi,bj) = sRef(K)           DO K=1,Nr
56              DO J=1,sNy
57               DO I=1,sNx
58                salt(I,J,K,bi,bj) = sRef(K)
59               ENDDO
60            ENDDO            ENDDO
61           ENDDO           ENDDO
62          ENDDO          ENDDO
63         ENDDO         ENDDO
64        ENDDO        ELSE
65           _BEGIN_MASTER( myThid )
66           CALL READ_FLD_XYZ_RL( hydrogSaltFile, ' ', salt, 0, myThid )
67           _END_MASTER(myThid)
68          ENDIF
69  C     Set initial tendency terms  C     Set initial tendency terms
70        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
71         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
72          DO K=1,Nz          DO K=1,Nr
73           DO J=1,sNy           DO J=1,sNy
74            DO I=1,sNx            DO I=1,sNx
75             gs   (I,J,K,bi,bj) = 0. _d 0             gs   (I,J,K,bi,bj) = 0. _d 0
76             gsNM1(I,J,K,bi,bj) = 0. _d 0             gsNM1(I,J,K,bi,bj) = 0. _d 0
77               IF (hFacC(I,J,K,bi,bj).EQ.0) salt(I,J,K,bi,bj) = 0.
78            ENDDO            ENDDO
79           ENDDO           ENDDO
80          ENDDO          ENDDO
# Line 74  C Line 84  C
84        _EXCH_XYZ_R8(salt  , myThid )        _EXCH_XYZ_R8(salt  , myThid )
85        _EXCH_XYZ_R8(gs , myThid )        _EXCH_XYZ_R8(gs , myThid )
86        _EXCH_XYZ_R8(gsNM1 , myThid )        _EXCH_XYZ_R8(gsNM1 , myThid )
87    
88          CALL PLOT_FIELD_XYZRL( salt, 'Initial Salinity' , Nr, 1, myThid )
89    
90        RETURN        RETURN
91        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22