/[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.12 by cnh, Sun Feb 4 14:38:47 2001 UTC revision 1.13 by adcroft, Tue May 29 14:01:37 2001 UTC
# Line 37  C     == Local variables == Line 37  C     == Local variables ==
37  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
38  C     I,J,K  C     I,J,K
39        INTEGER bi, bj        INTEGER bi, bj
40        INTEGER  I,  J, K        INTEGER  I,  J, K, localWarnings
41          CHARACTER*(MAX_LEN_MBUF) msgBuf
42    
43        _BARRIER  C--   Initialise salinity field to the vertical reference profile
44          DO bj = myByLo(myThid), myByHi(myThid)
45        IF ( hydrogSaltFile .EQ. ' ' ) THEN         DO bi = myBxLo(myThid), myBxHi(myThid)
46  C--    Initialise salinity field to the vertical reference profile          DO K=1,Nr
47         DO bj = myByLo(myThid), myByHi(myThid)           DO J=1-Oly,sNy+Oly
48          DO bi = myBxLo(myThid), myBxHi(myThid)            DO I=1-Olx,sNx+Olx
49           DO K=1,Nr             salt(I,J,K,bi,bj) = sRef(K)
           DO J=1,sNy  
            DO I=1,sNx  
             salt(I,J,K,bi,bj) = sRef(K)  
            ENDDO  
50            ENDDO            ENDDO
51           ENDDO           ENDDO
52          ENDDO          ENDDO
53         ENDDO         ENDDO
54        ELSE        ENDDO
55    
56          IF ( hydrogSaltFile .NE. ' ' ) THEN
57         _BEGIN_MASTER( myThid )         _BEGIN_MASTER( myThid )
58         CALL READ_FLD_XYZ_RL( hydrogSaltFile, ' ', salt, 0, myThid )         CALL READ_FLD_XYZ_RL( hydrogSaltFile, ' ', salt, 0, myThid )
59         _END_MASTER(myThid)         _END_MASTER(myThid)
60           _EXCH_XYZ_R8(salt  , myThid )
61        ENDIF        ENDIF
62  C     Set initial tendency terms  
63    C     Apply mask and test consistancy
64          localWarnings=0
65        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
66         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
67          DO K=1,Nr          DO K=1,Nr
68           DO J=1,sNy           DO J=1,sNy
69            DO I=1,sNx            DO I=1,sNx
            gs   (I,J,K,bi,bj) = 0. _d 0  
            gsNM1(I,J,K,bi,bj) = 0. _d 0  
70             IF (hFacC(I,J,K,bi,bj).EQ.0) salt(I,J,K,bi,bj) = 0.             IF (hFacC(I,J,K,bi,bj).EQ.0) salt(I,J,K,bi,bj) = 0.
71               IF (hFacC(I,J,K,bi,bj).NE.0.AND.salt(I,J,K,bi,bj).EQ.0.
72         &      .AND. sRef(k).NE.0.) THEN
73                 localWarnings=localWarnings+1
74                ENDIF
75            ENDDO            ENDDO
76           ENDDO           ENDDO
77          ENDDO          ENDDO
78         ENDDO         ENDDO
79        ENDDO        ENDDO
80  C        IF (localWarnings.NE.0) THEN
81        _EXCH_XYZ_R8(salt  , myThid )         WRITE(msgBuf,'(A,A)')
82        _EXCH_XYZ_R8(gs , myThid )       &  'S/R INI_SALT: salt = 0 identically. If this is intentional',
83        _EXCH_XYZ_R8(gsNM1 , myThid )       &  'you will need to edit ini_salt.F to avoid this safety check'
84           CALL PRINT_ERROR( msgBuf , myThid)
85           STOP 'ABNORMAL END: S/R INI_SALT'
86          ENDIF
87    
88        CALL PLOT_FIELD_XYZRL( salt, 'Initial Salinity' , Nr, 1, myThid )        CALL PLOT_FIELD_XYZRL( salt, 'Initial Salinity' , Nr, 1, myThid )
89    
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22