<div dir="ltr">Hello MITgcm Support experts,<br><div><br></div><div><div dir="ltr">    <font color="#000000">I work on idealized simulation on the curvilinear coordinates in MITgcm.  I am trying to prepare the <font face="arial, helvetica, sans-serif">horizGridFile for the </font>curvilinear coordinates.</font>   By trial and error I completed a matlab script below to generate <span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif">horizGridFile for the curvilinear coordinates.</span>  Do you think this file grid.bin as the attachment satisfies the need of 10N-90N <span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif">for the curvilinear coordinates </span>in MITgcm?  </div><div dir="ltr"><br></div><div dir="ltr">    I used grid.bin in data in execution of MITgcm and got almost the same temperature field as giving delX=320*0.25, and delY=320*0.25 in the data file.  In fact, I got almost the same temperature field with the same <span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif">horizGridFile</span><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif"> </span>in the following two conditions of f0 and ygOrigin after 750000 sec intergration.  Although I work on idealized simulations, I can't help but ask do I need to prepare different <span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif">horizGridFile</span><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif"> for 0N-80N and 10N-90N?  How to let MITgcm know the domain is 10N-90N instead of 1N-9N?</span></div><div dir="ltr"><br></div><div dir="ltr"><div dir="ltr">f0=2.531787E-5,,<br></div><div dir="ltr">ygOrigin=10.,</div><div dir="ltr"><br></div></div><div dir="ltr">f0=0.,<br></div><div dir="ltr">ygOrigin=0.<br></div><div dir="ltr"><br></div><div>    Thanks much,</div><div><br></div><div>                      Yi-Chih</div><div dir="ltr"><div><br></div><div>####################################################</div><div><div>addpath('C:\Users\huang\Documents\MITgcm\Matlab_scripts')</div><div>cd C:\Users\huang\Documents\MITgcm\vel_out\run_n_320_253t;</div><span class="gmail-im" style="color:rgb(80,0,80)"><div><br></div><div>xc = rdmds('XC');</div><div>yc = rdmds('YC');</div></span><div>lon = 10.25:0.25:90;</div><div>lat = 10.25:0.25:90; </div><span class="gmail-im" style="color:rgb(80,0,80)"><div>% 2d field fld2d on lon-lat grid:</div><div>% fld = interp2(lon,lat,fld2d,xc(:),yc(:),method); % method = 'linear','cubic', etc see help interp2</div></span><div>A=rdmds('T.0000000000.001.001');</div><div>size(A);</div><div>fld = interp2(lon,lat,A(:,:,1),xc(:),yc(:),'spline');</div><div>fld = reshape(fld,size(xc));</div><div>writebin('grid.bin', fld, 1, 'real*4');</div></div></div></div></div>