Parent Directory
| 
Revision Log
| 
Revision Graph
initial checkin of topography and hydrography interpolation scripts for the llc-grid, based on old matlab scripts by Alistair Adcroft Let's hope, they are useful.
| 1 | function b=mdsiocompact(a); |
| 2 | % hack to convert to mdsio compact format |
| 3 | |
| 4 | [nx,ny,nz,nt]=size(a); |
| 5 | if nx/4~=round(nx/4); |
| 6 | error('not a llc field'); |
| 7 | end |
| 8 | |
| 9 | n4=nx/4; |
| 10 | m=(ny-n4)/n4; |
| 11 | m4=m*n4; |
| 12 | for kt = 1:nt |
| 13 | for kz = 1:nz |
| 14 | for k=1:4 |
| 15 | sides{k} = a([1:n4]+(k-1)*n4,1:m4,kz,kt); |
| 16 | end |
| 17 | cap = a([1:n4]+1*n4,m4+1:end,kz,kt); |
| 18 | % reformat |
| 19 | btmp = [sides{1} ... |
| 20 | sides{2} ... |
| 21 | cap ... |
| 22 | reshape(rot90(sides{3},1),[n4 m4]) ... |
| 23 | reshape(rot90(sides{4},1),[n4 m4]) ... |
| 24 | cap*0]; |
| 25 | b(:,:,kz,kt)=btmp; |
| 26 | end |
| 27 | end |
| 28 | |
| 29 | return |
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |