Ocean Modeling Discussion

Hi all, I submitted a job, but (1) it shows the job is running, (2) no error is reported, (3) no file is created, such log file, history file, average file etc. After hours, I check the status, it gives: Job id Name User Time Use S Queue ...

Thanks very much! Now the problem has been solved and it is running smoothly. It turns out that the problem is caused by the wrong header file. I was using the header file I have been using from the old source code, which misses Ice boundary conditions . Again, Kate, thanks very much for your ...

Thanks very much for your reply. I turned debug on, still it is not working properly. log file with MPI: Time is Mon Jul 8 16:03:16 EDT 2013 This jobs runs on the following processors: iw-k30-24.pace.gatech.edu iw-k30-24.pace.gatech.edu iw-k30-24.pace.gatech.edu iw-k30-24.pace.gatech.edu iw-k30-24 ...

Hi all, When I was submitting a job. It stops immediately and the log file gives Caught error: Segmentation fault (signal 11) It seems there something with the mpi? But even I turn off the mpi, it gives error too: READ_PhyPar - Error while processing line: pu9 60848864 207 954289 59062278 40855 0 ...

Hi, I am struggling with the ICESHELF in my ROMS application. Here are what I did: 1. Turn on #define ICESHELF in the header file to compile. It works fine. 2. Include zice in the grid file. zice (negative values) is the thickness of iceshelf. Where there is no iceshelf, I manually set zice=0. 3 ...

I am running an application with iceshelf thickness ("zice" in the grid file), values being positive. But the first record is NaN and blows up. Without including "zice", the model runs well. So the problem is I am doing something wrong with the iceshelf thickness "zice"? Can anyone has done with ice ...

Thanks very much! Fan As a consistent data set for bedrock and ice shelf topography as well grounding line positions you might consider the following data set Timmermann, R, et al. (2010), " A consistent dataset of Antarctic ice sheet topography, cavity geometry, and global bathymetry ", Earth Syst ...

Thanks Kate! BTW, do you have any suggestions on where to obtain the iceshelf thickness data? I learned that some use this dataset: http://nsidc.org/api/metadata?id=gla12 Fan zice is on the rho grid. It is the thickness of ice above the ocean at each gridpoint. Since ice shelves can be hundreds of ...

I am trying to use ICESHELF for the model.

But it asks for the zice in grid file. Could someone tell what sizes should zice be? And which file requires specifying zice?

Thanks very much!

Fan

Yes. kate, thanks very much for you suggestion. I switched to the latest version of ROMS. Now it compiles smoothly. Fan Are you honestly using ROMS 3.0 and expecting us to remember anything about it? As far as I know, the most complete ice shelf code is from my branch. Check out the latest version ...

Hi, I am trying to include the ICESHELF into the cppdefs. It compiles well without ICESHLF. However, after I define ICESHLEF, it gave the following error about set_vbc.f90 mpif90 -c -O4 -fastsse -Mipa=fast,inline oi_update.f90 mpif90 -c -O4 -fastsse -Mipa=fast,inline set_avg.f90 mpif90 -c -O4 ...

Thanks Kate. Yes, the compiler is not configured properly. /usr/bin/cpp -P -traditional -I/usr/local/packages/netcdf/4.1.3/mvapich2-1.6/gcc-4.4.5/include -DMPI -DLINUX -DX86_64 -DPGI -IInclude -IAdjoint -INonlinear -IRepresenter -ITangent -IDrivers -ISeaIce Drivers/propagator.F > propagator.f90 Bin ...

When I was compiling ROMS, it gave me this error, even after I set in the makefile "MPI=off" 1. The problem is the "inline"? 2. Among f90 files, I could not find the mod_kinds.f90 3. In the propagator.f90, there are only two lines: MODULE propagator_mod END MODULE propagator_mod Your suggestions ...

Kate, thanks very much. Fan There's more than one kind of ice at sea. The ICESHELF option sounds like what you are looking for. It is a rigid ice cover for which you have to provide its thickness, which can be hundreds of meters. You may or may not also have sea ice, provided by the ICE_MODEL option ...

Hi, I am trying to run ROMS in the Ross Sea, which has sea ice shelf similar to the Weddell Sea. I find the Ice code by Kate from https://github.com/kshedstrom/roms/tree/kate_svn/ROMS/SeaIce I'd better use the ICESHELF in cppdefs or the ice code for the Ross Sea case? I suppose the Sea Ice module ...

We found out that it is due to the QCORRECTION in cppdefs file

Why is that? What factors may conflict with QCORRECTION?

Hi, I am running one configuration with all boundaries closed. But it blows up at the second time-step, gives NaN. I checked the pm and pn are all positive. Masks are OK. The input initial values are OK. Could anyone help guess what may be the potential causes? Thanks in advance for your help. I am ...

Kate, thanks very much for your reply.
Actually, the problem is in the cppdefs.h file; it is reading another application.

kate wrote:None of the files you show have a dimension of xi_psi=141. Are you sure that's not coming from the ocean.in file or some other file?

I am configuring a model. It seems that the dimensions of my input files are not consistent. After an examination of all 6 input files, I could not figure it out which file has the wrong dimension. I am posting the error and information about six input files. Hope someone could help me out! Thanks ...