Back to home page

MITgcm

 
 

    


Warning, /verification/global_ocean.cs32x15/input/mk_bathy4gcm.m is written in an unsupported language. File is not indexed.

view on githubraw file Latest commit a3685688 on 2003-10-03 15:54:46 UTC
a368568840 Jean*0001 setMinDepth=input('setMinDepth (>0) = ? ');
                0002 
                0003 h0=rdda('bathy_cs32.bin',[6*32 32],1,'real*8','b');
                0004 h1=min(0,h0);
                0005 [I J]=find(-setMinDepth<h1 & h1<0 );
                0006 fprintf('make deeper %i grid points\n',length(I));
                0007 h1(find(-setMinDepth<h1 & h1<0 ))=-setMinDepth;
                0008 
                0009 namfil=['bathy_Hmin',int2str(setMinDepth),'.bin'];
                0010 fid=fopen(namfil,'w','b'); fwrite(fid,h1,'real*8'); fclose(fid);
                0011 fprintf([' write bathy on file:',namfil,' \n']);
                0012 
                0013 return
                0014 
                0015 rac='grid_cs32/';
                0016 set_axis 
                0017 load_cs
                0018 var=min(h1,0); var(find(var==0))=NaN;
                0019 shift=-1; cbV=0; AxBx=[-180 180 -90 90];
                0020 ccB=[-5000 100];
                0021 grph_CS(var,xcs,ycs,xcg,ycg,ccB(1),ccB(2),shift,cbV,AxBx);
                0022