Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Thu Mar 21, 2019 4:23 pm




Post new topic Reply to topic  [ 2 posts ] 

All times are UTC

Author Message
PostPosted: Fri Mar 03, 2017 3:01 am 
Offline

Joined: Tue Mar 22, 2016 5:21 pm
Posts: 16
Location: Japan Weather Association
Hello.

I am a ROMSTOOLS user.
I want to analyze the region from 170E to 170W across the 180E.
Therefore I set lonmin=170 and lonmax=190 in the romstools_param.m.
But, when I ran make_clim after make_forcing, the following errors occurred.
Would you please tell me how to solve this problem ?

------------------------------------------------------------------------------------
Create the OA file...

Creating the file : /home/mathuser/work/roms/Roms_tools/Run/ROMS_FILES/roms_oa.nc


Horizontal extrapolations

Temperature...

Ext tracers: Roa = 0 km - default value = NaN
Ext tracers: horizontal interpolation of the seasonal data

time index: 1 of total: 12
Error: griddedInterpolant

Error: interp2>makegriddedinterp (line 228)
F = griddedInterpolant(varargin{:});

Error: interp2 (line 128)
F = makegriddedinterp({X, Y}, V, method,extrap);

Error: ext_tracers (line 162)
datazgrid(k,:,:)=interp2(x,y,data,lon,lat,'cubic');

Error: make_clim (line 134)
ext_tracers(oaname,temp_month_data,temp_ann_data,...

------------------------------------------------------------------------------------

Thanks for any help.


Top
 Profile  
Reply with quote  
PostPosted: Fri Mar 03, 2017 7:20 am 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3531
Location: IMS/UAF, USA
I am not a ROMSTOOLS user, nor is this the forum run by the developers of it. Still, I have run Arctic domains that need to interpolate over the full range of longitude. In general, you want to be consistent, say with a grid that has longitude values between 0 and 360, bathymetry with longitude values between 0 and 360, ditto the T and S fields you are interpolating from. It could just as well be -180 to 180, but still be consistent. I have Python or NCL scripts to convert my grid from one convention to the other, that being easier than converting SODA or whatever. Our MERRA scripts convert from the range (-180,180) to (0,360) as part of the download process.

Then you still need tools that can deal with the dateline or the prime meridian, whichever. For scrip, there's a periodic flag that needs to be set. It's probably true for any other interpolating tool as well. My branch of ROMS has a flag for dealing with periodic forcing when ROMS is interpolating internally, but I don't think the trunk code has it. Note to Hernan for his Arctic plans... :mrgreen:


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 2 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group