/[MITgcm]/MITgcm/pkg/gridalt/dyn2phys.F
ViewVC logotype

Diff of /MITgcm/pkg/gridalt/dyn2phys.F

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

revision 1.2 by molod, Tue Mar 30 18:53:12 2004 UTC revision 1.3 by molod, Tue Mar 30 23:44:43 2004 UTC
# Line 50  C*************************************** Line 50  C***************************************
50        integer Lbot(im1:im2,jm1:jm2,Nsx,Nsy)        integer Lbot(im1:im2,jm1:jm2,Nsx,Nsy)
51    
52        _RL weights(im1:im2,jm1:jm2,lmphy)        _RL weights(im1:im2,jm1:jm2,lmphy)
53        _RL dpkedyn, dpkephy, windsum        _RL dpkedyn, dpkephy, windsum, qd
54        integer  i,j,L,Lout1,Lout2,Lphy        integer  i,j,L,Lout1,Lout2,Lphy
55        _RL getcon, kappa        _RL getcon, kappa
56    
57        kappa = getcon('KAPPA')        kappa = getcon('KAPPA')
58    
 c do loop for all grid points  
       do j = jdim1,jdim2  
        do i = idim1,idim2  
59  c do loop for all dynamics (input) levels  c do loop for all dynamics (input) levels
60          do L = 1,lmdyn        do L = 1,lmdyn
61    c do loop for all grid points
62           do j = jdim1,jdim2
63            do i = idim1,idim2
64             qd = qdyn(i,j,L,bi,bj)
65  c Check to make sure we are above ground - if not, do nothing  c Check to make sure we are above ground - if not, do nothing
66           if(L.ge.Lbot(i,j,bi,bj))then           if(L.ge.Lbot(i,j,bi,bj))then
67            if(L.eq.Lbot(i,j,bi,bj)) then            if(L.eq.Lbot(i,j,bi,bj)) then
# Line 87  c do loop for all physics levels contain Line 88  c do loop for all physics levels contain
88             weights(i,j,Lphy) = 1.             weights(i,j,Lphy) = 1.
89             if( (flg.eq.1).and.(windsum.ne.0.) )             if( (flg.eq.1).and.(windsum.ne.0.) )
90       .                weights(i,j,Lphy)=windphy(i,j,Lphy,bi,bj)/windsum       .                weights(i,j,Lphy)=windphy(i,j,Lphy,bi,bj)/windsum
91             qphy(i,j,Lphy,bi,bj) = qdyn(i,j,L,bi,bj) * weights(i,j,Lphy)             qphy(i,j,Lphy,bi,bj) = qd * weights(i,j,Lphy)
92            enddo            enddo
93           endif           endif
94          enddo          enddo
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22