/[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.2 by cnh, Fri Apr 24 02:05:41 1998 UTC revision 1.3 by cnh, Tue Jun 9 18:57:38 1998 UTC
# Line 45  C     I,J,K Line 45  C     I,J,K
45        INTEGER  I,  J, K        INTEGER  I,  J, K
46        REAL rad, rD        REAL rad, rD
47    
48  C--   Initialise salinity field to the vertical reference profile        _BARRIER
49        DO bj = myByLo(myThid), myByHi(myThid)  
50         DO bi = myBxLo(myThid), myBxHi(myThid)        IF ( hydrogSaltFile .EQ. ' ' ) THEN
51          DO K=1,Nz  C--    Initialise salinity field to the vertical reference profile
52           DO J=1,sNy         DO bj = myByLo(myThid), myByHi(myThid)
53            DO I=1,sNx          DO bi = myBxLo(myThid), myBxHi(myThid)
54             salt(I,J,K,bi,bj) = sRef(K)           DO K=1,Nz
55              DO J=1,sNy
56               DO I=1,sNx
57                salt(I,J,K,bi,bj) = sRef(K)
58               ENDDO
59            ENDDO            ENDDO
60           ENDDO           ENDDO
61          ENDDO          ENDDO
62         ENDDO         ENDDO
63        ENDDO        ELSE
64           _BEGIN_MASTER( myThid )
65           CALL READ_FLD_XYZ_RS( hydrogSaltFile, ' ', salt, 0, myThid )
66           _END_MASTER(myThid)
67          ENDIF
68  C     Set initial tendency terms  C     Set initial tendency terms
69        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
70         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22