/[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.9 by ce107, Thu Jun 16 16:46:12 2005 UTC revision 1.10 by ce107, Fri Jun 17 16:51:24 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 _d 0, sqrt( max(0.0 _d 0,phis_var(i,j)) )/grav )           phis_std(i,j) = min( 400.0 _d 0, sqrt( max(0.0 _d 0,
129         $        phis_var(i,j)) )/grav )
130        enddo        enddo
131        enddo        enddo
132    
# Line 183  c -------------------------------------- Line 184  c --------------------------------------
184        do L = 1,Lm        do L = 1,Lm
185        do j = 1,jm        do j = 1,jm
186        do i = 1,im        do i = 1,im
187        dragu(i,j,L) = sign( min(0.006 _d 0,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
188        dragv(i,j,L) = sign( min(0.006 _d 0,abs(dragv(i,j,L))),dragv(i,j,L) )       $        ,j,L) )
189             dragv(i,j,L) = sign( min(0.006 _d 0,abs(dragv(i,j,L))), dragv(i
190         $        ,j,L) )
191        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) )
192       .                                                         *cpinv       .                                                         *cpinv
193         dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L)         dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22