39 |
& (smooth2Dsize(smoothOpNb).EQ.2)) then |
& (smooth2Dsize(smoothOpNb).EQ.2)) then |
40 |
write(fnamegeneric(1:80),'(1a,i3.3)') |
write(fnamegeneric(1:80),'(1a,i3.3)') |
41 |
& 'smooth2Dscales',smoothOpNb |
& 'smooth2Dscales',smoothOpNb |
42 |
call mdsreadfield(fnamegeneric,smoothprec,'RL',1, |
CALL READ_REC_3D_RL(fnamegeneric,smoothprec, |
43 |
& smooth2D_Lx,1,mythid) |
& 1, smooth2D_Lx,1,1,mythid) |
44 |
call mdsreadfield(fnamegeneric,smoothprec,'RL',1, |
CALL READ_REC_3D_RL(fnamegeneric,smoothprec, |
45 |
& smooth2D_Ly,2,mythid) |
& 1, smooth2D_Ly,2,1,mythid) |
46 |
_EXCH_XY_RL ( smooth2D_Lx, myThid ) |
_EXCH_XY_RL ( smooth2D_Lx, myThid ) |
47 |
_EXCH_XY_RL ( smooth2D_Ly, myThid ) |
_EXCH_XY_RL ( smooth2D_Ly, myThid ) |
48 |
else |
else |
82 |
c write diffusion operator to file |
c write diffusion operator to file |
83 |
write(fnamegeneric(1:80),'(1a,i3.3)') |
write(fnamegeneric(1:80),'(1a,i3.3)') |
84 |
& 'smooth2Doperator',smoothOpNb |
& 'smooth2Doperator',smoothOpNb |
85 |
call mdswritefield(fnamegeneric,smoothprec,.false.,'RL', |
CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec, |
86 |
& 1,smooth2D_Kux,1,1,mythid) |
& 1,smooth2D_Kux,1,1,mythid) |
87 |
call mdswritefield(fnamegeneric,smoothprec,.false.,'RL', |
CALL WRITE_REC_3D_RL(fnamegeneric,smoothprec, |
88 |
& 1,smooth2D_Kvy,2,1,mythid) |
& 1,smooth2D_Kvy,2,1,mythid) |
89 |
|
|
90 |
end |
end |