/[lmdze]/trunk/Sources/dyn3d/PPM3d/ppm3d.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/PPM3d/ppm3d.f

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

revision 168 by guez, Wed Sep 9 10:41:47 2015 UTC revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC
# Line 7  contains Line 7  contains
7    SUBROUTINE ppm3d(igd, q, ps1, ps2, u, v, w, ndt, iord, jord, kord, nc, imr, &    SUBROUTINE ppm3d(igd, q, ps1, ps2, u, v, w, ndt, iord, jord, kord, nc, imr, &
8         jnp, j1, nlay, ap, bp, pt, ae, fill, umax)         jnp, j1, nlay, ap, bp, pt, ae, fill, umax)
9    
     ! rajout de d\'eclarations  
     ! integer Jmax,kmax,ndt0,nstep,k,j,i,ic,l,js,jn,imh,iad,jad,krd  
     ! integer iu,iiu,j2,jmr,js0,jt  
     ! real dtdy,dtdy5,rcap,iml,jn0,imjm,pi,dl,dp  
     ! real dt,cr1,maxdt,ztc,d5,sum1,sum2,ru  
   
     ! ********************************************************************  
   
     ! =============  
10      ! INPUT:      ! INPUT:
11      ! =============      ! =============
12    
# Line 52  contains Line 43  contains
43      ! currently under development.      ! currently under development.
44      ! p(i,j,k) = (AP(k)*PT + BP(k)*PS(i,j))/(D(k)-C(k)*TE**(-1/kapa))      ! p(i,j,k) = (AP(k)*PT + BP(k)*PS(i,j))/(D(k)-C(k)*TE**(-1/kapa))
45    
46      ! /////////////////////////////////      ! Cf. ppm3d.txt.
     ! / \ ------------- PTOP --------------  AP(1), BP(1)  
     ! |  
     ! delp(1)    |  ........... Q(i,j,1) ............  
     ! |  
     ! W(1)    \ / ---------------------------------  AP(2), BP(2)  
   
   
   
     ! W(k-1)   / \ ---------------------------------  AP(k), BP(k)  
     ! |  
     ! delp(K)    |  ........... Q(i,j,k) ............  
     ! |  
     ! W(k)    \ / ---------------------------------  AP(k+1), BP(k+1)  
   
   
   
     ! / \ ---------------------------------  AP(NLAY), BP(NLAY)  
     ! |  
     ! delp(NLAY)   |  ........... Q(i,j,NLAY) .........  
     ! |  
     ! W(NLAY)=0  \ / ------------- surface ----------- AP(NLAY+1), BP(NLAY+1)  
     ! //////////////////////////////////  
47    
48      ! U(IMR,JNP,NLAY) & V(IMR,JNP,NLAY):winds (m/s) at mid-time-level (t+NDT/2)      ! U(IMR,JNP,NLAY) & V(IMR,JNP,NLAY):winds (m/s) at mid-time-level (t+NDT/2)
49      ! U and V may need to be polar filtered in advance in some cases.      ! U and V may need to be polar filtered in advance in some cases.
# Line 84  contains Line 53  contains
53      ! pole (j=1) to north pole (j=JNP) ]      ! pole (j=1) to north pole (j=JNP) ]
54    
55      ! IGD = 1  GEOS-GCM C-Grid      ! IGD = 1  GEOS-GCM C-Grid
56      ! [North]      ! Cf. ppm3d.txt.
   
     ! V(i,j)  
     ! |  
     ! |  
     ! |  
     ! U(i-1,j)---Q(i,j)---U(i,j) [EAST]  
     ! |  
     ! |  
     ! |  
     ! V(i,j-1)  
57    
58      ! U(i,  1) is defined at South Pole.      ! U(i,  1) is defined at South Pole.
59      ! V(i,  1) is half grid north of the South Pole.      ! V(i,  1) is half grid north of the South Pole.
Legend:
Removed from v.168  
changed lines
  Added in v.169
  ViewVC Help
Powered by ViewVC 1.1.21