Search found 108 matches

by jpringle
Tue Apr 30, 2024 3:01 pm
Forum: Job Opportunities
Topic: Post-doc in the physics and observation of coastal flooding events--University of New Hampshire
Replies: 0
Views: 4297

Post-doc in the physics and observation of coastal flooding events--University of New Hampshire

As coastal sea-levels rise, flooding is increasingly impacting coastal communities. The NSF funded EPSCOR project, "Community-Driven Coastal Climate Research & Solutions for the Resilience of New England Coastal Populations," seeks to help communities manage the increasing environmental risks caused ...
by jpringle
Fri Mar 01, 2024 4:27 pm
Forum: ROMS Discussion
Topic: CPU selection for self built servers
Replies: 4
Views: 3473

Re: CPU selection for self built servers

I would still benchmark, but unless you need all those cores for other reasons, I would go with a higher clockspeed, lower core count chip. A dual CPU machine will have twice as many memory channels. The memory size, if it is large enough for your model to run, will not affect the run speed, as long ...
by jpringle
Wed Feb 28, 2024 2:49 pm
Forum: ROMS Discussion
Topic: CPU selection for self built servers
Replies: 4
Views: 3473

Re: CPU selection for self built servers

Dear Enzo -- My first advice is to benchmark. Most server manufactures will allow you to bench mark systems remotely. How does your ROMS configuration scale with cores? My second advice is not to chase cores or threads. All my experience suggests that ROMS becomes quickly limited by bandwidth to ...
by jpringle
Wed Sep 20, 2023 2:19 pm
Forum: Job Opportunities
Topic: Postdoctoral position at University of New Hampshire in the physics of sea level and modeling coastal sea-level change.
Replies: 0
Views: 19109

Postdoctoral position at University of New Hampshire in the physics of sea level and modeling coastal sea-level change.

Postdoctoral position at University of New Hampshire in the physics of sea level and modeling coastal sea-level change. As coastal sea-levels rise, flooding will increasingly impact coastal communities. The NSF funded EPSCOR project, the Community-Driven Coastal Climate Research & Solutions for the ...
by jpringle
Sun Mar 19, 2023 3:00 am
Forum: Ocean News
Topic: Dale B. Haidvogel (1949-2023)
Replies: 18
Views: 62137

Re: Dale B. Haidvogel (1949-2023)

I remember meeting him at a conference at Rutgers, after which he invited everyone to his house. He was one of the first "big names" that I met, and realized that they were not just names on papers but actual people -- and idea that was still new to me, and a bit of a shock. He was very kind in ...
by jpringle
Tue Aug 30, 2022 12:09 pm
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 39033

Re: Consumer-grade desktop for toy-"operational" model running 24/7?

Just a quick callout to a past post of mine https://www.myroms.org/forum/viewtopic.php?t=4665 on performance versus memory bandwidth. I still find this to be true with the latest AMD hardware for ROMS -- The 16 core AMD 5950x is a beast for everyday work, but with its two channels of memory I can't ...
by jpringle
Mon Oct 25, 2021 3:37 pm
Forum: ROMS Discussion
Topic: relative/potential vorticity outputs
Replies: 7
Views: 19273

Re: relative/potential vorticity outputs

No problem. Note this is a straight translation from the code, accelerated via numba. Also note the many ways you can define Ertel's PV depending on how you define the conserved quantity... One note -- this code is accelerated with Numba; you need not do so, but it will be unbearably slow. If found ...
by jpringle
Tue Oct 12, 2021 1:42 pm
Forum: ROMS Discussion
Topic: relative/potential vorticity outputs
Replies: 7
Views: 19273

Re: relative/potential vorticity outputs

I have a python code to compute these from history file that is based on the ROMS Fortran code; ask if interested.
Jamie
by jpringle
Tue Oct 27, 2020 7:44 pm
Forum: Job Opportunities
Topic: Postdoctoral Fellowship at University of New Hampshire/Diversity emphasis
Replies: 0
Views: 2736

Postdoctoral Fellowship at University of New Hampshire/Diversity emphasis

The Department of Earth Sciences seeks a postdoctoral scholar hire in the field of Earth and Environmental Sciences. We will consider applicants from a broad range of disciplines in the geosciences (e.g., biogeochemistry, remote sensing, geochemistry, earth system modeling, surficial geology ...
by jpringle
Tue Jun 16, 2020 9:39 pm
Forum: ROMS Discussion
Topic: units of omega
Replies: 6
Views: 12518

Re: units of omega

Thanks
by jpringle
Tue Jun 16, 2020 8:59 pm
Forum: ROMS Discussion
Topic: units of omega
Replies: 6
Views: 12518

Re: units of omega

Hernan --

I agree with all that. All I am saying is that varinfo.dat should be changed to be consistent with the actual units written out.

Jamie
by jpringle
Tue Jun 16, 2020 7:37 pm
Forum: ROMS Discussion
Topic: units of omega
Replies: 6
Views: 12518

Re: units of omega

Yes, I agree. All that is left is to change the text in varinfo.dat to reflect this.
Jamie
by jpringle
Tue Jun 16, 2020 6:25 pm
Forum: ROMS Discussion
Topic: units of omega
Replies: 6
Views: 12518

units of omega

Dear all -- I, and some Europeans doing particle tracking with ROMS output, are confused by what the units of omega are supposed to be in the output. From equation 12 of https://www.myroms.org/wiki/Terrain-Following_Coordinate_Transformation , the units of omega should be w/Hz or 1/time since Hz ...
by jpringle
Thu Jun 11, 2020 1:58 pm
Forum: ROMS Bugs
Topic: bug in calculating average potential vorticity
Replies: 0
Views: 3081

bug in calculating average potential vorticity

Dear all -- I think I have found a bug in how the potential vorticity calculation is computed. I was trying to make the output match my expectations for a motionless, stratified f-plane ocean, and found this (I think) error. In set_avg.F, the function that calculates the potential vorticity is ...
by jpringle
Tue Jun 02, 2020 3:20 pm
Forum: ROMS Discussion
Topic: calculate potential vorticity flux J
Replies: 0
Views: 2716

calculate potential vorticity flux J

Dear all -- Does anyone know of a package to calculate the potential vorticity flux J either inline while the model is running, or from the saved output? I would prefer python, but would be happy with anything. I would prefer not to have to re-invent the wheel, especially if that wheel involves ...
by jpringle
Fri Apr 17, 2020 2:01 pm
Forum: ROMS Discussion
Topic: Request: Numerical modeling course syllabi
Replies: 11
Views: 11694

Re: Request: Numerical modeling course syllabi

Just to close this discussion. I ended up using Lars Roed's book, referenced above. I found it a very good fit to the class, and was very happy with it. It worked well with my approach of having students build a series of simple numerical models as various dynamics were introduced in Lars' book. The ...
by jpringle
Tue Jan 14, 2020 6:58 pm
Forum: ROMS Discussion
Topic: faster MPI on single computer / -bind-to core /Ubuntu 18.04
Replies: 1
Views: 3244

Re: faster MPI on single computer / -bind-to core /Ubuntu 18.04

With the same version of Ubuntu, but updated kernels (nothing special, just the default updates), this no longer appears to be true. A reminder to test your tiling and MPI options every so often.

Jamie
by jpringle
Tue Jan 07, 2020 7:39 pm
Forum: ROMS Discussion
Topic: Request: Numerical modeling course syllabi
Replies: 11
Views: 11694

Re: Request: Numerical modeling course syllabi

Lar's --

Thanks for bringing up your book. I like it very much.

Jamie
by jpringle
Mon Dec 16, 2019 2:35 pm
Forum: ROMS Discussion
Topic: Request: Numerical modeling course syllabi
Replies: 11
Views: 11694

Re: Request: Numerical modeling course syllabi

Thanks everyone!
Jamie
by jpringle
Mon Dec 16, 2019 2:34 pm
Forum: ROMS Installation
Topic: ROMS using Windows Subsystem for Linux
Replies: 3
Views: 3940

Re: ROMS using Windows Subsystem for Linux

I have used windows subsytem with linux with Ubuntu 18.04, and it worked just like ubuntu on bare metal. have you tried your problem on regular linux box?

Jamie
by jpringle
Wed Dec 11, 2019 8:50 pm
Forum: ROMS Discussion
Topic: Request: Numerical modeling course syllabi
Replies: 11
Views: 11694

Request: Numerical modeling course syllabi

Dear modelers -- I am teaching an advanced graduate course in numerical ocean modeling in the spring, and thus must call it into existence. Students will have had fluids and PDE's and introductory physical oceanography. Most will have had GFD. I would love to see your syllabi if you have taught such ...
by jpringle
Wed Dec 04, 2019 2:33 pm
Forum: ROMS Tools and Techniques
Topic: netCDF browser for remote work
Replies: 7
Views: 9331

Re: netCDF browser for remote work

A much less elegant solution, but one I use all the time, is XPRA -- it is a way to view X windows remotely that is MUCH MORE robust than forwarding a remote X11 connection via SSH. Essentially, it sets up a local X11 server, and then serves the windows that are created on that server via a screen ...
by jpringle
Mon Jul 08, 2019 8:23 am
Forum: ROMS Discussion
Topic: Alongshore currents with periodic alongshore boundaries
Replies: 9
Views: 8336

Re: Alongshore currents with periodic alongshore boundaries

Aryan -- I just want to amplify a point that Sasha made above, for it is important to understanding the dynamics of your case. An alongshore current in a periodic domain must be maintained by forcing; anything else is inconsistent. Why? what is an alongshore periodic domain? one in which there are ...
by jpringle
Fri Apr 26, 2019 7:04 pm
Forum: ROMS Discussion
Topic: example lightweight isosurface plotting code in python
Replies: 0
Views: 2280

example lightweight isosurface plotting code in python

Dear all -- I hope someone finds this usefull. I wanted to plot some iso-surfaces in 3D with ROMS data without re-gridding it to a rectangular grid or some such silliness. There are many heavy-weight environments that can do all sorts of cool things like that (paraview, vapor, yt) -- but getting ...
by jpringle
Tue May 15, 2018 11:57 pm
Forum: ROMS Discussion
Topic: cold water creeping upwards along topography
Replies: 3
Views: 3513

Re: cold water creeping upwards along topography

In your sections, the bathymetry has severe stair-stepping (discontinuities). Is that real?

Also, is there any along-isobath flow which could be driving a bottom Ekman transport?

Jamie
by jpringle
Fri Apr 27, 2018 5:37 pm
Forum: ROMS Discussion
Topic: faster MPI on single computer / -bind-to core /Ubuntu 18.04
Replies: 1
Views: 3244

faster MPI on single computer / -bind-to core /Ubuntu 18.04

Dear all -- I hope some of you find this useful. I am running ROMS on a dual processor 32 core AMD EPYC machine. I noticed on the latest version of Ubuntu my smaller machines showed good improvement with ROMS runs with openMP parallelization. So I upgraded the big machine to Ubuntu 18.04, and my ...
by jpringle
Wed Mar 14, 2018 3:20 pm
Forum: ROMS Tools and Techniques
Topic: Error of installing PyCNAL on conda
Replies: 11
Views: 12211

Re: Error of installing PyCNAL on conda

Dear All, The errors of the type Warning: distutils distribution has been initialized, it may be too late to add an extension _interp Warning: distutils distribution has been initialized, it may be too late to add an extension _obs_interp are a bit of a red herring. You can make the codes work by ...
by jpringle
Tue Feb 27, 2018 5:17 pm
Forum: ROMS Problems
Topic: excess mixing in k-epsilon GLS mixing; bug or misconfigured?
Replies: 8
Views: 11162

Re: excess mixing in k-epsilon GLS mixing; bug or misconfigu

Thanks to John Wilkin, John Warner, Nirnimesh Kumar, the excellent Warner et al. 2005 Ocean Modeling paper on mixing schemes in ROMS and the excellent note by Malcolm Scully (https://www.myroms.org/forum/viewtopic.php?f=17&t=2310&p=8441&hilit=malcolm+scully#p8458 ), I have run this issue down. In ...
by jpringle
Mon Feb 26, 2018 2:12 pm
Forum: ROMS Problems
Topic: excess mixing in k-epsilon GLS mixing; bug or misconfigured?
Replies: 8
Views: 11162

Re: excess mixing in k-epsilon GLS mixing; bug or misconfigu

Both in John Wilkin's reply and in some private emails, I have received some good hints related to minimum turbulence level parameters in the mixing scheme. I shall investigate and compare to the primary literature for this scheme in the next few days, and will post a summary here soon. Thanks for ...
by jpringle
Sat Feb 24, 2018 10:48 pm
Forum: ROMS Problems
Topic: excess mixing in k-epsilon GLS mixing; bug or misconfigured?
Replies: 8
Views: 11162

Re: excess mixing in k-epsilon GLS mixing; bug or misconfigu

ALT_BAK is 1.0d-6 for both T and S, so much less than observed... ¯\_(ツ)_/¯
by jpringle
Fri Feb 23, 2018 7:52 pm
Forum: ROMS Problems
Topic: excess mixing in k-epsilon GLS mixing; bug or misconfigured?
Replies: 8
Views: 11162

Re: excess mixing in k-epsilon GLS mixing; bug or misconfigu

I was asked offline why I expected the diffusion of temperature to be slower than ROMS's k-epsilon in the ocean -- why did I think the GOTM k-epsilon, ROMS MY25 and ROMS LMD was closer to the truth? In this stratified ocean, with no horizontal gradients of density and no forcing, there is no flow ...
by jpringle
Thu Feb 22, 2018 7:17 pm
Forum: ROMS Problems
Topic: excess mixing in k-epsilon GLS mixing; bug or misconfigured?
Replies: 8
Views: 11162

excess mixing in k-epsilon GLS mixing; bug or misconfigured?

Dear all -- I think I have stumbled across evidence of a bug in the GLS vertical mixing scheme in ROMS, at least for the k-epsilon scheme. It remains possible that I am just doing something stupid. I find that under the k-epsilon scheme vertical stratification diffuses away much more quickly than ...
by jpringle
Wed Jan 31, 2018 9:26 pm
Forum: ROMS Bugs
Topic: bug in ana_psource with openMP
Replies: 0
Views: 2943

bug in ana_psource with openMP

I am not well versed in openMP programming -- but I think I found a bug where ana_psource initializes the tracer input. I was having problems with different runs of my code producing wildly different answers from run to run under openMP. It did not do this when run in a single thread and when run ...
by jpringle
Wed Jan 31, 2018 4:37 pm
Forum: ROMS Bugs
Topic: bug in ana_psource.h with analytical initial conditions
Replies: 0
Views: 2746

bug in ana_psource.h with analytical initial conditions

I think I found a bug in the ana_psource.h when the model is initialized with ana_initial.h. In this case, the analytical forcing files are invoked in analytical.F before the first time step of integration as part of initializing the u, v etc fields. When this is done, the variable iic(ng) is 0. In ...
by jpringle
Fri Dec 15, 2017 5:45 pm
Forum: ROMS Problems
Topic: Compiling ROMs in Ubuntu
Replies: 7
Views: 9359

Re: Compiling ROMs in Ubuntu

Mark --

Ah -- I see -- I had not defined "USE_MY_LIBS" and so it ignored the changes I made to NF_CONFIG and NETCDF_INCDIR in my build.bash. I had made the changes, and never tested if it worked without them.


Cheers,
Jamie
by jpringle
Thu Dec 14, 2017 2:21 pm
Forum: ROMS Problems
Topic: Compiling ROMs in Ubuntu
Replies: 7
Views: 9359

Re: Compiling ROMs in Ubuntu

I have recently compiled ROMS on a bunch of ubuntu machines for benchmarking and production use, with versions 16.04, 17.04 and 17.10 (I recommend skipping 17.04). In all cases it is trivial to get ROMS going. Just install the following packages (which include the MPI packages, and you can compile a ...
by jpringle
Wed Dec 06, 2017 7:46 pm
Forum: ROMS Benchmarks
Topic: Benchmarking Epyc, Ryzen, and Xeon: Tyranny of Memory
Replies: 5
Views: 16864

Re: Benchmarking Epyc, Ryzen, and Xeon: Tyranny of Memory

Wayne -- A couple of things. To get the model to run on my Threadripper 1950x with only 8 threads, I used the environmental variable omp_num_threads. I suspect, based on my earlier post, that the speed of ROMS on this kind of system is limited by memory bandwidth (see Sasha's long thread on this ...
by jpringle
Mon Dec 04, 2017 12:44 am
Forum: ROMS Benchmarks
Topic: Benchmarking Epyc, Ryzen, and Xeon: Tyranny of Memory
Replies: 5
Views: 16864

Re: Benchmarking Epyc, Ryzen, and Xeon: Tyranny of Memory

Wayne-- This configuration seems reasonable -- but I have not experimented with the intel compilers yet. At some point I will have one of my students install the student version and let people know. If you would like to benchmark my threadripper system, let me know via email. On the 1920x (12 cores ...
by jpringle
Sat Oct 07, 2017 1:35 pm
Forum: ROMS Problems
Topic: Docker container for ROMS?
Replies: 7
Views: 11059

Re: Docker container for ROMS?

A naive question. For multi-threaded or multi-process programs like ROMS, is there any significant performance penalty for using something like docker?

Jamie
by jpringle
Tue Sep 12, 2017 6:14 pm
Forum: ROMS Benchmarks
Topic: Benchmarking Epyc, Ryzen, and Xeon: Tyranny of Memory
Replies: 5
Views: 16864

Benchmarking Epyc, Ryzen, and Xeon: Tyranny of Memory

Dear all -- I hope some of you find this useful; I have been researching purchasing a new modeling server, and the results should be interesting to those trying to optimize the performance on their own machines, or are thinking of purchasing new machines. There are a number of new architectures that ...
by jpringle
Wed Jul 26, 2017 2:37 pm
Forum: ROMS Benchmarks
Topic: Ryzen 1800x Benchmarking and memory bandwidth
Replies: 2
Views: 10670

Re: Ryzen 1800x Benchmarking and memory bandwidth

Sacha and All-- A couple of things. First, the AMD Ryzen has a dual memory channels, so is has less memory bandwidth than the x99 Intel boards. Compared to older AMD processors, the L1 and L2 cache's are private to each core, but the L3 cache is shared amoung 4 cores. But see below for more on this ...
by jpringle
Tue Jul 25, 2017 9:43 pm
Forum: ROMS Benchmarks
Topic: Ryzen 1800x Benchmarking and memory bandwidth
Replies: 2
Views: 10670

Ryzen 1800x Benchmarking and memory bandwidth

Dear all -- IF YOU ARE INTERESTED IN THIS POST, PLEASE ALSO SEE: https://www.myroms.org/forum/viewtopic.php?f=29&t=4665 I hope some of you find this helpful. I needed a computer for a new graduate student to run models and do data analysis. For 1500$ I built a Ryzen 1800x box. Since this is a new ...
by jpringle
Tue May 16, 2017 2:02 pm
Forum: ROMS Discussion
Topic: Cheap Ryzen modeling setup
Replies: 0
Views: 2277

Cheap Ryzen modeling setup

Dear all -- I need to build/acquire a linux box for a new graduate student for process oriented ROMS modeling and ocean model data analysis, with a side of embarrassingly parallel ecological modeling. I have a budget of 1500 (excluding the massive pile of disks for data storage). I am wondering if ...
by jpringle
Mon Nov 21, 2016 8:37 pm
Forum: Job Opportunities
Topic: Funding for graduate student
Replies: 0
Views: 4554

Funding for graduate student

I have funding available for graduate students at the University of New Hampshire to study the interaction of Western Boundary Currents and continental shelf flows. The applicant must posses strong quantiative skills, with a major in the sciences, mathematics or engineering. Applicants with either ...
by jpringle
Thu Sep 15, 2016 1:48 pm
Forum: ROMS Discussion
Topic: Full vorticity equations
Replies: 5
Views: 8615

Re: Full vorticity equations

I do not know if this helps, but attached is a calculation of relative and potential vorticity that exactly parallels the code in Vorticity.F, but implemented in python. It calculates both 3D and shallow water vorticity. Please look at the function defineVorticity() for the calculation. I have ...
by jpringle
Wed Dec 02, 2015 12:49 pm
Forum: ROMS Bugs
Topic: issue with coarse2fine.m
Replies: 4
Views: 10185

Re: issue with coarse2fine.m

It would be easier to debug if you printed out what the attribute name is. You could add an fprint by line 112 of nc_create, and print out what aname is. This might provide insight.

Cheers,
Jamie
by jpringle
Thu Oct 22, 2015 3:25 pm
Forum: ROMS Discussion
Topic: Rivers in nested grids
Replies: 2
Views: 3412

Re: Rivers in nested grids

While waiting for some other runs to complete, I did some simple test cases with nested grids, and can confirm that what John said is true.

Jamie
by jpringle
Mon Oct 12, 2015 12:37 am
Forum: ROMS Discussion
Topic: Rivers in nested grids
Replies: 2
Views: 3412

Rivers in nested grids

Dear all -- I am having some trouble with nesting and river inflows. I have justed started to implement rivers entering into my domain. Some rivers enter where the fine grid overlaps the coarse grid. When I run the model with nesting, it blows up quickly. When I run with rivers in just coarse grid ...
by jpringle
Wed Sep 09, 2015 8:51 pm
Forum: ROMS Problems
Topic: Bug in large spherical nested run?
Replies: 5
Views: 6654

Re: Bug in large spherical nested run?

Just to update -- the discrepancy in the area of grid cells between the fine and coarse grids was largely caused by an inconsistency between mat2roms_mw and coarse2fine. The former used sw_dist to compute distance between grid points, the later used gcircle. sw_dist uses the "plain sailing" method ...
by jpringle
Wed Sep 09, 2015 1:45 pm
Forum: ROMS Problems
Topic: Bug in large spherical nested run?
Replies: 5
Views: 6654

Re: Bug in large spherical nested run?

The area of the 9 rho cells of the fine grid is 0.144% greater than the area of the associated coarse grid rho cell, with a very weak dependency on latitude. This is of the same order of magnitude of the discrepancy of mass flux at the eastern boundary. I am investigating the origin of this issue. I ...
by jpringle
Tue Sep 08, 2015 7:29 pm
Forum: ROMS Problems
Topic: Bug in large spherical nested run?
Replies: 5
Views: 6654

Re: Bug in large spherical nested run?

Hernan -- I will compare the averaged metrics tomorrow, when I have time to code. A question and two points first. 1) It is a curvilinear spherical grids -- the delta-latitude and delta-longitudes are both constant, and the grid is north/south and east/west aligned. 2) There is no wind, so the wind ...
by jpringle
Tue Sep 08, 2015 4:06 pm
Forum: ROMS Bugs
Topic: bug in check_massflux in nesting.F
Replies: 2
Views: 4751

Re: bug in check_massflux in nesting.F

Hernan -- I have posted the information you want in viewtopic.php?f=17&t=3957 ; I did not put it as a bug yet since I have not isolated the exact cause.

Jamie
by jpringle
Tue Sep 08, 2015 4:04 pm
Forum: ROMS Problems
Topic: Bug in large spherical nested run?
Replies: 5
Views: 6654

Bug in large spherical nested run?

NEGLECT THIS MESSAGE -- THE PROBLEM WAS NOT WHERE I SUSPECTED IT WAS HERE. THERE IS SOME USEFUL INFORMATION BELOW ABOUT AREA CALCULATIONS WITH VARIOUS GRID PACKAGES, SO I AM NOT DELETING IT. I have been debugging on and off for several weeks what appears to be a bug (or my stupidity) in nesting in ...
by jpringle
Sat Aug 29, 2015 7:48 am
Forum: ROMS Bugs
Topic: bug in check_massflux in nesting.F
Replies: 2
Views: 4751

bug in check_massflux in nesting.F

In check_massflux in nesting.F, the variable ng in # ifdef NESTING_DEBUG IF (DOMAIN(ng)%SouthWest_Test(tile)) THEN IF (Master) THEN WRITE (300,10) 'Western Boundary Mass Fluxes: ', & & cr, dg, rg, iif(rg), iic(rg), INT(time(rg)) END IF END IF ! # endif should be "ngf". "ng" is not defined in this ...
by jpringle
Tue Aug 18, 2015 2:53 pm
Forum: ROMS Bugs
Topic: Bug in nesting when MASKING not defined
Replies: 1
Views: 3478

Bug in nesting when MASKING not defined

In the latest ROMS, the model does not compile when NESTING is defined but MASKING is not. The error is minor; mask_hweights is always defined at the top of nesting.F even when MASKING is on. The fix is straightforward: *** nesting.F_old 2015-08-18 10:49:28.000000000 -0400 --- nesting.F 2015-08-18 ...
by jpringle
Mon Aug 17, 2015 9:53 pm
Forum: ROMS Bugs
Topic: Bug in fine2coarse.m for nesting
Replies: 6
Views: 7472

Re: Bug in fine2coarse.m for nesting

My original edits to fine2coarse.m, and the "official" version on the Subversion repository, contain an error in which the curvilinear metrics dmde and dndx are not written to the coarse grid. Since those fields are defined in c_grid(), the grid file has fill values for those metrics. The fixed code ...
by jpringle
Mon Aug 17, 2015 5:46 pm
Forum: ROMS Discussion
Topic: mat2roms with modern matlab
Replies: 1
Views: 2627

mat2roms with modern matlab

Dear all -- I would like to create a ROMS grid file from a grid of h at longitude and latitude points. In the past this was possible with mat2roms_rps.m ( https://github.com/rsignell-usgs/seagrid/blob/master/seagrid/mat2roms_rps.m ) Unfortunately, the code I can find on the web for this uses Chuck ...
by jpringle
Thu Aug 13, 2015 5:40 pm
Forum: ROMS Bugs
Topic: contact.m issue with fine grid on border of coarse grid
Replies: 2
Views: 5214

Re: contact.m issue with fine grid on border of coarse grid

Ah... I noticed in an email from John Wilkins that he told me We are taking this a step at a time. As you say, presently the Matlab tools have been written assuming the nest is wholly within the parent, but the ROMS design anticipates the day that the contact file can allow a nest to cross the ...
by jpringle
Wed Aug 12, 2015 6:33 pm
Forum: ROMS Bugs
Topic: contact.m issue with fine grid on border of coarse grid
Replies: 2
Views: 5214

contact.m issue with fine grid on border of coarse grid

Dear all -- I just want to document a bug in contact.m, and to ask if anyone has succeed in putting a refined grid next to a model boundary? I am using the latest version of the matlab tools, from SVN revision 722. I create a coarse resolution grid, which I have attached to this post. I then create ...
by jpringle
Tue Jul 21, 2015 10:30 am
Forum: ROMS Bugs
Topic: Bug in fine2coarse.m for nesting
Replies: 6
Views: 7472

Re: Bug in fine2coarse.m for nesting

My fine2coarse code does not set either the coarse_factor or refine_factor attributes of either grid; see https://www.myroms.org/forum/viewtopic.php?f=1&t=3698&p=14962#p14962 ; I am not sure fine2coarse is a good place to do so, and will not be able to do so until I have a better understanding of ...
by jpringle
Tue Jul 21, 2015 7:50 am
Forum: ROMS Usage
Topic: How to generate and use composite grids?
Replies: 16
Views: 16714

Re: How to generate and use composite grids?

Hernan & Kate -- Kate -- to start to get things to work, I have had to add the refine_factor attributes by hand. My efforts have been interrupted by my vacation, but will pick up in two weeks or so. I am not sure if I will end up going fine2coarse or coarse2fine, then edit the bathymetry in the ...
by jpringle
Thu Jul 02, 2015 6:02 pm
Forum: ROMS Bugs
Topic: Bug in fine2coarse.m for nesting
Replies: 6
Views: 7472

Re: Bug in fine2coarse.m for nesting

Attached is a corrected fine2coarse.m, a testing code, a file with grid information for it to read, and a little routine to plot grids. The testing harness expects the rest of the ROMS Matlab tools to be in the path. Please make sure you are running the provided fine2coarse.m, and not the one that ...
by jpringle
Thu Jul 02, 2015 3:28 pm
Forum: ROMS Bugs
Topic: Bug in ana_grid.h for BENCHMARK
Replies: 0
Views: 3277

Bug in ana_grid.h for BENCHMARK

There is a bug in the definition of lonp,latp,lonr,latr,lonu,latu and lonv,latv in ana_grid.h for the benchmark case. It does not affect the solution, since the metric definitions are correct. However, it can make interpreting plots difficult, and will cause errors if this analytical grid is then ...
by jpringle
Thu Jul 02, 2015 3:25 pm
Forum: ROMS Bugs
Topic: Bug in fine2coarse.m for nesting
Replies: 6
Views: 7472

Re: Bug in fine2coarse.m for nesting

Before anyone else spends time on this, I am working on fixing this code. I will let people know when a working version exists and has been tested.

Jamie
by jpringle
Tue Jun 30, 2015 9:31 pm
Forum: ROMS Bugs
Topic: Bug in fine2coarse.m for nesting
Replies: 6
Views: 7472

Bug in fine2coarse.m for nesting

I am developing a new modeling setup with a simple nested grid configuration. Following Hernan's advice, I am starting from a full domain fine grid and extracting the coarse grid from it. However, the fine2coarse.m downloaded from the SVN repository as of this morning is completely broken (SVN ...
by jpringle
Fri May 22, 2015 3:52 pm
Forum: ROMS Bugs
Topic: bug in ana_dqdsst.F
Replies: 6
Views: 8347

Re: bug in ana_dqdsst.F

Agree.
by jpringle
Fri May 22, 2015 1:23 pm
Forum: ROMS Bugs
Topic: bug in ana_dqdsst.F
Replies: 6
Views: 8347

Re: bug in ana_dqdsst.F

John -- Thanks. You are right, the inclusion of the surface layer depth in the formulation is deeply problamatic. However, it is a science issue not a programming issue. In this code we prescribe the surface heat flux as a function of SST. Thus the time rate of change of temperature is set by the ...
by jpringle
Thu May 21, 2015 1:46 pm
Forum: ROMS Bugs
Topic: bug in ana_dqdsst.F
Replies: 6
Views: 8347

Re: bug in ana_dqdsst.F

Dear Zhang -- You are correct. My code is correct, but the comments are wrong. The comments should be changed to ! !----------------------------------------------------------------------- ! Set surface heat flux sensitivity to SST (m/s ...
by jpringle
Wed May 20, 2015 3:58 pm
Forum: ROMS Bugs
Topic: bug in ana_dqdsst.F
Replies: 6
Views: 8347

bug in ana_dqdsst.F

Dear Hernan et al -- I think I have found a bug in ana_dqdsst.F. The stub code from the most recent model version on subversion reads: ! !----------------------------------------------------------------------- ! Set surface heat flux sensitivity to SST (m/s/degC ...
by jpringle
Tue May 12, 2015 3:39 pm
Forum: ROMS Installation
Topic: OSX ./build.bash problems (upwelling test case).
Replies: 15
Views: 19015

Re: OSX ./build.bash problems (upwelling test case).

Old thread, but I hope this helps someone. You can run ROMS on OS X/mac using macports to install all the libraries. But was noted above, the netcdf is a bit tricky. To get it to work, you must specify the variant of netcdf-fortran to install that is consistent with your fortran compiler. To do this ...
by jpringle
Sat May 09, 2015 3:29 pm
Forum: ROMS Discussion
Topic: recent compiler benchmarks 2015
Replies: 0
Views: 2132

recent compiler benchmarks 2015

Dear ROMS runners -- I am starting a ROMS project after a hiatus. Has anyone benchmarked compilers for Intel architectures recently (i.e. gfortran 4.9 or 5, and recent Intel compilers)? On other codes of mine, the gfortran compiler has become much faster. I wonder if the same is true for ROMS? It ...
by jpringle
Mon Apr 20, 2015 1:14 pm
Forum: ROMS Discussion
Topic: Problems with unforced model -- unrealistic circulation
Replies: 9
Views: 8283

Re: Problems with unforced model -- unrealistic circulation

When you are testing for the effects of pressure gradient error and topography, I would set all the boundaries to closed. I am worried about how the currents are reacting to both the northern and southern boundaries. I would not be surprised if the currents were substantially less if the boundaries ...
by jpringle
Wed Jan 21, 2015 7:04 pm
Forum: Job Opportunities
Topic: Student funding available, University of New Hampshire
Replies: 0
Views: 4227

Student funding available, University of New Hampshire

I have funding available for graduate students at the University of New Hampshire to study the interaction of Western Boundary Currents and continental shelf flows. The applicant must posses strong quantiative skills, with a major in the sciences, mathematics or engineering. Applicants with either ...
by jpringle
Thu May 29, 2014 8:31 pm
Forum: ROMS Problems
Topic: ocean_frc.nc and KPP mixing
Replies: 1
Views: 2447

ocean_frc.nc and KPP mixing

I must be doing something silly, but it is escaping me. I have just turned on the LMD_* options to implement KPP mixing. I have ANA_STFLUX and ANA_SMFLUX defined, which should give me tracer fluxes and momentum fluxes. Temperature is the only tracer (there is no salinity). When I turned on the LMD ...
by jpringle
Fri May 24, 2013 3:09 pm
Forum: ROMS Discussion
Topic: Particle tracking, and differences between AGRIF and Rutgers
Replies: 2
Views: 4005

Re: Particle tracking, and differences between AGRIF and Rut

As always, thanks Kate! I will look at tracmass and see how easy it is to incorporate behavior and vertical mixing. LTRANS will require some recoding to recreate w or omega from u and v ARIANE has a series of weird coordinates to specify time and inititial condition that hurt my brain. If anyone ...
by jpringle
Thu May 23, 2013 8:32 pm
Forum: ROMS Discussion
Topic: Particle tracking, and differences between AGRIF and Rutgers
Replies: 2
Views: 4005

Particle tracking, and differences between AGRIF and Rutgers

Colleagues -- I am trying to do some particle tracking in ROMS output; the model runs were not done by me, and were made with ROMS-AGRIF V3. I have several questions here; I am also asking them on the AGRIF forum. 1) Has anyone used LTRANS v2 with AGRIF? 2) Is the vertical stretching and transform ...
by jpringle
Mon Dec 17, 2012 12:50 am
Forum: Job Opportunities
Topic: PhD student position at National University of Singapore
Replies: 1
Views: 3927

Re: PhD student position at National University of Singapore

This is the middle of the Christmas vacation season in most of Europe and the Americas -- perhaps reposting this in a month or two might get more of a response if you are interested in students from here...

cheers,
Jamie
by jpringle
Tue Oct 25, 2011 1:12 pm
Forum: ROMS Discussion
Topic: NONLIN_EOS in upwelling example creates velocities
Replies: 6
Views: 7161

Re: NONLIN_EOS in upwelling example creates velocities

This comment is orthogonal to Sasha's comment -- and his certainly the more interesting of the two. However, it is worth noting that if you have stratification over a sloping bottom, in the presence of vertical mixing, you will have an adjustment of the density near the bottom. If the mixing were ...
by jpringle
Sun Jul 31, 2011 6:23 pm
Forum: ROMS Installation
Topic: Netcdf 3.6.3 problem
Replies: 7
Views: 8521

Re: Netcdf 3.6.3 problem

Dear Linda--

It appears that it passes the fortran tests, and fails the C++ tests. Since you only need the fortran libraries for ROMS, it seems like you should be good to go once you install it.

Cheers,
Jamie Pringle
by jpringle
Wed May 11, 2011 2:28 pm
Forum: ROMS Problems
Topic: Poor-man's super-(mini?)-computer revisited in 2011? Advice?
Replies: 0
Views: 2255

Poor-man's super-(mini?)-computer revisited in 2011? Advice?

Dear Colleagues -- For the last 6 years, my ocean modeling needs have been well met by a 6 node cluster; each node a dual Opteron-250 with a single core for a total of 12 cores. This system has enough power to meet my modest process-oriented modeling needs; typical runs for me have a grid size of ...
by jpringle
Tue Mar 08, 2011 2:11 pm
Forum: ROMS Sediment
Topic: diagnosing vertical coordinate from history/grid files
Replies: 8
Views: 13517

Re: diagnosing vertical coordinate from history/grid files

Kate, Sasha & Hernan--

Thanks for your replies, they have been very helpful.

Cheers,
Jamie
by jpringle
Mon Mar 07, 2011 4:42 pm
Forum: ROMS Sediment
Topic: diagnosing vertical coordinate from history/grid files
Replies: 8
Views: 13517

diagnosing vertical coordinate from history/grid files

Hi all -- I have a simple question; I fear perhaps the answer is not so simple. I am involved in a couple of efforts to read in ROMS history and grid files, and either plot them or compute other products from them (e.g. Lagrangian particle paths). I have several questions 1) What is the best way to ...
by jpringle
Thu Aug 14, 2008 6:32 pm
Forum: ROMS Information
Topic: ROMS, the Cell Processor, and other exotic architecture
Replies: 1
Views: 3905

ROMS, the Cell Processor, and other exotic architecture

Hi all- Several people at my university are trying to purchase a computer cluster with either IBM's Cell Processor (same as in the Sony Playstation) or with many GPUs to use as co-processors. Has anyone tried to modify the ROMS code to take advantage of these exotic (for now) architectures, and what ...
by jpringle
Sat May 17, 2008 2:48 pm
Forum: ROMS Discussion
Topic: Left side of the upwelling example
Replies: 3
Views: 4960

Carl- The upwelling code is define to use a constant value of vertical mixing, as you can tell from this code in upwelling.h #if defined GLS_MIXING || defined MY25_MIXING # define KANTHA_CLAYSON # define N2S2_HORAVG #else # define ANA_VMIX #endif and by noting that neither GLS_MIXING or MY25_MIXING ...
by jpringle
Mon Jan 07, 2008 1:38 am
Forum: ROMS Tools and Techniques
Topic: A question about the gridgen
Replies: 7
Views: 9009

As evidenced by my silly answer above, I had not run into this problem with the simple domains of my application, so I had not ever worried how to solve this problem. I have thus never implemented a solution. When I have to, I will, but if someone beats me to it, I would be happy to incorporate ...
by jpringle
Sat Jan 05, 2008 3:05 am
Forum: ROMS Tools and Techniques
Topic: A question about the gridgen
Replies: 7
Views: 9009

I don't know if this is the issue, but when I was testing Pavel's gridgen, it occasionally threw up nan's if my domain was perfectly rectangular, but if I moved one of the corners by even an infinitesimal amount, the nan's went away. This was a flakey problem, that only occured for certain ...
by jpringle
Sun Sep 09, 2007 4:15 pm
Forum: ROMS Problems
Topic: How to increase ROMS parallel efficiency
Replies: 5
Views: 5878

Roughly, the efficiency of the parallelization is reduced by the following * Slow interconnects *Small tiles, so that there is a lot of communications (which scales as the perimeter of tile) for each bit of internal computation in the tile (which scales as the number of gridpoints in the tile). Your ...
by jpringle
Fri Sep 07, 2007 1:50 am
Forum: ROMS Problems
Topic: ROMS2.3 background running problem under PGI + MPI
Replies: 4
Views: 5365

Try putting the model run with a shell script, executing the shell script, and then logging out.

Our you can use the "screen" utility, which lets you disconnect from an active shell session, and reconnect to it later.

Cheers,
Jamie
by jpringle
Fri Sep 07, 2007 1:48 am
Forum: ROMS Problems
Topic: How to increase ROMS parallel efficiency
Replies: 5
Views: 5878

The efficiency of ROMS in parallel mode is a strong function of 1) the grid size 2) the interconnect speed 3) the configuration of MPI, and 4) the configuration of the model (parameters such as NDTFAST can make a big difference). It is hard to answer a question such as yours with no knowledge of the ...
by jpringle
Wed May 09, 2007 2:31 pm
Forum: ROMS Tools and Techniques
Topic: Grid generation
Replies: 61
Views: 210248

I have a matlab code that calls gridgen and allows you to specify multiple points along each edge, and to specify (roughly) a spatially variable resolution within the model domain. It works well for me, and is reasonably well documented. One caution, however -- if you have a very high-resolution ...
by jpringle
Tue Feb 27, 2007 2:52 pm
Forum: ROMS Discussion
Topic: Inflow-Outflow options for reservoirs
Replies: 12
Views: 11891

Dmitry-

I would be happy to, but I cannot find your email address. You can find my contact info at my web page, http://oxbow.sr.unh.edu .

Cheers,
Jamie
by jpringle
Sun Feb 25, 2007 3:07 pm
Forum: ROMS Discussion
Topic: Inflow-Outflow options for reservoirs
Replies: 12
Views: 11891

ROMS can handle multiple inflow points. This smells like a subtle coding error. My biggest errors in trying to specify point sources in ROMS have occurred because I calculated the positions in matlab, and did not properly take into account the fact that for array start in fortran at index 0 are ...
by jpringle
Sat Feb 24, 2007 5:06 pm
Forum: ROMS Discussion
Topic: Inflow-Outflow options for reservoirs
Replies: 12
Views: 11891

Have you turned on the UV_PSOURCE, TS_PSOURCE and ANA_PSOURCE in the cppdefs.h file for your application?

Jamie
by jpringle
Thu Feb 08, 2007 8:27 pm
Forum: ROMS Installation
Topic: ROMS on an Intel Mac with a Core Duo and free compilers
Replies: 1
Views: 4278

ROMS on an Intel Mac with a Core Duo and free compilers

Hi all- I have experimented with running ROMS 2.2 on a iMac 2Ghz Intel Core Duo processor, with 667 MHz DDR2 RAM. I used free compilers installed with fink. As you will see, I suspect that I would do better to use the Intel compiler, but I have not done so since for large runs I use my linux ...
by jpringle
Sat Dec 23, 2006 1:56 pm
Forum: ROMS Installation
Topic: ifort and Fedora, again
Replies: 1
Views: 3924

Craig (& others) I have gone through this on numerous systems. The main trick is figuring out how MPI-whatever runs its various sub-jobs on each system. Does it spawn a csh? a sh? Can you specify? Once you do, you can place the ulimit, etc commands in the appropriate .bashrc .cshrc or whatever. To ...
by jpringle
Fri Dec 01, 2006 1:44 pm
Forum: ROMS Tools and Techniques
Topic: Grid generation
Replies: 61
Views: 210248

Dear all- If you are using Pavel of CSIRO's gridgen, I have a little chunk of matlab code which allows you to specify the resolution of the grid along the edges of the model domain. If you are interested, mail me. I am not yet distributing it publicly because the rest of my evolving non-interactive ...
by jpringle
Wed Oct 11, 2006 5:04 pm
Forum: ROMS Usage
Topic: geostrophically adjusted initial condition.
Replies: 6
Views: 8832

Thanks a bunch yall. John Wilkin's code worked fine for my use, with minor modifications. It results in a cleaner initialization, and thus is easier for me to explain the initial evolution in my class.

Cheers,
Jamie
by jpringle
Mon Oct 09, 2006 6:31 pm
Forum: ROMS Usage
Topic: geostrophically adjusted initial condition.
Replies: 6
Views: 8832

I am starting a short, semi-idealized run with density features that are only 2-3 times larger than the radius of deformation. As part of the adjustment processes, gravity waves rattle around the domain, and bother the boundary condition. I could get away with a simpler model setup if I did not have ...
by jpringle
Mon Oct 09, 2006 4:46 pm
Forum: ROMS Usage
Topic: geostrophically adjusted initial condition.
Replies: 6
Views: 8832

geostrophically adjusted initial condition.

Dear hive mind- Does anyone have a code that takes an initial density field on a ROMS grid (or any other grid) and a level of no motion and uses it to return a geostrophically balanced flow field? It would be nice if it delt with curved grids... It would not be hard to do, but I am hoping to avoid ...
by jpringle
Tue Sep 26, 2006 1:41 pm
Forum: ROMS Discussion
Topic: Understanding T_cline
Replies: 1
Views: 3088

You can find scripts that allow you to experiment with the vertical grid parameters in part of the discussion thread "Stretching parameters - vertical coordinates" on this bulletin board. There is also some discussion of these parameters in the old manual for SCRUM, which can be found at http ...
by jpringle
Thu Sep 21, 2006 5:54 pm
Forum: ROMS Tools and Techniques
Topic: Grid generation
Replies: 61
Views: 210248

Hi all- I think Kate's problems are, indeed, due to her mexcdf install. Several other folks have used my code without an issue. I am not sure, however, what the issue is, which is less helpful than I would like to be. Also, Kate mentioned This works for ETOPO5, ETOPO2, and now GEBCO, .... If you've ...
by jpringle
Mon Sep 18, 2006 12:33 pm
Forum: ROMS Tools and Techniques
Topic: Grid generation
Replies: 61
Views: 210248

Hi everyone- A few weeks ago I had this problem, and I created a matlab code that uses Sakov's grid gen to make a matlab domain, useing Rich P's m_map routines and some code I got from Rich Signell to create ROMS grid files from matlab data. It is relatively primitive right now, but it is easy to ...
by jpringle
Mon Apr 24, 2006 6:12 pm
Forum: ROMS Problems
Topic: How to run ROMS with MPI + OpenPBS with every processors?
Replies: 2
Views: 6282

Hi- I do not recognize your particular error message. However, it smells like a problem in the configuration of your mpi system. You should certainly be able to run one part of your job on each processor. However, if you are having troubles running the code, you should try getting MPI to work on ...
by jpringle
Mon Nov 07, 2005 3:05 pm
Forum: ROMS Problems
Topic: Blow-up
Replies: 3
Views: 4727

John- You would be able to handle a freshwater wedge with a small vertical resolution in a layer model, like those from FSU, in which the depth of the layers can vary in time. However, in a model like ROMS where the level positions are fixed in time, you must have about six or more levels IN your ...
by jpringle
Fri Oct 07, 2005 8:23 pm
Forum: ROMS Benchmarks
Topic: Benchmarks
Replies: 31
Views: 81553

infiniband vrs. gigabit

Hi all-- I wondered if the varying results people get when they try to understand the relative importance of the interconnects could have to do with the relative importance of the 2D timestepping (lots of little messages) vrs. the 3D timesteps (fewer, larger messages). So I ran the benchmark1 test ...
by jpringle
Thu Jun 30, 2005 7:57 pm
Forum: ROMS Benchmarks
Topic: Benchmarks
Replies: 31
Views: 81553

Dual core/dual processor Opteron 275 benchmarks

Hi all-- Courtesy of Western Scientific Inc, I benchmarked ROMS 2.2 on a dual processor, dual core Opteron 275 rig. These chips have two processing cores per chip, so a workstation with two processor chips will appear to Linux to have four processors. The rig I tested had the following specs: 1) 8GB ...
by jpringle
Sun May 15, 2005 8:06 pm
Forum: ROMS Benchmarks
Topic: Benchmarks
Replies: 31
Views: 81553

dual core benchmarks

Has anyone used the standard benchmark cases (for examples, those in Rich Signell's benchmarks above) to exam the speed of the dual core operteron processors that AMD has recently released? For the problem I have in mind, I plan to be running on a 5 or 6 node linux cluster, with each node haveing ...
by jpringle
Sun May 15, 2005 8:04 pm
Forum: ROMS Discussion
Topic: Interconnect speed for MPI: Gigabit, Myrinet and Infiniband
Replies: 5
Views: 9867

Dual core benchmarks

Has anyone used the standard benchmark cases (for examples, those in Rich Signell's benchmarks above) to exam the speed of the dual core operteron processors that AMD has recently released? For the problem I have in mind, I plan to be running on a 5 or 6 node linux cluster, with each node haveing ...