/[MITgcm]/MITgcm_contrib/plumes/plume2dyn.F
ViewVC logotype

Diff of /MITgcm_contrib/plumes/plume2dyn.F

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

revision 1.2 by molod, Tue May 25 18:08:51 2004 UTC revision 1.3 by molod, Tue May 25 23:12:50 2004 UTC
# Line 76  C the angle given by uref and vref Line 76  C the angle given by uref and vref
76         do j = jdim1,jdim2         do j = jdim1,jdim2
77         do i = idim1,idim2         do i = idim1,idim2
78         do L = 1,Lmplume         do L = 1,Lmplume
79          sqrtarg = (qplumeav(i,j,L)*qplumeav(i,j,L)) /          if(vref(i,j,L,bi,bj).ne.0.) then
80             sqrtarg = (qplumeav(i,j,L)*qplumeav(i,j,L)) /
81       .  ( ( (uref(i,j,L,bi,bj)*uref(i,j,L,bi,bj)) /       .  ( ( (uref(i,j,L,bi,bj)*uref(i,j,L,bi,bj)) /
82       .      (vref(i,j,L,bi,bj)*vref(i,j,L,bi,bj)) ) + 1. )       .      (vref(i,j,L,bi,bj)*vref(i,j,L,bi,bj)) ) + 1. )
83          qdyn2(i,j,L,bi,bj) = sqrt(sqrtarg)           qdyn2(i,j,L,bi,bj) = sqrt(sqrtarg)
84          qdyn1(i,j,L,bi,bj) = qdyn2(i,j,L,bi,bj) *           qdyn1(i,j,L,bi,bj) = qdyn2(i,j,L,bi,bj) *
85       .                            (uref(i,j,L,bi,bj)/vref(i,j,L,bi,bj))       .                            (uref(i,j,L,bi,bj)/vref(i,j,L,bi,bj))
86            else
87             qdyn1(i,j,L,bi,bj) = qplumeav(i,j,L)
88             qdyn2(i,j,L,bi,bj) = 0.
89         enddo         enddo
90         enddo         enddo
91         enddo         enddo
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22