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

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

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

revision 1.8 by molod, Thu Jun 9 00:05:30 2005 UTC revision 1.9 by ce107, Thu Jun 16 16:46:12 2005 UTC
# Line 125  c        and constrain the Maximum Value Line 125  c        and constrain the Maximum Value
125  c ----------------------------------------------  c ----------------------------------------------
126        do j=1,jm        do j=1,jm
127        do i=1,im        do i=1,im
128        phis_std(i,j) = min( 400.0, sqrt( max(0.0,phis_var(i,j)) )/grav )        phis_std(i,j) = min( 400.0 _d 0, sqrt( max(0.0 _d 0,phis_var(i,j)) )/grav )
129        enddo        enddo
130        enddo        enddo
131    
# Line 183  c -------------------------------------- Line 183  c --------------------------------------
183        do L = 1,Lm        do L = 1,Lm
184        do j = 1,jm        do j = 1,jm
185        do i = 1,im        do i = 1,im
186        dragu(i,j,L) = sign( min(0.006,abs(dragu(i,j,L))),dragu(i,j,L) )        dragu(i,j,L) = sign( min(0.006 _d 0,abs(dragu(i,j,L))),dragu(i,j,L) )
187        dragv(i,j,L) = sign( min(0.006,abs(dragv(i,j,L))),dragv(i,j,L) )        dragv(i,j,L) = sign( min(0.006 _d 0,abs(dragv(i,j,L))),dragv(i,j,L) )
188        dragt(i,j,L) = -( uz(i,j,L)*dragu(i,j,L)+vz(i,j,L)*dragv(i,j,L) )        dragt(i,j,L) = -( uz(i,j,L)*dragu(i,j,L)+vz(i,j,L)*dragv(i,j,L) )
189       .                                                         *cpinv       .                                                         *cpinv
190         dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L)         dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L)
# Line 524  c exceeds the Critical Froude number Line 524  c exceeds the Critical Froude number
524  c ----------------------------------  c ----------------------------------
525          crifro = 1.0 - 0.25/richsn          crifro = 1.0 - 0.25/richsn
526          crif2 = crifro*crifro          crif2 = crifro*crifro
527          if( l.eq.2 ) crif2 = min(0.7,crif2)          if( l.eq.2 ) crif2 = min(0.7 _d 0,crif2)
528    
529          if( fro2.gt.crif2 ) then          if( fro2.gt.crif2 ) then
530           tensio(i,L) = crif2/fro2*tensio(i,L-1)           tensio(i,L) = crif2/fro2*tensio(i,L-1)
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22