Back to home page

MITgcm

 
 

    


Warning, /verification/aim.5l_cs/input.thSI/mk3regions_mask.m is written in an unsupported language. File is not indexed.

view on githubraw file Latest commit 60f2b34d on 2006-02-21 01:52:15 UTC
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