Search found 1139 matches
- Tue Mar 02, 2021 2:50 am
- Forum: ROMS Tools and Techniques
- Topic: Error to plot by the ROMS plotting package
- Replies: 17
- Views: 974
Re: Error to plot by the ROMS plotting package
Wow, are you really asking how Linux works? No wonder why you are having so many troubles at every step in your application. We have lots of posts in this forum from you. I am expecting users that select ocean modeling as a carrier or hobby to be more curious and familiar with computers. By the way...
- Tue Feb 02, 2021 4:09 pm
- Forum: ROMS Installation
- Topic: Compiling Error for COAWST
- Replies: 6
- Views: 681
Re: Compiling Error for COAWST
Hmm, the call to mp_scatter_state is correct. It is only used for the adjoint-based propagator analysis, which I doubt Users using those algorithms that compact the full ROMS state ( zeta , u , v , T , S ) into a 1D array for matrix operations with the ARPACK library. I think that the issue here is ...
- Thu Jan 28, 2021 3:00 am
- Forum: ROMS Installation
- Topic: Compilation error for build_roms.sh
- Replies: 17
- Views: 840
Re: Compilation error for build_roms.sh
Obviously, you need to work on the libraries available on your computer. There is plenty of information about Compiling and Linking in WikiROMS. If you are using the latest version of ROMS, you may benefit from using CMake . Check cbuild_roms.csh or cbuild_roms.sh . CMake is a different philosophy f...
- Wed Dec 16, 2020 4:25 pm
- Forum: ROMS Discussion
- Topic: Bottom flux for Fennel model
- Replies: 3
- Views: 419
Re: Bottom flux for Fennel model
Yes, good catch. Thank you. The stflx and btflx state variables are used in ROMS governing equations. Contrarily, the stflux and btflux are used exclusively to process data from NetCDF files or analytical functions. The problem here is that the assignment is missing in set_vbc.F : # if defined BIOLO...
- Fri Dec 11, 2020 1:40 am
- Forum: ROMS Discussion
- Topic: ROMS 3.7 nesting problems
- Replies: 6
- Views: 372
Re: ROMS 3.7 nesting problems
I don’t know what you are doing wrong. First, I highly recommend that you use the latest version of ROMS. Right now, we are in version 3.9. We have improved ROMS nesting. It is more efficient. The latest version of ROMS is consistent with all the files in the current test repository. That’s the reas...
- Tue Oct 20, 2020 3:26 pm
- Forum: ROMS Discussion
- Topic: [Edited] Vertical profiles of EKE budget terms
- Replies: 3
- Views: 567
Re: [Edited] Vertical profiles of EKE budget terms
Thank you, Lars. Those are great references. I added them to my reading list. A couple of decades ago, I tried to derive an energy conservation equation for ROMS, so the only way to add energy to the system is from external forcing. I was unable to eliminate the pressure gradient term because of the...
- Tue Oct 13, 2020 5:39 am
- Forum: ROMS Adjoint
- Topic: station file
- Replies: 3
- Views: 383
Re: station file
Nope, that bash script supports both SLURM (BATCH=1) and non SLURM (BATCH=0). You just need to change that flag. See the customizable section. I can run that script in my 14 CPUs Mac or in our SLURM supercomputer. The script is very generic, compact, and complex. It is well documented. You need to s...
- Mon Oct 12, 2020 5:22 pm
- Forum: ROMS Adjoint
- Topic: station file
- Replies: 3
- Views: 383
Re: station file
What is RPG? Did you mean RPCG option for RBL4D-Var minimization algorithm? I have never turned on the STATIONS option with the 4D-Var algorithms. I suppose that we can use the split 4D-Var algorithms that I released yesterday to only run the STATIONS in the background and analysis 4D-Varphases. And...
- Fri Sep 18, 2020 7:50 pm
- Forum: ROMS Discussion
- Topic: Evaporation + precipitation mass conservation
- Replies: 5
- Views: 1011
Re: Evaporation + precipitation mass conservation
I think that there is a confusion here with the forcing terminology. The surface and bottom heat ( Q ) and freshwater fluxes ( E-P ) are the boundary conditions to the vertical tracer (temperature and salinity, respectively) diffusion terms in the governing equations. We need to provide boundary con...
- Thu Aug 27, 2020 2:52 am
- Forum: ROMS Sediment
- Topic: Roms blowup at 1cell bays
- Replies: 2
- Views: 377
Re: Roms blowup at 1cell bays in Lake_jersey test case
The lake Jersey test case works for me. What you did is not LAKE_JERSEY :!: You changed the grid, land/sea mask, and added big and small islands. Stop calling it Lake Jersey because it is a different application and will confuse other users with misinformation! If doing nesting, it will require diff...
- Tue Jul 14, 2020 3:24 pm
- Forum: ROMS Bugs
- Topic: Bug in set_grid.F
- Replies: 1
- Views: 398
Re: Bug in set_grid.F
Yes, thank you. It just need:
Please, update.
Code: Select all
#ifdef NESTING
USE mod_nesting
#endif
- Thu Jun 25, 2020 2:12 am
- Forum: ROMS Discussion
- Topic: MPI - tiling on POWER8
- Replies: 2
- Views: 473
Re: MPI - tiling on POWER8
What do you mean by the same results? profiling time or state solution?
Usually, it is much efficient to do less partition in the I-direction (NtileI) and more in the J-direction (NtileJ) because of vectorization. Also, it depends on the computer architecture.
Usually, it is much efficient to do less partition in the I-direction (NtileI) and more in the J-direction (NtileJ) because of vectorization. Also, it depends on the computer architecture.
- Wed Jun 24, 2020 7:14 pm
- Forum: ROMS Bugs
- Topic: Small fortran bug in nl_ocean.h
- Replies: 1
- Views: 397
Re: Small fortran bug in nl_ocean.h
Yes, thank you!
- Wed Jun 24, 2020 6:13 pm
- Forum: ROMS Problems
- Topic: mp_collect bug when compiling with define PARALLEL_IO
- Replies: 18
- Views: 5816
Re: mp_collect bug when compiling with define PARALLEL_IO
Everybody has their own procedure for updating code from repositories. I have several research repositories, and I usually do the following steps: Copy the version that I want to update to a temporary one. So, I have a version of the code that works in case something goes wrong with the updating. I ...
- Tue Jun 23, 2020 4:21 pm
- Forum: ROMS Problems
- Topic: mp_collect bug when compiling with define PARALLEL_IO
- Replies: 18
- Views: 5816
Re: mp_collect bug when compiling with define PARALLEL_IO
Well, you cannot pick and choose what routines you need to update. It doesn't work that way. Every change done to ROMS is well documented in the trac system, and the changes sometimes depend on previous changes and may be modified in future changes. The last major update to the Parallel I/O was done...
- Tue Jun 23, 2020 3:07 am
- Forum: ROMS Problems
- Topic: mp_collect bug when compiling with define PARALLEL_IO
- Replies: 18
- Views: 5816
Re: mp_collect bug when compiling with define PARALLEL_IO
I can believe that this is done :shock: We are not responsible for corrections to the code by a third party. Your solution is straightforward. Update your code from the official ROMS repositories. The agreement to receive a freely distributed version of the ROMS framework from us is that Users need ...
- Sun Jun 21, 2020 5:17 pm
- Forum: ROMS Discussion
- Topic: inconsistent explicit time-step splitting in two-way nesting
- Replies: 3
- Views: 689
Re: inconsistent explicit time-step splitting in two-way nesting
Well, the error is telling your problem, if you think carefully about it. I assume that you previously defined SOLVE3D . Now, it is undefined , and you are solving the shallow-water equations of ROMS where you must have NDTFAST =1 for all nested grids, and you have NDTFAST = 30 :!: If you don't und...
- Tue Jun 16, 2020 9:27 pm
- Forum: ROMS Discussion
- Topic: units of omega
- Replies: 6
- Views: 714
Re: units of omega
The varinfo.dat is an user-changeable file. In the beginning, we were thinking of writing omega in the restart NetCDF file as m 3 /s, and true vertical velocity in history files. But that's no longer the case. Nowadays, omega is recomputed internally. I rarely activate the switch for writing omega i...
- Tue Jun 16, 2020 8:13 pm
- Forum: ROMS Discussion
- Topic: units of omega
- Replies: 6
- Views: 714
Re: units of omega
Well, the omega vertical velocity in ROMS is special and its orthogonality changes at each grid point because it depends on the bathymetry. It does not have a physical interpretation since it is only related to ROMS descrete, terrain-following coordinates, 3D continuity equation. Internally in ROMS,...
- Tue Jun 09, 2020 8:50 pm
- Forum: ROMS Problems
- Topic: forcing interpolation not working
- Replies: 7
- Views: 1313
Re: forcing interpolation not working
I also have been using this capability for years and frequently and I never have any problems. Our group at Rutgers uses the same strategy on a daily basis in our operational system. Something is wrong in your files. When I do nesting, I only specify the forcing files for the first grid, and ROMS wi...
- Tue Jun 09, 2020 3:01 pm
- Forum: ROMS Discussion
- Topic: two-way nested simulation blowing up with MPI
- Replies: 3
- Views: 565
Re: two-way nested simulation blowing up with MPI
The reason why nobody replayed to your post is that you didn't provide any concrete information about the ROMS or Fortran error. We are not magicians or have any potions to see the error written to your standard output. I have mentioned several times over the years in this forum to never never never...
- Mon Jun 08, 2020 3:48 pm
- Forum: ROMS Discussion
- Topic: the defination of density rho
- Replies: 16
- Views: 2220
Re: the defination of density rho
Are you trying to change ROMS equation of state? The definition of in situ density anomaly or potential density anomaly is that a 1000 kg m-3 value is subtracted to the density. The in situ density anomaly, as a field, is only used in ROMS in the pressure gradient algorithm where its value is not e...
- Mon Jun 01, 2020 12:45 am
- Forum: ROMS Bugs
- Topic: compiling 4DVar with WET_DRY + wclock_off problem
- Replies: 4
- Views: 805
Re: compiling 4DVar with WET_DRY + wclock_off problem
Yes, that's a small bug, thank you. I am running nowadays with the ROMS_STDOUT option that generates the output file log.roms instead of the standard Fortran output from unit 6. My bad, I forgot to run without the ROMS_STDOUT . Having this option activated is very convenient when running coupled sys...
- Thu May 28, 2020 4:01 pm
- Forum: ROMS Bugs
- Topic: compiling 4DVar with WET_DRY + wclock_off problem
- Replies: 4
- Views: 805
Re: compiling 4DVar with WET_DRY + wclock_off problem
Oh my goodness, no. We cannot use WET_DRY in any of the adjoint-based algorithms including any of the 4D-Var drivers. That's the reason that you will find WET_DRY_NOT_YET in the code. Every time that you see a multiplication of a state variable with a mask array, we will need another term, which is ...
- Wed May 27, 2020 5:16 pm
- Forum: ROMS Adjoint
- Topic: I4dvar compilation error about i4dvar.f90 and mod_fourdvar
- Replies: 3
- Views: 719
Re: I4dvar compilation error about i4dvar.f90 and mod_fourdvar
Yes, gfortran is a stricter compiler. Other compilers relax the Fortran standard rule that we cannot have a module, subroutine, or variable with the same name. I missed this one. Thank you for reporting this issue. Please, update your code.
- Mon May 11, 2020 10:22 pm
- Forum: ROMS Adjoint
- Topic: Adjoint model blow up
- Replies: 13
- Views: 2423
Re: Adjoint model blow up
In our experience, it is challenging for either the tangent linear ( TLM ) or adjoint models ( ADM ) to blow-up in the first iteration of the inner loops. Something must be wrong in your step-up. Both models are linearized from the nonlinear model ( NLM ) state trajectory (prior). The TLM and ADM st...
- Fri May 01, 2020 3:07 pm
- Forum: ROMS Bugs
- Topic: Possible bugs in def_state.F
- Replies: 1
- Views: 472
Re: Possible bugs in def_state.F
Yes, thank you. Please update your code.
- Tue Apr 28, 2020 3:07 pm
- Forum: ROMS Bugs
- Topic: Ltracing compiling error
- Replies: 2
- Views: 533
Re: Ltracing compiling error
Update your code again!
- Mon Apr 27, 2020 9:11 pm
- Forum: ROMS Bugs
- Topic: FORTRAN error - Probably related to January changes
- Replies: 4
- Views: 613
Re: FORTRAN error - Probably related to January changes
I actually removing that parameter in an update that I am releasing today. Remove the couple of ocurrences of Ltracing in read_phypar.F. I have several reseach repositories and sometimes unwanted features leaks into the trunk. Sorry about that.
- Wed Apr 22, 2020 4:42 pm
- Forum: ROMS Bugs
- Topic: Possible bug in wetdry.F for nested apps
- Replies: 14
- Views: 2015
Re: Possible bug in wetdry.F for nested apps
Great work, Stefan. I concur with john about your telescoping approach into meter scales. It will take me some time to digest the changes that you make to the code. If you provide us the modifications that you have to do for the shallow-water nested applications to work, we can incorporate such chan...
- Mon Apr 06, 2020 8:31 pm
- Forum: ROMS Discussion
- Topic: error in inp_par.F
- Replies: 4
- Views: 783
Re: error in inp_par.F
Okay, update your repository. It is now fixed. Here is the
trac ticket

- Thu Apr 02, 2020 6:59 pm
- Forum: ROMS Adjoint
- Topic: Where to introduce obs error
- Replies: 2
- Views: 550
Re: Where to introduce obs error
Of course, you need the observation error covariance (diagonal matrix) in the ROMS 4D-Var algorithm! The NetCDF variable is obs_error , which is read into ROMS variable ObsErr . You are searching for the wrong variable. Check ROMS/Modules/mod_fourdvar.F for the description of all variables.
- Wed Apr 01, 2020 9:35 pm
- Forum: ROMS Discussion
- Topic: error in inp_par.F
- Replies: 4
- Views: 783
Re: error in inp_par.F
Yes, it will be easy to fix. I will update the repository tonight. Thank you for reporting this issue.
- Sun Mar 22, 2020 3:56 pm
- Forum: ROMS Discussion
- Topic: Imbalance in Parallel writing in ROMS
- Replies: 3
- Views: 746
Re: Imbalance in Parallel writing in ROMS
Hi Koushik, great study, and thank you for the research that you continue doing about ROMS parallel I/O :!: Let me digest the information and come up with a strategy for more efficient parallel I/O. Are you aware that ROMS has the CPP option NO_WRITE_GRID , which suppresses the definition and writin...
- Tue Mar 10, 2020 3:34 pm
- Forum: ROMS Source
- Topic: Question about HSIMT in step3d_t.F (trunk version)
- Replies: 1
- Views: 778
Re: Question about HSIMT in step3d_t.F (trunk version)
Yes, good catch :!: Thank you. I have been looking for a bug in the HSIMT scheme on and off. We fixed the issue that it needs three ghost points a few days ago, but forgot to make changes to the range of oHz . We need have an additional switch, Lhsimt : logical :: LapplySrc, Lhsimt, Lmpdata ... Lhsi...
- Sun Mar 01, 2020 7:47 pm
- Forum: ROMS Source
- Topic: Typo in roms.in following overhaul of advection schemes
- Replies: 4
- Views: 1040
Re: Typo in roms.in following overhaul of advection schemes
Yes, thank you for reporting this typo in the documentation. However, why do you want to set a different tracer advection scheme for temperature and salinity across nested grids? I will never do such a thing in any of my nesting applications. It doesn't make sense to me. You need to understand how t...
- Tue Feb 25, 2020 7:07 pm
- Forum: ROMS Problems
- Topic: Error when compiling 2D test cases
- Replies: 2
- Views: 561
Re: Error when compiling 2D test cases
Yes, I already knew about this one and included in a repository update that will be released today.
- Tue Feb 18, 2020 4:06 pm
- Forum: ROMS Discussion
- Topic: Problem with surface heat flux
- Replies: 13
- Views: 1919
Re: Problem with surface heat flux
Hmmm, I noticed that you activated both BULK_FLUXES and CCSM_FLUXES . I am not familiar with that version of ROMS. However, my understanding is that you should use a single formulation. You need to choose one and not both. If the sea ice model is activated, we need BOTH freshwater (E-P) and saltwate...
- Tue Feb 18, 2020 1:12 am
- Forum: ROMS Bugs
- Topic: average file count logic
- Replies: 2
- Views: 516
Re: average file count logic
I revised this logic recently. My assumption is that ndefAVG is an exact multiple (integer factor) of ntimes . If it is not the case, the file creation definition is changed. What values are you using that gives you OutFiles=0 ? The time averaged files are special because of its delayed creation. Th...
- Mon Feb 17, 2020 4:36 pm
- Forum: ROMS Discussion
- Topic: Problem with surface heat flux
- Replies: 13
- Views: 1919
Re: Problem with surface heat flux
I don't know if this behaviour will be related, but I would turn off SPLINES_VDIFF and SPLINES_VVISC and see what happens. I also noticed that you are using a different version of ROMS that has sea ice modeling. This will complicate your analysis.
- Fri Feb 07, 2020 8:45 pm
- Forum: ROMS Bugs
- Topic: Possible bug in wetdry.F for nested apps
- Replies: 14
- Views: 2015
Re: Possible bug in wetdry.F for nested apps
Well, u2dbc_im.F and v2dbc_im.F need to be modified to include code relevant to 2D shallow-water nesting applications. I guess that we never accounted for a User needing nesting in a 2D application. It will need an equation to impose the transport from the donor coarser grid. We will have to figure ...
- Fri Feb 07, 2020 6:35 pm
- Forum: ROMS Bugs
- Topic: Possible bug in wetdry.F for nested apps
- Replies: 14
- Views: 2015
Re: Possible bug in wetdry.F for nested apps
Notice that the indices told and tnew are local (dummy) time indices: ! ! Set time snapshot indices for the donor grid data. ! told=3-RollingIndex(cr) tnew=RollingIndex(cr) However, if you are running a 2D shallow-water application ( SOLVE3D is off), why are you dealing with state variables ( DU_avg...
- Fri Feb 07, 2020 4:20 pm
- Forum: ROMS Bugs
- Topic: Possible bug in wetdry.F for nested apps
- Replies: 14
- Views: 2015
Re: Possible bug in wetdry.F for nested apps
Hmm. I will have to think about it but maybe it will be needed. It will need to be tested. There are a set of those indices for RHO-, PSI-, U-, and V-points. It will be only relevant it coastlines are present in the contact points.
- Thu Feb 06, 2020 11:44 pm
- Forum: ROMS Bugs
- Topic: ROMS 3.9 MPDATA bugs
- Replies: 2
- Views: 788
- Thu Feb 06, 2020 5:20 pm
- Forum: ROMS Adjoint
- Topic: 4DVAR assimilation of surface velocities destroy sst
- Replies: 3
- Views: 971
Re: 4DVAR assimilation of surface velocities destroy sst
It is all a matter of the assumed background error covariance hypothesis and its parameters. What do you trust more the model or the observations? HF radar velocity observations are very tricky since it is a very noisy data, and depends on the processing: averaging, detiding, super-observations, err...
- Sun Feb 02, 2020 3:36 pm
- Forum: ROMS Problems
- Topic: A ROMS nesting question
- Replies: 10
- Views: 1910
Re: A ROMS nesting question
You need to make sure that the variables in your forcing NetCDF files have the _FillValue = 1.e+35 attribute or whatever other value you used. Otherwise, ROMS doesn't know how to process such data without the _FillValue that have land/sea masks. However, the fact the you have Min = 1.0E+35 and Max =...
- Fri Jan 31, 2020 8:57 pm
- Forum: ROMS Problems
- Topic: A ROMS nesting question
- Replies: 10
- Views: 1910
Re: A ROMS nesting question
Check your standard input file, roms.in or your standard output log file, which contains lots of information. Notice that in roms.in we have: ! Logical switches (TRUE/FALSE) to nudge the desired climatology field(s). ! If not analytical climatology fields, users need to turn ON the logical ! switche...
- Sun Jan 26, 2020 8:33 pm
- Forum: ROMS Problems
- Topic: Three nested refined grids - Ngrids vs. Ncontact
- Replies: 2
- Views: 744
Re: Three nested refined grids - Ngrids vs. Ncontact
The nesting configuration is tricky since we support various types of grid connectivity. We only provide a Matlab interface to generate the nesting connectivity NetCDF file needed by ROMS. Unfortunately, we don't officially have any other software to create such a connectivity NetCDF file nor have t...
- Thu Jan 16, 2020 10:35 pm
- Forum: ROMS Discussion
- Topic: ROMS 3.7
- Replies: 1
- Views: 514
Re: ROMS 3.7
The frozen versions of ROMS are available only on SVN. Usually, one checkout the tags SVN repository: svn checkout https://www.myroms.org/svn/src/tags tags svn checkout https://www.myroms.org/svn/src/tags/roms-3.7 roms-3.7 The GIT repository works in a different way. If you need the version that was...
- Thu Jan 16, 2020 4:04 am
- Forum: ROMS Bugs
- Topic: Matlab nesting pre-processing program error.
- Replies: 2
- Views: 790
Re: Matlab nesting pre-processing program error.
Yes, your first correction is fine. Thank you :!: I updated the repository. We didn't catch that one before because all our grids always have the ( XI , ETA ) coordinates even if the applications grid is spherical. It turns out that in curvilinear grids is better to do the interpolation in ( XI , ET...
- Sat Jan 04, 2020 7:35 pm
- Forum: ROMS Bugs
- Topic: Segmentation fault
- Replies: 9
- Views: 1226
Re: Segmentation fault
This is not a bug. It always has to be equal or greater than one by design. If no input forcing files are needed, ROMS will ignore those values on standard input. I don't understand what it is the Fuzz about it. There are a lot of parameters in ROMS standard input that will be only relevant for a sp...
- Mon Dec 23, 2019 2:42 pm
- Forum: ROMS Discussion
- Topic: calculating ubar and vbar
- Replies: 10
- Views: 2691
Re: calculating ubar and vbar
Look at Nonlinear/ini_fields.F instead. It is very clear how it is computed there. It is how the inital ubar anf vbar is computed in ROMS.
- Tue Dec 17, 2019 6:16 am
- Forum: ROMS Discussion
- Topic: Request: Numerical modeling course syllabi
- Replies: 11
- Views: 2459
Re: Request: Numerical modeling course syllabi
I learned how to program and plot with the NCARG graphical library. I spent so much time reading the manuals and learning how to use that library. It was very satisfying plotting the results with computer-generated contours. I became a pro using that wonderful library. I still used it today. I wrote...
- Tue Dec 17, 2019 2:07 am
- Forum: ROMS Discussion
- Topic: Request: Numerical modeling course syllabi
- Replies: 11
- Views: 2459
Re: Request: Numerical modeling course syllabi
Indeed, I took that class from John Klinck when I was a graduate student at Texas A&M. I learned a lot about coding and ocean numerical modeling in that first class. It was a great introduction to a world that it is now my professional life. I recall that we, the students, and John, enjoyed that cla...
- Mon Dec 09, 2019 3:27 am
- Forum: ROMS Bugs
- Topic: typo in distribute.F
- Replies: 1
- Views: 549
Re: typo in distribute.F
Hi Frank, great catch. Thank you! I updated the repository. See the following
trac ticket for more information.

- Fri Dec 06, 2019 11:37 pm
- Forum: ROMS Discussion
- Topic: Parallel IO not scaling over Serial IO for ROMS
- Replies: 14
- Views: 1787
Re: Parallel IO not scaling over Serial IO for ROMS
An application of size 899x629x40 is TOO small to be run in 1350 PEs :!: No wonder why is not scaling after 240 PEs. You are killing scalability with excessive MPI collective communications and exchanges. Your 20x20 tile size is too small, and it is killing vectorization. The good news is that you f...
- Thu Dec 05, 2019 4:54 pm
- Forum: ROMS Discussion
- Topic: Parallel IO not scaling over Serial IO for ROMS
- Replies: 14
- Views: 1787
Re: Parallel IO not scaling over Serial IO for ROMS
Great study. We are missing important information about your simulations. We need the grid size and the parallel domain decomposition. I have mentioned frequently in this forum that ROMS scalability is a function of the application and depends on the domain size and parallel partition. There is alwa...
- Fri Nov 15, 2019 9:55 pm
- Forum: ROMS Installation
- Topic: Compiling error - undefined reference to ice routines
- Replies: 12
- Views: 2121
Re: Compiling error - undefined reference to ice routines
The version of the code that we distribute doesn't have a native sea ice model yet. I am sure that it will be included in the future. It has to be located in a subdirectory inside of the Nonlinear directory, like the biology and sediment models. The reason for it is quite logical and straightforward...
- Sun Aug 04, 2019 4:21 pm
- Forum: ROMS Installation
- Topic: upwelling: - Configuration error ..... exit_flag: 5
- Replies: 3
- Views: 913
Re: upwelling: - Configuration error ..... exit_flag: 5
The default upwelling example in ROMS test repository has UV_QDRAG activated in the build script. If you notice, there are various options in build_roms.sh : # can be used to write time-averaged fields. Notice that you can have as # many definitions as you want by appending values. setenv MY_CPP_FLA...
- Thu Jul 25, 2019 4:00 pm
- Forum: ROMS Releases
- Topic: ROMS 3.7 Released
- Replies: 0
- Views: 2313
ROMS 3.7 Released
ROMS/TOMS 3.7 Released The latest version of the :arrow: ROMS/TOMS svn repository , revision 974 , was frozen and tagged as ROMS 3.7 . We haven't tagged a frozen version of ROMS in several years. This version is quite stable, and it is the one to be used at the upcoming 2019 ROMS 4D-Var Workshop to...
- Mon Jul 22, 2019 3:03 pm
- Forum: ROMS Bugs
- Topic: Typo in ad_set_vbc.F line 1165 rev 973
- Replies: 1
- Views: 521
Re: Typo in ad_set_vbc.F line 1165 rev 973
Indeed, thank you! I didn't test with LIMIT_STFLX_COOLING activated.
- Tue Jul 09, 2019 8:33 pm
- Forum: ROMS Bugs
- Topic: Error regarding variable dP in prsgrd44.h or prsgrd42.h
- Replies: 1
- Views: 568
Re: Error regarding variable dP in prsgrd44.h or prsgrd42.h
Indeed. Thank you for reporting this issue. Please update.
- Thu Jun 27, 2019 7:48 pm
- Forum: ROMS Problems
- Topic: ROMS different src version leading to unalike results
- Replies: 7
- Views: 1657
Re: ROMS different src version leading to unalike results
The ROMS basic numerical kernel has not changed much in a few years. However, configuration, CPP options, and parameters in the standard input have been modified and enhanced. You need to be sure that your header file ( application.h ) and standard input script ( roms.in ) is equivalent between the ...
- Wed Jun 12, 2019 1:06 am
- Forum: ROMS Discussion
- Topic: Time type in nc file
- Replies: 5
- Views: 1040
Re: Time type in nc file
ROMS has nothing to do with the time convention and elapsed time. The user selects how the time is managed in an application. There is nothing that we can do about it. It is up to the person(s) that set-up that application to use a consistent time zone in all the input NetCDF files (initial conditio...
- Wed May 15, 2019 6:51 pm
- Forum: ROMS Installation
- Topic: compile ROMS error
- Replies: 6
- Views: 10330
Re: compile ROMS error
Your problem is that you haven't configured the build_roms.sh and build_roms.bash correctly. For example, notice that in the build_roms.sh script we have: #-------------------------------------------------------------------------- # If applicable, use my specified library paths. #-------------------...
- Sun May 12, 2019 1:44 pm
- Forum: ROMS Problems
- Topic: Segmentation fault with floats
- Replies: 6
- Views: 1369
Re: Segmentation fault with floats
Segmentation violations (SIGSEGV) are usually a sign of memory problems on the computer and available stack for the run. Perhaps, your computer doesn't have too much memory for your application, and you are deploying too many floats for it. I will try to run without floats to see if it works. Then a...
- Wed May 08, 2019 10:35 pm
- Forum: ROMS Installation
- Topic: Compile ROMS Error: Symbol 'eps' at (1) has no IMPLICIT type
- Replies: 4
- Views: 1003
Re: Compile ROMS Error: Symbol 'eps' at (1) has no IMPLICIT
You didn't understand at all what I was saying to you, or what WET_DRY does. The wetting and drying option ( WET_DRY ) is modeled in ROMS by changing the associated land/sea masking arrays at every barotropic time-step in step2d , which affect the ocean state variables. The wet/dry mask is a new arr...
- Wed May 08, 2019 6:31 pm
- Forum: ROMS Installation
- Topic: Compile ROMS Error: Symbol 'eps' at (1) has no IMPLICIT type
- Replies: 4
- Views: 1003
Re: Compile ROMS Error: Symbol 'eps' at (1) has no IMPLICIT
The WET_DRY option is not adjointed because it not fully differentiable, and because of every time-step of the forward trajectory need to be saved to have a consistent linearization of the tangent linear (TLM) and adjoint (ADM) models. You cannot use such time-depending land/sea masking in any ROMS ...
- Wed May 01, 2019 4:22 pm
- Forum: ROMS Discussion
- Topic: cppdefs for BULK_FLUX
- Replies: 11
- Views: 2435
Re: cppdefs for BULK_FLUX
My current strategy is to activate LONGWAVE_OUT and provide downward longwave radiation so bulk_flux will subtract the outgoing IR by using ROMS surface temperature ( K ), infrared emissivity, and Stefan-Boltzmann constant. I also use downward shortwave radiation and I activate COOL_SKIN to correct ...
- Tue Apr 16, 2019 9:50 pm
- Forum: ROMS Discussion
- Topic: CF compliant units attributes in Varinfo.dat
- Replies: 7
- Views: 1603
Re: CF compliant units attributes in Varinfo.dat
ROMS has nowadays a very robust calendar. See the routine ROMS/Utility/dateclock.F for more details. It is consistently internally and can handle different reference date ( epoch ) for each forcing field from NetCDF files regardless of what we choose in the ROMS input parameter TIME_REF . I have to ...
- Tue Apr 16, 2019 2:23 pm
- Forum: ROMS Discussion
- Topic: CF compliant units attributes in Varinfo.dat
- Replies: 7
- Views: 1603
Re: CF compliant units attributes in Varinfo.dat
I agree with John. This part is not desirable. We and others have been using this standard for years. There are lots of datasets out there that follow this nomenclature. I am sympathetic with the issue of the links but some of the changes suggested are radical.
- Thu Apr 11, 2019 9:33 pm
- Forum: ROMS Information
- Topic: units of specific humidity
- Replies: 10
- Views: 6475
Re: units of specific humidity
Indeed, I believe that unsavory piece of code was in the original COARE code. It is quite clever if you ask me. I recall putting enough comments in the code when I did some research about this threshold value many years ago. There is not a standardized unit for specific humidity in the atmospheric d...
- Thu Mar 21, 2019 9:55 pm
- Forum: ROMS Discussion
- Topic: CF compliant units attributes in Varinfo.dat
- Replies: 7
- Views: 1603
Re: CF compliant units attributes in Varinfo.dat
Thank you. The ROMS metadata file varinfo.dat is a user file and can be edited in any way that the user feels is necessary. You can even change the language. Some users in France wanted some of the attribute values to be in French. I haven't read the documentation about the UDUNITS-2 package. I look...
- Wed Mar 13, 2019 2:28 am
- Forum: ROMS Bugs
- Topic: A small bug in ana_srflux.h
- Replies: 12
- Views: 5143
Re: A small bug in ana_srflux.h
Finally, I have time to browse the literature, Googling, and write a Matlab script to check the issue mentioned above. Indeed, the 15.0 factor used to compute the zenith is incorrect :!: We need to have instead: zenith=cff1+cff2*COS(Hangle-lonr(i,j)*deg2rad) In my Matlab script, I compared against...
- Sun Mar 10, 2019 10:38 pm
- Forum: ROMS Bugs
- Topic: Pair and MPI version has problem with regridding
- Replies: 10
- Views: 2058
Re: Pair and MPI version has problem with regridding
I updated the code to correct the bug in regrid.F. Check the following trac ticket
src
808 for more details. The parallel bug was corrected. Good luck.

- Sat Mar 09, 2019 5:53 pm
- Forum: ROMS Bugs
- Topic: Pair and MPI version has problem with regridding
- Replies: 10
- Views: 2058
Re: Pair and MPI version has problem with regridding
It doesn't make sense to me. MyXout is a state, tiled variable and it is allocated as the others, and the pointer is passed correctly. It is the only way that this can be done. I bet that the problem is not in regrid . It seems like a memory leakage somewhere else. Yes, you can put the application i...
- Fri Mar 08, 2019 9:04 pm
- Forum: ROMS Bugs
- Topic: Pair and MPI version has problem with regridding
- Replies: 10
- Views: 2058
Re: Pair and MPI version has problem with regridding
I took a look in the debugger on our US East Coast application and I cannot find anything wrong. I activated both ATM_PRESS and PRESS_COMPENSATE . I am also using BULK_FLUXES , which also need Pair . I don't see a parallel bug. Pair.png shflux.png The Pair is kind of jagged but it is because of the ...
- Fri Mar 08, 2019 5:16 am
- Forum: ROMS Bugs
- Topic: Pair and MPI version has problem with regridding
- Replies: 10
- Views: 2058
Re: Pair and MPI version has problem with regridding
Yes, the Dmem is a diagnostic quantity to estimate the memory requirement for an application. It has nothing to do with the numerical kernel. It is in my TODO list to look at your problem in the debugger. The issues that you are talking about sound like a parallel bug. However, the regrid subroutine...
- Thu Mar 07, 2019 5:27 am
- Forum: ROMS Bugs
- Topic: A small bug in ana_srflux.h
- Replies: 12
- Views: 5143
Re: A small bug in ana_srflux.h
Both ana_srflux.h and ana_specir.h were coded longtime ago when we were testing the EcoSim bio-optical model in the late 90's. It was like twenty years ago. So I don't remember the details. At that time, both routines were intended for idealized toy problems for debugging purposes. The EcoSim model ...
- Thu Feb 28, 2019 3:36 am
- Forum: ROMS Source
- Topic: My ROMS
- Replies: 9
- Views: 4568
Re: My ROMS
Thank you, Mark. Your logic is correct. We don't need to add processing logic since it needs ATM_PRESS anyway. I loaded your request ticket to the repository.
- Mon Feb 11, 2019 10:34 pm
- Forum: ROMS Tools and Techniques
- Topic: Soundspeed for ROMS output
- Replies: 4
- Views: 1333
Re: Soundspeed for ROMS output
Just do it in Matlab.
- Mon Feb 11, 2019 5:12 am
- Forum: ROMS Tools and Techniques
- Topic: Soundspeed for ROMS output
- Replies: 4
- Views: 1333
Re: Soundspeed for ROMS output
It has been in ROMS plotting package that uses the NCAR's GKS library since the beginning. Check the file src/Utility/svel.F in the plot repository. It uses the Che and Millero (1976) formula. Chen, C.T. and F.J. Millero, 1976: re-evaluation of Willson sound-speed measurements in pure water. J. Acou...
- Sun Feb 10, 2019 1:14 am
- Forum: Ocean News
- Topic: Walter Munk (1917-2019)
- Replies: 0
- Views: 1934
Walter Munk (1917-2019)
The oceanography community and the world lost one of its most renowned scientists. Walter Munk was 101. I met Walter in the mid-'80s when I was a graduate student at Texas A&M University. He visited and gave us an excellent Seminar/Lecture. My thesis advisor Robert O. Reid introduced me to him. It w...
- Fri Feb 01, 2019 6:51 pm
- Forum: ROMS Results
- Topic: ROMS Granular model result in all variables
- Replies: 8
- Views: 2545
Re: ROMS Granular model result in all variables
:shock: Well, that's your problem. You are using biharmonic diffusivity (m4/s), which scales as 1/12 * ABS(U) * (dx ** 3) . A value for TS_DIF4 of 2 m4/s is too small :!: You need values of the order of 2.0d+8 m4/s . That is, you are using a value that is eight orders of magnitude too small :!: No w...
- Thu Jan 31, 2019 9:09 pm
- Forum: ROMS Results
- Topic: ROMS Granular model result in all variables
- Replies: 8
- Views: 2545
Re: ROMS Granular model result in all variables
I guess that the more fundamental question is what type of horizontal advection are you using? Are you specifying explicit diffusion and viscosity? What do you specify for vertical diffusion/viscosity parameterization closure? You should compute your Reynolds Number (U dx/nu), where U is the maximum...
- Tue Jan 22, 2019 10:35 pm
- Forum: ROMS Source
- Topic: My ROMS
- Replies: 9
- Views: 4568
Re: My ROMS
Hi Mark, Are you ready to submit the changes due to the PRESS_COMPENSATE?
Thank you.
Thank you.
- Tue Jan 22, 2019 9:22 pm
- Forum: ROMS Bugs
- Topic: Issue restarting solutions in fraction of a day
- Replies: 1
- Views: 924
Re: Issue restarting solutions in fraction of a day
Indeed, good catch
Thank you for reporting this bug. I will update the code soon.

- Tue Jan 22, 2019 7:24 pm
- Forum: ROMS Bugs
- Topic: tiny bugs? in metrics.F and nesting.F
- Replies: 1
- Views: 959
Re: tiny bugs? in metrics.F and nesting.F
Thank you for reporting this problem. I forgot to take care of the roundoff issues with the interpolation weights when I implemented the single precision ( SINGLE_PRECISION option) simulations. Your suggested solutions are incorrect :!: We cannot take the ABS value when computing the interpolation w...
- Sat Dec 15, 2018 8:53 pm
- Forum: ROMS Bugs
- Topic: Equation error in GLS model coded into gls_corstep.F
- Replies: 5
- Views: 1722
Re: Equation error in GLS model coded into gls_corstep.F
I haven't updated the code yet. I think that Ramsey Harcourt is correct. The shape functions used in mostly all vertical mixing parameterization should be smooth and continuous. We should stay away from sharp discontinuities that lead to eye bullet values in vertical diffusion and vertical viscosity...
- Fri Dec 14, 2018 8:32 pm
- Forum: ROMS Bugs
- Topic: Equation error in GLS model coded into gls_corstep.F
- Replies: 5
- Views: 1722
Re: Equation error in GLS model coded into gls_corstep.F
Thank you for looking this carefully and for providing a Matlab script to show the difference between the formulation currently coded in ROMS and the correction to the limiter (green curve). Obviously, the coded formulation has a discontinuity spike just before it remains constant. It is clearly wro...
- Mon Dec 03, 2018 2:14 am
- Forum: ROMS Problems
- Topic: 'nfast' and 'ndtfast' on Barotropic time step
- Replies: 5
- Views: 1829
Re: 'nfast' and 'ndtfast' on Barotropic time step
Again, Nope :!: The governing equations are timestep from n (right-hand-side terms) to n+1 (left-hand-side term: time rate of change). The vertically-integrated equations are integrated with smaller timestep to resolve fast dynamics due to gravity wave phenomena. For 3D total momentum coupling, the ...
- Sun Dec 02, 2018 4:19 pm
- Forum: ROMS Problems
- Topic: 'nfast' and 'ndtfast' on Barotropic time step
- Replies: 5
- Views: 1829
Re: 'nfast' and 'ndtfast' on Barotropic time step
Nope , there are actually M* barotropic timesteps :!: The reason for it is that the cosine-square shape filter needs to provide time-averaged values of ubar , vbar , and zeta centered at the n+1 baroclinic timestep, which gives us second-order temporal accuracy for time-averaged barotropic motions....
- Wed Nov 14, 2018 5:00 pm
- Forum: ROMS Bugs
- Topic: tiny little bug in get_state.F for time evolving bathymetry
- Replies: 2
- Views: 1065
Re: tiny little bug in get_state.F for time evolving bathyme
Done. Thank you for reporting this typo.
- Fri Nov 09, 2018 2:43 pm
- Forum: ROMS Bugs
- Topic: Error compiling ROMS with new flag SINGLE_PRECISION (svn928)
- Replies: 7
- Views: 2409
Re: Error compiling ROMS with new flag SINGLE_PRECISION (svn
Nope. What I said above is correct :!: If the last argument to the mp_bcastf routine is single- precision it will use mp_bcast_1d , and if it is double-precision it will use mp_bcastf_1dp . Just look the subroutines in distribute.F :idea: If you don't understard, you need to check a Fortran 1990, 19...
- Wed Nov 07, 2018 9:16 pm
- Forum: ROMS Bugs
- Topic: Error compiling ROMS with new flag SINGLE_PRECISION (svn928)
- Replies: 7
- Views: 2409
Re: Error compiling ROMS with new flag SINGLE_PRECISION (svn
You are doing something wrong or something is weird with the compiler. You need to read carefully the following :arrow: trac ticket . Notice that even if SINGLE_PRECISION is activated, there are several variables that are still in double precision ( dp 64-bit kind parameter) and the r8 kind paramete...
- Mon Oct 01, 2018 4:10 pm
- Forum: ROMS Discussion
- Topic: missing temp_time
- Replies: 2
- Views: 952
Re: missing temp_time
You don't need to do that. One of the great things about NetCDF is that it is self-described. I have mentioned in the past that we need to follow CF convections to guarantee the functionality of metadata and avoid missinterpretation issues in ROMS. I highly recommend users to add the units , time , ...
- Fri Sep 28, 2018 7:50 pm
- Forum: ROMS Problems
- Topic: severe (174): SIGSEGV, segmentation fault occurred. libpthre
- Replies: 13
- Views: 3391
Re: severe (174): SIGSEGV, segmentation fault occurred. libp
I am going to try again for the last time. Read carefully :arrow: trac ticket 747 . In the older version of the code, we choose either lower- or higher-level MPI function for exchanges. We no longer do that in the newer versions, you need to experiment and select which options are more efficient in ...
- Thu Sep 27, 2018 1:53 pm
- Forum: ROMS Problems
- Topic: severe (174): SIGSEGV, segmentation fault occurred. libpthre
- Replies: 13
- Views: 3391
Re: severe (174): SIGSEGV, segmentation fault occurred. libp
I think that you need to read the following :arrow: trac ticket and choose the MPI communication options that are more efficient in the computer environment that you are running. You should check the profiling information that ROMS reports to standard output to see in what region of the code are slo...
- Wed Sep 26, 2018 10:08 pm
- Forum: ROMS Discussion
- Topic: ROMS history file has suddenly included restart record
- Replies: 7
- Views: 1806
Re: ROMS history file has suddenly included restart record
Hmmm, I haven't had the time to examine the issue that you are having. I haven't modified this part of the code in years. I use the split of history files frequently and I have never observed the behavior that you are reporting. I wonder if it is a matter of configuration. Are you changing DSTART ev...
- Wed Sep 26, 2018 10:00 pm
- Forum: ROMS Problems
- Topic: severe (174): SIGSEGV, segmentation fault occurred. libpthre
- Replies: 13
- Views: 3391
Re: severe (174): SIGSEGV, segmentation fault occurred. libp
Yes, your problem is the stack size per CPU and it seems to be associated with the automatic arrays used in distributed-memory for I/O operations. This is not a direct ROMS problem, but a computer problem because of not enough memory to handle automatic arrays that are either allocated on stack or h...
- Tue Sep 25, 2018 1:48 pm
- Forum: ROMS Problems
- Topic: severe (174): SIGSEGV, segmentation fault occurred. libpthre
- Replies: 13
- Views: 3391
Re: severe (174): SIGSEGV, segmentation fault occurred. libp
It is very simple as I have mentioned several times before. You just need to edit your login script and add one of the lines below:
I wrote lots of information in previous
trac ticket.
Code: Select all
.cshrc, .tcshrc, etc.
limit stacksize unlimited
or .bashrc
ulimit -s unlimited

- Sat Sep 22, 2018 6:18 am
- Forum: ROMS Problems
- Topic: severe (174): SIGSEGV, segmentation fault occurred. libpthre
- Replies: 13
- Views: 3391
Re: severe (174): SIGSEGV, segmentation fault occurred. libp
Nowadays, severe segmentation errors are associated with stack size, which is used for allocating automatic arrays. They are allocated on stack or heap according to you choice of compiler options. I mentioned this in the last trac ticket.
- Thu Sep 13, 2018 9:12 pm
- Forum: ROMS Bugs
- Topic: Small bugs in CPP BIO_FENNEL with DIAG and DIAGNOSTICS_BIO
- Replies: 2
- Views: 1206
Re: Small bugs in CPP BIO_FENNEL with DIAG and DIAGNOSTICS_
Nope. That's not what I have in the repository for def_diags.F : # ifdef DIAGNOSTICS_BIO # if defined BIO_FENNEL || defined HYPOXIA_SRM DO ivar=1,NDbio2d ifield=iDbio2(ivar) IF (TRIM(var_name(i)).eq.TRIM(Vname(1,ifield))) THEN got_var(ifield)=.TRUE. DIA(ng)%Vid(ifield)=var_id(i) END IF END DO # endi...
- Tue Sep 11, 2018 9:23 pm
- Forum: ROMS Bugs
- Topic: bug in mod_tides.F
- Replies: 1
- Views: 996
Re: bug in mod_tides.F
Indeed, thank you. I have more stuff to update in the computation of ROMS memory requirements. I am missing arrays to account for. I am trying to come up with the best estimate.
- Mon Sep 10, 2018 7:25 pm
- Forum: ROMS Discussion
- Topic: ROMS prestep3d , Crack-Nicolson method, lambda
- Replies: 1
- Views: 792
Re: ROMS prestep3d , Crack-Nicolson method, lambda
Yes, we need to update the documentation.
- Fri Aug 31, 2018 9:00 pm
- Forum: ROMS Discussion
- Topic: Cannot find "coordinates" attribute for variable: Uwind
- Replies: 15
- Views: 3161
Re: Cannot find "coordinates" attribute for variable: Uwind
The ECMWF data has a minimal interval of 3 hours; you can easily create annual files instead monthly. It will be more efficient in ROMS to avoid inquiring and open so many NetCDF files. The time units for forcing file is usually days since the reference date ( epoch ). However, you can have any time...
- Sun Aug 26, 2018 3:44 pm
- Forum: ROMS Problems
- Topic: forrtl: severe (151): allocatable array is already allocated
- Replies: 17
- Views: 5074
Re: forrtl: severe (151): allocatable array is already alloc
Okay, that will do it. I put safeguards to the code when processing standard input file ocean.in . I made a couple of updates recently for it. It looks that your grid is not for a regional application, but it is a global grid. Is that correct? If that's the case, you have a problem with the design o...
- Wed Aug 22, 2018 2:07 pm
- Forum: ROMS Problems
- Topic: forrtl: severe (151): allocatable array is already allocated
- Replies: 17
- Views: 5074
Re: forrtl: severe (151): allocatable array is already alloc
It is a weird error if you didn't repeat the NtileJ parameter in ocean.in , which triggers the allocation of several modules. In yesterday update, I put safeguards for this to never happen. I don't know what to make of your case. I don't think that a corrupted ocean.in should trigger to process the ...
- Mon Aug 20, 2018 4:58 pm
- Forum: ROMS Discussion
- Topic: How to setup for rotating frame
- Replies: 10
- Views: 2418
Re: How to setup for rotating frame
It is weird that you are not activating periodic boundary conditions in the south and northern boundaries :!: It is the perfect and appropriate case for it because your grid is an annulus :roll: You can set the desired flow in the initial condition instead of using Flather boundary conditions. If yo...
- Mon Aug 13, 2018 6:08 pm
- Forum: ROMS Discussion
- Topic: blowing up on nested grid
- Replies: 5
- Views: 1702
Re: blowing up on nested grid
That can be also a problem. However, the main issue here is that there are one or more grid cells with a value of 1E+40 1/km for both inverse metric factors pm and pn . We cannot multiply those variables with the mask because of division by zero in ROMS. The grid spacing ( DX=1/pm and DY=1/pn ) is a...
- Sun Aug 12, 2018 3:36 pm
- Forum: ROMS Discussion
- Topic: blowing up on nested grid
- Replies: 5
- Views: 1702
Re: blowing up on nested grid
The information is in the front of you. If you don't understand that information, you have homework to do before running ROMS. Look the large and a small numbers. Check your grid files. ( DXmin , DYmin ) cannot be zero!!! It is a trillion times smaller than the scales of electrons... I am sure that ...
- Sun Jul 29, 2018 4:38 pm
- Forum: ROMS Problems
- Topic: huge memory per node when up scaling application
- Replies: 14
- Views: 3387
Re: huge memory per node when up scaling application
I usually put lots of information in trac when I update the code. We have noticed that this information is generally ignored and our emails are treated as spam. There is nothing that we can do about it, and we cannot force users to update and or read the detailed information provided. We have thousa...
- Fri Jul 27, 2018 8:00 pm
- Forum: ROMS Problems
- Topic: huge memory per node when up scaling application
- Replies: 14
- Views: 3387
Re: huge memory per node when up scaling application
PARALLEL_IO works for me. You need to use a recent version of the code. Also, make sure that all your input NetCDF files are NetCDF-4 format. You can check which format you have using ncdump : % ncdump -k filename.nc netCDF-4 To convert from NetCDF-3 to Netcdf-4 try: % nccopy -k netCDF-4 nc3name.nc...
- Wed Jul 25, 2018 4:30 pm
- Forum: ROMS Problems
- Topic: huge memory per node when up scaling application
- Replies: 14
- Views: 3387
Re: huge memory per node when up scaling application
First, I don't think that ROMS will run with the PNETCDF library. It will need special function calls in ROMS. I never tried that library, it is a third party library. If you use PARALLEL_IO, you need to compile with the NetCDF-4 library with parallel MPI I/O and HDF5 support and your computer need ...
- Wed Jul 18, 2018 9:43 pm
- Forum: ROMS Problems
- Topic: Problem generating contact file with three refined grids
- Replies: 5
- Views: 1861
Re: Problem generating contact file with three refined grids
Hi Marc and John,
I removed the restriction of having the noncoincident grid points in the indirect connectivity between grids for telescoping applications. See the following
trac ticket for more information.
Please update. Your particular set-up should work now.
I removed the restriction of having the noncoincident grid points in the indirect connectivity between grids for telescoping applications. See the following

Please update. Your particular set-up should work now.
- Tue May 22, 2018 7:05 am
- Forum: ROMS Bugs
- Topic: ifort compile (make) error
- Replies: 4
- Views: 1840
Re: ifort compile (make) error
Yes, it is Darwin configuration that needs it. As Kate said, just comment the line in the Linux-ifort configuration file. Linux will take the unlimited value from the environmental variable. It is the loader LD that it is having the problem. It depends on what you have in your computer.
- Thu Apr 26, 2018 7:04 pm
- Forum: ROMS Problems
- Topic: unable to open existing NetCDF file GET GRID
- Replies: 10
- Views: 2029
Re: unable to open existing NetCDF file GET GRID
By the way, you don't need to modify that routine. ROMS prints out the NetCDF error code at the end of the standard output file when it closes all the files. Notice that there are NetCDF-3 and NetCDF-4 input files for this application. You may be using the wrong one for the compiled code.
- Fri Mar 30, 2018 5:54 pm
- Forum: ROMS Problems
- Topic: Could not read in multiple boundary files correctly
- Replies: 2
- Views: 1162
Re: Could not read in multiple boundary files correctly
You need to provide more information. The new inquiry.F was modified to provide the correct transition between annual files.
- Wed Feb 28, 2018 5:30 pm
- Forum: ROMS Problems
- Topic: Contact point problem with salinity in Two Way Nesting run
- Replies: 6
- Views: 2455
Re: Contact point problem with salinity in Two Way Nesting r
I haven't noticed this in my applications that have run stable for years. Are you using those three refinement grids? If so, one of the refined grids is touching the boundary of the other. That will affect the order of the contact points between them. It is the chicken and egg question, what came fi...
- Tue Feb 20, 2018 3:54 pm
- Forum: ROMS Bugs
- Topic: List of analytical header files
- Replies: 3
- Views: 1413
Re: List of analytical header files
It cannot be done during initialization because the CPP include syntax <*****.h> allows to use the header files from user's directories other than ROMS/Functionals . Notice that in analytical.F , we have # ifdef SOLVE3D # if defined ANA_BIOLOGY && defined BIOLOGY # include <ana_biology.h> # endif # ...
- Tue Feb 20, 2018 5:30 am
- Forum: ROMS Bugs
- Topic: List of analytical header files
- Replies: 3
- Views: 1413
Re: List of analytical header files
There is a reason for it that I need to remember exactly. It is not a bug, but why does it matter? I think that it has to do with with the global attributes of output NetCDF files and the fact that the user can overwrite the location of the code header files instead of using the ones distributed by ...
- Tue Dec 12, 2017 6:35 pm
- Forum: ROMS Problems
- Topic: CPU time and MPI issues with nested grids
- Replies: 16
- Views: 6779
Re: CPU time and MPI issues with nested grids
Yes, you need to select a better option for MPI communications as I mentioned yesterday. Check the information on the svn trac ticket. If you are using too many processes for this application, the data exchanges due to nesting is a bottleneck. See the number for point data gathering . It explains th...
- Tue Dec 12, 2017 1:24 am
- Forum: ROMS Problems
- Topic: CPU time and MPI issues with nested grids
- Replies: 16
- Views: 6779
Re: CPU time and MPI issues with nested grids
Maybe your application is not big enough to justify 144 nodes, and your application is slowing down because of excessive MPI communications. It is a critical issue that some users ignore always. There is always an optimal parallel partition for each ROMS application. Now if you are using a recent ve...
- Sat Dec 09, 2017 5:27 pm
- Forum: ROMS Bugs
- Topic: A little bug in get_cycle.F of src 884 and the fixing
- Replies: 5
- Views: 2754
Re: A little bug in get_cycle.F of src 884 and the fixing
Do you have the cycle_length attribute in the time variable? For the forward time-stepping, we can also have the following logic in get_cycle.F: IF (job.gt.0) THEN ! forward: Tval(i) =< mday =< Tval(i+1) IF (Linside) THEN tstart=Tmin DO i=2,ntime IF ((tstart.le.mday).and.(mday.le.Tval(i))) THEN Tstr...
- Fri Dec 08, 2017 5:07 pm
- Forum: ROMS Bugs
- Topic: A little bug in get_cycle.F of src 884 and the fixing
- Replies: 5
- Views: 2754
Re: A little bug in get_cycle.F of src 884 and the fixing
Something is weird here. It doesn't make sense to me. Notice that the logic in get_cycle.F that set the switches is: ! ! Is the model time inside the data time range? If not, check if the ! data just has the LOWER- or the UPPER-snapshot interpolant. ! IF ((Tmin.le.mday).and.(mday.le.Tmax)) THEN Lins...
- Tue Nov 21, 2017 6:32 pm
- Forum: ROMS Bugs
- Topic: Bug in output time step in DEF_DIAGS after restart
- Replies: 5
- Views: 3309
Re: Bug in output time step in DEF_DIAGS after restart
Smaller values of NDIA are possible but its value needs to be an exact factor of NRST, say NDIA=NRST/2 should give you an integer. The time window (in seconds) should be exactly divided by DT.
- Fri Nov 17, 2017 6:26 pm
- Forum: ROMS Bugs
- Topic: Bug in output time step in DEF_DIAGS after restart
- Replies: 5
- Views: 3309
Re: Bug in output time step in DEF_DIAGS after restart
Well, this logic has worked for us for years. Your problem is that NAVG and NDIA are inconsistent with NRST . If you have AVERAGES and DIAGNOSTICS_TS and/or DIAGNOSTICS_UV activated, you need to set NAVG=NDIA . It doesn't make any sense to have different values because if you are going to compute ba...
- Tue Sep 26, 2017 1:24 am
- Forum: ROMS Problems
- Topic: mp_collect bug when compiling with define PARALLEL_IO
- Replies: 18
- Views: 5816
Re: mp_collect bug when compiling with define PARALLEL_IO
The PARALLEL_IO is slower than serial I/O, and it won't help you to improve I/O performance. The last time that I checked that option, it was not working well, and I haven't checked if the NetCDF4 library has made progress on it. I can fix the issue with nf_read2d.F but it will not help you much. Th...
- Tue Aug 29, 2017 4:27 pm
- Forum: ROMS Problems
- Topic: Error in ROMS execution
- Replies: 7
- Views: 2059
Re: Error in ROMS execution
You cannot change those values capriciously. They have to be exactly the same values used by the person that generated the initial conditions, boundary conditions, atmospheric forcing, and observations for data assimilation. The time clock has to be the same in all those files for consistency. The t...
- Mon Aug 28, 2017 7:35 pm
- Forum: ROMS Problems
- Topic: Error in ROMS execution
- Replies: 7
- Views: 2059
Re: Error in ROMS execution
There is a weird set-up of the ROMS time clock here. You need to ask the person that set-up this application. First, the input parameters: 13295.000 dstart Time-stamp assigned to model initialization (days). 19000101.00 time_ref Reference time for units attribute (yyyymmdd.dd) It implies that the RO...
- Mon Aug 28, 2017 4:36 pm
- Forum: ROMS Problems
- Topic: Error in ROMS execution
- Replies: 7
- Views: 2059
Re: Error in ROMS execution
Provide more information from the log file, not just the end. If you are using the latest version of ROMS, there is detailed information in the log file with source code and line number where the error occurred.
- Tue Jul 11, 2017 10:10 pm
- Forum: ROMS Benchmarks
- Topic: Understanding ROMS output
- Replies: 5
- Views: 2917
Re: Understanding ROMS output
Yes, the primal formulation I4DVAR is very expensive because the minimization is done in the physical space (i.e. full application grid for the state vector, 10^5 to 10^7 points) so we need to do a lot of I/O manipulations in cgradient to compute and orthogonalize the Lanczos vectors. We need to rea...
- Mon Jul 10, 2017 2:27 pm
- Forum: ROMS Benchmarks
- Topic: Understanding ROMS output
- Replies: 5
- Views: 2917
Re: Understanding ROMS output
What do you mean by the Lanczos algorithm? Are you talking about the minimization routine, rpcg_lanzos (RPCG)? The Lanczos vectors are stored in a NetCDF and use the adjoint model iterations over all outer and inner loops. I recently profiled WC13 with W4DPSAS : Nouter=2 and Ninner=50 on 4 CPUs on m...
- Wed Jun 07, 2017 6:14 pm
- Forum: Job Opportunities
- Topic: Postdoctoral Research Position at UCLA
- Replies: 0
- Views: 1966
Postdoctoral Research Position at UCLA
Postdoctoral Researcher Position in Continental Shelf Currents and Their Surface Wave Interactions Departments of Atmospheric and Oceanic Sciences University of California, Los Angeles (UCLA) We are seeking a researcher to work on problems of oceanic circulation dynamics using the Regional Oceanic ...
- Fri May 26, 2017 8:17 pm
- Forum: ROMS Bugs
- Topic: A new Y2K bug!
- Replies: 2
- Views: 1927
Re: A new Y2K bug!
Fixed the bug. It was not a Y2K bug. I coded datevec incorrectly in module dateclock.F . I was able to reproduce your problem with a simple test. Now, I get: 1 2000-02-28 00:02:30.00 9.193018E-09 1.954312E+04 1.954312E+04 1.659839E+17 (001,32,30) 1.799562E-05 8.159871E-10 0.000000E+00 4.667957E-03 2...
- Thu May 25, 2017 11:04 pm
- Forum: ROMS Bugs
- Topic: A new Y2K bug!
- Replies: 2
- Views: 1927
Re: A new Y2K bug!
Weird. What calendar are you using? The behavior of your time reporting indicates to me that you are using the Proleptic Gregorian Calendar ( time_ref=0 or time_ref>0 ). The datenum and datevec routines in ROMS starts in in March 1, 0000 to facilitate the management of leap years: MyMonth=MOD(month+...
- Sat Apr 29, 2017 2:00 pm
- Forum: ROMS Installation
- Topic: ran_state compilation error in ROMS
- Replies: 5
- Views: 2442
Re: ran_state compilation error in ROMS
Weird, I compiled with ifort , gfortran , and pgi and didn't get those errors. It turns out that I missed to transfer the updated files gasdev.F , nrutil.F , and ran1.F when I committed the ticket in another computer. We need these files to have i8b instead i4b . This is the reason why you are havin...
- Wed Apr 26, 2017 8:22 pm
- Forum: ROMS Ecosystem
- Topic: Compiling error: dateclock?
- Replies: 5
- Views: 2610
Re: Compiling error: dateclock?
Yes, thank you for reporting this one. The weird thing is that I ran and debugged ( TotalView ) this version several times with the ifort compiler and didn't report any errors or warnings. It seems to be happy with the dot syntax. Of course, the correct syntax in Fortran is to use Rclock % . This ha...
- Fri Apr 14, 2017 8:29 pm
- Forum: ROMS Adjoint
- Topic: nonlinear cost function not calculated in PSAS
- Replies: 1
- Views: 1591
Re: nonlinear cost function not calculated in PSAS
What are you reporting does not make sense to me. The NL_iDataPenalty and NL_fDataPenalty is computed and any of the 4D-Var drivers; w4dpsas_ocean.h in your case. In none of those drivers the VERIFICATION option is used. You just need to use the grep command to confirn this. The definition of these ...
- Fri Apr 14, 2017 6:46 pm
- Forum: ROMS Discussion
- Topic: ROMS not bit-reproducible: northern boundary?
- Replies: 32
- Views: 9556
Re: ROMS not bit-reproducible: northern boundary?
Well, it seems that there is a problem with the latest version of OpenMPI . I wonder what did they change. This kind of problem is difficult to track and frustrating. Compilers and libraries are known to have bugs. One thing to check is if the OpenMPI library was compiled correctly with the same com...
- Wed Jan 25, 2017 11:17 pm
- Forum: ROMS Bugs
- Topic: History file output fails using WRITE_WATER and river inflow
- Replies: 2
- Views: 2645
Re: History file output fails using WRITE_WATER and river in
Good catch, thank you
I updated the repository. Check the following
trac ticket for more information.


- Tue Jan 24, 2017 3:18 am
- Forum: ROMS Bugs
- Topic: Error in Flux Integrals in Correct_Tracer - two-way nesting
- Replies: 1
- Views: 5771
Re: Error in Flux Integrals in Correct_Tracer - two-way nest
Yes, thank you very much for reporting this problem with great detail and providing a solution :!: I know that this kind of problems are kind of frustrating and takes a lot of patience to find them. It is refreshing when users get deep into the code and try to understand and find solutions to the pr...
- Tue Jan 24, 2017 2:51 am
- Forum: ROMS Bugs
- Topic: Error with Multifiles
- Replies: 9
- Views: 7647
Re: Error with Multifiles
Yes, thank you for reporting this. It is difficult to guess the length of that string. I set it to 4096. See the following
trac ticket for more information.

- Tue Jan 24, 2017 2:44 am
- Forum: ROMS Bugs
- Topic: Quicksave bug
- Replies: 1
- Views: 1649
Re: Quicksave bug
Yes, thank you. It is always a good strategy to initialize all the switches and variables. I updated the repository. For more information check the following
trac ticket.

- Tue Jan 24, 2017 1:50 am
- Forum: ROMS Bugs
- Topic: A possible bug in get_data.F
- Replies: 1
- Views: 7261
Re: A possible bug in get_data.F
Yes, I can see why some compilers may have problems with this type of syntax. Thank you for reporting this issue. This is now fixed. See following
trac ticket.

- Mon Oct 10, 2016 7:48 pm
- Forum: Meetings/Workshops
- Topic: 2016 ROMS Asia-Pacific Workshop, 4D-Var Tutorial Slides
- Replies: 0
- Views: 7384
2016 ROMS Asia-Pacific Workshop, 4D-Var Tutorial Slides
2016 ROMS/TOMS Asia-Pacific Workshop Institute for Marine and Antarctic Studies Hobart Waterfront, Tasmania, Australia October 17 - 20, 2016 Andrew M. Moore will be giving a half day tutorial on ROMS 4D-Var. We are posting the slides for this presentation :arrow: here ahead of time so attendees can...
- Tue Oct 04, 2016 5:03 pm
- Forum: Meetings/Workshops
- Topic: 2016 ROMS Asia-Pacific Workshop, Tasmania - Final Notice
- Replies: 0
- Views: 5207
2016 ROMS Asia-Pacific Workshop, Tasmania - Final Notice
2016 ROMS/TOMS Asia-Pacific Workshop Institute for Marine and Antarctic Studies Hobart Waterfront, Tasmania, Australia October 17 - 20, 2016 We have had several cancellations, making it necessary to rework the agenda. If you are presenting , please check :arrow: the agenda again as your presentatio...
- Fri Sep 23, 2016 4:41 pm
- Forum: Meetings/Workshops
- Topic: 2016 ROMS/TOMS Workshop T-shirt
- Replies: 0
- Views: 5443
2016 ROMS/TOMS Workshop T-shirt
This is our design for the ROMS T-shirt for the upcoming Asia-Pacific Workshop. The front is the an artistic rendition of the Hobart skyline with Mt. Wellington as a backdrop.
- Fri Sep 23, 2016 3:00 pm
- Forum: Meetings/Workshops
- Topic: 2016 ROMS Asia-Pacific Workshop, Tasmania - 5th Notice
- Replies: 0
- Views: 4966
2016 ROMS Asia-Pacific Workshop, Tasmania - 5th Notice
2016 ROMS/TOMS Asia-Pacific Workshop Institute for Marine and Antarctic Studies Hobart Waterfront, Tasmania, Australia October 17 - 20, 2016 :!: NOTE: :!: Due to scheduling conflicts with other meetings the workshop has been compressed into 4 days. We apologize for the late notice on this. Please r...
- Wed Aug 17, 2016 2:11 am
- Forum: Meetings/Workshops
- Topic: 2016 ROMS Asia-Pacific Workshop, Tasmania - 4th Notice
- Replies: 0
- Views: 6138
2016 ROMS Asia-Pacific Workshop, Tasmania - 4th Notice
2016 ROMS/TOMS Asia-Pacific Workshop Institute for Marine and Antarctic Studies Hobart Waterfront, Tasmania, Australia October 17 - 21, 2016 This is a reminder that the abstract deadline has passed. Please check the :arrow: participants list to confirm that your abstract has been submitted. If your...
- Wed Aug 10, 2016 2:11 pm
- Forum: ROMS Bugs
- Topic: ROMS attempts to close un-named files
- Replies: 3
- Views: 4157
Re: ROMS attempts to close un-named files
Yes, please update. The error is due to the lack of the QCKNAME keyword in ocean.in . Notice that in :arrow: ticket 706 the field for NetCDF file ID, XXX%ncid=-1 , is initialized to close state in all the I/O structure arrays. This will allow to run newest ROMS version with old input ocean.in scripts.
- Wed Jul 27, 2016 3:56 pm
- Forum: Meetings/Workshops
- Topic: 2016 ROMS Asia-Pacific Workshop, Tasmania - 3rd Notice
- Replies: 0
- Views: 6817
2016 ROMS Asia-Pacific Workshop, Tasmania - 3rd Notice
2016 ROMS/TOMS Asia-Pacific Workshop Institute for Marine and Antarctic Studies Hobart Waterfront, Tasmania, Australia October 17 - 21, 2016 The extended registration deadline for the 2016 ROMS/TOMS Asia-Pacific Workshop is less than a week away. This is a reminder that registration deadline was ex...
- Fri Jul 15, 2016 2:12 pm
- Forum: ROMS Tools and Techniques
- Topic: standard_name for U,V current components
- Replies: 6
- Views: 4849
Re: standard_name for U,V current components
Nope. It is not a vertical sum divided by the number of levels. It is a matter of semantics here. Perhaps, depth-averaged will be more appropriate. The shallow water equations are obtained by vertically integrating momentum and continuity equations in flux form. In curvilinear coordinates (xi, eta) ...
- Thu Jul 14, 2016 5:06 pm
- Forum: ROMS Tools and Techniques
- Topic: standard_name for U,V current components
- Replies: 6
- Views: 4849
Re: standard_name for U,V current components
Nope, this cannot be done in varinfo.dat . I have thought about this in the past but we have so many variables. However, I can see how this is important for velocity components and third-party potting tools. It is very easy to add the standard_name to the variables in def_avg.F and def_his.F . We ju...
- Thu Jun 09, 2016 5:55 pm
- Forum: Meetings/Workshops
- Topic: 2016 ROMS Asia-Pacific Workshop, Tasmania - 2nd Notice
- Replies: 0
- Views: 7852
2016 ROMS Asia-Pacific Workshop, Tasmania - 2nd Notice
2016 ROMS/TOMS Asia-Pacific Workshop Institute for Marine and Antarctic Studies Hobart Waterfront, Tasmania, Australia October 17 - 21, 2016 The registration deadline for the 2016 ROMS/TOMS Asia-Pacific Workshop is approaching. This is a reminder that registration is open until June 30, 2016 . The ...
- Thu Apr 21, 2016 3:38 pm
- Forum: ROMS Bugs
- Topic: Bug with grid coordinates for station file
- Replies: 2
- Views: 2495
Re: Bug with grid coordinates for station file
Thank you. I have to think about this. I need to reproduce your problem so I can look in the debugger.
- Fri Feb 19, 2016 10:15 pm
- Forum: ROMS Problems
- Topic: Unrealisitic thickness of mixed layer
- Replies: 15
- Views: 12081
Re: Unrealisitic thickness of mixed layer
The MY2.5 scheme is coded correctly in ROMS and has several corrections that are not available in the original published formulation of the Mellor-Yamada 2.5. There are several papers in the literature about this. Also some people argue that the MY2.5 was not originally formulated correctly and the ...
- Fri Feb 19, 2016 3:38 am
- Forum: Meetings/Workshops
- Topic: 2016 ROMS Asia-Pacific Workshop, Tasmania - 1st Notice
- Replies: 0
- Views: 4725
2016 ROMS Asia-Pacific Workshop, Tasmania - 1st Notice
2016 ROMS/TOMS Asia-Pacific Workshop Institute for Marine and Antarctic Studies Hobart Waterfront, Tasmania, Australia October 17 - 21, 2016 We are inviting you to participate at the 2016 ROMS Asia-Pacific Workshop to be held at the Institute for Marine and Antarctic Studies building on the Hobart ...
- Thu Oct 29, 2015 2:54 pm
- Forum: ROMS Problems
- Topic: MP_EXCHANGE2D - communication buffer too small
- Replies: 6
- Views: 4150
Re: MP_EXCHANGE2D - communication buffer too small
This is not related to earlier posts in this forum thread. Anyway, this is weird. If you added biology, only mp_exchange4d is used in just two routines: pre_step3d.F: # ifdef DISTRIBUTE CALL mp_exchange4d (ng, tile, iNLM, 1, & & LBi, UBi, LBj, UBj, 1, N(ng), 1, NT(ng), & & NghostPoints, & & EWperiod...
- Fri Oct 02, 2015 4:19 pm
- Forum: ROMS Problems
- Topic: CPU time and MPI issues with nested grids
- Replies: 16
- Views: 6779
Re: CPU time and MPI issues with nested grids
That's good that you guys use the full stencil at the interface to compute the mass fluxes entering the refined grid. I haven't observed smoothness problems in my testing applications with one-way nesting and 1:5 ratios. Anyway, we should always use two-way nesting since the one-way interaction has ...
- Mon Sep 21, 2015 6:53 pm
- Forum: ROMS Problems
- Topic: CPU time and MPI issues with nested grids
- Replies: 16
- Views: 6779
Re: CPU time and MPI issues with nested grids
We need to take into account that ROMS_AGRIF has a reduced and simpler barotropic stepping engine ( step2d.F ), so it is more efficient. Therefore, we are not comparing apples with apples. Rutgers version of step2d.F includes additional terms so it is slower. In ROMS_AGRIF, all of those terms are in...
- Fri Sep 18, 2015 1:15 pm
- Forum: ROMS Problems
- Topic: Online Nesting slow
- Replies: 5
- Views: 3242
Re: Online Nesting slow
What is the size of your grids? Also, what are the partitions that you are using? I want to be sure that you are not using excessive parallel nodes. There is always an optimal partition in ROMS. This is computer engineering and we should abstain to put a lot of processes just because they are availa...
- Wed Sep 16, 2015 2:15 am
- Forum: ROMS Installation
- Topic: Can ROMS run just in Matlab?
- Replies: 2
- Views: 2654
Re: Can ROMS run just in Matlab?
Absolutely not. Why do you want to do such a thing? ROMS has over a million lines of code. There are ways to run simple and small Fortran programs in Matlab but the capability is very limited and you need to write an interface.
- Fri Sep 11, 2015 10:47 pm
- Forum: ROMS Problems
- Topic: two-way nesting only in the first refined grid
- Replies: 12
- Views: 5347
Re: two-way nesting only in the first refined grid
Please update. I fixed the problem that you were having. See the following
trac ticket for details. Thank you for reporting this problem 


- Thu Sep 10, 2015 7:57 pm
- Forum: ROMS Problems
- Topic: Online Nesting slow
- Replies: 5
- Views: 3242
Re: Online Nesting slow
Yes, we are aware of this and it is in our TO-DO list. The issue here is that in two-way nesting (default) we are currently making global gathering of data (full 2D and 3D arrays) between all the MPI nodes to compute easily the fine-to-coarse averaging. This causes a bottle neck because of the exces...
- Wed Sep 09, 2015 2:59 pm
- Forum: ROMS Problems
- Topic: Bug in large spherical nested run?
- Replies: 5
- Views: 2738
Re: Bug in large spherical nested run?
I have a new version of the Matlab scripts that I haven't released yet. I am still thinking about them. They are experimental and I will send them to your e-mail latter today. It provides alternatives to computing the metrics via interpolation instead of the great circle formula. Please sent me an e...
- Tue Sep 08, 2015 5:09 pm
- Forum: ROMS Problems
- Topic: Bug in large spherical nested run?
- Replies: 5
- Views: 2738
Re: Bug in large spherical nested run?
Thanks Jamie for the detailed information. I have to think about this one. Obviously, the spherical grid is doing something to it. It makes me suspicious of the area/volume conservation. I wonder if the grid spacing (area) between coarse and fine grids is conserved in the sherical case. I assume tha...
- Sat Sep 05, 2015 3:17 pm
- Forum: ROMS Problems
- Topic: two-way nesting only in the first refined grid
- Replies: 12
- Views: 5347
Re: two-way nesting only in the first refined grid
I will have to reproduce this problem with a test case to see what it is going on. Your fix to main3d.F is not correct.
- Fri Sep 04, 2015 10:10 pm
- Forum: ROMS Problems
- Topic: two-way nesting only in the first refined grid
- Replies: 12
- Views: 5347
Re: two-way nesting only in the first refined grid
Refined grid by a factor of 4 is illegal. I don't know why contact.m allowed you to work by such a factor. You should get an error warning and stop. You only can have refimnement ratios of 3 , 5 , 7 , 9 , etc in ROMS because it is discritized on a C-grid . As a matter of dynamics and because of nume...
- Fri Sep 04, 2015 5:47 pm
- Forum: ROMS Problems
- Topic: two-way nesting only in the first refined grid
- Replies: 12
- Views: 5347
Re: two-way nesting only in the first refined grid
Perhaps, you should show us a diagram (plot) of all grids and their refinement ratio with respect the donor grid to check if you are setting the grids correctly.
- Fri Sep 04, 2015 4:24 pm
- Forum: ROMS Problems
- Topic: two-way nesting only in the first refined grid
- Replies: 12
- Views: 5347
Re: two-way nesting only in the first refined grid
Nope, this is not a good idea! because of the internal manipulations in several routine inside module nesting . It seems to me that you didn't specified the parameters correctly. The order of the grids are crucial. You need to think about the logical order when generating the contact NetCDF file wi...
- Fri Sep 04, 2015 1:31 am
- Forum: ROMS Bugs
- Topic: bug in check_massflux in nesting.F
- Replies: 2
- Views: 2504
Re: bug in check_massflux in nesting.F
Yes, I did moved this code after I tested. Thank you for report this. I updated the repository. Nope, the values of BOUNDARY_CONTACT(:,:)%Mflux(:) are only used during debugging and testing to write several diagnostics to Fortran fort.100 , fort.200 , and fort.300 output files. Yes, they are loaded ...
- Wed Aug 19, 2015 5:09 pm
- Forum: ROMS Bugs
- Topic: Bug in nesting when MASKING not defined
- Replies: 1
- Views: 1786
Re: Bug in nesting when MASKING not defined
Yep, good catch. Thank yTo be precise we need instead:
I will update the code today.
Code: Select all
# if defined MASKING || defined WET_DRY
PRIVATE :: mask_hweights
# endif
- Tue Jul 21, 2015 8:41 pm
- Forum: ROMS Bugs
- Topic: Bug in fine2coarse.m for nesting
- Replies: 6
- Views: 3043
Re: Bug in fine2coarse.m for nesting
There are two ways to construct refined nested grids: (1) coarse to fine grid refinement or (2) fine to coarse grid extraction. I prefer to do option (2) since all the quantities needed in the contact area are already available in the donor grid. If you don't see this, think about refinement and com...
- Tue Jul 21, 2015 12:12 am
- Forum: ROMS Usage
- Topic: How to generate and use composite grids?
- Replies: 16
- Views: 8821
Re: How to generate and use composite grids?
Yes, it needs that information from the global attributes. Otherwise, contact.m will not work. These scripts are complex and there is a lot of internal details. They need to be used as a full package. I won't guarantee that a mixture of third party scripts will work together. There is a lot of subtl...
- Mon Jul 20, 2015 7:06 pm
- Forum: ROMS Usage
- Topic: How to generate and use composite grids?
- Replies: 16
- Views: 8821
Re: How to generate and use composite grids?
That's your problem! You need to use Lmask=false when generating the contact points with the Matlab script. Please read the instructions! Notice that the default value is false . :idea: I have said several times in this forum that if we don't have experience in nesting applications for a particular ...
- Thu Jul 16, 2015 7:02 pm
- Forum: ROMS Bugs
- Topic: Possible IO bug in close_io.F
- Replies: 1
- Views: 1894
Re: Possible IO bug in close_io.F
This is weird because close_inp is used in the adjoint-based algorithms to facilitate the iterations in the propagators and 4D-Var algorithms. This means that these files were already processed when reading input forcing fields. In the example that you provide above, which is exactly what I wrote in...
- Thu Jul 16, 2015 6:29 pm
- Forum: ROMS Bugs
- Topic: Small bug in ad_v2dbc_im.F
- Replies: 4
- Views: 2723
Re: Small bug in ad_v2dbc_im.F
I don't know what are you talking about. I checked the code and it already have that correction! I think that you need to update. As you see below, ad_v2dbc_im.F already have the change as shown on
trac for Revision 767:

- Thu Jul 16, 2015 5:15 pm
- Forum: ROMS Problems
- Topic: Create a sponge layer
- Replies: 5
- Views: 3465
Re: Create a sponge layer
I did a major overhaul of the sponge areas in ROMS sometime ago. We can have now visc_factor and diff_factor in the Grid NetCDF file. This is the way all application (simple or complex) that require areas of enhanced viscosity and diffusion. See the following :arrow: trac ticket for details. Notice ...
- Thu Jul 16, 2015 4:35 pm
- Forum: ROMS Wish List
- Topic: Wish list
- Replies: 2
- Views: 2521
Re: Wish list
The nesting in ROMS (one- or two-way) is designed for all tracers active and passive. The nesting interface is generic for all tracers. The only constraint is that the bio-geochemical or sediment model needs to be activated in all grids in order to compute the contact point data between all grids (d...
- Wed Apr 01, 2015 12:36 am
- Forum: ROMS Bugs
- Topic: Possible Issue with White_noise
- Replies: 4
- Views: 2508
Re: Possible Issue with White_noise
The routines in the white noise module white_noise.F are specific to the adjoint algorithms and we pass the Imin , Imax , Jmin , Jmax dummy variables as arguments. These are tiled parallel routines and need to be used by passing the appropriate values for these I- and J-range parameters. Using the a...
- Thu Mar 12, 2015 5:10 pm
- Forum: ROMS Bugs
- Topic: Small bug in ad_v2dbc_im.F
- Replies: 4
- Views: 2723
Re: Small bug in ad_v2dbc_im.F
Yes, thank you.
- Tue Mar 10, 2015 6:26 pm
- Forum: ROMS Tools and Techniques
- Topic: the depth of rho points for ROMS
- Replies: 2
- Views: 58285
Re: the depth of rho points for ROMS
The function zlevs is not part of our official distribution code and we are not responsible for third party functions coded by others. In the ROMS Matlab repository we offer: set_depth.m . This is the official Matlab version of such a function for computing depths at any C-grid location, if you use ...
- Thu Mar 05, 2015 1:50 am
- Forum: ROMS Discussion
- Topic: DISCUSSION: Parallel Global Reduction & Volume Conservation
- Replies: 0
- Views: 1303
DISCUSSION: Parallel Global Reduction & Volume Conservation
Users can use the following thread to discuss the
Parallel Global Reduction and Volume Conservation post from the ROMS Messages forum.

- Thu Mar 05, 2015 1:28 am
- Forum: ROMS Messages
- Topic: IMPORTANT: Parallel Global Reduction and Volume Conservation
- Replies: 0
- Views: 2647
IMPORTANT: Parallel Global Reduction and Volume Conservation
WARNING: You cannot replay or post comments in this forum thread. If you have any comments, use the following :arrow: post instead. In ROMS, we have the option to impose a volume conservation in applications with open boundaries if tidal forcing is not anabled: ! Set lateral open boundary edge volu...
- Wed Mar 04, 2015 1:13 am
- Forum: ROMS Usage
- Topic: How to generate and use composite grids?
- Replies: 16
- Views: 8821
Re: How to generate and use composite grids?
I already fixed this problem in contact.m . I pretty much rewrote this logic so it is more robust. However, I am not ready to release it yet. There are other changes that we need to think carefully about and discuss between us before we can release the updated scripts. I don't have the time to check...
- Fri Feb 27, 2015 6:26 pm
- Forum: ROMS Usage
- Topic: Is OUT_DOUBLE strictly necessary when using PERFECT_RESTART?
- Replies: 2
- Views: 1751
Re: Is OUT_DOUBLE strictly necessary when using PERFECT_REST
Yes, the ROMS restart file is always, by default, in double precision since all the computations in ROMS are in double precision. The option RST_SINGLE was added because some users were complaining that the restart file was getting too big in very large applications. The restart file has only two ro...
- Mon Feb 16, 2015 6:40 pm
- Forum: ROMS Problems
- Topic: Assimilating SSH – high res model blowup! Stability issue?
- Replies: 5
- Views: 3406
Re: Assimilating SSH – high res model blowup! Stability issu
The error covariance balanced operator is kind of tricky. It is based on Weaver et al. (2005) method. It imposes a multivariate constraint on the background error covariance such that the unobserved variable(s) information is extracted from the observed (assimilated) data using T-S empirical relatio...
- Mon Feb 16, 2015 5:29 pm
- Forum: ROMS Ecosystem
- Topic: light units for ec_Fobar
- Replies: 1
- Views: 2073
Re: light units for ec_Fobar
I don't know what are the units. This is used in the downwelling spectral irradiance routine ana_specir.h . This was adapted from a routine provided by Paul Bissett for the EcoSim model. That routine says that it is based on Gregg and Carder (1990) model and few units were provided. It is used to co...
- Wed Feb 11, 2015 6:15 pm
- Forum: ROMS Discussion
- Topic: How can i get the depth from sima coordinate to Z- ??
- Replies: 2
- Views: 1371
Re: How can i get the depth from sima coordinate to Z- ??
In the ROMS Matlab repository you can find the function depths.m or set_depth.m . More importantly, the above plot was produced with scoord.m , which outputs the quantity that you seek: [z,s,C]=scoord(h, x, y, Vtransform, Vstretching, theta_s, theta_b, ... hc, N, kgrid, column, index, plt, Zzoom) Yo...
- Wed Feb 11, 2015 6:02 pm
- Forum: ROMS Bugs
- Topic: Bug in wet_dry- the lake_jersey case has the same issue
- Replies: 8
- Views: 3243
Re: Bug in wet_dry- the lake_jersey case has the same issue
You need to start paying more attention to the updates to the code, specially if you are using nesting. You need to read carefully the following :arrow: post . It contains information about the quadratic interpolation. It is not ready yet :!: Also, the information about the wet/dry update was provid...
- Fri Feb 06, 2015 4:09 pm
- Forum: ROMS Problems
- Topic: Assimilating SSH – high res model blowup! Stability issue?
- Replies: 5
- Views: 3406
Re: Assimilating SSH – high res model blowup! Stability issu
I hope that you are plotting the cost function to see if 4D-Var is converging. It is not clear to me if you are using 5 or 30 inner loops. If you are using 5 inner loops, I doubt that the 4D-Var algorithm converged and this may explain your blow-up when running the nonlinear model. If your prior SSH...
- Thu Jan 29, 2015 5:49 pm
- Forum: ROMS Tools and Techniques
- Topic: Coastline file for the ROMS Plotting Package
- Replies: 4
- Views: 4291
Re: Coastline file for the ROMS Plotting Package
I pretty much doubt that you will be able to do this. These coastlines files are very special and extremely difficult to reproduce. We need closed polygons and the order of values matters to distinguish between continents, islands, and lakes. The filling algorithm in the plotting package requires th...
- Wed Jan 28, 2015 4:18 pm
- Forum: ROMS Tools and Techniques
- Topic: Coastline file for the ROMS Plotting Package
- Replies: 4
- Views: 4291
Re: Coastline file for the ROMS Plotting Package
Yes, that's correct. There are several coastline files at various resolutions including global datasets: global_low.cst, global_crude.cst, and global_inter.cst. Western longitudes and southern latitudes are negative with ranges -180 to 180 and -90 to 90, respectively.
- Thu Jan 22, 2015 5:14 pm
- Forum: ROMS Usage
- Topic: How to generate and use composite grids?
- Replies: 16
- Views: 8821
Re: How to generate and use composite grids?
Although I have not tried passive tracers yet in nesting, the algorithm is designed to work on all tracers (active and passive). We operate on the full 5D tracer array. If works on temperature and salinity, it will work also on passive tracers. As designed, you need passive tracers (like biology) in...
- Wed Jan 21, 2015 6:00 pm
- Forum: ROMS Bugs
- Topic: bug report in BIO_fennel module
- Replies: 3
- Views: 2228
Re: bug report in BIO_fennel module
Nope, what I am saying is that in a model like BIO_FENNEL with so many options for additional biological components, you need to know what the order of biological tracer is in idbio(:) and that's the order that you need follow in LBC(isTvar) . If you don't know what that order is, you need to look a...
- Wed Jan 21, 2015 5:21 pm
- Forum: ROMS Usage
- Topic: How to generate and use composite grids?
- Replies: 16
- Views: 8821
Re: How to generate and use composite grids?
Okay, I need to track this problem in contact.m . I was able to reproduce you error. I will look at this soon to see what is the solution for this out-of-range values when computing the interpolation weights. I also have in my to do list to have a realistic application with composite grids. However,...
- Wed Jan 21, 2015 5:12 pm
- Forum: ROMS Bugs
- Topic: bug report in BIO_fennel module
- Replies: 3
- Views: 2228
Re: bug report in BIO_fennel module
Your first panel is completely illegal and a numerical absurdity; the western boundary is periodic and the eastern boundary closed for zooplankton and others. You better understand what periodicity is. Periodic boundary conditions works in pairs: west/east boundaries and south/north boundaries. ROMS...
- Thu Jan 15, 2015 2:07 pm
- Forum: ROMS Problems
- Topic: 4DVAR observation impact driver problem WC13
- Replies: 3
- Views: 2185
Re: 4DVAR observation impact driver problem WC13
Yes, that fixed the problem.
- Wed Jan 14, 2015 4:32 pm
- Forum: ROMS Problems
- Topic: make_ clim error
- Replies: 5
- Views: 3638
Re: make_ clim error
The only Matlab scripts that we support for processing our version of ROMS are located in the Matlab svn repository: https://www.myroms.org/svn/src/matlab :!: These are the official version of these scripts containing the correct functions to process vertical coordinate schemes, nesting, land/sea ma...
- Thu Jan 08, 2015 2:55 pm
- Forum: ROMS Problems
- Topic: TS_MPDATA and TS_U3HADVECTION_TL... Dsrc not found? I4DVAR
- Replies: 5
- Views: 2911
Re: TS_MPDATA and TS_U3HADVECTION_TL... Dsrc not found? I4DV
I changed the TS_MPDATA to TS_MPDATA_NOT_YET in several TLM , RPM , and ADM routines. It think that this will solve your compilation problem. See :arrow: trac ticket for more information. You need to update. By the way, your header file angola_3km_4dvar.h needs some cleaning, specially when you use ...
- Tue Jan 06, 2015 5:23 pm
- Forum: ROMS Problems
- Topic: TS_MPDATA and TS_U3HADVECTION_TL... Dsrc not found? I4DVAR
- Replies: 5
- Views: 2911
Re: TS_MPDATA and TS_U3HADVECTION_TL... Dsrc not found? I4DV
Yes, it needs to be SOURCES(ng)%Dsrc(is) , thank you :!: That's was easy an obvious. The ARPACK library is needed in the propagators drivers for Generalized Stability Analysis that uses the tangent linear and adjoint models. What are the undefined references in ad_step3d_t ? We don't need the ARPACK...
- Wed Dec 17, 2014 2:49 pm
- Forum: ROMS Bugs
- Topic: nesting_bug
- Replies: 3
- Views: 2300
Re: nesting_bug
I have checked the new released codes carefully for the last two months (countless hours in the debugger) with idealized and a realistic applications in the South China Sea and the US East Coast and they work fine. The one that you were using before have several serious bugs. I don't know what is wr...
- Tue Dec 16, 2014 2:26 pm
- Forum: ROMS Installation
- Topic: Compiling with TS_MPDATA and I4DVAR... Dsrc not found?
- Replies: 3
- Views: 2650
Re: Compiling with TS_MPDATA and I4DVAR... Dsrc not found?
The option TS_MPDATA is not adjointed. It is highly nonlinear. You can still use the TS_MPDATA for the forward nonlinear model and any other option for the adjoint and tangent linear models. In cppdefs.h , we have: ** The flags below are optional. By default, the same options chosen ** ** for the no...
- Tue Dec 16, 2014 12:02 pm
- Forum: ROMS Releases
- Topic: Major update to the Nesting Algorithms
- Replies: 0
- Views: 11794
Major update to the Nesting Algorithms
Recommendations: :idea: This post has lots of detailed information. Please carefully read and study the provided information before loading this forum with lots of questions that are probably answered below. The nesting algorithms are complex and require experience and patience. I highly recommend ...
- Mon Dec 15, 2014 5:18 pm
- Forum: ROMS Information
- Topic: Question about Sea level
- Replies: 5
- Views: 3685
Re: Question about Sea level
This is very simple. You are gaining volume and this is due to the open boundaries. You may plot the domain volume for your solution as a function of time. It seems that using climatology at the boundary is not working for you. Nowadays, we get boundary condition from global models. Imposing climato...
- Tue Dec 09, 2014 3:59 am
- Forum: ROMS Bugs
- Topic: WET_DRY and NESTING (refined) possible bug
- Replies: 6
- Views: 3189
Re: WET_DRY and NESTING (refined) possible bug
I changed the logic how is this done. Please be patient. An upgrade to the nesting algorithms will be released in the next two days. I am working on the documentation.
- Wed Dec 03, 2014 2:14 am
- Forum: ROMS Discussion
- Topic: one problem in run time
- Replies: 7
- Views: 1994
Re: one problem in run time
One person reported much better performance from the ROMS_Agrif code over both the Rutgers ROMS and the COAWST ROMS when using the nesting. You might try that. Actually, this has a good explanation that mostly everybody is not aware. The barotropic engine in either the UCLA ROMS or ROMS_Agrid (I be...
- Tue Nov 25, 2014 2:49 pm
- Forum: ROMS Discussion
- Topic: Latent Heat Flux Problem with ROMS 741
- Replies: 3
- Views: 3023
Re: Latent Heat Flux Problem with ROMS 741
It is up to you what values to use according to the dataset that you have. Notice that in ocean.in you need to specify the height (m) of the atmospheric fluxes: ! Height (m) of atmospheric measurements for Bulk fluxes parameterization. BLK_ZQ == 10.0d0 ! air humidity BLK_ZT == 10.0d0 ! air temperatu...
- Thu Nov 20, 2014 12:07 am
- Forum: ROMS Discussion
- Topic: inert passive tracer boundary not being read
- Replies: 8
- Views: 2500
Re: inert passive tracer boundary not being read
Please update your codes. I corrected this bug in the repository. See following
trac ticket for details.

- Wed Nov 19, 2014 11:03 pm
- Forum: ROMS Tools and Techniques
- Topic: TNU2 and VISC2 values
- Replies: 4
- Views: 2982
Re: TNU2 and VISC2 values
Your question cannot be answered in this forum because pretty much you need to be aware of the full history of Physical Oceanography and Ocean Modeling. There are decades of excellent literature about the subject starting with Henry Stommel, Walter Munk, and others. Any Physical Oceanographer or Oce...
- Fri Nov 14, 2014 9:54 pm
- Forum: ROMS Discussion
- Topic: inert passive tracer boundary not being read
- Replies: 8
- Views: 2500
Re: inert passive tracer boundary not being read
Yes, I know. I has an update for it. Mark Hadfield put a post about this and the solution.
- Thu Nov 13, 2014 6:39 pm
- Forum: ROMS Discussion
- Topic: way of reducing vertical mixing ...
- Replies: 3
- Views: 1755
Re: way of reducing vertical mixing ...
The turbulent schemes like MY-2.5 or GLS tend to give high vertical mixing coefficients if you have strong winds and buoyancy forcing. Try the KPP scheme to see what you get. The mixing is more controlled in LMD_MIXING.
- Fri Nov 07, 2014 7:33 pm
- Forum: ROMS Problems
- Topic: tidal forcing file: inconsistent restart and harmonic time
- Replies: 2
- Views: 1931
Re: tidal forcing file: inconsistent restart and harmonic ti
I don't know why you have to do this. I far as recall, we don't need to change the input script except for the restart flag and the initial conditions NetCDF: NRREC == -1 ... ININAME == ocean_rst.nc You shouldn't change: DSTART = xxxxx ! days TIDE_START = xxxxx ! days TIME_REF = xxxxx ! yyyymmdd.dd ...
- Thu Nov 06, 2014 4:04 pm
- Forum: ROMS Discussion
- Topic: Compiling 'AVERAGES_TIDE'
- Replies: 13
- Views: 4322
Re: Compiling 'AVERAGES_TIDE'
Come on people... With a little of curiosity and using " grep -r Hcount * | grep -v .svn " in the top of the ROMS tree directory you will find where Hcount is used. You will noticed that it is defined in def_tides.F and written in wrt_tides.F . If you do another search for def_tides you will find th...
- Thu Nov 06, 2014 3:42 pm
- Forum: ROMS Bugs
- Topic: PROFILE report
- Replies: 4
- Views: 3036
Re: PROFILE report
This is in my to do list. I had noticed it too. The logic is delicate because of shared-memory. I am busy with other codes now.
- Sat Nov 01, 2014 4:55 pm
- Forum: ROMS Bugs
- Topic: Problem with compact tracer Rad + Nud boundaries
- Replies: 2
- Views: 2502
Re: Problem with compact tracer Rad + Nud boundaries
Yes Mark, good catch
Thank you.

- Wed Oct 29, 2014 2:30 am
- Forum: ROMS Discussion
- Topic: how to active ADJUST_BOUNDARY
- Replies: 2
- Views: 946
Re: how to active ADJUST_BOUNDARY
This is only used in the 4D-Var algorithms to include the boundary conditions in the control vector. You just need to activate ADJUST_BOUNDARY. Additional input files are required for this option. See tutorials.
- Wed Oct 29, 2014 2:25 am
- Forum: ROMS Discussion
- Topic: ROMS Boundary Condition: Clamped vs. Flather
- Replies: 4
- Views: 1721
Re: ROMS Boundary Condition: Clamped vs. Flather
Apparently, you are using a very old version the ROMS. All the CPP for open boundary conditions were removed and replaced with logical switches in ocean.in to facilitate nesting.
- Wed Oct 29, 2014 2:22 am
- Forum: ROMS Bugs
- Topic: a small bug in inp_par.F
- Replies: 2
- Views: 2080
Re: a small bug in inp_par.F
Yes, good catch. Thank you.
- Sat Oct 25, 2014 11:51 pm
- Forum: ROMS Problems
- Topic: Unexpected(?) blow up
- Replies: 7
- Views: 3793
Re: Unexpected(?) blow up
My experience with blows-up of this type after several years running is due to atmospheric forcing, specially wind stress due to strong storms. By the way, this happen to us in the Mediterranean part of our grids too. In the Aegean, there are strong winds sometimes. What we usually do is restart wit...
- Wed Oct 22, 2014 4:35 pm
- Forum: ROMS Problems
- Topic: wetting and drying problem
- Replies: 22
- Views: 12115
Re: wetting and drying problem
Applying the wetting and drying mask in step_t or to any tracer anywhere is really bad :!: Lyon, this was actually the problem that you were having with the perfect restart in your Cook Inlet application. We were multiplying with such mask during I/O. I mentioned this before. If a particular cell is...
- Mon Oct 20, 2014 4:57 pm
- Forum: ROMS Bugs
- Topic: WET_DRY and NESTING (refined) possible bug
- Replies: 6
- Views: 3189
Re: WET_DRY and NESTING (refined) possible bug
Well, this is a hard problem. The interpolation weights need to be adjusted if you get wet and dry cells at every time step because of the masking. I need to look for a solution for this case. Bringing refinement grids all the way to the coast with wetting and drying requires special considerations....
- Mon Oct 20, 2014 4:46 pm
- Forum: ROMS Tools and Techniques
- Topic: nested grid layer and min depth
- Replies: 7
- Views: 3912
Re: nested grid layer and min depth
It has been mentioned several times in this forum that more than two nesting layers for refinement grids still need some debugging for NestLayers > 2 . Your configuration with four grids is very poor and ill-posed in my opinion. Obviously, you need to gain more experience with nested grids and ocean...
- Mon Oct 20, 2014 4:30 pm
- Forum: ROMS Bugs
- Topic: latest ROMS version problem
- Replies: 6
- Views: 3931
Re: latest ROMS version problem
You are getting this error because you are using an older version of the metadata file varinfo.dat . You need to use the updated file from the latest version of the code. Notice that I added to varinfo.dat the metadata for the wetting and drying mask at PSi-points that it is needed for restart: 'wet...
- Mon Oct 20, 2014 4:22 pm
- Forum: ROMS Bugs
- Topic: tiny nesting/mp_aggregate2d bug
- Replies: 1
- Views: 1517
Re: tiny nesting/mp_aggregate2d bug
Yes, good catch. Thank you. I am still debugging this routines for various configurations.
- Mon Oct 13, 2014 5:35 pm
- Forum: ROMS Bugs
- Topic: WET_DRY and Dcrit
- Replies: 13
- Views: 7838
Re: WET_DRY and Dcrit
I can see how the rotated tensor for tracer diffusion can go unstable near wetting and drying areas. We cannot numerically develop a wetting and drying algorithm from first principles in a model like ROMS, but we can use an engineering approach. I think that if you need horizontal diffusion in your ...
- Mon Oct 13, 2014 4:50 pm
- Forum: ROMS Problems
- Topic: wetting and drying problem
- Replies: 22
- Views: 12115
Re: wetting and drying problem
Hi Kate, I think that my big update to the wetting and drying will be released this week and the restart problem is now solved :D I am still running NOAA's Cook Inlet application but very slowly :( The grid is huge ( 1044x724x30 ) and I am running in my 8-CPUs desktop in case that I need to examine ...
- Wed Oct 08, 2014 8:58 pm
- Forum: ROMS Discussion
- Topic: initial-nest
- Replies: 13
- Views: 4557
Re: initial-nest
The answer for both questions is that you need separate files. Nesting is complex and users need to explore first with the larger grid before trying to run with nested grid. It requires expertise to set the nested grids correctly. Users need lots of experience with ROMS, ocean modeling, and the ocea...
- Tue Oct 07, 2014 1:09 am
- Forum: ROMS Installation
- Topic: How can I write concise global attributes into netcdf files
- Replies: 2
- Views: 2128
Re: How can I write concise global attributes into netcdf fi
This is absurd. In never heard anything like it. The buffer for the global attributes is extremely smaaaaaaall compared with the data. Obviously, you need to read more about the NetCDF files and its unique format. You cannot remove the global attributes in any of the ROMS history file because the mo...
- Mon Sep 01, 2014 5:34 pm
- Forum: ROMS Problems
- Topic: 4DVAR observation impact driver problem WC13
- Replies: 3
- Views: 2185
Re: 4DVAR observation impact driver problem WC13
I don't know what it is the problem here. The CLIMA(:) structure has a -1 index. That's is kind of weird. I will need to reproduce the problem. I added it to my to-do list.
- Tue Aug 19, 2014 5:59 pm
- Forum: ROMS Bugs
- Topic: Test cases failing to compile in ROMS 3.7
- Replies: 2
- Views: 2168
Re: Test cases failing to compile in ROMS 3.7
Indeed, thank you for reporting this.
- Mon Jun 30, 2014 4:17 pm
- Forum: ROMS Discussion
- Topic: tidal analysis of large output files
- Replies: 22
- Views: 8345
Re: tidal analysis of large output files
It is possible that your solution is not long enough to separate all the 9 tidal frequencies in your application. Recall that we have in ROMS the nonlinear interaction between all the frequencies and resonance. If you have an application with just M 2 and S 2 , it will take around 30 days of solutio...
- Fri May 16, 2014 3:54 pm
- Forum: Meetings/Workshops
- Topic: 2014 ROMS/TOMS User Workshop, Rovinj, Croatia - Final Notice
- Replies: 0
- Views: 1872
2014 ROMS/TOMS User Workshop, Rovinj, Croatia - Final Notice
2014 ROMS/TOMS User Workshop Island Hotel Istra, Rovinj, Croatia, May 26 - 29, 2014 Please check the following link for the workshop agenda: http://www.myroms.org/index.php?page=ROMS_2014_agenda Notice that it includes oral talks, lectures, and poster sessions. The oral presentations are 30 minutes...
- Tue May 06, 2014 3:25 pm
- Forum: ROMS Discussion
- Topic: 3-level nesting
- Replies: 3
- Views: 1479
Re: 3-level nesting
Yes, thank you. It is in my list of things to do. I am updating the script contact.m . In the mean time, you may experiment with the two-level nesting. It takes awhile to get the grids just right. So you need to run and play with the configuration that works the best. The placement of nested grids i...
- Tue May 06, 2014 12:29 pm
- Forum: Meetings/Workshops
- Topic: 2014 ROMS/TOMS Workshop T-shirt
- Replies: 0
- Views: 1845
2014 ROMS/TOMS Workshop T-shirt
This is out current design of the ROMS T-shirt for the upcoming User Workshop. The symbol on the front is the letter M (mislite) from the Glagolitic alphabet and means thought or think. The alternating red and white symbolizes the well-known red and white checkerboard from the Croatian coat of arms ...
- Tue Apr 29, 2014 2:32 pm
- Forum: ROMS Discussion
- Topic: different restart files
- Replies: 7
- Views: 2466
Re: different restart files
You have LcycleRST = F . This implies that you will have a single restart file with multiple time records (in your case, each record is written every 1920 time-steps). The restart records are not recycled to have the latest two records in the restart file. Therefore, the restart file will grow with ...
- Fri Apr 18, 2014 10:22 pm
- Forum: ROMS Bugs
- Topic: Potential bug in ROMS perfect restart
- Replies: 16
- Views: 6173
Re: Potential bug in ROMS perfect restart
I don't know. I will need to reproduce the problem to check if this is the case.
- Fri Apr 18, 2014 8:00 pm
- Forum: ROMS Bugs
- Topic: Potential bug in ROMS perfect restart
- Replies: 16
- Views: 6173
Re: Potential bug in ROMS perfect restart
The correction suggested by Kate above is incorrect. We ALWAYS need to use the rmask_io, umask_io, and vmask_io in the output NetCDF files
This is even more crucial when we have rivers and wetting and drying...

- Fri Apr 11, 2014 4:08 pm
- Forum: ROMS Usage
- Topic: Construction of forcing file (stationary forcing)
- Replies: 5
- Views: 3200
Re: Construction of forcing file (stationary forcing)
The cycle_length attribute is for other purposes. There is even a better solution for perpetual forcing. You don't need two-time records and avoid all that unnecessary time interpolation at every time-step :!: You guys really need to read the information that I put in the ticket changes to the repos...
- Thu Apr 10, 2014 10:33 pm
- Forum: ROMS Bugs
- Topic: Possibly 3 nesting-related bugs
- Replies: 16
- Views: 6792
Re: Possibly 3 nesting-related bugs
Of course, this makes a lot of sense to me. We cannot apply volume conservation in a regional tidally forced application :!: Tides do not conserve volume in such applications. If you do so, you will alter the tidal signal. Therefore, this explains very well your pink curve. So now you know what happ...
- Thu Apr 03, 2014 9:01 pm
- Forum: ROMS Bugs
- Topic: Beware of netCDF STRINGs
- Replies: 6
- Views: 5112
Re: Beware of netCDF STRINGs
Yes, I read your postings in the NetCDF forum. This is an annoying one. I always use the CHAR flags. Characters and strings are awkward in NetCDF. I am pretty conservative about the new features in the NetCDF library. I struggled a lot with them when designing the NetCDF file for the nesting contact...
- Thu Apr 03, 2014 6:35 pm
- Forum: ROMS Discussion
- Topic: 3-level nesting
- Replies: 3
- Views: 1479
Re: 3-level nesting
Okey, I haven't tried the three-level case yet. I did tried 3-grids and two-levels. I will need to reproduce this problem. I will put in my list of things to do. I hope that it is simple to fix. The order of how you specify the grid file names is extremely important (coarser to finer). Thank you for...
- Thu Apr 03, 2014 5:48 pm
- Forum: ROMS Bugs
- Topic: Possibly 3 nesting-related bugs
- Replies: 16
- Views: 6792
Re: Possibly 3 nesting-related bugs
Yes, the calls to the exchange routines exchange_v2d_tile and mp_exchange3d need to be corrected as you mentioned above. However, your corrections to get_data and set_data are incorrect :!: We still can use climatology fields in refined grids for other purposes :!: Notice that the processing of clim...
- Thu Apr 03, 2014 5:15 pm
- Forum: ROMS Bugs
- Topic: missing routine
- Replies: 7
- Views: 2820
Re: missing routine
Okey, all my refinement test cases have a ratio of 1:3 . I will need to have a test case with a 1:5 ratio to reproduce this problem and check in the debugger for a fix to the out-of-bounds error. It shouldn't that difficult to fix. This is only used in coupled routines. I think that something is wro...
- Thu Apr 03, 2014 5:05 pm
- Forum: ROMS Bugs
- Topic: missing routine
- Replies: 7
- Views: 2820
Re: missing routine
Yes, the uniform field in the grid structure was added recently. The fix in coarse2fine.m and fine2coarse.m is very simple. ITS value should be available and inherited from the donor grid.
Thank you for bringing this to my attention. Please update.
Thank you for bringing this to my attention. Please update.
- Tue Apr 01, 2014 12:13 pm
- Forum: ROMS Discussion
- Topic: How to define the input forcing files when having many?
- Replies: 2
- Views: 1969
Re: How to define the input forcing files when having many?
:idea: The instructions are very clear in all input scripts. You just have to read :!: ! Input parameters can be entered in ANY order, provided that the parameter ! ! KEYWORD (usually, upper case) is typed correctly followed by "=" or "==" ! ! symbols. Any comment lines are allowed and must begin wi...
- Tue Apr 01, 2014 2:55 am
- Forum: ROMS Bugs
- Topic: missing routine
- Replies: 7
- Views: 2820
Re: missing routine
The gcircle.m script is in the matlab/utility directory. The G.uniform should be set-up in get_roms_grid.m . Do you have the latest version of this script? or perhaps it is shadowed by an old version in another disrectory in your path. If not, it missed something about this. I always get the uniform...
- Fri Mar 28, 2014 10:03 pm
- Forum: ROMS Bugs
- Topic: 2 possible bugs in Shc boundary conditions
- Replies: 6
- Views: 4203
Re: 2 possible bugs in Shc boundary conditions
Indeed, I fixed those bugs. Please update. Thank you for bringing this to my attention. I also corrected the tangent linear, representer, and adjoint versions of u2dbc_im.F and v2dbc_im.F.
- Fri Mar 28, 2014 6:26 pm
- Forum: Meetings/Workshops
- Topic: 2014 ROMS/TOMS User Workshop, Rovinj, Croatia - 4th Notice
- Replies: 0
- Views: 1664
2014 ROMS/TOMS User Workshop, Rovinj, Croatia - 4th Notice
2014 ROMS/TOMS User Workshop Island Hotel Istra, Rovinj, Croatia, May 26 - 29, 2014 The deadline to submit abstracts for oral presentations and posters is only a week away on April 4, 2014 . Please submit your abstract by that date. The registration deadline for the 2014 ROMS/TOMS User Workshop was...