/[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.32 by jmc, Tue Apr 6 00:23:08 2004 UTC revision 1.41 by heimbach, Mon Jan 2 21:17:02 2006 UTC
# Line 35  C     == Local variables == Line 35  C     == Local variables ==
35  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
36  C     I,J  C     I,J
37        INTEGER bi, bj        INTEGER bi, bj
38        INTEGER  I,  J        INTEGER  I, J, K
39  CEOP  CEOP
40    
41        _BARRIER        _BARRIER
# Line 54  CEOP Line 54  CEOP
54            Qsw             (i,j,bi,bj) = 0. _d 0            Qsw             (i,j,bi,bj) = 0. _d 0
55  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
56            pload           (i,j,bi,bj) = 0. _d 0            pload           (i,j,bi,bj) = 0. _d 0
57              sIceLoad        (i,j,bi,bj) = 0. _d 0
58  #endif  #endif
59            surfaceTendencyU(i,j,bi,bj) = 0. _d 0            surfaceForcingU(i,j,bi,bj) = 0. _d 0
60            surfaceTendencyV(i,j,bi,bj) = 0. _d 0            surfaceForcingV(i,j,bi,bj) = 0. _d 0
61            surfaceTendencyT(i,j,bi,bj) = 0. _d 0            surfaceForcingT(i,j,bi,bj) = 0. _d 0
62            surfaceTendencyS(i,j,bi,bj) = 0. _d 0            surfaceForcingS(i,j,bi,bj) = 0. _d 0
63            surfaceTendencyTice(i,j,bi,bj) = 0. _d 0            surfaceForcingTice(i,j,bi,bj) = 0. _d 0
64  #ifndef ALLOW_EXF  #ifndef ALLOW_EXF
65            taux0           (i,j,bi,bj) = 0. _d 0            taux0           (i,j,bi,bj) = 0. _d 0
66            taux1           (i,j,bi,bj) = 0. _d 0            taux1           (i,j,bi,bj) = 0. _d 0
# Line 69  CEOP Line 70  CEOP
70            Qnet1           (i,j,bi,bj) = 0. _d 0            Qnet1           (i,j,bi,bj) = 0. _d 0
71            EmPmR0          (i,j,bi,bj) = 0. _d 0            EmPmR0          (i,j,bi,bj) = 0. _d 0
72            EmPmR1          (i,j,bi,bj) = 0. _d 0            EmPmR1          (i,j,bi,bj) = 0. _d 0
73              saltFlux0       (i,j,bi,bj) = 0. _d 0
74              saltFlux1       (i,j,bi,bj) = 0. _d 0
75            SST0            (i,j,bi,bj) = 0. _d 0            SST0            (i,j,bi,bj) = 0. _d 0
76            SST1            (i,j,bi,bj) = 0. _d 0            SST1            (i,j,bi,bj) = 0. _d 0
77            SSS0            (i,j,bi,bj) = 0. _d 0            SSS0            (i,j,bi,bj) = 0. _d 0
# Line 87  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          _BARRIER
115        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
116        IF ( zonalWindFile .NE. ' '  ) THEN        IF ( zonalWindFile .NE. ' '  ) THEN
117         CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )         CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
# Line 96  C Line 121  C
121        ENDIF        ENDIF
122        IF ( surfQFile .NE. ' '  ) THEN        IF ( surfQFile .NE. ' '  ) THEN
123         CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )         CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )
124          ELSEIF ( surfQnetFile .NE. ' '  ) THEN
125           CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid )
126        ENDIF        ENDIF
127        IF ( EmPmRfile .NE. ' '  ) THEN        IF ( EmPmRfile .NE. ' '  ) THEN
128         CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )         CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
129        ENDIF        ENDIF
130          IF ( saltFluxFile .NE. ' '  ) THEN
131           CALL READ_FLD_XY_RS( saltFluxFile, ' ', saltFlux, 0, myThid )
132          ENDIF
133        IF ( thetaClimFile .NE. ' '  ) THEN        IF ( thetaClimFile .NE. ' '  ) THEN
134         CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )         CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )
135        ENDIF        ENDIF
136        IF ( saltClimFile .NE. ' '  ) THEN        IF ( saltClimFile .NE. ' '  ) THEN
137         CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )         CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
138        ENDIF        ENDIF
139          IF ( lambdaThetaFile .NE. ' '  ) THEN
140           CALL READ_FLD_XY_RS( lambdaThetaFile, ' ',
141         &  lambdaThetaClimRelax, 0, myThid )
142          ENDIF
143          IF ( lambdaSaltFile .NE. ' '  ) THEN
144           CALL READ_FLD_XY_RS( lambdaSaltFile, ' ',
145         &  lambdaSaltClimRelax, 0, myThid )
146          ENDIF
147  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
148        IF ( surfQswFile .NE. ' '  ) THEN        IF ( surfQswFile .NE. ' '  ) THEN
149         CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )         CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
150           IF ( surfQFile .NE. ' '  ) THEN
151    C-     Qnet is now (after c54) the net Heat Flux (including SW)
152            DO bj = 1,nSy
153             DO bi = 1,nSx
154              DO j=1-OLy,sNy+OLy
155               DO i=1-OLx,sNx+OLx
156                Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
157               ENDDO
158              ENDDO
159             ENDDO
160            ENDDO
161           ENDIF
162        ENDIF        ENDIF
163  #endif  #endif
164  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
# Line 116  C Line 166  C
166         CALL READ_FLD_XY_RS( pLoadFile, ' ', pload, 0, myThid )         CALL READ_FLD_XY_RS( pLoadFile, ' ', pload, 0, myThid )
167        ENDIF        ENDIF
168  #endif  #endif
169    
170        _END_MASTER(myThid)        _END_MASTER(myThid)
171  C  C
172        _EXCH_XY_R4(fu   , myThid )        _EXCH_XY_R4(fu   , myThid )
173        _EXCH_XY_R4(fv   , myThid )        _EXCH_XY_R4(fv   , myThid )
174        _EXCH_XY_R4(Qnet , myThid )        _EXCH_XY_R4(Qnet , myThid )
175        _EXCH_XY_R4(EmPmR, myThid )        _EXCH_XY_R4(EmPmR, myThid )
176          _EXCH_XY_R4( saltFlux, myThid )
177        _EXCH_XY_R4(SST  , myThid )        _EXCH_XY_R4(SST  , myThid )
178        _EXCH_XY_R4(SSS  , myThid )        _EXCH_XY_R4(SSS  , myThid )
179          _EXCH_XY_R4(lambdaThetaClimRelax , myThid )
180          _EXCH_XY_R4(lambdaSaltClimRelax , myThid )
181  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
182        _EXCH_XY_R4(Qsw  , myThid )        _EXCH_XY_R4(Qsw  , myThid )
183  #endif  #endif
# Line 131  C Line 185  C
185        _EXCH_XY_R4(pload  , myThid )        _EXCH_XY_R4(pload  , myThid )
186  C     CALL PLOT_FIELD_XYRS( pload, 'S/R INI_FORCING pload',1,myThid)  C     CALL PLOT_FIELD_XYRS( pload, 'S/R INI_FORCING pload',1,myThid)
187  #endif  #endif
   
188  C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)  C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
189  C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)  C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
190    
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22