Search found 79 matches

by hetland
Fri Jun 27, 2014 2:11 pm
Forum: ROMS Tools and Techniques
Topic: GRIDGEN by Pavel Sakov url is wrong
Replies: 4
Views: 2731

Re: GRIDGEN by Pavel Sakov url is wrong

You can find a very recent version in the octant tree. They are basically mirrored here:

https://github.com/hetland/octant/tree/master/external

I assume since you got octant, you can also get this.
by hetland
Fri Jun 27, 2014 1:59 pm
Forum: ROMS Tools and Techniques
Topic: GRIDGEN by Pavel Sakov url is wrong
Replies: 4
Views: 2731

Re: GRIDGEN by Pavel Sakov url is wrong

You can find it here: https://code.google.com/p/gridgen-c/ But you need a few more things to compile it. All of the relevant links are: https://code.google.com/p/gridgen-c/ https://code.google.com/p/gridutils-c/ https://code.google.com/p/nn-c/ https://code.google.com/p/csa-c/ You're welcome, -Rob
by hetland
Thu Jun 26, 2014 2:26 pm
Forum: ROMS Tools and Techniques
Topic: octan install problem
Replies: 4
Views: 1850

Re: octan install problem

Frist of all, the version of octant on googlecode is very old, and should be considered obsolete. I have migrated the codebase to github. You can find it here: https://github.com/hetland/octant Some recent changes include pulling some of the compiled code out so that it is easier to install. Notably...
by hetland
Mon Sep 30, 2013 12:54 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS installation
Replies: 40
Views: 16729

Re: PYROMS installation

Did you make either package before using it? You need to do something like this: python setup.py build sudo python setup.py install In order for the C-libraries to be built, and to install everything in the right place. If you did this, then your build is failing somewhere, and you need to look at t...
by hetland
Fri Sep 27, 2013 1:10 am
Forum: ROMS Tools and Techniques
Topic: PYROMS installation
Replies: 40
Views: 16729

Re: PYROMS installation

Which pyroms?
by hetland
Mon Sep 02, 2013 5:34 pm
Forum: Job Opportunities
Topic: Postdoc/Research Associate position open at Texas A&M
Replies: 0
Views: 2600

Postdoc/Research Associate position open at Texas A&M

Texas A&M University, Department of Oceanography seeks to hire a Postdoctoral Research Associate in Coastal Ocean Hydrodynamic Modeling, starting as soon as possible . The appointment is for 12 months, with an additional funding of at least one year available upon satisfactory performance. The posit...
by hetland
Sun Jul 14, 2013 2:29 am
Forum: ROMS Tools and Techniques
Topic: Problems on import octant
Replies: 2
Views: 1772

Re: Problems on import octant

The googlecode version of octant is deprecated. You should use the one on github: https://github.com/hetland/octant You can download it using the command: git clone https://github.com/hetland/octant.git For this version, I separated out the gridgen stuff, so it should be much easier to install. -Rob...
by hetland
Sat Feb 02, 2013 4:41 pm
Forum: ROMS Discussion
Topic: Problem with the river discharge
Replies: 6
Views: 2110

Re: Problem with the river discharge

There are a number of reasons for this discrepancy. First of all, what does the data represent? A climatology? A synoptic cruise? At what level were the salinities measured (usually a few meters depth for standard hydrography)? As for errors in the model, the differences you show are too far away fo...
by hetland
Sat Nov 17, 2012 4:59 pm
Forum: ROMS Tools and Techniques
Topic: Pyroms or Octant and Enthought
Replies: 2
Views: 1746

Re: Pyroms or Octant and Enthought

Can you give a more specific error? I think this is a very easy fix -- just changing one line of an import. But we need to know which file is the problem.
by hetland
Mon Nov 05, 2012 5:34 pm
Forum: ROMS Usage
Topic: River discharge issues
Replies: 10
Views: 5786

Re: River discharge issues

I am not sure how to specify the data in the forcing file. I have the following but I am not sure if it is correct, do I need to have river_Vshape? And should I have a transport value for each time and river? [/code] Yes, and yes. I am not sure if there is a 'default' Vshape, but I think it is best...
by hetland
Fri Nov 02, 2012 2:00 am
Forum: ROMS Usage
Topic: River discharge issues
Replies: 10
Views: 5786

Re: River discharge issues

Here is a (slightly old) python script that will give you some guidance in making river forcing files. Also, if you don't see the rivers at all, you need to back up some steps. Make sure you can get the RIVERPLUME test case to work, and then work from there. Make sure you are compiling the code you ...
by hetland
Thu Mar 08, 2012 1:01 am
Forum: ROMS Information
Topic: Central East Florida Shelf Realistic Application Setup Help
Replies: 7
Views: 7175

Re: Central East Florida Shelf Realistic Application Setup H

All ROMS applications can take forcing and bathymetry (and other) information from either an analytical file (look in the Analytical directory for many, many examples) or from a netCDF file. The variables, in the end, are the same to ROMS, it is just a matter of whether they are defined by FORTRAN c...
by hetland
Thu Sep 15, 2011 9:54 pm
Forum: Job Opportunities
Topic: Postdoc position at Texas A&M University
Replies: 0
Views: 1677

Postdoc position at Texas A&M University

Texas A&M University, Department of Oceanography seeks to hire a Postdoctoral Research Associate in Coastal Ocean Hydrodynamic Modeling, starting as soon as possible . The appointment is for 12 months, with an additional funding of at least one year available upon satisfactory performance. This proj...
by hetland
Tue Mar 29, 2011 2:10 pm
Forum: ROMS Tools and Techniques
Topic: pyroms problem running/installing
Replies: 8
Views: 3112

Re: pyroms problem running/installing

I'm not sure what is going on. Since you can import netCDF4 by hand, it seems like that is all fine. You can try to find where pypurene is loaded in pyroms (it will be in a try/except block), and hardwire in a netcdf4 import. I don't think that the latest version of octant does this, though. You may...
by hetland
Tue Mar 29, 2011 2:06 pm
Forum: ROMS Discussion
Topic: vertical mixing & noise in a channel flow
Replies: 11
Views: 2736

Re: vertical mixing & noise in a channel flow

Your noise looks like Rossby waves (note the pattern that the noise is least at the equator, and seems to be propagating faster to the south). I suspect that you have some strangeness at your open boundaries that is causing this noise.

What are your OBCs?
by hetland
Tue Mar 22, 2011 1:55 pm
Forum: ROMS Tools and Techniques
Topic: netcdf4-python
Replies: 7
Views: 4260

Re: netcdf4-python

I suspect that your netcdf build actually built netcdf3 libraries. Double check the include and library dirs. Look at the output of configure and make.

Also, DAP support is built into the netCDF4 libraries, not HDF. So there should be no problem there.
by hetland
Wed Dec 22, 2010 5:25 pm
Forum: ROMS Tools and Techniques
Topic: A Grid value problem by Octant
Replies: 9
Views: 2833

Re: A Grid value problem by Octant

I usually do something like this to fill in undefined values (assuming you use masked arrays for your variables): grd.dx = grd.dx.filled( grd.dx.mean() ) grd.dy = grd.dy.filled( grd.dy.mean() ) grd.dndx = grd.dndx.filled( grd.dndx.mean() ) grd.dmde = grd.dmde.filled( grd.dmde.mean() ) grd.angle = gr...
by hetland
Fri Dec 17, 2010 5:09 pm
Forum: ROMS Tools and Techniques
Topic: gridgen-c warning: F(): NaN detected or Segmentation fault
Replies: 11
Views: 4450

Re: gridgen-c warning: F(): NaN detected or Segmentation fau

Right... After looking at the octant code - it seems that it does not save sigmas to disk (Rob?). If you are using a later version of octant that calls the gridgen libraries directly (instead of a straight command line call to gridgen), no files are saved. I think this is what you are using, so you...
by hetland
Mon Dec 13, 2010 12:58 am
Forum: ROMS Tools and Techniques
Topic: gridgen-c warning: F(): NaN detected or Segmentation fault
Replies: 11
Views: 4450

Re: gridgen-c warning: F(): NaN detected or Segmentation fau

Well, since you can get octant.gridgen to work (I assume the output of the square looked sensible?), then it must be the gridgen library. I still suspect this might be a 32/64 bit issue. There are many options to put into griggen to tell how many times to itterate, for example. You can try to modify...
by hetland
Sun Dec 12, 2010 2:05 am
Forum: ROMS Tools and Techniques
Topic: gridgen-c warning: F(): NaN detected or Segmentation fault
Replies: 11
Views: 4450

Re: gridgen-c warning: F(): NaN detected or Segmentation fau

Can you see if your code works with a square, that is x = [1, 0, 0, 1] y = [1.1, 1.0, 0., 0.] beta = [1, 1, 1, 1] That way you can isolate the problem to the gridgen library. Also, I seem to remember some issues with compiling the gridgen library with 64 bit. I have found that 32 bit python, numpy, ...
by hetland
Mon May 03, 2010 12:47 am
Forum: ROMS Ecosystem
Topic: Quick question for vertical mixing
Replies: 6
Views: 3748

Re: Quick question for vertical mixing

The mixing value for salinity is used for all passive tracers, as this represents the mixing of actual molecules. Heat can be transferred without moving molecules around.
by hetland
Mon Apr 12, 2010 10:53 pm
Forum: Job Opportunities
Topic: Postdoc at Texas A&M University
Replies: 0
Views: 1604

Postdoc at Texas A&M University

Texas A&M University, Department of Oceanography seeks to hire a Postdoctoral Research Associate in Coastal Ocean Hydrodynamic Modeling, starting September 1, 2010 or sooner if possible. The appointment is for 12 months , with an additional 12 months of funding available upon satisfactory performanc...
by hetland
Fri Apr 09, 2010 2:35 pm
Forum: ROMS Discussion
Topic: do i have to add the UV_TIDES for a real case?
Replies: 7
Views: 3180

Re: do i have to add the UV_TIDES for a real case?

I think it is important to understand if you have a propagating or standing wave. Usually, this boils down to large vs. small domains, respectively. I usually just specify sea-level on my boundaries for small domains, since the tides at my boundaries have both the incoming and outgoing waves in appr...
by hetland
Sat Oct 31, 2009 12:48 am
Forum: ROMS Problems
Topic: MPI process terminated unexpectedly
Replies: 7
Views: 3997

Re: MPI process terminated unexpectedly

There is a known issue with openmpi that gives an error like this. I used version openmpi-1.2.8 and got this error (this is the version that ships with the latest version of rocks, so you probably have this, or older..). This flag helped: -mca mpi_leave_pinned 0 as in: /usr/mpi/gcc/openmpi-1.2.8/bin...
by hetland
Sat Oct 10, 2009 7:23 pm
Forum: ROMS Bugs
Topic: writing to nonexisting history files
Replies: 8
Views: 4522

Re: writing to nonexisting history files

I think that the current logic is flexible enough to allow you to manipulate the I/O in a satisfactory way. The code is not flexible, it is broken. It does not do what is expected for the obvious cases. There is absolutely no reason that this simple subroutine should be so complicated that it depen...
by hetland
Wed May 27, 2009 9:07 pm
Forum: Job Opportunities
Topic: Postdoc research associate at Texas A&M University
Replies: 0
Views: 1527

Postdoc research associate at Texas A&M University

Texas A&M University, Department of Oceanography seeks to hire a Postdoctoral Research Associate in Coastal Ocean Hydrodynamic Modeling, starting September 1, 2009. The appointment is for 24 months . Primary responsibilities include assessing and improving a hydrodynamic model of surface currents ov...
by hetland
Mon Mar 30, 2009 5:10 pm
Forum: ROMS Tools and Techniques
Topic: Grid generation
Replies: 61
Views: 34372

Re: Grid generation

Pavel moved to Norway, and his web site was taken down by CSIRO. You can find his codes here:

http://code.google.com/p/octant/source/ ... ant/extern

Get it with a command like this:

svn checkout http://octant.googlecode.com/svn/trunk/octant/extern sakov

-Rob
by hetland
Tue Feb 10, 2009 1:20 am
Forum: ROMS Problems
Topic: UV_U3ADV_SPLIT in 2D applications broken
Replies: 1
Views: 897

UV_U3ADV_SPLIT in 2D applications broken

It seems the new(?) UV_U3ADV_SPLIT does not work in 2D applications. It expects a tracer variable to be set, which of course does not exist.
by hetland
Tue Feb 10, 2009 12:56 am
Forum: ROMS Benchmarks
Topic: Parallel I/O benchmarks?
Replies: 13
Views: 5476

Re: Parallel I/O benchmarks?

I can't see why it would depend on application. I guess it would depend on hardware, compiler, etc. The previous benchmark explicitly left out i/o, so as only to test the computational engine. For an i/o benchmark, we would need a very very large domain that did as little computation as possible (di...
by hetland
Fri Nov 07, 2008 4:55 pm
Forum: ROMS Discussion
Topic: timestep musings
Replies: 8
Views: 2697

Re: timestep musings

For what it's worth, I have seen this exact sort of behavior. Sometimes, I need to decrease the timestep, but sometimes, I just need to restart the run, and it sails through the previous crash without any problems. I usually don't even see any problems in the history file that suggests that a proble...
by hetland
Thu Oct 23, 2008 2:14 am
Forum: ROMS Discussion
Topic: "missing_value" for masked regions in NetCDF output files?
Replies: 16
Views: 3888

Re: "missing_value" for masked regions in NetCDF output files?

I can't think of a single good reason not to use _FillValue. I've wondered for years why ROMS does not do this (but have been too lazy to write the wrapper). Other models, like GETM, do this, and suffer no ill effects. And it makes quick looks in ncview much nicer. Note, the values in the model woul...
by hetland
Sun Sep 28, 2008 3:53 pm
Forum: ROMS Installation
Topic: ROMS3, g95, OS X 10.5.4, netcdf-3.6.3
Replies: 4
Views: 2846

Re: ROMS3, g95, OS X 10.5.4, netcdf-3.6.3

I would use gfortran (found here: http://r.research.att.com/tools/ ) instead of g95. I experimented with a number of gnu fortran variants on Mac OS X a while ago, and gfortran is definately the way to go. Here are some instructions for compiling all the stuff you need using ifort and gfortran: http:...
by hetland
Mon Jul 21, 2008 8:12 pm
Forum: ROMS Tools and Techniques
Topic: pyroms -- a suite of tools for working with ROMS in Python
Replies: 5
Views: 4579

Re: pyroms -- a suite of tools for working with ROMS in Python

This problem is due to the changes in the netCDF4 API. Basically, the multifile support has been folded into the regular netCDF4 object. In other words, there is no more MFNetCDF4. I would recommend that you use the newer package (that has already fixed this bug) called octant. It is very similar to...
by hetland
Wed Apr 02, 2008 6:30 am
Forum: ROMS Installation
Topic: ROMS 3.1 - OSX 10.5.2 Intel - ifort 10.1.006 - netcdf 3.6.2:
Replies: 17
Views: 6100

Usually temporary environmental variable definitions go before the command (e.g., ./configure). Look at the script to make sure it is using the right fortran compiler. Also, try a which ifort command to make sure you are using the one you think (you will have to follow a few softlinks to get to the ...
by hetland
Tue Apr 01, 2008 6:39 pm
Forum: ROMS Installation
Topic: ROMS 3.1 - OSX 10.5.2 Intel - ifort 10.1.006 - netcdf 3.6.2:
Replies: 17
Views: 6100

Versions? Which ifort, which netcdf?
by hetland
Mon Mar 31, 2008 5:48 pm
Forum: ROMS Installation
Topic: ROMS 3.1 - OSX 10.5.2 Intel - ifort 10.1.006 - netcdf 3.6.2:
Replies: 17
Views: 6100

I mean 10.1.012, of course...

I have edited the above posts so that no one will ever know I was wrong. Sneaky, eh?
by hetland
Mon Mar 31, 2008 3:31 pm
Forum: ROMS Tools and Techniques
Topic: Python build script
Replies: 1
Views: 1884

Python build script

Inspired by some of the other build scripts I have seen using either the C shell or bash, I have made a python script that does the same. I feel the advantage of this script is that the options are more flexible -- it is easier to have more logic in a real programing language than a shell script. #!...
by hetland
Mon Mar 31, 2008 3:26 pm
Forum: ROMS Installation
Topic: ROMS 3.1 - OSX 10.5.2 Intel - ifort 10.1.006 - netcdf 3.6.2:
Replies: 17
Views: 6100

For the record, here are some build notes for netcdf and openmp on Leopard. I used ifort 10.1.012 and gfortran 4.3.0. Nothing too surprising here, but I finally decided to write all this down after I kept on forgetting which options needed to be specified. Notes for making netcdf3 and netcdf4: -----...
by hetland
Mon Mar 31, 2008 8:02 am
Forum: ROMS Installation
Topic: ROMS 3.1 - OSX 10.5.2 Intel - ifort 10.1.006 - netcdf 3.6.2:
Replies: 17
Views: 6100

I had exactly the same problem. The solution appears to be to get the latest version of the compiler (10.1.012). Intel has had some pretty sever problems with Leopard, and I'm not sure whether to blame Apple or Intel. Either way, it has not been good. See: http://softwarecommunity.intel.com/isn/Comm...
by hetland
Wed Jan 30, 2008 3:39 pm
Forum: ROMS Information
Topic: Movie on installing roms (via svn) and making a project
Replies: 9
Views: 4382

OK, so Jacapo sent me the movie (that I made... but lost... now have again..). I put it back on my website; get it here:

http://pong.tamu.edu/~rob/movies/roms_project.mov
by hetland
Fri Jan 18, 2008 7:34 am
Forum: ROMS Information
Topic: Movie on installing roms (via svn) and making a project
Replies: 9
Views: 4382

Indeed it is not there. At all. I checked... This is a bit embarrassing, but I hope that somebody downloaded a copy and can send it back to me again. Or, better, put it on the myroms website. However, in the worst case, I will redo the movie. The first was done in one take, and (as you can see by He...
by hetland
Tue Jan 08, 2008 7:24 am
Forum: ROMS Tools and Techniques
Topic: A question about the gridgen
Replies: 7
Views: 3518

Running a number of separate appendage grids may have advantages in some particular cases -- in particular, a long estuary attached to a coastal sea. However, in the long run, I would rather see ROMS do something more like GETM. GETM allows you to ignore tiles that are completely land. Presently, RO...
by hetland
Sun Jan 06, 2008 9:32 pm
Forum: ROMS Tools and Techniques
Topic: A question about the gridgen
Replies: 7
Views: 3518

You need to specify values for everything where there is a NaN. Gridgen simply does not define the grid here, and they are not water points, so you just need to fill them in with values that do no harm. First, you need to set your mask to include the NaNs (and possibly other points if you have a mor...
by hetland
Wed Dec 19, 2007 9:16 am
Forum: ROMS Usage
Topic: Including a river discharge
Replies: 1
Views: 1679

The first thing to do is take a look in the analytical section of the code (now, Functionals/ana_psource.h). All of the variables that are described there need to be set. It will be more clear what each variable does when you see how it is specified in the analytical code. Perhaps the least clear fl...
by hetland
Wed Aug 29, 2007 5:48 pm
Forum: ROMS Information
Topic: Movie on installing roms (via svn) and making a project
Replies: 9
Views: 4382

Movie on installing roms (via svn) and making a project

I recently made a movie for a student on how to download roms via svn and create a project. I would like to say it is 'little' but it is 28 min and 200 megs... but perhaps you might still find it useful:

http://pong.tamu.edu/~rob/?p=98
by hetland
Tue Jul 31, 2007 5:58 pm
Forum: ROMS FAQ
Topic: Tidal forcing at the open boundary
Replies: 23
Views: 42297

The way to include tides in a painless way is to use a FLATHER/CHAPMAN combination, then specify the climatological sea surface height with tides in it. Thus your cpp options for an eastern open boundary with tides would look like this: #define EAST_FSCHAPMAN #define EAST_M2FLATHER #define EAST_M3RA...
by hetland
Sat Apr 21, 2007 3:09 am
Forum: ROMS Information
Topic: ROMS new version
Replies: 6
Views: 2960

There is not much to think about, as the release of tagged versions is done often for many open source software projects. This is how it goes: Official release versions (in the future, of course) are nothing more than snapshots of the svn repository. These are usually tagged in the svn repository wi...
by hetland
Wed Feb 28, 2007 3:47 pm
Forum: ROMS Benchmarks
Topic: ROMS 3.0 on the Altix 4700, benchmarks
Replies: 1
Views: 2786

This short suite of benchmarks was done on a cluster of 2.0 GHz Duel Opterons. The cluster has eight nodes, two processors per node, and two cores per processor for a total of 16 actual CPUs and 32 kind-of CPUs. The nodes are connected with Infiniband. The total cost of the cluster was about 40K, wi...
by hetland
Sun Feb 25, 2007 3:12 pm
Forum: ROMS Discussion
Topic: Inflow-Outflow options for reservoirs
Replies: 12
Views: 3871

Another thing to try is to move the points one further into the interior in analytical.F. Since this specifies the transport on a cell face, you will see your flux as a source of bouyancy on one side of the source, and a sink on the other. This will help you diagnose if the points are in the right p...
by hetland
Sat Feb 24, 2007 5:54 pm
Forum: ROMS Discussion
Topic: Inflow-Outflow options for reservoirs
Replies: 12
Views: 3871

Perhaps your source point indicies are one off? Make sure they are on the edge of the domain interior, and not a point off..

-r
by hetland
Fri Feb 09, 2007 3:10 pm
Forum: ROMS Installation
Topic: ROMS on an Intel Mac with a Core Duo and free compilers
Replies: 1
Views: 1850

If you hate fink as much as I do (I _really_ hate fink), you can also install gfortran very easily from the binary distribution found at

http://hpc.sourceforge.net

I have used these versions successfully for a few years now.

-r
by hetland
Thu Feb 08, 2007 3:18 pm
Forum: ROMS Tools and Techniques
Topic: pyroms -- a suite of tools for working with ROMS in Python
Replies: 5
Views: 4579

pyroms -- a suite of tools for working with ROMS in Python

NOTE: THIS SECTION USED TO BE ABOUT PYROMS, WHICH HAS CEASED DEVELOPMENT. PLEASE USE OCTANT INSTEAD. I have started to create a tool suite for working with ROMS in python, called octant. You can get the source and see other project resources at http://octant.google.com/ This package is based on pyr...
by hetland
Thu Feb 08, 2007 3:09 pm
Forum: ROMS Documentation
Topic: svn resources and documentation
Replies: 0
Views: 2364

svn resources and documentation

Given that roms_sed is maintained with subversion (as is all the stuff on google code, and newer sourceforge projects), I thought I would share the best subversion book I have found yet: http://svnbook.red-bean.com There are others that are almost as good, but this one does a good job describing wor...
by hetland
Thu Jan 25, 2007 3:24 pm
Forum: ROMS Installation
Topic: Intel Core 2 Duo processor
Replies: 14
Views: 7365

Regarding MPI on the Mac, I got this note from Enrique Rob, I can't seem to be able to log on to the forum. Anyway, I have ROMS compiled on my new Mac Quad PRO using mpich2 and g95. e. ____________________________ Enrique Curchitser Institute of Marine and Coastal Sci. 71 Dudley Rd New Brunswick, NJ...
by hetland
Wed Jan 24, 2007 10:17 pm
Forum: ROMS Problems
Topic: Necessary grid components.
Replies: 2
Views: 1378

Well, I already know that it ..works.. if you fill in pm and pn with non-zeros. I don't know if the solution is sensitive to the fill value. I suspect that it can really only come in through the viscous/diffusive terms, if at all. Rather than try the empirical approach, and see what happens with dif...
by hetland
Wed Jan 24, 2007 3:29 pm
Forum: ROMS Problems
Topic: Necessary grid components.
Replies: 2
Views: 1378

Necessary grid components.

I am working more on using Sakov's gridgen utility for ROMS grids. The output of this utility only returns 'wet' cells, that is, there is no information about the grid at all in some locations. Generally, this is not a problem -- these points would be masked anyways. However, it leads me to the ques...
by hetland
Wed Jan 24, 2007 3:23 pm
Forum: ROMS Installation
Topic: Intel Core 2 Duo processor
Replies: 14
Views: 7365

I have worked on this (a bit), and Brian Powel has worked on this more. Right now, we have both been unsuccessful at getting this to work with MPI. If you figure it out, let us know. I would start with LAMMPI, as they say they now have Mac compatible libraries available.

Good luck,
by hetland
Tue Nov 28, 2006 4:28 pm
Forum: ROMS Bugs
Topic: Parallel bug in set_data.F for point winds.
Replies: 3
Views: 2423

Parallel bug in set_data.F for point winds.

There is a parallel bug in roms-2.2 (and beyond, I think) which does not correctly read in point winds for parallel cases. Here is the fix to Nonlinear/set_data.F: ! ! If input point surface winds, rotate to curvilinear grid. ! IF (.not.Linfo(1,idUair,ng)) THEN DO j=MIN(Jstr-1,JstrR),JendR !!! chang...
by hetland
Tue Nov 21, 2006 4:37 pm
Forum: ROMS Discussion
Topic: River point source direction
Replies: 2
Views: 1796

I have also see people hack in an addition to the momentum equation one cell in from the source. That way, if you have a straight coastline, you can still have a river that comes in at an angle. However, I would not recommend this. If you really want to resolve the details of the circulation field a...
by hetland
Thu Oct 12, 2006 2:01 pm
Forum: ROMS Installation
Topic: compile error on Mac with gfortran
Replies: 6
Views: 3319

Did you install XTools? Did you get the latest version from the Developer Connection? As you probably know, gcc does not come installed by default on a Mac. Also, there have been pretty significant changes recently between gcc4.0 and gcc4.1. Make sure you are running gcc 4.1: $ gcc --version i686-ap...
by hetland
Wed Oct 11, 2006 1:54 pm
Forum: ROMS Installation
Topic: compile error on Mac with gfortran
Replies: 6
Views: 3319

I've found that I need to copy libc or libc++ to a place the linker can find. This was a particular problem when the new macs updated to gcc4. See if you can find libc, and if it's not already in /usr/lib, put a copy there.
by hetland
Tue Oct 10, 2006 1:45 pm
Forum: ROMS Usage
Topic: geostrophically adjusted initial condition.
Replies: 6
Views: 3470

If your features are only 2-3 times the deformation radius, your flow field might not be all that geostrophic. I would think about spinning up diagnostically, to get a cyclostrophic/geostrophic/QG (or whatever) balance that better matches your density field *and* the physics at the scales you are in...
by hetland
Mon Oct 09, 2006 3:59 pm
Forum: ROMS Information
Topic: ROMS Grid: what is essential?
Replies: 3
Views: 3267

Below is an ncdump -h of a minimalist grid, created with gridgen so land points are not defined, that works in roms. I did linear interpolation on a retangular grid (i.e., using gridpoint values for x and y) of all the values, so that all of the values were defined everywhere, but just not used by t...
by hetland
Thu Oct 05, 2006 3:43 pm
Forum: ROMS Problems
Topic: MPI 1 by 1 tile configuration --segmentation error
Replies: 11
Views: 4536

I have had much success on linux machines, after increasing the stack size by applying ulimit -s unlimited I think the heap is typically unlimited (or very large) by default on most linux distributions. The BENCHMARK cases are good to see how well you can deal with very large domains. Unfortunately ...
by hetland
Tue Oct 03, 2006 1:51 pm
Forum: ROMS Tools and Techniques
Topic: Grid generation
Replies: 61
Views: 34372

Gridgen with python

All I have made a gridgen wrapper for python. This is not too terribly different from using gridgen from the command line, but it does keep everything in one place. You can get it with subversion using the command svn co http://pong.tamu.edu/svn/python/trunk/gridgen gridgen You can browse the code w...
by hetland
Fri Oct 28, 2005 2:52 pm
Forum: ROMS Tools and Techniques
Topic: Python as a platform for ROMS analysis
Replies: 0
Views: 2432

Python as a platform for ROMS analysis

MATLAB is presently the defacto standard in analysis of ROMS model output. This is in part because MATLAB is the standard analysis tool for observational oceanogrpahy (although not always), but mostly because of the hurculean efforts of Rich Signell (with much help from Chuck Denham and John Evans)....
by hetland
Mon Sep 05, 2005 3:10 pm
Forum: ROMS Discussion
Topic: Simulating tides with ROMS
Replies: 2
Views: 2486

If this still works the same way it did when we first put tides in, all the _TIDES routines do is add the appropriate high-frequency signals to the low-frequency climatology signals. So, what you need to specify depends on how you set up your open boundary conditions. I would reccoment the CHAPMAN/F...
by hetland
Mon Jul 25, 2005 2:53 pm
Forum: ROMS Problems
Topic: River biology
Replies: 7
Views: 3476

You also need to make sure that you specify *all* of the other tracer variables in your river inflow values. Like, if you are running an NPZD model, you need to specify phytoplankton (et al.) as well. There may be a way to make ROMS default to a zero tracer value, but this is not the way it is set u...
by hetland
Wed Jun 22, 2005 7:31 pm
Forum: ROMS Tools and Techniques
Topic: Simple matlab tools for working with structured grids
Replies: 0
Views: 1874

Simple matlab tools for working with structured grids

I know that we have all developed little hacks in matlab to deal with structured grids, in particular, roms output. I have tried to make a few very simple, minimalist functions that do the things I find myself doing all the time. The goal is to have each function be simple enough that you do not hav...
by hetland
Wed Apr 27, 2005 3:45 pm
Forum: ROMS Source
Topic: Using a version control system for ROMS.
Replies: 4
Views: 4685

2.1 has subroutines that have been fixed (for example the wieghts for the 2D stepping) and typos (e.g., #ifdef CANUTO_A || CANUTO_B). Why not fix these things in the standard releases so that everybody doesn't need to discover these things for themselves (or reapply them each time they start from sc...
by hetland
Wed Apr 27, 2005 1:45 pm
Forum: ROMS Source
Topic: Using a version control system for ROMS.
Replies: 4
Views: 4685

Monotone vs. subversion.

Subversion seems to be the hippest CVS replacement, and many thought it should be used for the Linux source code (which used BitKeeper, until recently...). In the discussions on slashdot about replacing BitKeeper, I found some comments about monotone that seemed positive, so I checked it out and li...
by hetland
Mon Apr 25, 2005 7:55 pm
Forum: ROMS Tools and Techniques
Topic: Parameter space driver
Replies: 0
Views: 1592

Parameter space driver

Often, I find that I want to run roms over some set of parameters. The best way to do this is to make critical parameters dependent on the USER parameters in the roms.in file. I have written a program in python which does just this, automatically. The program, pspace.py , creates a series of input f...
by hetland
Mon Apr 25, 2005 7:50 pm
Forum: ROMS Source
Topic: Using a version control system for ROMS.
Replies: 4
Views: 4685

Using a version control system for ROMS.

I think we sould consider using some sort of version control system for distributing ROMS. This would be useful, because users could fix bugs in the code, and the code that we download could also be bug free. The second thing is a big deal for me. Every time I want to start over with ROMS, I need to...
by hetland
Mon Jan 31, 2005 5:01 pm
Forum: ROMS Problems
Topic: running ROMS on an OSX cluster
Replies: 5
Views: 4033

Yes, that did it. Youv'e saved me an afternoon with the debugger.

Thanks,

-Rob.
by hetland
Wed Jan 26, 2005 11:38 pm
Forum: ROMS Problems
Topic: running ROMS on an OSX cluster
Replies: 5
Views: 4033

Trouble with Mac xlf and netcdf

The trick above worked fine to get the input files read. I rebuilt netcdf with xlf (various different ways and versions), but I still get errors that DEF_VAR can't write. E.g., [...] DEF_VAR - unable to define variable: xl in NetCDF file: ocean_his.nc [...] Obviously the code compiles error free, an...
by hetland
Thu Nov 18, 2004 8:17 pm
Forum: Meetings/Workshops
Topic: Where would you like the ROMS/TOMS workshop to be held?
Replies: 1
Views: 2292

Vote for Australia

How about Australia? (Or is that covered by 'Asia'?)
by hetland
Tue Apr 06, 2004 5:02 pm
Forum: ROMS Discussion
Topic: Treatment of rivers
Replies: 20
Views: 17099

First of all, I was unclear when I said unphysical numbers. I meant local values of the outflow speed and Froude number not matching measured values but otherwise physically consistent, rather than numerical dispersion (i.e., the very high salinity). These slightly elevated values of salinity near t...
by hetland
Fri Apr 02, 2004 4:33 pm
Forum: ROMS Discussion
Topic: Treatment of rivers
Replies: 20
Views: 17099

The Mississippi Delta is like a leaky sponge, with a few cuts in it. Southwest Pass (the major outflow) is resolved with three horizontal points, (slightly less than 1 km wide each), and it carries about 1/3 of the total flow. Maximum transports are about 4000 m3/s, and velocities in the vicinity of...
by hetland
Fri Apr 02, 2004 4:13 pm
Forum: ROMS Discussion
Topic: Treatment of rivers
Replies: 20
Views: 17099

Rivers are a boundary condition within the model. As such, they should be formed just like any other boundary condition -- in a way that is numerically sound and physically consistent. That said, a river boundary condition should look like river flow, especially if this is put at the head of some so...