1 |
|
2 |
%- write a 3 regions mask (based on latitude) for DiagStat output: |
3 |
% region 1 = South , yLat <= -24 ; |
4 |
% region 2 = Tropics , -24 < yLat < 24 ; |
5 |
% region 2 = North , 24 <= yLat ; |
6 |
|
7 |
yc=rdmds('YC'); |
8 |
regMsk=ones(size(yc)); |
9 |
|
10 |
regMsk(find( -24 < yc))=2; |
11 |
regMsk(find( 24 <= yc ))=3; |
12 |
|
13 |
namfil='regMask_lat24.bin'; |
14 |
fid=fopen(namfil,'w','b'); fwrite(fid,regMsk,'real*8'); fclose(fid); |
15 |
fprintf([' write mask to file: ',namfil,'\n']); |
16 |
|
17 |
return |
18 |
|
19 |
%- a plot to check: |
20 |
%rac='ref/'; |
21 |
%load_cs; |
22 |
%ccB=[0 4]; shift=0; cbV=0; AxBx=[-180 180 -90 90]; kEnv=0; |
23 |
%var=regMsk; |
24 |
%figure(1);clf; |
25 |
%grph_CS(var,xcs,ycs,xcg,ycg,ccB(1),ccB(2),shift,cbV,AxBx,kEnv); |
26 |
|