150 |
& (4*k1AtC_str(i,j,bi,bj)*exx+2*k1AtC_str(i,j,bi,bj)*eyy+ |
& (4*k1AtC_str(i,j,bi,bj)*exx+2*k1AtC_str(i,j,bi,bj)*eyy+ |
151 |
& 4*0.5*k2AtC_str(i,j,bi,bj)*exy) |
& 4*0.5*k2AtC_str(i,j,bi,bj)*exy) |
152 |
|
|
|
! IF (bi.eq.2.and.bj.eq.2.and.i.eq.15.and. |
|
|
! & (exx.ne.0.0 .or. eyy.ne.0.0 .or. exy.ne.0.0)) THEN |
|
|
! PRINT *, "CG_FUNCTION", j, v(i,j,bi,bj),v(i+1,j,bi,bj), |
|
|
! & v(i,j+1,bi,bj),v(i+1,j+1,bi,bj) |
|
|
! ENDIF |
|
153 |
|
|
154 |
|
|
155 |
uret(i-1+inode,j-1+jnode,bi,bj) = |
uret(i-1+inode,j-1+jnode,bi,bj) = |
205 |
|
|
206 |
#ifdef ALLOW_STREAMICE |
#ifdef ALLOW_STREAMICE |
207 |
|
|
208 |
|
#ifdef STREAMICE_CONSTRUCT_MATRIX |
209 |
|
|
210 |
C the linear action of the matrix on (u,v) with triangular finite elements |
C the linear action of the matrix on (u,v) with triangular finite elements |
211 |
C as of now everything is passed in so no grid pointers or anything of the sort have to be dereferenced, |
C as of now everything is passed in so no grid pointers or anything of the sort have to be dereferenced, |
212 |
C but this may change pursuant to conversations with others |
C but this may change pursuant to conversations with others |
460 |
enddo |
enddo |
461 |
|
|
462 |
#endif |
#endif |
463 |
|
#endif |
464 |
RETURN |
RETURN |
465 |
END SUBROUTINE |
END SUBROUTINE |
466 |
|
|