source: lmdz_wrf/WRFV3/external/io_netcdf/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: 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
22else
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
34endif
35
36  return
Note: See TracBrowser for help on using the repository browser.