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 |
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. |
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) |