| /phasta/phSolver/common/ |
| H A D | renum.f | 150 do jj = 1, nfint
152 xts1 = x(nrint(jj),1) !*cos(xcyl(nrint(jj),2)+tolerence)
153 xts2 = x(nrint(jj),2) !*sin(xcyl(nrint(jj),2)+tolerence)
154 xts3 = x(nrint(jj),3) + plandist
157 xsinfin(jj,1,:) = xtmp(:)
158 elcnfin(jj,1) = etmp
159 imax(jj) = (thmax-thmin)*xcyl(nrint(jj),1)/ds
160 do i=1,imax(jj)
161 if ( xcyl(nrint(jj),1) .eq. radcyl) then
162 xts1 = (xcyl(nrint(jj),1)-tolerence)
[all …]
|
| /phasta/phSolver/compressible/ |
| H A D | timedata.f90 | 113 do jj=1,ntspts ! read coordinate data where solution desired
114 read(626,*) ptts(jj,1), ptts(jj,2), ptts(jj,3)
159 do jj = 1, ntspts
160 write(1001, "(I5, T12, 3(F16.12))") jj, ptts(jj,1:3)
199 do jj = 1, ntspts
200 vartssoln((jj-1)*ndof+1:jj*ndof)=varts(jj,:)
221 do jj = 1, ntspts
222 indxvarts = (jj-1)*ndof
225 varts(jj,k) =
241 do jj = 1, ntspts
[all …]
|
| H A D | timeseries.f | 20 do jj = 1, ntspts
22 if(statptts(jj,1).gt.0) then
23 if(statptts(jj,1).eq.iblkts) then
25 call shphex (ipord, parptts(jj,:),shape(:),
28 call shptet (ipord, parptts(jj,:),shape(:),
31 founde=statptts(jj,2)
34 xts1 = ptts(jj,1)
35 xts2 = ptts(jj,2)
36 xts3 = ptts(jj,3)
212 statptts(jj,1)=iblkts
[all …]
|
| H A D | e3b.f | 502 do jj = 1, nshl ! first find the value of levelset at point ii
505 & + shape(ii,jj) * ycl(ii,jj,6)
|
| H A D | e3source.f | 439 do jj = 1, nshl ! first find the value of levelset at point ii
442 & + shape(ii,jj) * ycl(ii,jj,6)
|
| /phasta/phSolver/incompressible/ |
| H A D | timeseries.f | 20 do jj = 1, ntspts
22 if(statptts(jj,1).gt.0) then
23 if(statptts(jj,1).eq.iblkts) then
25 call shphex (ipord, parptts(jj,:),shape(:),
28 call shptet (ipord, parptts(jj,:),shape(:),
31 founde=statptts(jj,2)
34 xts1 = ptts(jj,1)
35 xts2 = ptts(jj,2)
36 xts3 = ptts(jj,3)
212 statptts(jj,1)=iblkts
[all …]
|
| H A D | itrdrv.f | 1374 do jj=1,ntspts ! read coordinate data where solution desired
1375 read(626,*) ptts(jj,1),ptts(jj,2),ptts(jj,3)
1397 do jj=1,ntspts
1400 jjm1 = jj-1
1406 iv_rank(jj) = iv_node*iv_rankpernode
1411 write(*,*) ' Probe', jj, 'handled by rank',
1412 & iv_rank(jj), ' on node', iv_node
1416 if(iv_rank(jj) .lt.0 .or. iv_rank(jj) .ge. numpe) then
1417 write(*,*) 'WARNING: iv_rank(',jj,') is ', iv_rank(jj),
1419 iv_rank(jj) = numpe-1
[all …]
|
| H A D | e3source.f | 483 do jj = 1, nshl
484 kay(ii) = kay(ii) + shape_funct(ii,jj) * yl(ii,jj,6)
486 & + shape_funct(ii,jj) * yl(ii,jj,7)
527 do jj = 1, nshl ! first find the value of levelset at point ii
530 & + shape_funct(ii,jj) * yl(ii,jj,6)
|
| /phasta/phSolver/AMG/ |
| H A D | ramg_coarse.f | 37 integer :: I_nnz,ki,kj,jj local
260 jj = amg_A_rowp(level1)%p(kj)
261 m = aLoc(jj)
264 aLoc(jj) = m
265 amg_Fn(m) = jj
271 if ( (amg_F(jj).eq.1) .and. (mod(amg_S(kj),2).eq.1) ) then
|
| H A D | ramg_tools.f | 946 integer ii,jj local
967 jj = ncorp_map(j)!ipmap(j))
969 write(264,mformat)ii,jj,(alhs(k,p),k=1,redun)
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 459 do jj = 1, nshl
460 kay(ii) = kay(ii) + shape_funct(ii,jj) * yl(ii,jj,6)
462 & + shape_funct(ii,jj) * yl(ii,jj,7)
503 do jj = 1, nshl ! first find the value of levelset at point ii
506 & + shape_funct(ii,jj) * yl(ii,jj,6)
|
| H A D | e3source.fRHSonly | 445 do jj = 1, nshl
446 kay(ii) = kay(ii) + shape_funct(ii,jj) * yl(ii,jj,6)
448 & + shape_funct(ii,jj) * yl(ii,jj,7)
489 do jj = 1, nshl ! first find the value of levelset at point ii
492 & + shape_funct(ii,jj) * yl(ii,jj,6)
|