Search found 3883 matches

by kate
Wed Feb 17, 2021 5:53 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

All those Dout flags control the output, not the tides. I have this section in my .h file: /* tides */ #define LTIDES #ifdef LTIDES # if defined AVERAGES && !defined USE_DEBUG # define FILTERED # endif # define SSH_TIDES # define UV_TIDES # define ADD_FSOBC # define ADD_M2OBC # undef RAMP_TIDES # de...
by kate
Wed Feb 17, 2021 5:48 pm
Forum: ROMS Discussion
Topic: Discharge a chemical substance
Replies: 3
Views: 363

Re: Discharge a chemical substance

There are people who study oil spills, maybe not on this forum though. What I know of it is that they simulate a spill for all those cases and see where it goes. For some hypothetical spill, you want to be prepared for all the cases, all the sea levels, all the wind directions. Is it from a leaky ba...
by kate
Sat Feb 13, 2021 2:29 pm
Forum: ROMS Tools and Techniques
Topic: Error to plot by the ROMS plotting package
Replies: 10
Views: 756

Re: Error to plot by the ROMS plotting package

You need to clean out the old build and start fresh. Getting that error in debug mode should give you more useful information than that.
by kate
Fri Feb 12, 2021 7:39 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

Simply add statements such as these marked with the plus signs: diff --git a/ROMS/Nonlinear/initial.F b/ROMS/Nonlinear/initial.F index 6e20ae7d9..25fd9c4e6 100644 --- a/ROMS/Nonlinear/initial.F +++ b/ROMS/Nonlinear/initial.F @@ -677,6 +677,7 @@ END DO END IF #endif + print, *, "after get_state" #ifd...
by kate
Thu Feb 11, 2021 6:36 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

Once the model gets the initial conditions, it should start reading in your surface forcing fields. I see you have BULK_FLUXES, so it should read in shortwave, longwave, winds, etc. There's no mention of that in your output so I have no idea what went wrong. There is a listing of your input forcing ...
by kate
Wed Feb 10, 2021 7:13 pm
Forum: ROMS Discussion
Topic: How to setup the river discharges with only salinity but no temperature?
Replies: 6
Views: 521

Re: How to setup the river discharges with only salinity but no temperature?

For numerical stability you must give it something other than the water temperature at that grid point. Using the water temperature there is an unconditionally unstable numerical scheme - downwind advection. In MOM6, the default is to use air temperature instead, something you'd have anyway if you w...
by kate
Wed Feb 10, 2021 6:35 pm
Forum: ROMS Discussion
Topic: How to setup the river discharges with only salinity but no temperature?
Replies: 6
Views: 521

Re: How to setup the river discharges with only salinity but no temperature?

It doesn't make sense numerically to do what you want. As you see, the thing will become unstable. Lacking river temperatures, some have used air temperature instead.
by kate
Mon Feb 08, 2021 7:14 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

I would use the NCO package for this job, specifically the "ncatted" command. Something like:

Code: Select all

ncatted -a time,,m,c,'time'
Note that the missing value between two commas means apply to all variables.
by kate
Mon Feb 08, 2021 5:20 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

You just need to change some attributes in your frc file. You now have: double time(time) ; time:long_name = "time since intialization" ; time:units = "days since 2011-01-01 00:00:00" ; time:calendar = "gregorian" ; float Uwind(time, lat, lon) ; Uwind:long_name = "surface u-wind component (east) at ...
by kate
Sun Feb 07, 2021 10:48 pm
Forum: ROMS Installation
Topic: Compilation error for plotting package
Replies: 3
Views: 219

Re: Compilation error for plotting package

make: /opt/gfortransoft/serial/netcdf4/bin/nc-config: Command not found
If you have a modern netcdf-c and netcdf-fortran correctly installed, nc-config should be on your search path. Did you install them? What operating system is this? Your laptop or a major computing platform?
by kate
Thu Feb 04, 2021 6:01 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

From your frc file:

Code: Select all

	double time(time) ;
		time:long_name = "time since intialization" ;
		time:units = "days since 2011-01-01 00:00:00" ;
		time:calendar = "gregorian" ;
Can you use 20110101 as your time origin for everything else?
by kate
Tue Feb 02, 2021 6:09 am
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

I always have a TIME_REF of 19000101, which is year 1900, day 1. What year is your time origin? My DSTART is the offset in days from that time, 29225 for 1980, something larger for more recent times. Again, how does the time look in the forcing and what are its attributes?
by kate
Tue Feb 02, 2021 4:00 am
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

The forcing times need to be along the same time axis as everything else. What are the times in your files? It might be easier to make ROMS consistent with them rather than vice versa.
by kate
Sat Jan 30, 2021 3:48 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 23
Views: 1142

Re: Error while executing ROMS for a self created grid region

NFFILES == 1 ! number of unique forcing files FRCNAME == frc_TPXO9_alderneyrace_09Dec2068.nc | ! forcing file 1, grid 1 (tidal forcing) frc_watl_ERA5_bulkflux_new_201805.nc !(wind forcing) This should instead be: NFFILES == 2 ! number of unique forcing files FRCNAME == frc_TPXO9_alderneyrace_09Dec2...
by kate
Thu Jan 28, 2021 5:11 pm
Forum: ROMS Discussion
Topic: How to run COAWST model with real-time data
Replies: 1
Views: 175

Re: How to run COAWST model with real-time data

I went through this exercise in pyroms, now in the pyroms/examples/Arctic_HYCOM_GLBy directory.
by kate
Wed Jan 27, 2021 4:59 pm
Forum: ROMS Discussion
Topic: Use Bulk Formula or Net Flux Form
Replies: 4
Views: 304

Re: Use Bulk Formula or Net Flux Form

I personally always use bulk_fluxes, in part because the sea ice model requires it.

If you don't use bulk_fluxes, it is advised to also use the dQdSST term to account for your SST differing from that used to compute the fluxes.
by kate
Tue Jan 26, 2021 5:54 pm
Forum: ROMS Discussion
Topic: Use Bulk Formula or Net Flux Form
Replies: 4
Views: 304

Re: Use Bulk Formula or Net Flux Form

Think about say the Gulf Stream or the Kuroshio. If your fluxes were computed when it was wobbling south, colder water would be present and the fluxes would reflect that. If your Gulf Stream was instead wobbling north, the fluxes would therefore be off. Instead, using BULK_FLUXES, you get the fluxes...
by kate
Mon Jan 25, 2021 6:29 pm
Forum: ROMS Discussion
Topic: problem in svn and nf_config
Replies: 13
Views: 530

Re: problem in svn and nf_config

Code: Select all

NETCDF_CREATE - unable to create output NetCDF file:
roms_his.nc
call from: ROMS/Utility/def_his.F
NetCDF: Malformed URL
Found Error: 03 Line: 105 Source: ROMS/Utility/def_his.F
I don't think your problem is in varinfo.dat.
by kate
Mon Jan 25, 2021 5:29 pm
Forum: ROMS Discussion
Topic: problem in svn and nf_config
Replies: 13
Views: 530

Re: problem in svn and nf_config

Wow, I've never seen that one before. I would check for hidden characters in the ocean.in file where you provide the file names. Or is that file somehow in DOS mode?
by kate
Sun Jan 24, 2021 8:25 pm
Forum: ROMS Discussion
Topic: problem in svn and nf_config
Replies: 13
Views: 530

Re: problem in svn and nf_config

Code: Select all

chinook02.rcs.alaska.edu 203% grep id3dPVI var*
varinfo.dat:  'id3dPVI'
varinfo.dat.cobalt:  'id3dPVI'
varinfo_jra.dat:  'id3dPVI'
I ran this in the Apps directory.
by kate
Sun Jan 24, 2021 7:32 pm
Forum: ROMS Discussion
Topic: problem in svn and nf_config
Replies: 13
Views: 530

Re: problem in svn and nf_config

This is due to a change in the varinfo.dat file. Search through mine for that string idStr and you'll find it.
by kate
Sun Jan 24, 2021 7:42 am
Forum: ROMS Discussion
Topic: problem in svn and nf_config
Replies: 13
Views: 530

Re: problem in svn and nf_config

Just change that path to my nf-config to your path to your nf-config.
by kate
Sat Jan 23, 2021 7:21 pm
Forum: ROMS Discussion
Topic: problem in svn and nf_config
Replies: 13
Views: 530

Re: problem in svn and nf_config

The version of netcdf-fortran doesn't matter, but it has to be compiled with the same version of gfortran as you are using. Just search for "netcdf-fortran download" and you find it. Then you need to follow the instructions for compiling it. Or if you have control of the computer, something like "su...
by kate
Sat Jan 23, 2021 7:03 pm
Forum: ROMS Discussion
Topic: problem in svn and nf_config
Replies: 13
Views: 530

Re: problem in svn and nf_config

I don't know that the svn error matters.
/import/home/kshedstrom/bin/nf-config
This is a hard-coded path on a computer here to my home-built netcdf-fortran. You'll have to change it.
by kate
Wed Jan 20, 2021 6:13 pm
Forum: ROMS Discussion
Topic: cmpiling error when using roms-kate_svn
Replies: 1
Views: 154

Re: cmpiling error when using roms-kate_svn

It is failing to find the netcdf library, part of the build system. Since I don't use Hernan's build script, you'll have to figure out what's different between the two and do the appropriate thing.
by kate
Tue Jan 19, 2021 4:53 pm
Forum: ROMS Ecosystem
Topic: Problems with river input forcing when using NEMURO
Replies: 7
Views: 788

Re: Problems with river input forcing when using NEMURO

Do you have coastal river sources? Are you providing values of the BGC tracers on them? They can easy go unstable if you don't provide values for all incoming tracers. Look at the global NEWS product for nutrient values. Don't worry about varinfo.dat. Best practice is to have the name of the time va...
by kate
Mon Jan 18, 2021 6:03 pm
Forum: ROMS Ice
Topic: How add ice
Replies: 6
Views: 699

Re: How add ice

You should read the manual. It has sections on the ice and also on how to use git.
by kate
Sun Jan 17, 2021 10:30 pm
Forum: ROMS Ice
Topic: How add ice
Replies: 6
Views: 699

Re: How add ice

My code is at: https://github.com/kshedstrom/roms. I'm not now finding an svn option, but you can download a tarball of it from that site. It has been over ten years since I've used svn...
by kate
Sun Jan 17, 2021 7:50 pm
Forum: ROMS Ice
Topic: How add ice
Replies: 6
Views: 699

Re: How add ice

There is no ice in the ROMS trunk code. There is ice in my code on github, also in COAWST from John Warner. There is also ice in the metroms code on github, where they couple to a more complete ice model (CICE). As for forcing, if you are reading winds, air temp, etc. and using BULK_FLUXES, you are ...
by kate
Thu Jan 14, 2021 5:55 pm
Forum: ROMS Discussion
Topic: Curvilinear grid for ROMS
Replies: 4
Views: 628

Re: Curvilinear grid for ROMS

Yes, I can make a grid like that. However, you want your corners to be right angles. Also, a boundary across the pack ice will be hard to make well-behaved. We had a grid with a boundary at Bering Strait and we didn't like how that behaved either. I guess it depends on what your goals are. How to do...
by kate
Wed Jan 13, 2021 6:06 pm
Forum: ROMS Ecosystem
Topic: Problems with river input forcing when using NEMURO
Replies: 7
Views: 788

Re: Problems with river input forcing when using NEMURO

This is when you need to build up some debugging skills. You need to find out if the bio fields are sensible in the model, getting corrupted on output, or if they're just nuts everywhere. Feel free to add print statements at will.
by kate
Tue Jan 12, 2021 5:18 pm
Forum: ROMS Ecosystem
Topic: Problems with river input forcing when using NEMURO
Replies: 7
Views: 788

Re: Problems with river input forcing when using NEMURO

It has to be something in the varinfo.dat. I seem to have both: 'river_fed' 'fed river source' 'mol.m-2.s-1' 'river_fed, scalar, series' 'runoff_time' 'idriver_fed' 'r2dvar' 1.0d0 and 'river_fed' 'river runoff Dissolved Iron' 'mol/kg' 'fed, scalar, series' 'ocean_time' 'idRtrc(ifed)' 'r3dvar' 1.0d0 ...
by kate
Sat Jan 09, 2021 2:27 pm
Forum: ROMS Tools and Techniques
Topic: NCL
Replies: 16
Views: 6905

Re: NCL

This is not the time to start learning NCL. Support for NCL development has been dropped and the focus is now on developing Python tools instead. While we have some plotting capability in pyroms, it too is now out of date, what with the rapid development in new Python tools. I am hardly the expert o...
by kate
Sat Jan 09, 2021 2:17 pm
Forum: ROMS Tools and Techniques
Topic: plot package problem
Replies: 21
Views: 5607

Re: plot package problem

by kate
Fri Jan 08, 2021 5:09 pm
Forum: ROMS Tools and Techniques
Topic: plot package problem
Replies: 21
Views: 5607

Re: plot package problem

unknown option: --flib Usage: nc-config
Try --flibs instead of --flib
by kate
Fri Jan 08, 2021 5:36 am
Forum: ROMS Tools and Techniques
Topic: plot package problem
Replies: 21
Views: 5607

Re: plot package problem

There's a freetype library too, I would guess. Perhaps a quick web search is what you need.

The warnings are a little concerning too, like maybe there's a version inconsistency.
by kate
Thu Jan 07, 2021 5:51 pm
Forum: ROMS Tools and Techniques
Topic: plot package problem
Replies: 21
Views: 5607

Re: plot package problem

cairo is a library which you'll need to link to as well. Try adding -lcairo to your link command.
by kate
Wed Jan 06, 2021 6:46 pm
Forum: ROMS Discussion
Topic: Erorr when applying rivers
Replies: 8
Views: 660

Re: Erorr when applying rivers

Here's what one of my rivers files looks like: netcdf JRA-1.4_Arctic_rivers_1980 { dimensions: river_time = UNLIMITED ; // (366 currently) river = 18430 ; s_rho = 50 ; variables: double river_time(river_time) ; river_time:units = "days since 1900-01-01 00:00:00" ; river_time:long_name = "river runof...
by kate
Tue Jan 05, 2021 9:23 pm
Forum: ROMS Discussion
Topic: Erorr when applying rivers
Replies: 8
Views: 660

Re: Erorr when applying rivers

Yes, exactly. How much do you want shallow and how much deep? For the coastal Gulf of Alaska, I use a linear profile, max at the surface, zero at the very bottom.
by kate
Tue Jan 05, 2021 8:03 pm
Forum: ROMS Discussion
Topic: Erorr when applying rivers
Replies: 8
Views: 660

Re: Erorr when applying rivers

You could well be dividing by zero because it doesn't have the Vshape. Yes, you need it to be non-zero.
by kate
Thu Dec 31, 2020 3:01 am
Forum: ROMS Problems
Topic: Issues when running ROMS coupled with bio_UMAINE15 model
Replies: 14
Views: 1223

Re: Issues when running ROMS coupled with bio_UMAINE15 model

If you don't have a working debugger, I would add print statements in there to make sure everything is as you expect it to be. You can add a print of min(pm), max(pm), min(pn), max(pn). Whatever it takes.
by kate
Wed Dec 30, 2020 7:17 pm
Forum: ROMS Problems
Topic: Issues when running ROMS coupled with bio_UMAINE15 model
Replies: 14
Views: 1223

Re: Issues when running ROMS coupled with bio_UMAINE15 model

758 CALL mp_exchange2d (ng, tile, model, 2, & 759 & LBi, UBi, LBj, UBj, & 760 & NghostPoints, & 761 & EWperiodic(ng), NSperiodic(ng), & 762 & pm, pn) It is unhappy with your pm, pn values.
by kate
Tue Dec 29, 2020 5:45 pm
Forum: ROMS Problems
Topic: Issues when running ROMS coupled with bio_UMAINE15 model
Replies: 14
Views: 1223

Re: Issues when running ROMS coupled with bio_UMAINE15 model

It means it's trying to pass a bad value in an MPI call. You need to look up a level to:
758 analytical.f90
by kate
Mon Dec 28, 2020 5:32 pm
Forum: ROMS Problems
Topic: Issues when running ROMS coupled with bio_UMAINE15 model
Replies: 14
Views: 1223

Re: Issues when running ROMS coupled with bio_UMAINE15 model

Now you need to make sure idpthR is in your varinfo.dat.
by kate
Sun Dec 27, 2020 6:09 pm
Forum: ROMS Problems
Topic: Issues when running ROMS coupled with bio_UMAINE15 model
Replies: 14
Views: 1223

Re: Issues when running ROMS coupled with bio_UMAINE15 model

So now what's on this line?
read_phypar_ 2611 read_phypar.f90
If you're getting your code from Professor Peng, shouldn't you be asking him how to get it to run?
by kate
Fri Dec 25, 2020 5:43 pm
Forum: ROMS Problems
Topic: Issues when running ROMS coupled with bio_UMAINE15 model
Replies: 14
Views: 1223

Re: Issues when running ROMS coupled with bio_UMAINE15 model

I did run CoSine successfully, but it has been some years since I've tried. Claudine Hauri and I added an iron limitation to it.

As for your line 2002, you still need to make sure your varinfo.dat has everything you need in it. Did professor Peng give you a working one?
by kate
Thu Dec 24, 2020 7:44 pm
Forum: ROMS Problems
Topic: Issues when running ROMS coupled with bio_UMAINE15 model
Replies: 14
Views: 1223

Re: Issues when running ROMS coupled with bio_UMAINE15 model

2002 mod_ncparam.f90 What is on line 2002 of your mod_ncparam.f90? Mine has: CASE ('idRtrc(icadet_calc)') You then need to see if your equivalent to idRtrc(icadet_calc) is in your varinfo.dat. The only varinfo.dat I have that might contain the CoSine stuff is this one: https://github.com/kshedstrom...
by kate
Wed Dec 23, 2020 6:58 pm
Forum: ROMS Discussion
Topic: blowing up error
Replies: 15
Views: 1056

Re: blowing up error

Do you want to have ice in the model? There are numerous options to do so, including my branch, COAWST, Metroms, and Scott Durski's branch. Metroms (on github) has a version coupled to the community CICE model, the most complete/mature of the ice models. On the other hand, before I had ice in the mo...
by kate
Wed Dec 23, 2020 5:45 pm
Forum: ROMS Discussion
Topic: Velocity is too large near open boundary
Replies: 8
Views: 645

Re: Velocity is too large near open boundary

There are reanalyses out there at better than monthly, starting before 1979. You might have to contact the authors and mail them some harddrives to load it onto to get it, though. A colleague of mine has done that to get both SODA and Mercator.
by kate
Tue Dec 22, 2020 5:55 pm
Forum: ROMS Discussion
Topic: Velocity is too large near open boundary
Replies: 8
Views: 645

Re: Velocity is too large near open boundary

You haven't told us anything about your boundary conditions except which options you have. What is the frequency of providing new information and what is the source of it? I remember seeing a vast improvement when going from monthly climatology of just T and S to 5-daily SODA fields. The climatology...
by kate
Tue Dec 22, 2020 5:34 pm
Forum: ROMS Discussion
Topic: blowing up error
Replies: 15
Views: 1056

Re: blowing up error

Yes, exactly. However, you know how you have the #ifdef GLS_MIXING... section? Here's what I have for LMD: # define LMD_MIXING # ifdef LMD_MIXING # define LMD_RIMIX # define LMD_CONVEC # define LMD_SKPP # define LI_FOX_KEMPER # undef LMD_BKPP # define LMD_NONLOCAL # define LMD_SHAPIRO # define LMD_D...
by kate
Tue Dec 22, 2020 7:48 am
Forum: ROMS Discussion
Topic: Zeta data in boundary files can't be read
Replies: 2
Views: 296

Re: Zeta data in boundary files can't be read

What is the exact error you are getting?
by kate
Tue Dec 22, 2020 7:44 am
Forum: ROMS Discussion
Topic: blowing up error
Replies: 15
Views: 1056

Re: blowing up error

I bet you have some vertically unstable T and S profiles. If they are what you want, you should turn on a vertical mixing scheme like GLS or LMD.
by kate
Sun Dec 20, 2020 9:55 pm
Forum: ROMS Discussion
Topic: blowing up error
Replies: 15
Views: 1056

Re: blowing up error

You don't have an open boundary. Where is the velocity always 1.4 m/s? Do you expect it to be less? Is that where the blow-up is happening or not? Actually, it reports: REASON: MaxSpeed = 6.2E+03 Rather a larger velocity. So where is that and why? I usually use U3H advection and some fourth-order ve...
by kate
Fri Dec 18, 2020 11:36 pm
Forum: ROMS Discussion
Topic: blowing up error
Replies: 15
Views: 1056

Re: blowing up error

This is when you need to start looking at the model's netCDF output as well as just the text. Where is the model blowing up and in what way? Is it a boundary problem or not?
by kate
Wed Dec 16, 2020 9:39 pm
Forum: ROMS Discussion
Topic: DYNAMICS OPTION: Eulerian Mass Coordinate forrtl: severe (174): SIGSEGV, segmentation fault occurred
Replies: 18
Views: 1115

Re: DYNAMICS OPTION: Eulerian Mass Coordinate forrtl: severe (174): SIGSEGV, segmentation fault occurred

3082 read_phypar.f90
What is on this line in your code? You say you are using COAWST, but where did you get your roms.in file? Does it have lines for how many boundary and climatology files you have? NCLMFILES and NBCFILES, also NSSFFILES. These need to be non-zero.
by kate
Wed Dec 16, 2020 9:35 pm
Forum: ROMS Installation
Topic: git clone
Replies: 11
Views: 788

Re: git clone

Sorry, maybe Dave can answer? I only use git via github.
by kate
Tue Dec 15, 2020 5:50 pm
Forum: ROMS Installation
Topic: git clone
Replies: 11
Views: 788

Re: git clone

Did you do a "git clone" first? You can checkout different branches/tags once you have the clone.
by kate
Mon Dec 14, 2020 12:41 am
Forum: ROMS Documentation
Topic: Abnormal in boundary
Replies: 2
Views: 249

Re: Abnormal in boundary

The people who run climate models at one degree resolution have spent *years* working on their solutions. ROMS was designed to be run at more eddy-resolving scales while the climate models have parameterizations to try to take care of what the eddies would be doing. Is one of your ovals on the equat...
by kate
Fri Dec 11, 2020 6:09 pm
Forum: ROMS FAQ
Topic: Segmentation Faults
Replies: 9
Views: 4680

Re: Segmentation Faults

Oh gosh, that was two and a half years ago! I have no memory of it whatsoever.

As for metroms, I haven't played with that lately either. I might have to go back to it if I can't get this other monster (CESM) working on our supercomputer.
by kate
Fri Dec 11, 2020 5:57 pm
Forum: ROMS Discussion
Topic: Running error
Replies: 1
Views: 152

Re: Running error

All of your times have to be consistent, however you choose to do that. Pick a time origin and stick with it, in both your input files and ocean.in.
by kate
Thu Dec 10, 2020 5:08 pm
Forum: ROMS Discussion
Topic: Vertical eddy viscosity
Replies: 4
Views: 313

Re: Vertical eddy viscosity

It is the vertical viscosity used by the model, the sum of all the bits that go into it. The details depend on which vertical mixing scheme you are using.
by kate
Wed Dec 09, 2020 4:49 pm
Forum: ROMS Discussion
Topic: Vertical eddy viscosity
Replies: 4
Views: 313

Re: Vertical eddy viscosity

I believe it is called Akv in the ROMS output.
by kate
Tue Dec 08, 2020 6:06 pm
Forum: ROMS Discussion
Topic: Temperature blowing up near boundary while use nudge
Replies: 4
Views: 452

Re: Temperature blowing up near boundary while use nudge

So this is only in the bottom temperature and salinity? What are the velocities doing? Rather than nudging to cover it up, it would be better to understand where it is coming from.
by kate
Mon Dec 07, 2020 12:08 am
Forum: ROMS Discussion
Topic: Temperature blowing up near boundary while use nudge
Replies: 4
Views: 452

Re: Temperature blowing up near boundary while use nudge

If the boundary is at the edge of the plot, you want to mask out those two purple points. Otherwise you should unmask the point just north of them. At the boundary, the algorithm is using values a few gridpoints in from the boundary so they all have to be wet for it to be well-behaved.
by kate
Wed Dec 02, 2020 6:23 pm
Forum: ROMS Usage
Topic: Differences with CROCO
Replies: 6
Views: 941

Re: Differences with CROCO

I don't know what the capabilities of PISCES is, but my ROMS "master" branch has the COBALT model. The person who put it into ROMS has moved on, taken a different job, and I'm hardly an expert. The files for a Gulf of Alaska domain are in my Apps_master repository (time to rename some of these thing...
by kate
Mon Nov 30, 2020 6:38 pm
Forum: ROMS Discussion
Topic: forrtl: severe (174): SIGSEGV, segmentation fault occurred
Replies: 9
Views: 806

Re: forrtl: severe (174): SIGSEGV, segmentation fault occurred

When you get a seg fault, you can recompile with USE_DEBUG to get line numbers in the output. You then need to look at the relevant lines in *your* .f90 files.
by kate
Mon Nov 23, 2020 7:41 pm
Forum: ROMS Problems
Topic: upwelling case READ_PHYPAR - Variable index not yet loaded, idStr
Replies: 1
Views: 206

Re: upwelling case READ_PHYPAR - Variable index not yet loaded, idStr

You need to use a varinfo.dat file which is consistent with your ROMS version. idStr is in some varinfo.dat files: 'strain' ! Output 'strain' 'second-1' ! [1/s] 'strain, scalar, series' 'ocean_time' 'idStr' 'p3dvar' 1.0d0
by kate
Wed Nov 11, 2020 6:27 pm
Forum: ROMS Discussion
Topic: Compiling error
Replies: 1
Views: 223

Re: Compiling error

'use_mpif90=on' means that the executable name to use is mpif90. However, which: no mpif90 in (/root/.local/bin:/root/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/var/lib/snapd/snap/bin) On the supercomputer I use, openmpi is in a module. Loading the module adds mpif90 to my pat...
by kate
Wed Nov 11, 2020 3:24 am
Forum: ROMS Problems
Topic: ROMS<kshedstrom version> question about the "zeta <free surface elevation>"
Replies: 4
Views: 321

Re: ROMS<kshedstrom version> question about the "zeta <free surface elevation>"

What is it that you hope to gain by switching from POM to ROMS? I haven't modeled such a small lake with ROMS, but I expect it can handle it. It smells like a setup problem or who knows what. Anyway, I am in the process of leaving ROMS for yet another model - MOM6. It too should be able to handle a ...
by kate
Tue Nov 10, 2020 10:56 pm
Forum: ROMS Problems
Topic: ROMS<kshedstrom version> question about the "zeta <free surface elevation>"
Replies: 4
Views: 321

Re: ROMS<kshedstrom version> question about the "zeta <free surface elevation>"

Are the red areas elevated land, as in for wetting and drying? I too would have expected the currents to be more closely related to the ssh contours - it looks like they align right at the end?
by kate
Sat Nov 07, 2020 3:07 am
Forum: ROMS Problems
Topic: ROMS<kshedstrom version>blowing up problems with the "specific heat for seawater" "
Replies: 3
Views: 298

Re: ROMS<kshedstrom version>blowing up problems with the "specific heat for seawater" "

I meant you should #define LIMIT_BSTRESS, though I haven't tried it with LOGDRAG.
by kate
Fri Nov 06, 2020 3:52 pm
Forum: ROMS Problems
Topic: ROMS<kshedstrom version>blowing up problems with the "specific heat for seawater" "
Replies: 3
Views: 298

Re: ROMS<kshedstrom version>blowing up problems with the "specific heat for seawater" "

You have gone for many steps before it blows up. You might as well try a shorter step and see if that helps. You also know that it blew up on the bottom: point 28, 7, 1. What do you have for bottom drag options? Are you limiting the bottom drag?
by kate
Tue Nov 03, 2020 12:54 am
Forum: ROMS Discussion
Topic: Problem with ssh response to uniform winds.
Replies: 4
Views: 382

Re: Problem with ssh response to uniform winds.

If you want a solution that is invariant along-channel, making it periodic that way is certainly the way to go. You can save some cycles by making the channel quite short, say only 3-4 points in that direction. Did you try a wall at the shallow side and open offshore?
by kate
Mon Nov 02, 2020 8:59 pm
Forum: ROMS Discussion
Topic: running error
Replies: 5
Views: 317

Re: running error

You are fine, your builds are all clean. As for how to figure out what's going on, you say you are running in serial mode? Can you use the gdb debugger? I would run roms in a debugger to see exactly what is going on, but I know most don't have access to totalview or ddt (costly debuggers that work i...
by kate
Mon Nov 02, 2020 8:09 pm
Forum: ROMS Discussion
Topic: running error
Replies: 5
Views: 317

Re: running error

I know the default with the build scripts is to always do clean builds - "make clean" before "make". It looks like you didn't change that.
by kate
Mon Nov 02, 2020 7:27 pm
Forum: ROMS Discussion
Topic: running error
Replies: 5
Views: 317

Re: running error

To clarify, you changed to your grid file and it ran. Now you've changed so that ANA_INITIAL is no longer on and that's when it failed?
Did you do a clean build in debug mode? I would have expected more of a backtrace.
by kate
Sat Oct 31, 2020 9:54 pm
Forum: ROMS Discussion
Topic: run error when applying forcing file
Replies: 7
Views: 523

Re: run error when applying forcing file

You can see my various .h and .in files here: https://github.com/kshedstrom/Apps Of course, they go with my ROMS branch, so you have to know what you are doing to pick and choose features. I've described things as best I can here: https://github.com/kshedstrom/roms_manual/blob/master/roms_manual.pdf
by kate
Sat Oct 31, 2020 8:53 pm
Forum: ROMS Discussion
Topic: run error when applying forcing file
Replies: 7
Views: 523

Re: run error when applying forcing file

Maybe? I don't see BULK_FLUXES in your .h file, so without ANA_SMFLUX, the model will only read your sms_era file, not many of the rest of them. There may be other errors in your roms.in file though. The line: allocate ( FRC(max_Ffiles,Ngrids) ) happens once in read_phypar.F, so I assume that's what...
by kate
Sat Oct 31, 2020 8:01 pm
Forum: ROMS Discussion
Topic: run error when applying forcing file
Replies: 7
Views: 523

Re: run error when applying forcing file

I see that you turned off ANA_SMFLUX. Did you turn on anything else to compensate? I use BULK_FLUXES with: # if defined BULK_FLUXES || defined CCSM_FLUXES # define LONGWAVE_OUT # undef DIURNAL_SRFLUX # define SOLAR_SOURCE # define EMINUSP # undef ALBEDO_CLOUD # define ALBEDO_CURVE /* for water */ # ...
by kate
Sat Oct 24, 2020 3:33 pm
Forum: User Applications
Topic: Forcing file
Replies: 1
Views: 401

Re: Forcing file

Once the forcing is on a regular lat, lon grid, ROMS should be able to interpolate correctly internally. Just make sure your winds are east (u) and north (v) and let ROMS rotate to its grid.
by kate
Tue Oct 20, 2020 7:33 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS MERRA-2 Example
Replies: 3
Views: 454

Re: PYROMS MERRA-2 Example

I haven't run the MERRA scripts since we switched to using JRA55-do some time back. The "lat" on line 56 is creating the variable and setting it based on what's in the first file. The use on line 85 is using the existing "lat" to write it into a second file. You can search on the error message and f...
by kate
Wed Oct 07, 2020 5:16 pm
Forum: ROMS Problems
Topic: ROMS<kshedstrom version> blowing-up with ' REASON: KEchar = NaN, PEchar = NaN' report
Replies: 2
Views: 491

Re: ROMS<kshedstrom version> blowing-up with ' REASON: KEchar = NaN, PEchar = NaN' report

That's a very intensified surface layer. Did you try a shorter timestep? There are numerical limitations on the vertical operations which can cause blow-ups like that.
by kate
Wed Sep 30, 2020 6:13 pm
Forum: ROMS Information
Topic: How to link an external library (CFIO) to ROMS
Replies: 1
Views: 369

Re: How to link an external library (CFIO) to ROMS

Simply linking to it may not be enough, but adding a library to link to is done in the Compilers file for your system. Just add to the LIBS macro - there are plenty of examples where it gets added to with a "+=". To get ROMS to call your library is another story unless the library has routines with ...
by kate
Wed Sep 23, 2020 4:25 pm
Forum: ROMS Bugs
Topic: Error with wclock_on subroutine
Replies: 1
Views: 515

Re: Error with wclock_on subroutine

If you recompile with USE_DEBUG you should get some line numbers in the error report. Those line numbers refer to your files in the SCRATCH_DIR, so it helps if you go looking there to find out say which call to wclock_on is failing.
by kate
Tue Sep 22, 2020 9:29 pm
Forum: ROMS Discussion
Topic: mpirun has exited due to process rank 28 with PID 108752 on
Replies: 2
Views: 357

Re: mpirun has exited due to process rank 28 with PID 108752 on

Because each of us runs with different cpp options, we get differing line numbers. You should look at your lines here: coawstG 00000000009A1EFC load_s2d_ 1722 inp_par.f90 coawstG 0000000000A83CE5 read_phypar_ 1456 read_phypar.f90 However, I suspect that you might be using a roms.in file which is inc...
by kate
Fri Aug 21, 2020 5:08 pm
Forum: ROMS Discussion
Topic: ROMS-CICE test case
Replies: 2
Views: 529

Re: ROMS-CICE test case

Some config files are in my https://github.com/kshedstrom/Apps repository, though most were run with just ROMS ice. In terms of all the big netcdf files you'd need, those are not available online anywhere. Maybe you should ask the Norwegians who are actually running that coupled code.
by kate
Fri Aug 14, 2020 6:01 am
Forum: ROMS Discussion
Topic: Input Initial Conditions from External sources
Replies: 2
Views: 503

Re: Input Initial Conditions from External sources

There are tools in various languages for creating these files from various input sources. Sometimes we have to get our hands dirty and create the exact tool we need from what's out there. I'm working with the pyroms package which can read the old HYCOM files but not the new HYCOM files. Also, the py...
by kate
Wed Aug 12, 2020 6:30 am
Forum: ROMS Discussion
Topic: NBCFILES' unread or misssing from input script 'roms.in'
Replies: 2
Views: 431

Re: NBCFILES' unread or misssing from input script 'roms.in'

You must use roms.in files which are compatible with whatever version of ROMS you are using. I added the ability to use multiple BRY, etc. files to my branch of ROMS. If you are using my branch, you must use an ocean.in that provides everything my branch wants.
by kate
Wed Aug 05, 2020 10:22 pm
Forum: ROMS Discussion
Topic: River Sinks
Replies: 7
Views: 1022

Re: River Sinks

Do you only have river sinks or are there sources as well? You always want the tracer advection to be "upwind", using the tracer value of the upwind point. For sources, that means specifying the river temperature. For sinks, that means using the ambient temperature.
by kate
Wed Jul 29, 2020 6:16 pm
Forum: ROMS Discussion
Topic: bottom freshwater flux
Replies: 2
Views: 427

Re: bottom freshwater flux

In the model, the variable is called "btflx". Most of us set it to zero and #define ANA_BTFLUX and #define ANA_BSFLUX so that no files are used for the zeros.
by kate
Thu Jul 16, 2020 6:34 am
Forum: ROMS Discussion
Topic: River input file in ROMS.
Replies: 3
Views: 654

Re: River input file in ROMS.

I never have enough data to make use of the full power of the model to import depth-dependent information. To make smaller river files, I hacked ROMS to read in river tracers either as a function of time only or as a function of time and river, but not depth.
by kate
Wed Jul 15, 2020 4:04 pm
Forum: ROMS Discussion
Topic: River input file in ROMS.
Replies: 3
Views: 654

Re: River input file in ROMS.

Yes, this is a problem. For salinity, we just use 0.0. For the temperature, some climate modelers use the air temperature. If you were already reading air temperature for the bulk fluxes, you could modify the code to just use it. For Gulf of Alaska waters, Seth Danielson gave me a year of daily temp...
by kate
Mon Jul 13, 2020 6:48 pm
Forum: ROMS Discussion
Topic: error compiling ana_lrflux
Replies: 2
Views: 373

Re: error compiling ana_lrflux

In my code, I added the option to have an analytic longwave flux for the ICE_BOX test case. See if you can copy what I did into COAWST to make it work if that's what you need. The sea ice needs the separate longwave and shortwave fluxes to compute its own bulk fluxes (if you are using sea ice).
by kate
Thu Jul 09, 2020 4:29 pm
Forum: ROMS Discussion
Topic: Budgell vs. CICE sea ice model with ROMS
Replies: 7
Views: 1231

Re: Budgell vs. CICE sea ice model with ROMS

There's a CICE forum at https://bb.cgd.ucar.edu/cesm/forums/cic ... rtium.146/ which would be another place to ask this sort of question.
by kate
Wed Jul 08, 2020 5:29 pm
Forum: ROMS Discussion
Topic: Running error: unable to open existing NetCDF file
Replies: 6
Views: 884

Re: Running error: unable to open existing NetCDF file

There is a dying art known as "debugging". The lucky few have access to interactive debuggers such as totalview, ddt, ddd. Many have access to gdb, without a gui. If you compile ROMS with USE_DEBUG then you can run it within one of these debuggers and watch what's going on line by line. The fancy de...
by kate
Thu Jul 02, 2020 5:44 pm
Forum: ROMS Discussion
Topic: Running error: unable to open existing NetCDF file
Replies: 6
Views: 884

Re: Running error: unable to open existing NetCDF file

Weird, you've got ANA_INITIAL so it shouldn't have a reason to read a file instead.
by kate
Wed Jul 01, 2020 8:02 pm
Forum: ROMS Discussion
Topic: Roms with ice module in Coawst
Replies: 13
Views: 1329

Re: Roms with ice module in Coawst

FILTERED is the five-day running mean tidal filter on AVERAGE output, one way of detiding things.
by kate
Wed Jul 01, 2020 6:45 am
Forum: ROMS Discussion
Topic: Add the path of a new Fortran90 file
Replies: 5
Views: 613

Re: Add the path of a new Fortran90 file

If I were to overhaul this now, I would check to see if the compilers now in use can all invoke the C-preprocessor automatically for .F90 files. If so, take out the cpp step and the cpp_clean step and just have compile-time errors point to the line number of the .F90 file. I did once meet a compiler...
by kate
Wed Jul 01, 2020 5:55 am
Forum: ROMS Discussion
Topic: Add the path of a new Fortran90 file
Replies: 5
Views: 613

Re: Add the path of a new Fortran90 file

I don't think so. It instead uses:

Code: Select all

define compile-rules
  $(foreach f, $(local_src),       \
    $(call one-compile-rule,$(call source-to-object,$f), \
    $(call f90-source,$f),$f))
endef
and friends. The compile-rules function is called in the Module.mk files.
by kate
Wed Jul 01, 2020 5:01 am
Forum: ROMS Discussion
Topic: Add the path of a new Fortran90 file
Replies: 5
Views: 613

Re: Add the path of a new Fortran90 file

The makefile has rules for going from .F files to .f90 files to compile them into .o files. You would have to add rules for .F90 files if you refuse to rename your file.
by kate
Wed Jun 24, 2020 5:44 pm
Forum: ROMS Discussion
Topic: Budgell vs. CICE sea ice model with ROMS
Replies: 7
Views: 1231

Re: Budgell vs. CICE sea ice model with ROMS

CICE is a much fancier model, with compile time choices of the number of ice categories, choices for melt ponds, you name it. I don't know your domain well enough to know if you need the extra features. For our needs, we ran the ROMS ice because it's cheaper. Now we're moving away from that for the ...
by kate
Wed Jun 24, 2020 5:32 pm
Forum: ROMS Problems
Topic: mp_collect bug when compiling with define PARALLEL_IO
Replies: 18
Views: 5784

Re: mp_collect bug when compiling with define PARALLEL_IO

I suggest learning about git. You can have a git repo with (a) the old ROMS you started from (b) the old ROMS with your modifications and (c) new ROMS. Or the crude thing to do is diff old ROMS vs. your code and save those diffs to a file. Then apply those changes to the new ROMS by hand. I assume t...
by kate
Tue Jun 23, 2020 5:25 am
Forum: ROMS Discussion
Topic: Roms with ice module in Coawst
Replies: 13
Views: 1329

Re: Roms with ice module in Coawst

What is it trying to read at line 222 of your get_data.f90?
by kate
Sat Jun 20, 2020 11:44 pm
Forum: ROMS Tools and Techniques
Topic: pyroms with mercator dataset
Replies: 4
Views: 864

Re: pyroms with mercator dataset

Here is the ncdump of the thing: chinook02.rcs.alaska.edu 207% ncdump -h GL2V1_mesh_mask_new.nc netcdf GL2V1_mesh_mask_new { dimensions: y = 1021 ; x = 1442 ; z = 75 ; variables: int y(y) ; y:axis = "Y" ; y:standard_name = "projection_y_coordinate" ; y:units = " " ; int x(x) ; x:axis = "X" ; x:stand...
by kate
Sat Jun 20, 2020 6:25 am
Forum: ROMS Discussion
Topic: Roms with ice module in Coawst
Replies: 13
Views: 1329

Re: Roms with ice module in Coawst

You can attack this with print statements or by doing the big diff between the two codes. At least you know it's going bad very, very quickly. That sort of problem is easier to debug than those which go bad after days or years of run time.
by kate
Fri Jun 19, 2020 5:28 pm
Forum: ROMS Discussion
Topic: Roms with ice module in Coawst
Replies: 13
Views: 1329

Re: Roms with ice module in Coawst

trouble in ice_frazil 40 31 9 21.255285453705273 -3705.3862141090412
If you check this print statement, you'll see that it is telling you the i,j,k location of where it thinks there is trouble. Next are the T and S values there. That S value is not good.
by kate
Fri Jun 19, 2020 5:23 pm
Forum: ROMS Installation
Topic: Shared libraries error: libnetcdff.so.5
Replies: 2
Views: 502

Re: Shared libraries error: libnetcdff.so.5

Is it possible for you to recompile the program so that it uses the new library instead of the old library? It is assumed that you will be compiling this code to avoid these problems.
by kate
Fri Jun 19, 2020 5:21 pm
Forum: ROMS Tools and Techniques
Topic: pyroms with mercator dataset
Replies: 4
Views: 864

Re: pyroms with mercator dataset

The thing is half a GB in size. Can you ask for it from your source of GLORYS files?
by kate
Tue Jun 16, 2020 5:12 pm
Forum: ROMS Discussion
Topic: Boundary condition problems
Replies: 4
Views: 685

Re: Boundary condition problems

I'm not sure about the boundary conditions, but for the initial conditions, the ubar/vbar fields are set in ini_fields.F.
by kate
Fri Jun 05, 2020 4:20 pm
Forum: ROMS Discussion
Topic: the defination of density rho
Replies: 16
Views: 2181

Re: the defination of density rho

OK, but I'm surprised.
by kate
Fri Jun 05, 2020 4:15 am
Forum: ROMS Discussion
Topic: the defination of density rho
Replies: 16
Views: 2181

Re: the defination of density rho

The variable "rho" in the model is insitu density anomaly. You add rho0 to it to get full density. Note that rho0 may or may not be 1000.
by kate
Sat May 16, 2020 7:35 pm
Forum: ROMS Discussion
Topic: Blow up problem
Replies: 3
Views: 800

Re: Blow up problem

For once I can't say try a shorter timestep because you are already in trouble before you've taken any timesteps: 0 0001-01-01 00:00:00.00 NaN NaN NaN NaN (000,000,0) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 I'm not sure what has gone wrong.
by kate
Thu Apr 30, 2020 4:47 pm
Forum: ROMS Discussion
Topic: Could I repeat the simulation forced by bulk flux by using its output?
Replies: 3
Views: 704

Re: Could I repeat the simulation forced by bulk flux by using its output?

Also, it depends on the frequency of the forcing. I'm using three-hourly forcing files and if I were to apply them as daily averages from the ROMS output, that would be different.
by kate
Fri Apr 24, 2020 4:41 pm
Forum: ROMS Tools and Techniques
Topic: pyroms error during initial condition remapping process
Replies: 2
Views: 674

Re: pyroms error during initial condition remapping process

I think perhaps you need to change this in make_ic_file.py: xrange=(400, 500), yrange=(180, 280) To save time, it doesn't search the whole world, but just the part of the world where your grid is (or in this case where the CCS1 grid is). In my case, num_links is rather large: netcdf remap_weights_SO...
by kate
Sat Apr 18, 2020 3:08 pm
Forum: ROMS Discussion
Topic: grid generation for idealized geometry
Replies: 4
Views: 867

Re: grid generation for idealized geometry

I just want to point out that your two western corners are not right angles and will lead to weirdness for an orthogonal grid.
by kate
Thu Apr 16, 2020 3:17 am
Forum: ROMS Discussion
Topic: Test of open boundary conditions
Replies: 12
Views: 1660

Re: Test of open boundary conditions

Flather does allow introduction of outside information.
by kate
Wed Apr 15, 2020 5:50 pm
Forum: ROMS Discussion
Topic: Writing out wT (time average)
Replies: 2
Views: 384

Re: Writing out wT (time average)

One usually simply does a simple average between T points to get values at w points. For the top and bottom, just use the nearest value.
by kate
Mon Apr 13, 2020 4:34 pm
Forum: ROMS Discussion
Topic: Test of open boundary conditions
Replies: 12
Views: 1660

Re: Test of open boundary conditions

So should I keep Chapman+Flather or Chapman+Shc for 2D variables and use radiation+nudging for 3D u,v, t, and s?
Yes, exactly.
by kate
Mon Apr 13, 2020 12:54 am
Forum: ROMS Discussion
Topic: Test of open boundary conditions
Replies: 12
Views: 1660

Re: Test of open boundary conditions

Because you don't have stratification, you are mostly testing the barotropic OBCs. You should try both the Chapman/Flather combo and the Sch option (described in a Mason et al., 2010 paper). Flather is known to work pretty well for the barotropic mode. What's more challenging are the baroclinic sign...
by kate
Sat Apr 04, 2020 8:17 pm
Forum: ROMS Problems
Topic: An error for two-way nesting with restart
Replies: 2
Views: 538

Re: An error for two-way nesting with restart

Two-way nesting is slow in ROMS. That's the way it is. The French croco code may be faster.
by kate
Fri Apr 03, 2020 4:56 pm
Forum: ROMS Discussion
Topic: nested grid with pyroms: tutorial and python routines
Replies: 7
Views: 1339

Re: nested grid with pyroms: tutorial and python routines

The question is what are your needs? Then you can ask how to do that in Python. There is a Python community for geoscientists called Pangeo and they have other resources. My use of pyroms can be seen through the examples. This is not how things would be done today by the Pangeo community because new...
by kate
Mon Mar 30, 2020 4:20 pm
Forum: ROMS Discussion
Topic: Question about radiation boundary condition
Replies: 7
Views: 1032

Re: Question about radiation boundary condition

How will you know until you try it?
by kate
Sun Mar 29, 2020 10:08 pm
Forum: ROMS Discussion
Topic: Question about radiation boundary condition
Replies: 7
Views: 1032

Re: Question about radiation boundary condition

I guess it depends on how well you trust your climatology vs your boundary fields. Things have come a long way since that paper was written in terms of the global model reanalysis products that are available as boundary conditions.
by kate
Wed Mar 25, 2020 9:52 pm
Forum: ROMS Tools and Techniques
Topic: How to cite pyroms ?
Replies: 4
Views: 754

Re: How to cite pyroms ?

It took a little fussing, but here you go: Image
by kate
Wed Mar 25, 2020 4:14 pm
Forum: ROMS Tools and Techniques
Topic: How to cite pyroms ?
Replies: 4
Views: 754

Re: How to cite pyroms ?

Do you need a DOI? Is that what you're asking for? I can generate one from the current version.
by kate
Tue Mar 17, 2020 9:22 pm
Forum: ROMS Problems
Topic: Blowing up problem
Replies: 9
Views: 1673

Re: Blowing up problem

I assume you are showing us bottom density, which really only tells something about the bottom depth. It looks like you still have lakes. Note that the blue square here is a lake, not a bay (red being land):
bitmap2.png
bitmap2.png (1.76 KiB) Viewed 1206 times
by kate
Tue Mar 17, 2020 5:24 pm
Forum: ROMS Problems
Topic: Blowing up problem
Replies: 9
Views: 1673

Re: Blowing up problem

I can't tell from that picture if the mask is a problem or not. You need to zoom in and edit the mask so that there are no enclosed lakes. Also, near the boundary, you want no bits of land leaving just one grid of water at the boundary, like: bitmap.png The edge of this grid is the edge of the compu...
by kate
Tue Mar 17, 2020 5:09 am
Forum: ROMS Problems
Topic: Blowing up problem
Replies: 9
Views: 1673

Re: Blowing up problem

That Courant number could still be high. Does the model run longer with a shorter step?

Did you smooth the bathymetry? Did you look at where the blow-up happens?
by kate
Wed Mar 04, 2020 4:56 pm
Forum: ROMS Discussion
Topic: compiling ROMS with DIAGNOSTICS_TS
Replies: 3
Views: 613

Re: compiling ROMS with DIAGNOSTICS_TS

This is a matter of personal preference. I have a friend here to swears by never updating his code such that if he now had to update, it would be a matter of months of work to merge his stuff with modern ROMS and he's hoping to retire before it comes to that. I go the other way, trying to merge earl...
by kate
Mon Mar 02, 2020 4:14 pm
Forum: ROMS Discussion
Topic: Inundation Blowing Up
Replies: 2
Views: 558

Re: Inundation Blowing Up

Have you tried a shorter timestep? It is blowing up so fast that is likely to be a contributing factor. After that, you need to be looking at how and where it is going bad. Wetting and drying problems can be tricky. I have gone so far as to manually edit the bathymetry in the wet-dry places to make ...
by kate
Wed Feb 26, 2020 6:19 pm
Forum: ROMS Discussion
Topic: Problem with surface heat flux
Replies: 13
Views: 1888

Re: Problem with surface heat flux

CCSM_FLUXES want kg/kg for specific humidity while BULK_FLUXES has some weird check to see if you are providing relative humidity vs. specific humidity. ROMS wants Celsius for air temp. I hacked the reading of it to accept Kelvin if that's what you have. ROMS wants millibars, hence the 0.01 factor i...
by kate
Wed Feb 26, 2020 5:22 pm
Forum: ROMS Discussion
Topic: Problem with surface heat flux
Replies: 13
Views: 1888

Re: Problem with surface heat flux

I am using surface forcing from JRA55-do, which has units of Pascals for pressure, kg/kg for specific humidity, and Kelvin for air temperature. That's how my setup works.
by kate
Tue Feb 25, 2020 6:45 pm
Forum: ROMS Problems
Topic: Error when compiling 2D test cases
Replies: 2
Views: 561

Re: Error when compiling 2D test cases

You should be able to surround all instances of T_ADV with #ifdef SOLVE3D. There is no tracer advection in the shallow water equations.
by kate
Tue Feb 25, 2020 6:23 pm
Forum: ROMS Discussion
Topic: Problem with surface heat flux
Replies: 13
Views: 1888

Re: Problem with surface heat flux

Simply turn off BULK_FLUXES. However, if you have ice where the driving model does not, things could get weird.
by kate
Wed Feb 19, 2020 6:00 am
Forum: ROMS Discussion
Topic: Problem with surface heat flux
Replies: 13
Views: 1888

Re: Problem with surface heat flux

What are you doing for the ice albedo? I have: # define ALBEDO_CURVE /* for water */ # undef ICE_ALB_EC92 /* for ice */ # define ALBEDO_CSIM /* for ice */ # undef ALBEDO_FILE /* for both */ Someone pointed out that ICE_ALB_EC92 can be unstable, so maybe I should delete it. ETA: Oh, I see you are usi...
by kate
Wed Feb 19, 2020 5:56 am
Forum: ROMS Discussion
Topic: Problem with surface heat flux
Replies: 13
Views: 1888

Re: Problem with surface heat flux

You could always try turning off CCSM_FLUXES to get the other flux formulation to see if it's any better. Beyond that, I'd worry about units of all the inputs.
by kate
Wed Jan 29, 2020 11:55 pm
Forum: ROMS Problems
Topic: A ROMS nesting question
Replies: 10
Views: 1896

Re: A ROMS nesting question

The climatology files are meant to be provided on your grid. Yes, if you need climatology nudging, you will need to make files for all grids that need it, though one might hope that online nesting will mean the children don't need climatology nudging.
by kate
Wed Jan 15, 2020 6:41 pm
Forum: ROMS Discussion
Topic: Wave absorption boundary condition
Replies: 6
Views: 1086

Re: Wave absorption boundary condition

If your barotropic solution is working for you, don't change that.

Using Rad alone will not require external fields. You may or may not do better with RADIATION_2D. Just try it and see if it helps.

If you want more specific advice, you should probably tell us more about what you are trying to do.
by kate
Wed Jan 15, 2020 5:01 pm
Forum: ROMS Discussion
Topic: some problems about bulk_fluxes
Replies: 16
Views: 1731

Re: some problems about bulk_fluxes

No, the .h files are included in the fortran files, compiled by the Fortran compiler. Linking to C code is a different beast. Not impossible, just different.
by kate
Wed Jan 15, 2020 4:50 pm
Forum: ROMS Bugs
Topic: WET_DRY and Gridbuilder Problem
Replies: 21
Views: 3149

Re: WET_DRY and Gridbuilder Problem

You mean it died inside the land mask? Do you know where in the code this happened?

Dcrit is the depth of the water in the "dry" places. You can never remove that last bit of water from a cell. That's what keeps things stable (if all is working right).
by kate
Tue Jan 14, 2020 7:40 pm
Forum: ROMS Discussion
Topic: Wave absorption boundary condition
Replies: 6
Views: 1086

Re: Wave absorption boundary condition

This is a classic problem with regional ocean models. You can try the radiation boundary conditions instead, if you like. However, what is working for us in practice is using values from a global ocean model on the boundaries to keep things behaving better (if not entirely well).
by kate
Tue Jan 14, 2020 7:36 pm
Forum: ROMS Discussion
Topic: some problems about bulk_fluxes
Replies: 16
Views: 1731

Re: some problems about bulk_fluxes

The answer might be in ROMS/Functionals//ana_smflux.h, which is the file with the analytic surface stresses for some of the standard test problems. Some of them have things which vary in time as examples. If you want to change things by the timestep number, that is iic in the code.
by kate
Mon Jan 13, 2020 10:43 pm
Forum: ROMS Ice
Topic: extreme salinity values if ICE MODEL used in COAWST branch
Replies: 2
Views: 955

Re: extreme salinity values if ICE MODEL used in COAWST branch

I don't know where those salinities come from. The ice_frazil routine makes the coupled ice-ocean model more sensitive to some issues than the stand-alone ocean model. I know that ice_frazil complains when the level thickness Hz goes negative too. What I would do is add a check to diag.F to complain...
by kate
Mon Jan 13, 2020 1:59 pm
Forum: ROMS Discussion
Topic: some problems about bulk_fluxes
Replies: 16
Views: 1731

Re: some problems about bulk_fluxes

It is called every time step. Changes in time happen in your code.
by kate
Mon Jan 13, 2020 11:17 am
Forum: ROMS Discussion
Topic: some problems about bulk_fluxes
Replies: 16
Views: 1731

Re: some problems about bulk_fluxes

It sounded like you wanted a linear progression in time for your stress. If so, ROMS does a linear interpolation in time from one record to the next and could be doing that work for you. Or for a very simple pattern you can use the ana_xxx.h code which gets called every timestep.
by kate
Mon Jan 13, 2020 1:05 am
Forum: ROMS Discussion
Topic: some problems about bulk_fluxes
Replies: 16
Views: 1731

Re: some problems about bulk_fluxes

GET_2DFLD - surface u-momentum stress, t = 0 02:24:00 (Rec=0000002, Index=2, File: palu_stress.nc) (Tmin= 0.0000 Tmax= 3001.0000) (Min = 0.00000000E+00 Max = 0.00000000E+00) It is reading the second record here. What are the units on time in the stress file? What units do you think it should be? Pe...
by kate
Thu Jan 09, 2020 11:48 am
Forum: ROMS Bugs
Topic: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Replies: 12
Views: 1467

Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281

With WET_DRY, zeta in the land area is always some DCRIT above the land, nothing to worry about. It should not be changing in time.

Have you tried fussing with your timestep?
by kate
Thu Jan 09, 2020 3:04 am
Forum: ROMS Bugs
Topic: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Replies: 12
Views: 1467

Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281

The Flather condition requires values from outside the domain for both ubar/vbar and zeta. You can provide them via ANA_xxx options or you can specify a file: NBCFILES == 1 ! number of boundary files BRYNAME == /center1/AKWATERS/kshedstrom/Arctic/Files/Arctic4_bdry_1980_SODA3.3.1.nc | /center1/AKWAT...
by kate
Thu Jan 09, 2020 2:47 am
Forum: ROMS Bugs
Topic: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Replies: 12
Views: 1467

Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281

Since you are using COAWST, you need to set the number of boundary and climatology files, not just the number of forcing files. See the coawst roms.in for an example. Same for the river inputs.
by kate
Wed Jan 08, 2020 12:53 pm
Forum: ROMS Bugs
Topic: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Replies: 12
Views: 1467

Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281

So you think you have already remapped the thing to your grid. You have to have the dimensions *exactly* right or else ROMS will decide to remap for you. Some words on grid sizes. If you want powers of two or three in the grid for your decompositions without needing leftovers, you need Lm, Mm to hav...
by kate
Wed Jan 08, 2020 11:40 am
Forum: ROMS Bugs
Topic: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Replies: 12
Views: 1467

Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281

at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:104 I suspect you are reading an atmospheric field here? The model wouldn't try to interpolate from your grid file, assuming all those fields are already on your grid. So, what atmospheric fields are you reading and are you expecting ROMS to do the...
by kate
Mon Jan 06, 2020 8:43 pm
Forum: ROMS Discussion
Topic: Fail to active T and S boundary conditons
Replies: 2
Views: 536

Re: Fail to active T and S boundary conditons

You only have two tracers, right? LBC(isTvar) == Clo Clo RadNud Clo \ ! temperature Clo Clo RadNud Clo \ ! salinity Clo Clo Clo Clo \ ! inert(1) Clo Clo Clo Clo \ ! inert(2) Clo Clo Clo Clo \ ! inert(3) Clo Clo Clo Clo \ ! inert(4) Clo Clo Clo Clo \ ! inert(5) Clo Clo Clo Clo ! inert(6) You need to ...
by kate
Sun Jan 05, 2020 1:35 pm
Forum: ROMS Source
Topic: how to choose cpp?
Replies: 5
Views: 922

Re: how to choose cpp?

What advection scheme are you using? There's this block in globaldefs.h: #ifdef UV_U3ADV_SPLIT # define VISC_3DCOEF # ifndef UV_C4ADVECTION # define UV_C4ADVECTION # endif # ifndef UV_VIS4 # define UV_VIS4 # endif # ifdef UV_VIS2 # undef UV_VIS2 # endif # ifdef UV_SMAGORINSKY # undef UV_SMAGORINSKY ...
by kate
Sat Jan 04, 2020 7:17 pm
Forum: ROMS Source
Topic: how to choose cpp?
Replies: 5
Views: 922

Re: how to choose cpp?

The code I have calls it UV_VIS2, not VISC2. Asking for a flag ROMS doesn't know about is like asking for #define MARS, not going to do anything. Which ROMS are you using?
by kate
Sat Jan 04, 2020 7:07 pm
Forum: ROMS Bugs
Topic: Segmentation fault
Replies: 9
Views: 1220

Re: Segmentation fault

If you log into: https://www.myroms.org/projects/src you can submit a bug report.
by kate
Sat Jan 04, 2020 7:38 am
Forum: ROMS Discussion
Topic: some problems about bulk_fluxes
Replies: 16
Views: 1731

Re: some problems about bulk_fluxes

1. Undef BULK_FLUXES in .h file. and define ANA_STFLUX, ANA_SSFLUX, ANA_BSFLUX and ANA_BTFLUX (totally 4 ?) in .h file (I wonder how ana function set these 4 variables, just as zero automatically?). 2. Set NFFFILE=1 and input test_stress .nc file Yes, exactly. The way to tell what the default stres...
by kate
Sat Jan 04, 2020 7:34 am
Forum: ROMS Bugs
Topic: Segmentation fault
Replies: 9
Views: 1220

Re: Segmentation fault

Well, you could try submitting a bug report for this case of no forcing files, but most of us just set NFFILES to 1 and move on.
by kate
Fri Jan 03, 2020 5:37 pm
Forum: ROMS Discussion
Topic: some problems about bulk_fluxes
Replies: 16
Views: 1731

Re: some problems about bulk_fluxes

I wonder that if I just input sustrs and svstrs, must I set ANA_STFLUX, ANA_SSFLUX, ANA_BSFLUX and ANA_BTFLUX?
Yes, if you provide the stresses, you need to provide the tracer fluxes too, and turn off BULK_FLUXES.
by kate
Fri Jan 03, 2020 5:32 pm
Forum: ROMS Bugs
Topic: Segmentation fault
Replies: 9
Views: 1220

Re: Segmentation fault

What is on this line?
at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:101
Oops, didn't see this:
Fortran runtime error: Index '1' of dimension 2 of array 'frc' above upper bound of 0
Did you tell it NFFILES before the list of forcing files?
by kate
Fri Jan 03, 2020 5:09 am
Forum: ROMS Discussion
Topic: d_ecmwf2roms.m
Replies: 9
Views: 2080

Re: d_ecmwf2roms.m

Then you have to either use NCO to merge the files, use something else to merge the files, or hack that ROMS to be able to read multiple files.
by kate
Thu Jan 02, 2020 5:03 pm
Forum: ROMS Discussion
Topic: d_ecmwf2roms.m
Replies: 9
Views: 2080

Re: d_ecmwf2roms.m

You have two choices. One is to leave the files as they are and just list them all in the FRCNAME. Here is an example: NFFILES == 14 ! number of forcing files FRCNAME == /center1/AKWATERS/kate/JRA55-do/JRA55DO_1.4_Pair_2015.nc | /center1/AKWATERS/kate/JRA55-do/JRA55DO_1.4_Pair_2016.nc \ /center1/AKW...
by kate
Wed Jan 01, 2020 5:42 pm
Forum: ROMS Discussion
Topic: SET_CYCLE ERROR: non-cycling regime, but model time exeeds time of the last available data record for
Replies: 1
Views: 358

Re: SET_CYCLE ERROR: non-cycling regime, but model time exeeds time of the last available data record for

The model checks the current model time and the time of all the forcing records. It has found an inconsistency and reported it. Check all your times.
by kate
Thu Dec 26, 2019 6:06 pm
Forum: ROMS Discussion
Topic: Best practices for ROMS bathymetric values
Replies: 10
Views: 1567

Re: Best practices for ROMS bathymetric values

I don't think those points that are MASKED should cause trouble. I know wet-dry points are very problematic in say coastal fjords with one point at -1 m and the next point at 30 m depth. This is in the same domain with tidal flats with a much more gradual bottom depth. It's all about resolving the f...
by kate
Tue Dec 24, 2019 8:42 pm
Forum: ROMS Discussion
Topic: Best practices for ROMS bathymetric values
Replies: 10
Views: 1567

Re: Best practices for ROMS bathymetric values

Really, you should be looking at the fields as it goes bad. Now that you know it goes bad quickly, I would set ninfo to one to catch it going bad ASAP.

Also, did you mean to have ANA_SSFLUX? Bulk_flux should provide salt fluxes based on precipitation.
by kate
Tue Dec 24, 2019 8:15 am
Forum: ROMS Discussion
Topic: Best practices for ROMS bathymetric values
Replies: 10
Views: 1567

Re: Best practices for ROMS bathymetric values

Are all those depths of -100 always dry? Fully masked? Then they shouldn't be causing any trouble.

If it all behaved before you used ALBEDO_CLOUD, maybe you should check that out. I tend not to use it and it might have bugs. Have you looked at your surface fluxes, before and after the change?
by kate
Tue Dec 24, 2019 8:09 am
Forum: ROMS Discussion
Topic: calculating ubar and vbar
Replies: 10
Views: 2678

Re: calculating ubar and vbar

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME C => (i,j,k) Cu Cv Cw Max Speed 144 1 00:00:00 NaN NaN NaN 5.758336E+16 (163,187,19) 2.195484E-02 1.405542E-02 0.000000E+00 1.002493E+00 Have you tried a shorter timestep? Have you tried smoothing your bathymetry? You need to get ROMS ...
by kate
Wed Dec 18, 2019 6:43 pm
Forum: ROMS Bugs
Topic: Error in ROMS execution
Replies: 12
Views: 1572

Re: Error in ROMS execution

Code: Select all

GET_GRID - unable to open grid NetCDF file: wtype_grid
I am adding this field to the grid file and it is saying you are missing the field. You can turn off having a spatially variable water type or you can provide the field. I set it as a function of depth using the attached script.
by kate
Tue Dec 17, 2019 6:52 pm
Forum: ROMS Bugs
Topic: Error in ROMS execution
Replies: 12
Views: 1572

Re: Error in ROMS execution

READ_IcePar - variable info not yet loaded, idUiceE 'uice_eastward' ! Output 'ice eastward momentum component at RHO-points' 'meter second-1' ! [m/s] 'uice_eastward, scalar, series' 'ocean_time' 'idUiceE' 'r2dvar' 1.0d0 'vice_northward' ! Output 'ice northward momentum component at RHO-points' 'met...
by kate
Tue Dec 17, 2019 5:46 am
Forum: ROMS Discussion
Topic: Request: Numerical modeling course syllabi
Replies: 11
Views: 2449

Re: Request: Numerical modeling course syllabi

Harper and I were talking today about how he learned Matlab in grad school. I said I did not, but I did have Bob Parker's
code for creating contour plots back in the 1980's.
by kate
Mon Dec 16, 2019 7:36 pm
Forum: ROMS Bugs
Topic: Error in ROMS execution
Replies: 12
Views: 1572

Re: Error in ROMS execution

You need to set NSSFFILES in your roms.in file, such as: NSSFFILES == 2 ! number of boundary files SSFNAME == /center1/AKWATERS/kshedstrom/Chukchi/Files/JRA_Chukchi_rivers_1993.nc | /center1/AKWATERS/kshedstrom/Chukchi/Files/JRA_Chukchi_rivers_1994.nc | /center1/AKWATERS/kshedstrom/Chukchi/Files/JRA...
by kate
Mon Dec 16, 2019 7:25 pm
Forum: ROMS Bugs
Topic: Error in ROMS execution
Replies: 12
Views: 1572

Re: Error in ROMS execution

Probaby, if you tell me what is on this line:

Code: Select all

 1914  read_phypar.f90
Note that this file is in the build directory.
by kate
Mon Dec 16, 2019 6:09 pm
Forum: ROMS Bugs
Topic: Error in ROMS execution
Replies: 12
Views: 1572

Re: Error in ROMS execution

Oh, yeah, sorry. There needs to be a chunk in varinfo.dat with: 'Uwind_eastward' ! Output 'eastward surface wind component at RHO-points' 'meter second-1' ! [m/s] 'Uwind_eastward, scalar, series' 'ocean_time' 'idUairE' 'r2dvar' 1.0d0 'Vwind_northward' ! Output 'northward surface wind component at RH...
by kate
Mon Dec 09, 2019 3:18 am
Forum: ROMS Discussion
Topic: Very early blowing up problem
Replies: 13
Views: 4360

Re: Very early blowing up problem

When doing wetting and drying, it is possible to have huge rx0 numbers in the wet-dry zone and still be OK. It is also possible to be in trouble there, so you need to really know what you are doing. I know that a "dry" cell 1 m above sea level next to a fjord can spell trouble. There might still be ...
by kate
Fri Dec 06, 2019 7:19 pm
Forum: ROMS Discussion
Topic: nested grid with pyroms: tutorial and python routines
Replies: 7
Views: 1339

Re: nested grid with pyroms: tutorial and python routines

I use the best bathymetry I can on each grid, smoothed appropriately. I don't worry about matching them, though some do (like the Mason et al. paper). When interpolating fields, pyroms puts the parent fields onto a z-grid, then remaps laterally, then puts them onto the child grid's vertical grid. Th...
by kate
Fri Dec 06, 2019 7:11 pm
Forum: ROMS Discussion
Topic: Parallel IO not scaling over Serial IO for ROMS
Replies: 14
Views: 1774

Re: Parallel IO not scaling over Serial IO for ROMS

Where did you get your ROMS code from? That's what I was asking. Hernan's code has PARALLEL_IO, but you say turning that on isn't doing anything. My code does not recognize PARALLEL_IO but instead recognizes PARALLEL_IN and PARALLEL_OUT. A quick check of cppdefs.h should tell you what the options ar...
by kate
Fri Dec 06, 2019 6:06 pm
Forum: ROMS Discussion
Topic: Parallel IO not scaling over Serial IO for ROMS
Replies: 14
Views: 1774

Re: Parallel IO not scaling over Serial IO for ROMS

As for parallel I/O, what code are you using? I tried splitting PARALLEL_IO into PARALLEL_IN and PARALLEL_OUT so that I wouldn't have to convert all my inputs to HDF5, yet still get the benefits of parallel output. It's been a while since I checked, but it didn't speed things up. I heard that at th...
by kate
Fri Dec 06, 2019 6:02 pm
Forum: ROMS Discussion
Topic: nested grid with pyroms: tutorial and python routines
Replies: 7
Views: 1339

Re: nested grid with pyroms: tutorial and python routines

Is this online or offline nesting? If online, the nest needs to be a very specific subgrid of the parent, at an odd ratio such a 1:3, 1:5, etc. If offline, just make another grid as you normally would. It can be at any orientation as long as it is entirely within the parent domain. If online, I thin...
by kate
Thu Dec 05, 2019 4:46 pm
Forum: ROMS Discussion
Topic: Parallel IO not scaling over Serial IO for ROMS
Replies: 14
Views: 1774

Re: Parallel IO not scaling over Serial IO for ROMS

Many years ago our supercomputer center bought a Cray and an IBM. The timings I did without I/O led me to believe that the Cray was the machine for me. However, my real jobs have I/O and that meant that the IBM was actually faster for my problems. At that time, I didn't have a vectorized netCDF libr...
by kate
Thu Dec 05, 2019 3:08 am
Forum: ROMS Discussion
Topic: doubts about BULK_FLUX
Replies: 3
Views: 715

Re: doubts about BULK_FLUX

If you want bulk fluxes, you don't want any of those ANA_SXFLUX options on.

I believe you need air pressure. I don't recall what all it is used for, but I think it goes into converting relative humidity to specific humidity. In a pinch, you can try to set it to 1000 everywhere.
by kate
Sat Nov 30, 2019 11:46 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS: SCRIP installation error
Replies: 6
Views: 1308

Re: PYROMS: SCRIP installation error

However you try installing things, try installing the jpeg library:

Code: Select all

ld: library not found for -ljpeg
The pyroms message about not finding scrip is just a warning. If you don't want to be remapping things with scrip, it is safe to ignore it.
by kate
Fri Nov 29, 2019 4:36 pm
Forum: ROMS Trivia
Topic: u,v high resolution database
Replies: 7
Views: 1796

Re: u,v high resolution database

Again, when and where? We have almost 40 years of Arctic hindcasts, but we'd have to load up a hard drive to get it to you. Or you can run your own regional domain using the global HYCOM boundaries, then run a nest within that. I do that in offline mode quite routinely.
by kate
Fri Nov 29, 2019 3:47 am
Forum: ROMS Trivia
Topic: u,v high resolution database
Replies: 7
Views: 1796

Re: u,v high resolution database

Is there a time and a place you are thinking of? There are daily global HYCOM fields out there. https://www.hycom.org/dataserver
by kate
Wed Nov 27, 2019 5:12 pm
Forum: ROMS Discussion
Topic: restart file error
Replies: 4
Views: 701

Re: restart file error

Trial and error? My runs can go for months, then suddenly blow up. It is often an instability in momentum, in the middle of the water column, just at one point. In that case, trying a shorter timestep gets me through it.
by kate
Tue Nov 26, 2019 6:22 pm
Forum: ROMS Discussion
Topic: restart file error
Replies: 4
Views: 701

Re: restart file error

I think this is actually a ROMS issue. I had a recent run fail the same way. I "fixed" it by changing DSTART to an earlier date. I was reading a restart file and had changed DSTART to match the restart time - that's the job that failed. The job that ran didn't need to read an averages file (I didn't...
by kate
Mon Nov 25, 2019 9:52 pm
Forum: ROMS Problems
Topic: variables don't change at the dimension of time
Replies: 2
Views: 576

Re: variables don't change at the dimension of time

If the model is blowing up, it has to be changing in time.

Looking at the surface is looking at the highest k level. To look at other depths, there are codes for computing that in matlab, python, etc. What are you using now?
by kate
Sat Nov 23, 2019 4:44 am
Forum: ROMS Installation
Topic: Compiling error - undefined reference to ice routines
Replies: 12
Views: 2113

Re: Compiling error - undefined reference to ice routines

The code is such that one could make a simpler version of ice_mk.h and include that instead. It's all thermodynamics, not dynamics.
by kate
Thu Nov 21, 2019 8:31 pm
Forum: ROMS Installation
Topic: Compiling error - undefined reference to ice routines
Replies: 12
Views: 2113

Re: Compiling error - undefined reference to ice routines

Switching on ICE_THERMO should be independent of the ice dynamics (rheology and such). However, you are delving into uncharted territory and should be prepared to look at code to see what is really going on. Perhaps even look at the .f90 files to see what is left without the ice dynamics.
by kate
Thu Nov 21, 2019 7:28 am
Forum: ROMS Benchmarks
Topic: Speed-up advice?
Replies: 2
Views: 970

Re: Speed-up advice?

Thanks, John, I'd forgotten about that TS_HSIMT option. I'll give it a whirl.

I haven't tried other tilings. I heard from their consultants saying the two systems are remarkably similar so they would expect the performance I got. Cheyenne has 36 cores/node.
by kate
Tue Nov 19, 2019 9:12 pm
Forum: ROMS Benchmarks
Topic: Speed-up advice?
Replies: 2
Views: 970

Speed-up advice?

We are hoping to have a long simulation of a Pan-Arctic grid with a BGC model. My colleagues have therefore obtained a computing grant at NCAR. I have done a few timings of a day-long simulation on Cheyenne. I'm finding it no faster than a not-so-new Penguin Supercomputing cluster at home: On 216 co...
by kate
Tue Nov 19, 2019 6:50 pm
Forum: ROMS Installation
Topic: Compiling error - undefined reference to ice routines
Replies: 12
Views: 2113

Re: Compiling error - undefined reference to ice routines

Actually, there is an AInudgcof for nudging to ice concentration climatology. See the Apps_master/Chukchi or the Apps/Beaufort versions of ana_nudgcoef.h.
by kate
Sun Nov 17, 2019 4:13 pm
Forum: ROMS Discussion
Topic: unmatch surface lateral boundary condition data
Replies: 2
Views: 632

Re: unmatch surface lateral boundary condition data

I would start by upping the temporal frequency of your atmospheric fields. Monthly is going wash out all the storms. I use three-hourly to resolve diurnal cycles as well. I haven't tried monthly boundary information since we stopped using monthly Levitus many, many moons ago, but going to more frequ...
by kate
Sun Nov 17, 2019 4:57 am
Forum: ROMS Installation
Topic: Compiling error - undefined reference to ice routines
Replies: 12
Views: 2113

Re: Compiling error - undefined reference to ice routines

COAWST has my sea ice in it now, a one-layer dynamic, thermodynamic sea ice model as part of ROMS. This is a rather simpler model than CICE, the community sea ice model. There are several efforts to couple ROMS with CICE, but that may not be what you want now. There's an effort doing sea ice in a mo...
by kate
Fri Nov 15, 2019 8:00 pm
Forum: ROMS Installation
Topic: Compiling error - undefined reference to ice routines
Replies: 12
Views: 2113

Re: Compiling error - undefined reference to ice routines

There are long-standing disagreements between me and Hernan on a number of things, including where the ice code should live and how to compile. I maintain my own Makefiles which point to the ice code where I put it. Since he's not actually distributing the ice code I feel no pressure to comply. Call...
by kate
Thu Nov 14, 2019 10:37 pm
Forum: ROMS Tools and Techniques
Topic: pyroms grid missing libgridgen, unable to solve
Replies: 5
Views: 1203

Re: pyroms grid missing libgridgen, unable to solve

Most people separate the vertical grid from the horizontal grid. The ROMS grid file really only has the horizontal grid information while the vertical grid comes from the stretching parameters, THETA_S and so forth. You should be able to create a horizontal grid file with those other tools, then rea...
by kate
Thu Nov 14, 2019 8:57 pm
Forum: ROMS Tools and Techniques
Topic: pyroms grid missing libgridgen, unable to solve
Replies: 5
Views: 1203

Re: pyroms grid missing libgridgen, unable to solve

If it is grid generation you are after, I suggest you check out pygridgen and pygridtools . I'm sorry this happened, but pyroms installation has always been troublesome and Mark Hadfield has offered to work on cleaning things up. There has been progress, but there are bumps in the road, as you can s...
by kate
Mon Nov 11, 2019 5:34 pm
Forum: ROMS Discussion
Topic: Is it a must to use VolCons if RadNud activated?
Replies: 8
Views: 1373

Re: Is it a must to use VolCons if RadNud activated?

It's blowing up at point 2, 13, 8. Your Courant numbers are getting rather large, so you might start by trying a shorter timestep. You could also look at the fields it has saved for you in the restart file. 122 2018-01-01 20:20:00.00 3.461952E-03 7.113528E+02 7.113563E+02 5.699591E+12 (002,013,10) 1...
by kate
Thu Nov 07, 2019 6:14 pm
Forum: ROMS Discussion
Topic: How to set up outgoing nudging in nud netCDF
Replies: 2
Views: 519

Re: How to set up outgoing nudging in nud netCDF

The answer to this question is in the source code, say in t3dbc_im.F. It looks like it is getting the outgoing value from the nudging field you have specified and the incoming value is obcfac times that. Remember that it has been scaled to be an inverse timescale before being used, so an obcfac of 1...
by kate
Thu Nov 07, 2019 4:36 am
Forum: ROMS Bugs
Topic: DIAG speed trouble
Replies: 4
Views: 1111

Re: DIAG speed trouble

GET_NGFLD - river runoff mass transport, 2000-01-08 00:00:00.00 (Grid= 01, Rec=0000008, Index=1, File: create-river_daily2017_cycle_add9rivers_30layers-15-change-huanghe-1996coastl_dye.nc) (Tmin= 0.0000 Tmax= 364.0000) t = 7.0000 (Min = -1.27772940E+23 Max = 2.81650723E-37) GET_NGFLD - river runoff...
by kate
Sun Nov 03, 2019 5:08 pm
Forum: ROMS Discussion
Topic: error in comile a file
Replies: 2
Views: 520

Re: error in comile a file

Actually the error is that it can't find the makefile. Do you have the correct path to your ROMS sources? The makefile should be in the top directory there.
by kate
Thu Oct 31, 2019 5:09 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 28
Views: 2762

Re: Problem with setting up inlet_test

And yet:
INP_PAR - Unable to open ROMS/TOMS input script file.
coupling_inlet_test.in is in your current run directory? Plus any files it refers to?
by kate
Wed Oct 30, 2019 5:14 pm
Forum: ROMS Problems
Topic: problems with the addition of tide forcing
Replies: 9
Views: 3363

Re: problems with the addition of tide forcing

There are errors and warnings coming from SWAN. Is this from COAWST or not?
by kate
Tue Oct 29, 2019 5:48 pm
Forum: ROMS Problems
Topic: problems with the addition of tide forcing
Replies: 9
Views: 3363

Re: problems with the addition of tide forcing

You're missing parts of that file - where did you get it from? How did that even compile, since it is missing an ending #endif or two? For the problem you are seeing, before the #ifdef SSW_BBL bit, there should be: /* define only one of the following 5 */ #undef UV_LOGDRAG #undef UV_QDRAG #undef MB_...
by kate
Mon Oct 28, 2019 6:16 pm
Forum: ROMS Problems
Topic: problems with the addition of tide forcing
Replies: 9
Views: 3363

Re: problems with the addition of tide forcing

inlet_test.h has SSW_BBL defined, which should satisfy that bottom drag condition.
by kate
Mon Oct 28, 2019 5:21 pm
Forum: ROMS Problems
Topic: problems with the addition of tide forcing
Replies: 9
Views: 3363

Re: problems with the addition of tide forcing

Which case are you trying to run? I suggest starting with a proven case like UPWELLING. There is an upwelling.h file containing the cppdefs for that case.
by kate
Sun Oct 27, 2019 6:12 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 28
Views: 2762

Re: Problem with setting up inlet_test

Error: Invalid character in name at (1) make: *** [/Users/venuscheung/ocean_model/inlet_test/Build_roms/ocpmix.o] Error 1 I would look at these files in an editor which shows hidden characters, including the MakeDepend file. The vi(m) mode to see hidden characters is ":set list" - others can say ho...
by kate
Thu Oct 24, 2019 7:14 pm
Forum: ROMS Bugs
Topic: Error while compiling ROMS
Replies: 4
Views: 870

Re: Error while compiling ROMS

Since my Perl is dated 2015 and yours is rather older, it might help if you could try a newer Perl.
by kate
Thu Oct 24, 2019 6:49 pm
Forum: ROMS Bugs
Topic: Error while compiling ROMS
Replies: 4
Views: 870

Re: Error while compiling ROMS

This is a problem with your Perl. What version of Perl is this? My Perl works:
chinook03.rcs.alaska.edu 435% perl -v

This is perl 5, version 22, subversion 1 (v5.22.1) built for x86_64-linux-thread-multi
by kate
Wed Oct 23, 2019 3:45 pm
Forum: ROMS Source
Topic: Huon_eastward and Hvom_northward
Replies: 9
Views: 1586

Re: Huon_eastward and Hvom_northward

If you want transport across some section not aligned with the grid, the proper thing to do is described here.
by kate
Tue Oct 22, 2019 2:02 pm
Forum: ROMS Discussion
Topic: Tracers don't radiate out at boundary
Replies: 4
Views: 881

Re: Tracers don't radiate out at boundary

No, that's not the problem. I use those numbers too. I haven't tried monthly HYCOM, but back in the day we would try monthly Levitus and have a terrible time with the boundaries. I use daily HYCOM and I also nudge to monthly HYCOM in a band around the edge of the domain.
by kate
Mon Oct 21, 2019 6:03 pm
Forum: ROMS Source
Topic: Huon_eastward and Hvom_northward
Replies: 9
Views: 1586

Re: Huon_eastward and Hvom_northward

Rotating a velocity is easy while rotating flux through a cell face is less so. You'd have to think carefully about exactly what it is that you want to output, especially if dx and dy are not the same. As for the averages, adding new things to the averages output typically depends on adding new vari...
by kate
Mon Oct 21, 2019 5:54 pm
Forum: ROMS Discussion
Topic: Shchepetkin (Shc) boundary condition reference?
Replies: 1
Views: 442

Re: Shchepetkin (Shc) boundary condition reference?

I believe it to be:
@article(Mason2010,
author = "E. Mason and J. Molemaker and A.~F. Shcheptkin and F. Colas
and J.~C. McWilliams and P. Sangr{\`{a}}",
title = "Procedures for offline grid nesting in regional ocean models",
year = "2010",
journal = om,
volume = "35",
pages = "1--15"
by kate
Mon Oct 21, 2019 5:51 pm
Forum: ROMS Discussion
Topic: Tracers don't radiate out at boundary
Replies: 4
Views: 881

Re: Tracers don't radiate out at boundary

Turning off the temperature of the river source is numerically unstable.

Do you mean that you use monthly means from HYCOM for the boundary files? I have always had better luck with more frequent boundary information and even then it's not perfect.
by kate
Sat Oct 12, 2019 3:35 am
Forum: ROMS Installation
Topic: Finished compiling COAWST Installation - Ubuntu 14.04
Replies: 26
Views: 16358

Re: Finished compiling COAWST Installation - Ubuntu 14.04

The modern way to install netcdf-fortran on Ubuntu is:
sudo apt install netcdf-bin libnetcdf-dev libnetcdff-dev
Then you will have to check to see if nf-config is working because some versions gave a weird error about cmake.
by kate
Wed Oct 09, 2019 5:23 am
Forum: ROMS Discussion
Topic: how to improve vertical profile velocity
Replies: 10
Views: 1958

Re: how to improve vertical profile velocity

1. Probably easiest in Matlab or something with similar capabilities, something that knows the vertical coordinate system. 2. If you only want to nudge to surface T/S, there is separate code for that as part of the surface boundary condition. Check out for instance SCORRECTION and SRELAXATION, using...
by kate
Wed Oct 02, 2019 4:00 pm
Forum: ROMS Installation
Topic: Compile error
Replies: 1
Views: 582

Re: Compile error

cp: /opt/gfortransoft/serial/netcdf3/include/netcdf.mod: No such file or directory
make: *** No rule to make target
You need to have netcdf and netcdff installed and you need to know where they are. Best to get netcdf4 and set USE_NETCDF4 which has an nc-config to tell you where things are.
by kate
Mon Sep 23, 2019 10:45 pm
Forum: ROMS Discussion
Topic: how to improve vertical profile velocity
Replies: 10
Views: 1958

Re: how to improve vertical profile velocity

If you turn on climatology nudging, it will change the OBC time scales. If you don't want that, change the open boundary code to do what you want it to do - that's what I did in my branch.
by kate
Mon Sep 23, 2019 10:42 pm
Forum: ROMS Discussion
Topic: Asking for cpp and control files
Replies: 1
Views: 416

Re: Asking for cpp and control files

At line 960 of file inp_par.f90
You need to look at your inp_par.f90 and see what is in there at line 960.
by kate
Sun Sep 22, 2019 6:20 pm
Forum: ROMS Discussion
Topic: Amplification in zeta after adding river forcing
Replies: 2
Views: 589

Re: Amplification in zeta after adding river forcing

What about the surface salinity and velocity? Are you getting a run-away estuary where the surface is too fresh and fails to vertically mix down? I would expect there to be enough fresh water input there to change the results, but are they now more sensible or not? I can't tell from zeta.
by kate
Sun Sep 22, 2019 6:08 pm
Forum: ROMS Discussion
Topic: PYROMS error
Replies: 2
Views: 604

Re: PYROMS error

For #1, I get the same error, I'm not using lpsolve, I don't care enough to figure it out. For #2, I don't get that error. It is finding _remapping.cpython-37m-x86_64-linux-gnu.so in my $HOME/python/lib/python3.7/site-packages/pyroms which is where I told it to put stuff. Do you have such a _remappi...
by kate
Fri Sep 20, 2019 6:00 pm
Forum: ROMS Discussion
Topic: Some confusion about boundary condition
Replies: 4
Views: 760

Re: Some confusion about boundary condition

If any part of a boundary is open, you have to set the whole boundary to open. The land masking will take care of the rest. Your north can be closed. If you set it to open, no harm done, but you'll be asking ROMS to do extra work of reading fields (or computing analytic fields).
by kate
Thu Sep 19, 2019 5:15 pm
Forum: ROMS Discussion
Topic: Some confusion about boundary condition
Replies: 4
Views: 760

Re: Some confusion about boundary condition

Yes, ROMS is expecting ubar, vbar in the initial and boundary files. It will ensure consistency between baroclinic and barotropic flows during initialization and after every timestep, but the boundary files ought to be generated in such a way that the u and ubar fields are consistent. If you are doi...
by kate
Thu Sep 19, 2019 5:07 pm
Forum: ROMS Bugs
Topic: Some questions about river input
Replies: 4
Views: 845

Re: Some questions about river input

You want river inputs to be from land to water as in your left arrow. The right arrow appears to be all within the land and will therefore do nothing at all. How much transport are you asking for? What is the corresponding flow rate to go with it? Are you violating the ROMS speed limit? I have done ...
by kate
Thu Sep 19, 2019 3:19 am
Forum: ROMS Bugs
Topic: Some questions about river input
Replies: 4
Views: 845

Re: Some questions about river input

Are we looking at depth there? Where is the land mask? Are you sure you are putting the rivers where you think you are? Blows up how? Are you providing the river tracers (T and S)?
by kate
Wed Sep 18, 2019 5:07 pm
Forum: ROMS Discussion
Topic: WET_DRY option results in too much land_mask?
Replies: 2
Views: 482

Re: WET_DRY option results in too much land_mask?

You are telling the model that it can never have less than 10 m water depth in the "wet" parts. You might try something rather smaller for DCRIT, like 0.5 m. Then you'd have to have 9.5 m of water leaving due to the storm surge before a spot became "dry".
by kate
Tue Sep 17, 2019 4:22 pm
Forum: ROMS Discussion
Topic: erro about MPDATA
Replies: 2
Views: 604

Re: erro about MPDADA

Any of the other advection schemes should work. Is there any reason you can't use MPI on that shared-memory cluster?
by kate
Sat Sep 14, 2019 3:25 am
Forum: ROMS Discussion
Topic: open boundary for future projection of CMIP6 downscaling
Replies: 11
Views: 1716

Re: open boundary for future projection of CMIP6 downscaling

Salt looks like you have boundary weirdness, low values right at the edge. Zeta seems to be in balance with a strong boundary current which I'm not seeing in the u, v plots. Have you looked at the boundary files to see if the fields there look reasonable?

What values did you pick for LBC all around?
by kate
Sat Sep 14, 2019 1:22 am
Forum: ROMS Discussion
Topic: open boundary for future projection of CMIP6 downscaling
Replies: 11
Views: 1716

Re: open boundary for future projection of CMIP6 downscaling

The blow-up messages are not helpful. You need to look at the fields in the restart record that it saves.
by kate
Fri Sep 13, 2019 5:44 pm
Forum: ROMS Discussion
Topic: open boundary for future projection of CMIP6 downscaling
Replies: 11
Views: 1716

Re: open boundary for future projection of CMIP6 downscaling

Well, I also hope to use CMIP for open boundaries in the future. Which of the models did you try and what software did you use to create the OBC files? Did the blow-up happen at the boundary and what did it look like?
by kate
Thu Sep 12, 2019 5:53 am
Forum: ROMS Discussion
Topic: Is it a must to use VolCons if RadNud activated?
Replies: 8
Views: 1373

Re: Is it a must to use VolCons if RadNud activated?

I never, ever use VolCons. If you have the Flather boundary condition, the model should not fill or drain excessively - unless of course the zeta boundary conditions you give it suggest that it should fill or drain excessively.
by kate
Wed Sep 11, 2019 4:56 pm
Forum: ROMS Discussion
Topic: if I can set hc=Tcline>hmin?
Replies: 7
Views: 1124

Re: if I can set hc=Tcline>hmin?

hmin is no longer used in computing z. If you keep reading in the Warning, it says the restriction was removed.
by kate
Wed Sep 11, 2019 4:52 pm
Forum: ROMS Bugs
Topic: DIAG speed trouble / gridbuilder
Replies: 3
Views: 914

Re: DIAG speed trouble / gridbuilder

You might still need to reduce the timestep. You ran long enough that you should be getting close.
by kate
Wed Sep 11, 2019 4:21 pm
Forum: ROMS Installation
Topic: Error during compilation of build file
Replies: 1
Views: 560

Re: Error during compilation of build file

I think you should have: NETCEDF_INCDIR=/opt/gfortransoft/serial/netcdf3/include if that's where the netcdf.mod and friends live. The FIRST_2D_STEP is a cpp macro defined in globaldefs.h, not a regular Fortran variable. If the compile goes correctly, the compiler should never see it. Step one is a "...
by kate
Wed Sep 11, 2019 3:08 am
Forum: ROMS Discussion
Topic: ROMS processing water sluice in upstream
Replies: 1
Views: 379

Re: ROMS processing water sluice in upstream

I think that should work. What would be problematic is a source turning into a sink sometimes, because the tracers need to come from upstream to be stable.
by kate
Tue Sep 10, 2019 11:23 pm
Forum: ROMS Discussion
Topic: if I can set hc=Tcline>hmin?
Replies: 7
Views: 1124

Re: if I can set hc=Tcline>hmin?

Actually, no, you want hc = Tcline = 10 m. It was only with the bad old stretching that the code enforced hc = Tcline = hmin. The stretching is meant to be "z-like" above Tcline (hc). Having the old way meant that the z-like part of your domain could only be above hmin. Meanwhile, I was getting away...
by kate
Mon Sep 09, 2019 4:44 pm
Forum: ROMS Discussion
Topic: if I can set hc=Tcline>hmin?
Replies: 7
Views: 1124

Re: if I can set hc=Tcline>hmin?

It depends on which vertical "transform" function you use. The old one needed Tcline <= hmin, but the new one does not. I'm sure there's plenty of code still around from the bad old days. I use: Vtransform == 2 ! transformation equation Vstretching == 4 ! stretching function THETA_S == 7.0d0 ! surfa...
by kate
Fri Sep 06, 2019 8:58 pm
Forum: ROMS Discussion
Topic: The xi_position and eta_position of rivers.nc input
Replies: 4
Views: 770

Re: The xi_position and eta_position of rivers.nc input

It refers to i, j. The western-most u-point is at i=1.
by kate
Mon Sep 02, 2019 4:06 pm
Forum: ROMS Discussion
Topic: surface boundary layer depth
Replies: 2
Views: 600

Re: surface boundary layer depth

You can have positive Hsbl where zeta is positive. The boundary layer thickness is (zeta-Hsbl).
by kate
Sun Sep 01, 2019 3:13 pm
Forum: ROMS Problems
Topic: Interannual simulation
Replies: 8
Views: 1732

Re: Interannual simulation

It depends on your ROMS version. Older versions required you to change NTIMES to get a second year while newer ones don't.
by kate
Fri Aug 30, 2019 3:05 pm
Forum: ROMS Problems
Topic: Interannual simulation
Replies: 8
Views: 1732

Re: Interannual simulation

I don't see an error in the output you show. That's just the normal ROMS report at the end of a run. How many timesteps did you ask it to do?
by kate
Thu Aug 29, 2019 11:03 pm
Forum: ROMS Tools and Techniques
Topic: problem with pyroms.grid.Gridgen
Replies: 8
Views: 3950

Re: problem with pyroms.grid.Gridgen

For generating grids, you might look into pygridtools. They seem more advanced in their Python support.
by kate
Thu Aug 29, 2019 5:12 am
Forum: ROMS Discussion
Topic: i have a question
Replies: 11
Views: 3390

Re: i have a question

I don't see what the problem is. Yes, it is reading the forcing files for time 0 and for time 6 hours, then it linearly time-interpolates between them to the current model time. It always has two time records in memory for this interpolation. Once six hours pass, it will read hour 12, and so on.
by kate
Wed Aug 28, 2019 3:29 pm
Forum: ROMS Problems
Topic: Interannual simulation
Replies: 8
Views: 1732

Re: Interannual simulation

#define TCLM_NUDGING #define TCLIMATOLOGY This is from an older ROMS version. #define BULF_FLUXES /*bulk formulation for surface heat fluxes*/ You have misspelled BULK_FLUXES. If you spell it correctly, it will ask for winds instead of surface stress. ROMS will ignore random strings that it doesn't...
by kate
Tue Aug 27, 2019 4:43 pm
Forum: ROMS Problems
Topic: Interannual simulation
Replies: 8
Views: 1732

Re: Interannual simulation

It blew up in the first step - this could be any number of things. The first thing I would try is a shorter barotropic timestep - i.e., a larger NDTFAST. Your surface inputs appear to be well behaved. What other inputs do you have?
by kate
Mon Aug 26, 2019 6:00 pm
Forum: ROMS Problems
Topic: nudging/relaxation time scale field
Replies: 1
Views: 613

Re: nudging/relaxation time scale field

Answers 1, 2 and 3: In t3dbc_im.F, the nudging code has a different nudging constant depending on whether or not LnudgeTCLM is true. You can search for LnudgeTCLM there to see. That is not the behavior I want, so I have modified my code accordingly.
by kate
Fri Aug 23, 2019 4:18 pm
Forum: ROMS Discussion
Topic: Barotropic model Blowing-up with high courant and dt = 1.
Replies: 6
Views: 921

Re: Barotropic model Blowing-up with high courant and dt = 1

The OP said it was a barotropic model which requires NDTFAST = 1.
by kate
Fri Aug 23, 2019 1:58 pm
Forum: ROMS Discussion
Topic: Barotropic model Blowing-up with high courant and dt = 1.
Replies: 6
Views: 921

Re: Barotropic model Blowing-up with high courant and dt = 1

Why not try dt=0.5 or shorter? Your horizontal grid spacing is quite a bit smaller than anything I've tried outside of test problems, so it makes sense that you'd need a short timestep.
by kate
Fri Aug 23, 2019 5:02 am
Forum: ROMS Discussion
Topic: blow up the model
Replies: 3
Views: 639

Re: blow up the model

So save a record right before it blows up and look at it. I can't do it from here.
by kate
Thu Aug 22, 2019 4:06 pm
Forum: ROMS Discussion
Topic: Shortwave and Longwave radiation.
Replies: 3
Views: 938

Re: Shortwave and Longwave radiation.

I always let ROMS compute its own upward longwave flux. There's an option for that.

For the shortwave, it might be worth trying what you say. I use an albedo option because the ice model wants to have its own albedo.