% svd analysis of temperature and precipitation during winter
clear all
load r160
load t160
load sta160
kk = find( lon160 >= 100 );
lon = lon160( kk );
lat = lat160( kk );
r160 = r160( kk,1:end-1,12 ) + r160( kk,2:end,1 ) + ...
r160( kk,2:end,2 ) ;
t160 = t160( kk,1:end-1,12 ) + t160( kk,2:end,1 ) + ...
t160( kk,2:end,2 ) ;
t160 = t160./3./10;
r160 = zscore( r160',1 ); % using n
t160 = zscore( t160',1 );
clear ans
[ UL,VL,TU,TV,expl ] = svdgdy( t160,r160,10 );
xi = 100:130;
yi = 20:55;
[ glon,glat ] = meshgrid( xi,yi );
gz = griddata( lon,lat,UL(:,1),glon,glat );
clf
subplot(2,2,1);
domain = [ 70 140 10 60 ];
set( gcf,'color','w' );
hh = worldmap( [ domain(3) domain(4) ],...
[ domain(1) domain(2) ] );
setm( hh,'mapprojection','mercator' );
[ cs,hh ] = contourfm( glat,glon,gz,'linestyle','none' );
load coast
linem( lat,long,'color',[ .5 .5 .5 ] );
load china_lon_lat
linem( lat_china,lon_china,'linestyle','-', ...
'color',[ .5 .5 .5 ] );
subplot(2,2,2);
lat = lat160( kk );
gz2 = griddata( lon,lat,VL(:,1),glon,glat );
domain = [ 70 140 10 60 ];
set( gcf,'color','w' );
hh = worldmap( [ domain(3) domain(4) ],...
[ domain(1) domain(2) ] );
setm( hh,'mapprojection','mercator' );
[ cs,hh ] = contourfm( glat,glon,gz2,'linestyle','none' );
load coast
linem( lat,long,'color',[ .5 .5 .5 ] );
load china_lon_lat
linem( lat_china,lon_china,'linestyle','-', ...
'color',[ .5 .5 .5 ] );
subplot(2,2,3);
lat = lat160( kk );
gz3 = griddata( lon,lat,UL(:,2),glon,glat );
domain = [ 70 140 10 60 ];
set( gcf,'color','w' );
hh = worldmap( [ domain(3) domain(4) ],...
[ domain(1) domain(2) ] );
setm( hh,'mapprojection','mercator' );
[ cs,hh ] = contourfm( glat,glon,gz3,'linestyle','none' );
load coast
linem( lat,long,'color',[ .5 .5 .5 ] );
load china_lon_lat
linem( lat_china,lon_china,'linestyle','-', ...
'color',[ .5 .5 .5 ] );
subplot(2,2,4);
lat = lat160( kk );
gz4 = griddata( lon,lat,VL(:,2),glon,glat );
domain = [ 70 140 10 60 ];
set( gcf,'color','w' );
hh = worldmap( [ domain(3) domain(4) ],...
[ domain(1) domain(2) ] );
setm( hh,'mapprojection','mercator' );
[ cs,hh ] = contourfm( glat,glon,gz4,'linestyle','none' );
load coast
linem( lat,long,'color',[ .5 .5 .5 ] );
load china_lon_lat
linem( lat_china,lon_china,'linestyle','-', ...
'color',[ .5 .5 .5 ] );
