/[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.2 by molod, Mon Feb 23 20:33:59 2004 UTC revision 1.3 by molod, Tue Mar 30 18:53:12 2004 UTC
# Line 7  C               grid to the higher resol Line 7  C               grid to the higher resol
7  C Algorithm:  C Algorithm:
8  C   Routine works one layer (edge to edge pressure) at a time.  C   Routine works one layer (edge to edge pressure) at a time.
9  C   Physics -> Dynamics computes the physics layer mean value,  C   Physics -> Dynamics computes the physics layer mean value,
10  C   weighted by dp**kappa  C   weighted by dp**kappa (interp1) or by dp.
11  C  C
12  C Input:  C Input:
13  C   qphy..... [im,jm,lmphy] Arbitrary Quantity on Input Grid  C   qphy..... [im,jm,lmphy] Arbitrary Quantity on Input Grid
# Line 45  C*************************************** Line 45  C***************************************
45        integer Lbot(im1:im2,jm1:jm2,Nsx,Nsy)        integer Lbot(im1:im2,jm1:jm2,Nsx,Nsy)
46    
47        integer  i,j,L,Lout1,Lout1p1,Lout2,Lphy        integer  i,j,L,Lout1,Lout1p1,Lout2,Lphy
48        _RL getcon, kappa, dpkephy, dpkedyn, dpksum, dpsum        _RL getcon, kappa, dpkephy, dpkedyn, dpksum
49    
50        kappa = getcon('KAPPA')        kappa = getcon('KAPPA')
51    
# Line 64  c Check to make sure we are above ground Line 64  c Check to make sure we are above ground
64            endif            endif
65            Lout2 = nlperdyn(i,j,L,bi,bj)            Lout2 = nlperdyn(i,j,L,bi,bj)
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)-  cinterp1  dpkedyn = (pedyn(i,j,L,bi,bj)**kappa)-
68       .                                     (pedyn(i,j,L+1,bi,bj)**kappa)  cinterp1                                   (pedyn(i,j,L+1,bi,bj)**kappa)
69              dpkedyn = pedyn(i,j,L,bi,bj)-pedyn(i,j,L+1,bi,bj)
70            dpksum = 0.            dpksum = 0.
           dpsum = 0.  
71            Lout1p1 = Lout1+1            Lout1p1 = Lout1+1
72            do Lphy = Lout1p1,Lout2            do Lphy = Lout1p1,Lout2
73             dpkephy = (pephy(i,j,Lphy,bi,bj)**kappa)-  cinterp1   dpkephy = (pephy(i,j,Lphy,bi,bj)**kappa)-
74       .                                  (pephy(i,j,Lphy+1,bi,bj)**kappa)  cinterp1                                (pephy(i,j,Lphy+1,bi,bj)**kappa)
75               dpkephy = pephy(i,j,Lphy,bi,bj)-pephy(i,j,Lphy+1,bi,bj)
76             dpksum=dpksum+dpkephy             dpksum=dpksum+dpkephy
            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.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22