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

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

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 10 by guez, Fri Apr 18 14:45:53 2008 UTC
# Line 80  C INPUT: Line 80  C INPUT:
80  C =============  C =============
81  C  C
82  C Q(IMR,JNP,NLAY,NC): mixing ratios at current time (t)  C Q(IMR,JNP,NLAY,NC): mixing ratios at current time (t)
83  C NC: total # of constituents  C NC: total number of constituents
84  C IMR: first dimension (E-W); # of Grid intervals in E-W is IMR  C IMR: first dimension (E-W); number of Grid intervals in E-W is IMR
85  C JNP: 2nd dimension (N-S); # of Grid intervals in N-S is JNP-1  C JNP: 2nd dimension (N-S); number of Grid intervals in N-S is JNP-1
86  C NLAY: 3rd dimension (# of layers); vertical index increases from 1 at  C NLAY: 3rd dimension (number of layers); vertical index increases from 1 at
87  C       the model top to NLAY near the surface (see fig. below).  C       the model top to NLAY near the surface (see fig. below).
88  C       It is assumed that 6 <= NLAY <= JNP (for dynamic memory allocation)  C       It is assumed that 6 <= NLAY <= JNP (for dynamic memory allocation)
89  C  C
# Line 240  C A larger polar cap is used if j1=3 (re Line 240  C A larger polar cap is used if j1=3 (re
240  C winds are noisy near poles).  C winds are noisy near poles).
241  C  C
242  C Flux-Form Semi-Lagrangian transport in the East-West direction is used  C Flux-Form Semi-Lagrangian transport in the East-West direction is used
243  C when and where Courant # is greater than one.  C when and where Courant number is greater than one.
244  C  C
245  C The user needs to change the parameter Jmax or Kmax if the resolution  C The user needs to change the parameter Jmax or Kmax if the resolution
246  C is greater than 0.5 deg in N-S or 150 layers in the vertical direction.  C is greater than 0.5 deg in N-S or 150 layers in the vertical direction.
# Line 441  C Line 441  C
441        do 1500 k=1,NLAY        do 1500 k=1,NLAY
442  C  C
443        if(IGD.eq.0) then        if(IGD.eq.0) then
444  C Convert winds on A-Grid to Courant # on C-Grid.  C Convert winds on A-Grid to Courant number on C-Grid.
445        call A2C(U(1,1,k),V(1,1,k),IMR,JMR,j1,j2,CRX,CRY,dtdx5,DTDY5)        call A2C(U(1,1,k),V(1,1,k),IMR,JMR,j1,j2,CRX,CRY,dtdx5,DTDY5)
446        else        else
447  C Convert winds on C-grid to Courant #  C Convert winds on C-grid to Courant number
448        do 45 j=j1,j2        do 45 j=j1,j2
449        do 45 i=2,IMR        do 45 i=2,IMR
450  45    CRX(i,J) = dtdx(j)*U(i-1,j,k)  45    CRX(i,J) = dtdx(j)*U(i-1,j,k)
Legend:
Removed from v.3  
changed lines
  Added in v.10
  ViewVC Help
Powered by ViewVC 1.1.21