/[MITgcm]/MITgcm/pkg/fizhi/fizhi_tendency_apply.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_tendency_apply.F

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

revision 1.8 by molod, Mon Feb 14 22:56:58 2005 UTC revision 1.9 by molod, Fri Feb 18 19:43:04 2005 UTC
# Line 135  C======================================= Line 135  C=======================================
135        if(klev.eq.Nr .or. rC(klev).lt.7000.)then        if(klev.eq.Nr .or. rC(klev).lt.7000.)then
136         cp = getcon('CP')         cp = getcon('CP')
137         kappa = getcon('KAPPA')         kappa = getcon('KAPPA')
138         pNrkappa = (rC(Nr)/100000.)**kappa         pNrkappa = (rC(klev)/100000.)**kappa
139         rayleighdrag = 1./((31.*86400.*2.)*(pNrkappa*cp))         rayleighdrag = 1./((31.*86400.*2.)*(pNrkappa*cp))
140        else        else
141         rayleighdrag = 0.         rayleighdrag = 0.
# Line 145  C======================================= Line 145  C=======================================
145         do i=iMin,iMax         do i=iMin,iMax
146          gT(i,j,kLev,bi,bj) = maskC(i,j,kLev,bi,bj)          gT(i,j,kLev,bi,bj) = maskC(i,j,kLev,bi,bj)
147       .       *( gT(i,j,kLev,bi,bj) + gthphy(i,j,kLev,bi,bj) )       .       *( gT(i,j,kLev,bi,bj) + gthphy(i,j,kLev,bi,bj) )
148       . + rayleighdrag *       . + rayleighdrag * 0.5 *
149       . (maskW(i,j,kLev,bi,bj)*uVel(i,j,kLev,bi,bj)*uVel(i,j,kLev,bi,bj)+       . (maskW(i,j,kLev,bi,bj)*uVel(i,j,kLev,bi,bj)*uVel(i,j,kLev,bi,bj)+
150       .  maskS(i,j,kLev,bi,bj)*vVel(i,j,kLev,bi,bj)*vVel(i,j,kLev,bi,bj))       .  maskW(i+1,j,kLev,bi,bj)*
151         .                uVel(i+1,j,kLev,bi,bj)*uVel(i+1,j,kLev,bi,bj)+
152         .  maskS(i,j,kLev,bi,bj)*vVel(i,j,kLev,bi,bj)*vVel(i,j,kLev,bi,bj)+
153         .  maskS(i,j+1,kLev,bi,bj)*
154         .                vVel(i,j+1,kLev,bi,bj)*vVel(i,j+1,kLev,bi,bj))
155         enddo         enddo
156        enddo        enddo
157    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22