Ocean Modeling Discussion

It does seem to be a parallelization problem. I reset NAVG to 100 for diagnostic, and tried with a single processor - no crash. Then 2 - they're dual-cpu nodes, no crash. 4 nodes at 2x2 worked fine. 8 nodes at 2x4 or 1x8 crashed right at timestep 100. I'm rebooting nodes and doing a disk check for ...

The baroclinic timestep is very small due to complex topography in a shallow region where rivers are coming in. I was running with 15-20 second dt, but river inputs of ~100 m3/sec during the spring rains was causing blowups. To that end, I do a fair amount of restarting since I hate to run with a dt ...

Hi all, I just updated to the newest release and tried PERFECT_RESTART with similar issues. Here's the application info from the outfile: Operating system : Linux CPU/hardware : x86_64 Compiler system : ifort Compiler command : /opt/roms/mpich-1.2.7p1/bin/mpif90 Compiler flags : -i-static -ip -O2 ...

I am trying something similar - I would like to do spatially uniform winds over my curvilinear domain that vary in time. I have done this with BULK_FLUXES and full wind vectors (Uwind with dimensions (time, eta_rho, xi_rho), but this is space-intensive and seems silly to specify winds at every point ...

HI all, I'm modeling Rhode Island Sound and Narragansett Bay and we have been using ADCIRC harmonics from 2004 for the southern boundary tidal forcing. I'm trying to model summer 2006 now and things are not going as well. I tried just adding two years with tide_start but haven't had much luck. Would ...

To hopefully settle this - Kate, you are right, editmask is failing while reading the mask variables, but for some reason it tries to read hraw as well. This could easily be edited out of the code.

line 123 of read_mask.m in my version:
if (got.hraw),
bath=nc_read(Gname,Vname.hraw,1);

Hi Sankar,

hraw is tricky... it needs a time parameter. For example,
ncwrite('grid.nc','hraw',MyH,1)

I'm not certain what hraw is actually used for, 'h' is the bathymetry variable in the gridfile that is active in my setup. Editmask will likely work once this is assigned.

Thank you both - a combination seemed to work just fine. I commented the old version checks out and changed mexcdf.m to read %if eval(v(1)) > 4 % fcn = 'mexcdf53'; % Matlab-5 or 6. %elseif eval(v(1)) == 4 % fcn = 'mexcdf4'; % Matlab-4 only. %elseif eval(v(1)) == '7' fcn = 'mexnc'; %else % error ...

We are trying to upgrade to matlab 7 but editmask does not seem to want to work. It looks like the change from mexcdf60 to mexnc is the issue: >> editmask spherical is a scalar and has a value of T mask_rho has the following dimensions (input order): 1) eta_rho = 150 2) xi_rho = 300 mask_rho loaded ...