/[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.1 by molod, Thu Jan 29 14:23:09 2004 UTC revision 1.2 by molod, Tue Mar 30 18:53:12 2004 UTC
# Line 71  c Check to make sure we are above ground Line 71  c Check to make sure we are above ground
71            Lout2 = nlperdyn(i,j,L,bi,bj)            Lout2 = nlperdyn(i,j,L,bi,bj)
72  c for U and V fields, need to compute for the weights:  c for U and V fields, need to compute for the weights:
73            if(flg.eq.1)then            if(flg.eq.1)then
74             dpkedyn = (pedyn(i,j,L,bi,bj)**kappa)-  cinterp1   dpkedyn = (pedyn(i,j,L,bi,bj)**kappa)-
75       .                                     (pedyn(i,j,L+1,bi,bj)**kappa)  cinterp1                                   (pedyn(i,j,L+1,bi,bj)**kappa)
76               dpkedyn = pedyn(i,j,L,bi,bj)-pedyn(i,j,L+1,bi,bj)
77             windsum = 0.             windsum = 0.
78             do Lphy = Lout1+1,Lout2             do Lphy = Lout1+1,Lout2
79              dpkephy = (pephy(i,j,Lphy,bi,bj)**kappa)-  cinterp1    dpkephy = (pephy(i,j,Lphy,bi,bj)**kappa)-
80       .                                  (pephy(i,j,Lphy+1,bi,bj)**kappa)  cinterp1                                (pephy(i,j,Lphy+1,bi,bj)**kappa)
81                dpkephy = pephy(i,j,Lphy,bi,bj)-pephy(i,j,Lphy+1,bi,bj)
82              windsum = windsum+(windphy(i,j,Lphy,bi,bj)*dpkephy)/dpkedyn              windsum = windsum+(windphy(i,j,Lphy,bi,bj)*dpkephy)/dpkedyn
83             enddo             enddo
84            endif            endif
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22