/[MITgcm]/MITgcm/verification/hs94.1x64x5/code/external_forcing.F
ViewVC logotype

Diff of /MITgcm/verification/hs94.1x64x5/code/external_forcing.F

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

revision 1.5 by adcroft, Mon Jun 4 20:20:18 2001 UTC revision 1.6 by adcroft, Thu Jun 7 18:01:19 2001 UTC
# Line 49  C     _RL velDragHeightFac Line 49  C     _RL velDragHeightFac
49        _RL termP,kV,kF        _RL termP,kV,kF
50    
51  C--   Forcing term(s)  C--   Forcing term(s)
52        kF=1./86400.        kF=1. _d 0/86400. _d 0
53        DO J=jMin,jMax        DO J=jMin,jMax
54         DO I=iMin,iMax         DO I=iMin,iMax
55          IF ( HFacW(i,j,kLev,bi,bj) .GT. 0. ) THEN          IF ( HFacW(i,j,kLev,bi,bj) .GT. 0. ) THEN
56  C        termP=0.5*( rF(kLev) + min( rF(kLev+1) ,  C        termP=0.5*( rF(kLev) + min( rF(kLev+1) ,
57  C    &           min(H(I,J,bi,bj),H(I,J-1,bi,bj))            ) )  C    &           min(H(I,J,bi,bj),H(I,J-1,bi,bj))            ) )
58           termP=0.5*( rF(kLev) + rF(kLev+1) )           termP=0.5 _d 0*( rF(kLev) + rF(kLev+1) )
59           kV=kF*MAX(0. _d 0,           kV=kF*MAX(0. _d 0,
60       &  (termP*recip_Rcol(I,J,bi,bj)-0.7 _d 0)/(1. _d 0-0.7 _d 0) )       &  (termP*recip_Rcol(I,J,bi,bj)-0.7 _d 0)/(1. _d 0-0.7 _d 0) )
61           gU(i,j,kLev,bi,bj)=gU(i,j,kLev,bi,bj)           gU(i,j,kLev,bi,bj)=gU(i,j,kLev,bi,bj)
# Line 112  C     _RL velDragHeightFac Line 112  C     _RL velDragHeightFac
112        _RL termP,kV,kF        _RL termP,kV,kF
113    
114  C--   Forcing term(s)  C--   Forcing term(s)
115        kF=1./86400.        kF=1. _d 0/86400. _d 0
116        DO J=jMin,jMax        DO J=jMin,jMax
117         DO I=iMin,iMax         DO I=iMin,iMax
118          IF ( HFacS(i,j,kLev,bi,bj) .GT. 0. ) THEN          IF ( HFacS(i,j,kLev,bi,bj) .GT. 0. ) THEN
119  C        termP=0.5*( rF(kLev) + min( rF(kLev+1) ,  C        termP=0.5*( rF(kLev) + min( rF(kLev+1) ,
120  C    &           min(H(I,J,bi,bj),H(I,J-1,bi,bj))            ) )  C    &           min(H(I,J,bi,bj),H(I,J-1,bi,bj))            ) )
121           termP=0.5*( rF(kLev) + rF(kLev+1) )           termP=0.5 _d 0*( rF(kLev) + rF(kLev+1) )
122           kV=kF*MAX(0. _d 0,           kV=kF*MAX(0. _d 0,
123       &  (termP*recip_Rcol(I,J,bi,bj)-0.7 _d 0)/(1. _d 0-0.7 _d 0) )       &  (termP*recip_Rcol(I,J,bi,bj)-0.7 _d 0)/(1. _d 0-0.7 _d 0) )
124           gV(i,j,kLev,bi,bj)=gV(i,j,kLev,bi,bj)           gV(i,j,kLev,bi,bj)=gV(i,j,kLev,bi,bj)
# Line 177  C      termP=0.5*( rF(kLev) + min( rF(kL Line 177  C      termP=0.5*( rF(kLev) + min( rF(kL
177         term2=10. _d 0*log(termP/rSurf)         term2=10. _d 0*log(termP/rSurf)
178       &          *(cos(yC(1,J,bi,bj)*deg2rad)**2)       &          *(cos(yC(1,J,bi,bj)*deg2rad)**2)
179         thetaLim = 200. _d 0/ ((termP/rSurf)**(2. _d 0/7. _d 0))         thetaLim = 200. _d 0/ ((termP/rSurf)**(2. _d 0/7. _d 0))
180         thetaEq=315.-term1-term2         thetaEq=315. _d 0-term1-term2
181         thetaEq=MAX(thetaLim,thetaEq)         thetaEq=MAX(thetaLim,thetaEq)
182         DO I=iMin,iMax         DO I=iMin,iMax
183          kT=ka+(ks-ka)          kT=ka+(ks-ka)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22