aim.5l_cs/input.thSI/mk3regions_mask.m

60f2b34db4 Jean*0001
                0002 %- write a 3 regions mask (based on latitude) for DiagStat output:
                0003 %   region 1 = South   ,     yLat <= -24 ;
                0004 %   region 2 = Tropics , -24 < yLat < 24 ;
                0005 %   region 2 = North   ,  24 <= yLat ;
                0006
                0007 yc=rdmds('YC');
                0008 regMsk=ones(size(yc));
                0009
                0010 regMsk(find( -24 < yc))=2;
                0011 regMsk(find( 24 <= yc ))=3;
                0012
                0013 namfil='regMask_lat24.bin';
                0014 fid=fopen(namfil,'w','b'); fwrite(fid,regMsk,'real*8'); fclose(fid);
                0015 fprintf([' write mask to file: ',namfil,'\n']);
                0016
                0017 return
                0018
                0019 %- a plot to check:
                0020 %rac='ref/';
                0021 %load_cs;
                0022 %ccB=[0 4]; shift=0; cbV=0; AxBx=[-180 180 -90 90]; kEnv=0;
                0023 %var=regMsk;
                0024 %figure(1);clf;
                0025 %grph_CS(var,xcs,ycs,xcg,ycg,ccB(1),ccB(2),shift,cbV,AxBx,kEnv);
                0026