Warning, /utils/matlab/cs_grid/latloncap/quikplot_llc.m is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit 370f350e on 2013-10-28 16:09:40 UTC
0524fa3d09 Dimi*0001 function quikplot_llc(fld);
cb1c045303 Dimi*0002
0524fa3d09 Dimi*0003 % Function quikplot_llc(fld)
cb1c045303 Dimi*0004 % plot lat-lon-cap field
0005 %
0006 % INPUTS
0007 % fld input array of dimension nx*nx*13
0008
0524fa3d09 Dimi*0009 if nargin < 1, error('please specify field to plot'); end
cb1c045303 Dimi*0010
0011 [nx ny]=size(fld);
0012 ny=nx*3;
0013
0014 % read face 1, nx*nx*3
0015 f1=fld(:,1:(nx*3));
0016
0017 % read face 2, nx*nx*3
0018 f2=fld(:,(nx*3+1):(nx*6));
0019
0020 % read face 3, nx*nx
0021 f3=fld(:,(nx*6+1):(nx*7));
0022
0023 % read face 4, nx*3*nx
0024 f4=f1';
0025 for f=8:10
0026 i1=(1:nx)+(f-8)*nx;
0027 i2=(1:3:(nx*3))+7*nx+f-8;
0028 f4(i1,:)=fld(:,i2);
0029 end
0030
0031 % read face 5, nx*3*nx
0032 f5=f1';
0033 for f=11:13
0034 i1=(1:nx)+(f-11)*nx;
0035 i2=(1:3:(nx*3))+10*nx+f-11;
0036 f5(i1,:)=fld(:,i2);
0037 end
0038
0039 % plot field
0040 f=nan*ones(4*nx,ny+nx/2);
0041 f(1:nx,1:ny)=f1;
0042 f((nx+1):(2*nx),1:ny)=f2;
0043 f(1:nx,(ny+1):(ny+nx/2))=rot90(f3(1:(nx/2),:),1);
0044 f((2*nx+1):(3*nx),(ny+1):(ny+nx/2))=rot90(f3((nx/2+1):nx,:),3);
0045 f((2*nx+1):(3*nx),1:ny)=rot90(f4,3);
0046 f((3*nx+1):(4*nx),1:ny)=rot90(f5,3);
370f350e58 Dimi*0047 quikpcolor(f');
cb1c045303 Dimi*0048 set(gca,'xtick',[],'ytick',[])
