|
Last change
on this file since 3567 was
11,
checked in by aslmd, 14 years ago
|
|
spiga@svn-planeto:ajoute le modele meso-echelle martien
|
|
File size:
676 bytes
|
| Line | |
|---|
| 1 | ix=0 |
|---|
| 2 | jx=0 |
|---|
| 3 | kx=0 |
|---|
| 4 | call reorder(MemoryOrder,MemO) |
|---|
| 5 | if(IACHAR(MemO(1:1)) > MaxUpperCase) ix=i2+i1 |
|---|
| 6 | if(IACHAR(MemO(2:2)) > MaxUpperCase) jx=j2+j1 |
|---|
| 7 | if(IACHAR(MemO(3:3)) > MaxUpperCase) kx=k2+k1 |
|---|
| 8 | |
|---|
| 9 | ! pjj/cray |
|---|
| 10 | if(IO == 'write') then |
|---|
| 11 | !dir$ concurrent |
|---|
| 12 | do k=k1,k2 |
|---|
| 13 | do j=j1,j2 |
|---|
| 14 | !dir$ prefervector |
|---|
| 15 | !dir$ concurrent |
|---|
| 16 | !cdir select(vector) |
|---|
| 17 | do i=i1,i2 |
|---|
| 18 | DFIELD = Field(1:di,abs(ix-i),abs(jx-j),abs(kx-k)) |
|---|
| 19 | enddo |
|---|
| 20 | enddo |
|---|
| 21 | enddo |
|---|
| 22 | else |
|---|
| 23 | !dir$ concurrent |
|---|
| 24 | do k=k1,k2 |
|---|
| 25 | do j=j1,j2 |
|---|
| 26 | !dir$ prefervector |
|---|
| 27 | !dir$ concurrent |
|---|
| 28 | !cdir select(vector) |
|---|
| 29 | do i=i1,i2 |
|---|
| 30 | Field(1:di,abs(ix-i),abs(jx-j),abs(kx-k)) = DFIELD |
|---|
| 31 | enddo |
|---|
| 32 | enddo |
|---|
| 33 | enddo |
|---|
| 34 | endif |
|---|
| 35 | |
|---|
| 36 | return |
|---|
Note: See
TracBrowser
for help on using the repository browser.
