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

Diff of /MITgcm/pkg/gridalt/phys2dyn.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, Mon Feb 23 20:33:59 2004 UTC
# Line 44  C*************************************** Line 44  C***************************************
44        _RL qdyn(im1:im2,jm1:jm2,lmdyn,Nsx,Nsy)        _RL qdyn(im1:im2,jm1:jm2,lmdyn,Nsx,Nsy)
45        integer Lbot(im1:im2,jm1:jm2,Nsx,Nsy)        integer Lbot(im1:im2,jm1:jm2,Nsx,Nsy)
46    
47        integer  i,j,L,Lout1,Lout2,Lphy        integer  i,j,L,Lout1,Lout1p1,Lout2,Lphy
48        _RL getcon, kappa, dpkephy, dpkedyn        _RL getcon, kappa, dpkephy, dpkedyn, dpksum, dpsum
49    
50        kappa = getcon('KAPPA')        kappa = getcon('KAPPA')
51    
# Line 66  c Check to make sure we are above ground Line 66  c Check to make sure we are above ground
66  c do loop for all physics levels contained in this dynamics level  c do loop for all physics levels contained in this dynamics level
67            dpkedyn = (pedyn(i,j,L,bi,bj)**kappa)-            dpkedyn = (pedyn(i,j,L,bi,bj)**kappa)-
68       .                                     (pedyn(i,j,L+1,bi,bj)**kappa)       .                                     (pedyn(i,j,L+1,bi,bj)**kappa)
69            do Lphy = Lout1+1,Lout2            dpksum = 0.
70              dpsum = 0.
71              Lout1p1 = Lout1+1
72              do Lphy = Lout1p1,Lout2
73             dpkephy = (pephy(i,j,Lphy,bi,bj)**kappa)-             dpkephy = (pephy(i,j,Lphy,bi,bj)**kappa)-
74       .                                  (pephy(i,j,Lphy+1,bi,bj)**kappa)       .                                  (pephy(i,j,Lphy+1,bi,bj)**kappa)
75               dpksum=dpksum+dpkephy
76               dpsum=dpsum+(pephy(i,j,Lphy,bi,bj)-pephy(i,j,Lphy+1,bi,bj))
77             qdyn(i,j,L,bi,bj) = qdyn(i,j,L,bi,bj) +             qdyn(i,j,L,bi,bj) = qdyn(i,j,L,bi,bj) +
78       .                            qphy(i,j,Lphy,bi,bj)*(dpkephy/dpkedyn)       .                            qphy(i,j,Lphy,bi,bj)*(dpkephy/dpkedyn)
79            enddo            enddo
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22