[MITgcm-support] Increasing the number of grid points

Matthew Mazloff mmazloff at ucsd.edu
Fri Oct 1 15:44:44 EDT 2010


check that the eedata file is in your run directory?  If not copy it  
over from your previous set-up

-Matt




On Oct 1, 2010, at 12:29 PM, Joseph Schoonover wrote:

> So, I've successfully changed the binary files 'bathyPolR.bin' and  
> 'thetaPolR.bin' and have changed the proper parameters in 'data' and  
> 'size.h'. Upon execution, I get the following error:
>
> " Unable to open execution environment parameter file "eedata"
> STOP ABNORMAL END: S/R EESET_PARMS statement executed "
>
> I'm not quite sure what this means. After looking under the  
> rotating_tank portion of the manual, I do not see anything on eedata  
> other than the fact that this files contains no customizations for  
> this experiment.
>
> Is there a particular portion of the manual I should look over?
>
> Thanks,
> Joe Schoonover
>
> ----- Original Message -----
> From: "Christoph Völker" <Christoph.Voelker at awi.de>
> Date: Wednesday, September 22, 2010 12:57 am
> Subject: Re: [MITgcm-support] Increasing the number of grid points
> To: mitgcm-support at mitgcm.org
>
> <div style="font-family: 'Times New Roman'; font-size: 16px;  
> "><div>Dear Joe,</div><div><br /></div><div>the file that you  
> mention contains the topography for the experiment, i.e. one water  
> depth for each grid point. No wonder it does not fit anymore after  
> changing the number of grid points. You need to provide a new one  
> with the right number of grid points. </div><div>Do you have access  
> to matlab? Then I'd read in the old file like</div><div><br /></ 
> div><div>fid = fopen('bathyPolR.bin','r','ieee-be'); % files are  
> written in big-endian byte ordering</div><div>a =  
> fread(fid,'real*8');   % read in file contents either as double  
> precision binary (real*8) </ 
> div><div>                                     % or as single  
> precision (real*4); check with which way you</ 
> div><div>                                     % get the right number  
> of grid points</div><div>fclose(fid);</div><div>a = reshape(a, 
> [nx,ny]); % the grid is written as one long vector with nx*ny  
> elements;</div><div>                                     % you have  
> to make is a nx by ny matrix</div><div><br /></div><div>and then  
> think how to modify it, and write a new file using the same  
> methods.</div><div><br /></div><div>Cheers, Christoph</div><div><br / 
> ></div><div>On 21.09.10, <b class="name">Joseph Schoonover </b> &lt;js08s at fsu.edu 
> &gt; wrote:</div><blockquote  
> cite="mid:f749ae1322e01.4c98e106 at fsu.edu" class="iwcQuote"  
> style="border-left: 1px solid #00F; padding-left: 13px; margin-left:  
> 0;" type="cite"><div class="mimepart text plain">To whom it may  
> concern,<br /><br />I am working off of the rotating_tank experiment  
> under ../verification.  I have been able to run the model with  
> default parameters with no problem. Now, I would like to increase  
> the tank size without losing resolution, i.e., I need to increase  
> the number of grid points. I've gone into '../input/size.h' and  
> changed sNy to 60. Upon compilation I receive an error referring to  
> the binary 'bathyPolR.bin'<br /><br /> I'm not too familiar with how  
> binary files work, in general, but I know that I need to modify the  
> file 'bathyPolR.bin' in order to increase the number of grid  
> points.<br /><br /> Given that I want to keep the same precision,  
> how can I go about modifying 'bathyPolR.bin' ?<br /><br />Thanks in  
> advance,<br />Joe Schoonover<br /><br />Undergraduate Researcher,  
> GFDI<br /><br />_______________________________________________<br / 
> >MITgcm-support mailing list<br />MITgcm-support at mitgcm.org<br /><a  
> href="http://mitgcm.org/mailman/listinfo/mitgcm-support" target="1">http://mitgcm.org/mailman/listinfo/mitgcm-support 
> </a><br /></div></blockquote></div>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list