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

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

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

revision 1.37 by jmc, Tue Jan 4 15:57:22 2005 UTC revision 1.38 by heimbach, Wed Apr 6 20:18:19 2005 UTC
# Line 90  CEOP Line 90  CEOP
90         ENDDO         ENDDO
91        ENDDO        ENDDO
92  C  C
93          DO bj = myByLo(myThid), myByHi(myThid)
94           DO bi = myBxLo(myThid), myBxHi(myThid)
95            DO J=1-Oly,sNy+Oly
96             DO I=1-Olx,sNx+Olx
97              IF ( doThetaClimRelax .AND.
98         &         abs(yC(i,j,bi,bj)).LE.latBandClimRelax ) THEN
99               lambdaThetaClimRelax(I,J,bi,bj) = 1./tauThetaClimRelax
100              ELSE
101               lambdaThetaClimRelax(I,J,bi,bj) = 0.D0
102              ENDIF
103              IF ( doSaltClimRelax .AND.
104         &         abs(yC(i,j,bi,bj)).LE.latBandClimRelax ) THEN
105               lambdaSaltClimRelax(I,J,bi,bj) = 1./tauSaltClimRelax
106              ELSE
107               lambdaSaltClimRelax(I,J,bi,bj) = 0.D0
108              ENDIF
109             ENDDO
110            ENDDO
111           ENDDO
112          ENDDO
113    C
114        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
115        IF ( zonalWindFile .NE. ' '  ) THEN        IF ( zonalWindFile .NE. ' '  ) THEN
116         CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )         CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
# Line 114  C Line 135  C
135        IF ( saltClimFile .NE. ' '  ) THEN        IF ( saltClimFile .NE. ' '  ) THEN
136         CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )         CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
137        ENDIF        ENDIF
138          IF ( lambdaThetaFile .NE. ' '  ) THEN
139           CALL READ_FLD_XY_RS( lambdaThetaFile, ' ',
140         &  lambdaThetaClimRelax, 0, myThid )
141          ENDIF
142          IF ( lambdaSaltFile .NE. ' '  ) THEN
143           CALL READ_FLD_XY_RS( lambdaSaltFile, ' ',
144         &  lambdaSaltClimRelax, 0, myThid )
145          ENDIF
146  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
147        IF ( surfQswFile .NE. ' '  ) THEN        IF ( surfQswFile .NE. ' '  ) THEN
148         CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )         CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
# Line 145  C Line 174  C
174        _EXCH_XY_R4( saltFlux, myThid )        _EXCH_XY_R4( saltFlux, myThid )
175        _EXCH_XY_R4(SST  , myThid )        _EXCH_XY_R4(SST  , myThid )
176        _EXCH_XY_R4(SSS  , myThid )        _EXCH_XY_R4(SSS  , myThid )
177          _EXCH_XY_R4(lambdaThetaClimRelax , myThid )
178          _EXCH_XY_R4(lambdaSaltClimRelax , myThid )
179  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
180        _EXCH_XY_R4(Qsw  , myThid )        _EXCH_XY_R4(Qsw  , myThid )
181  #endif  #endif
Legend:
Removed from v.1.37  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22