/[MITgcm]/MITgcm/verification/hs94.cs-32x32x5/code/external_forcing.F
ViewVC logotype

Diff of /MITgcm/verification/hs94.cs-32x32x5/code/external_forcing.F

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

revision 1.3 by adcroft, Wed Jun 6 19:46:43 2001 UTC revision 1.4 by adcroft, Wed Jun 6 20:31:48 2001 UTC
# Line 177  C      termP=0.5*( rF(kLev) + min( rF(kL Line 177  C      termP=0.5*( rF(kLev) + min( rF(kL
177         termP=0.5 _d 0*( rF(kLev) + rF(kLev+1) )         termP=0.5 _d 0*( rF(kLev) + rF(kLev+1) )
178         term2=10. _d 0*log(termP/rSurf)         term2=10. _d 0*log(termP/rSurf)
179       &          *(cos(yC(I,J,bi,bj)*deg2rad)**2)       &          *(cos(yC(I,J,bi,bj)*deg2rad)**2)
180         thetaLim = 200. _d 0/ ((termP/rSurf)**(2./7.))         thetaLim = 200. _d 0/ ((termP/rSurf)**(2. _d 0/7. _d 0))
181         thetaEq=315. _d 0-term1-term2         thetaEq=315. _d 0-term1-term2
182         thetaEq=MAX(thetaLim,thetaEq)         thetaEq=MAX(thetaLim,thetaEq)
183          kT=ka+(ks-ka)          kT=ka+(ks-ka)

Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4

  ViewVC Help
Powered by ViewVC 1.1.22