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

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

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

revision 1.11 by heimbach, Tue Apr 10 22:35:25 2001 UTC revision 1.12 by adcroft, Tue May 29 14:01:37 2001 UTC
# Line 1  Line 1 
1  C     $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
# Line 105  C     Add windstress momentum impulse in Line 105  C     Add windstress momentum impulse in
105  CStartOfInterface  CStartOfInterface
106        SUBROUTINE EXTERNAL_FORCING_T(        SUBROUTINE EXTERNAL_FORCING_T(
107       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin, iMax, jMin, jMax,bi,bj,kLev,
      I           maskC,  
108       I           myCurrentTime,myThid)       I           myCurrentTime,myThid)
109  C     /==========================================================\  C     /==========================================================\
110  C     | S/R EXTERNAL_FORCING_T                                   |  C     | S/R EXTERNAL_FORCING_T                                   |
# Line 137  C     iMax Line 136  C     iMax
136  C     jMin  C     jMin
137  C     jMax  C     jMax
138  C     kLev  C     kLev
       _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
139        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
140        _RL myCurrentTime        _RL myCurrentTime
141        INTEGER myThid        INTEGER myThid
# Line 153  C     Add heat in top-layer Line 151  C     Add heat in top-layer
151         DO j=jMin,jMax         DO j=jMin,jMax
152          DO i=iMin,iMax          DO i=iMin,iMax
153           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
154       &     +maskC(i,j)*surfaceTendencyT(i,j,bi,bj)       &     +maskC(i,j,kLev,bi,bj)*surfaceTendencyT(i,j,bi,bj)
155          ENDDO          ENDDO
156         ENDDO         ENDDO
157        ENDIF        ENDIF
# Line 168  C Penetrating SW radiation Line 166  C Penetrating SW radiation
166       O     swfracb)       O     swfracb)
167        DO j=jMin,jMax        DO j=jMin,jMax
168         DO i=iMin,iMax         DO i=iMin,iMax
169          gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj)          gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj)
170       &    -maskC(i,j)*Qsw(i,j,bi,bj)*(swfracb(1)-swfracb(2))       &   -maskC(i,j,klev,bi,bj)*Qsw(i,j,bi,bj)*(swfracb(1)-swfracb(2))
171       &    *recip_Cp*recip_rhoNil*recip_dRf(klev)       &    *recip_Cp*recip_rhoNil*recip_drF(klev)
172         ENDDO         ENDDO
173        ENDDO        ENDDO
174  #endif  #endif
# Line 179  C Penetrating SW radiation Line 177  C Penetrating SW radiation
177  CStartOfInterface  CStartOfInterface
178        SUBROUTINE EXTERNAL_FORCING_S(        SUBROUTINE EXTERNAL_FORCING_S(
179       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin, iMax, jMin, jMax,bi,bj,kLev,
      I           maskC,  
180       I           myCurrentTime,myThid)       I           myCurrentTime,myThid)
181  C     /==========================================================\  C     /==========================================================\
182  C     | S/R EXTERNAL_FORCING_S                                   |  C     | S/R EXTERNAL_FORCING_S                                   |
# Line 204  C     iMax Line 201  C     iMax
201  C     jMin  C     jMin
202  C     jMax  C     jMax
203  C     kLev  C     kLev
       _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
204        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
205        _RL myCurrentTime        _RL myCurrentTime
206        INTEGER myThid        INTEGER myThid
# Line 220  C     Add fresh-water in top-layer Line 216  C     Add fresh-water in top-layer
216         DO j=jMin,jMax         DO j=jMin,jMax
217          DO i=iMin,iMax          DO i=iMin,iMax
218           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)
219       &   +maskC(i,j)*surfaceTendencyS(i,j,bi,bj)       &   +maskC(i,j,kLev,bi,bj)*surfaceTendencyS(i,j,bi,bj)
220          ENDDO          ENDDO
221         ENDDO         ENDDO
222        ENDIF        ENDIF
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22