/[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.2 by cnh, Fri Nov 6 22:44:46 1998 UTC revision 1.3 by adcroft, Tue Dec 15 00:20:34 1998 UTC
# Line 43  C     Add windstress momentum impulse in Line 43  C     Add windstress momentum impulse in
43         DO j=jMin,jMax         DO j=jMin,jMax
44          DO i=iMin,iMax          DO i=iMin,iMax
45           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)
46       &    + foFacMom*fu(i,j,bi,bj)*_maskW(i,j,kLev,bi,bj)       &   +foFacMom*fu(i,j,bi,bj)
47         &   *horiVertRatio*recip_rhoNil*recip_dRf(kLev)
48         &   *_maskW(i,j,kLev,bi,bj)
49          ENDDO          ENDDO
50         ENDDO         ENDDO
51        ENDIF        ENDIF
# Line 91  C     Add windstress momentum impulse in Line 93  C     Add windstress momentum impulse in
93         DO j=jMin,jMax         DO j=jMin,jMax
94          DO i=iMin,iMax          DO i=iMin,iMax
95           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)
96       &    + foFacMom*fv(i,j,bi,bj)*_maskS(i,j,kLev,bi,bj)       &   +foFacMom*fv(i,j,bi,bj)
97         &   *horiVertRatio*recip_rhoNil*recip_dRf(kLev)
98         &   *_maskS(i,j,kLev,bi,bj)
99          ENDDO          ENDDO
100         ENDDO         ENDDO
101        ENDIF        ENDIF
# Line 143  C     Add heat in top-layer Line 147  C     Add heat in top-layer
147           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)           gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
148       &  +maskC(i,j)*(       &  +maskC(i,j)*(
149       &   -lambdaThetaClimRelax*(theta(i,j,kLev,bi,bj)-SST(i,j,bi,bj))       &   -lambdaThetaClimRelax*(theta(i,j,kLev,bi,bj)-SST(i,j,bi,bj))
150       &   -Qnet(i,j,bi,bj) )       &   -Qnet(i,j,bi,bj)*recip_Cp*recip_rhoNil*recip_dRf(kLev) )
151          ENDDO          ENDDO
152         ENDDO         ENDDO
153        ENDIF        ENDIF
# Line 195  C     Add fresh-water in top-layer Line 199  C     Add fresh-water in top-layer
199           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)           gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)
200       &   +maskC(i,j)*(       &   +maskC(i,j)*(
201       &   -lambdaSaltClimRelax*(salt(i,j,kLev,bi,bj)-SSS(i,j,bi,bj))       &   -lambdaSaltClimRelax*(salt(i,j,kLev,bi,bj)-SSS(i,j,bi,bj))
202       &   +EmPmR(i,j,bi,bj) )  #ifndef USE_NATURAL_BCS
203         &   +EmPmR(i,j,bi,bj)*recip_dRf(1)*35.
204    #endif
205         &   )
206          ENDDO          ENDDO
207         ENDDO         ENDDO
208        ENDIF        ENDIF
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22