source: lmdz_wrf/WRFV3/external/io_pnetcdf/transpose.code @ 1

Last change on this file since 1 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 737 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
24!dir$ concurrent
25     do k=k1,k2
26        do j=j1,j2
27!dir$ prefervector
28!dir$ concurrent
29!cdir select(vector)
30           do i=i1,i2
31              Field(1:di,abs(ix-i),abs(jx-j),abs(kx-k)) = DFIELD
32           enddo
33        enddo
34     enddo
35  endif
36
37  return
Note: See TracBrowser for help on using the repository browser.