665 |
real VT1(IM,JM) |
real VT1(IM,JM) |
666 |
integer i |
integer i |
667 |
|
|
668 |
|
print *,' In msc2grd ',(igrd(i),i=1,nchp) |
669 |
|
print *,' In msc2grd ',(chfr(i),i=1,nchp) |
670 |
|
|
671 |
IF(NCHP.GE.0) THEN |
IF(NCHP.GE.0) THEN |
672 |
DO I=1,IM*JM |
DO I=1,IM*JM |
673 |
VT1(I,1) = ZERO |
VT1(I,1) = ZERO |
677 |
VT1(IGRD(I),1) = VT1(IGRD(I),1) + B(I)*CHFR(I) |
VT1(IGRD(I),1) = VT1(IGRD(I),1) + B(I)*CHFR(I) |
678 |
ENDDO |
ENDDO |
679 |
|
|
680 |
|
stop |
681 |
|
|
682 |
DO I=1,IM*JM |
DO I=1,IM*JM |
683 |
A(I,1) = A(I,1)*(ONE-FRACG(I,1)) + VT1(I,1) |
A(I,1) = A(I,1)*(ONE-FRACG(I,1)) + VT1(I,1) |
684 |
ENDDO |
ENDDO |