Search found 3930 matches

by kate
Fri May 07, 2021 5:11 pm
Forum: ROMS Ice
Topic: Landfast ice
Replies: 2
Views: 2044

Re: Landfast ice

You may well be right, but as I said elsewhere, I'm no longer maintaining that code.
by kate
Fri May 07, 2021 5:08 pm
Forum: ROMS Ice
Topic: How to use SEAICE
Replies: 3
Views: 103

Re: How to use SEAICE

It is true that Paul Budgell put the ICE_SHOREFAST in there for the Caspian, but that was a long, long time ago. The ICE_LANDFAST is what I would use now, it's the Lemieux landfast ice parameterization which is also in CICE6. If you need it in CICE5, I can help with that. Why would you read the manu...
by kate
Fri May 07, 2021 4:36 pm
Forum: ROMS Ice
Topic: How to use SEAICE
Replies: 3
Views: 103

Re: How to use SEAICE

I'm sorry, I'm no longer supporting that ice code. I'm five years into a plan to switch to using MOM6 instead, finally getting somewhere.

For ice with ROMS, the Norwegians are supporting their metroms package (but don't try to learn git from their example).
by kate
Wed May 05, 2021 6:03 pm
Forum: ROMS Discussion
Topic: How to split the advection and diffusion term
Replies: 4
Views: 378

Re: How to split the advection and diffusion term

The vertical diffusivity is in Akt, which is dimensioned by three spatial dimensions and the NT (number of tracers) dimension. You can set some to zero. There will still be implicit changes which I don't think you can do anything about other than try different advection schemes.
by kate
Tue May 04, 2021 6:28 pm
Forum: ROMS Discussion
Topic: How to split the advection and diffusion term
Replies: 4
Views: 378

Re: How to split the advection and diffusion term

In modeling, there is always implicit diffusion from your advection operator. You can reduce this by going to a fourth-order advection scheme. I use the third-order one with its implicit diffusion, then turn off explicit diffusion by turning off the cppdefs for it. Vertical diffusion is trickier. Pe...
by kate
Mon May 03, 2021 12:23 am
Forum: ROMS Discussion
Topic: Blow up ( River run off)
Replies: 7
Views: 770

Re: Blow up ( River run off)

Did you provide river tracer values?
by kate
Mon Apr 26, 2021 5:27 pm
Forum: ROMS Discussion
Topic: What's the difference between online and offline running?
Replies: 2
Views: 246

Re: What's the difference between online and offline running?

As I posted elsewhere, I've given up on offline biology. The point of it was that you run the physics once, saving daily outputs. Then you read those fields in as "climalotogy" for running the biology. The dream is that you can maybe take a longer timestep when you only have the biology going on and...
by kate
Mon Apr 26, 2021 5:21 pm
Forum: ROMS Bugs
Topic: main3d_offline.F and step3d_t.F in ROMS 3.9
Replies: 2
Views: 309

Re: main3d_offline.F and step3d_t.F in ROMS 3.9

If this is my code, I've given up on OFFLINE_BIOLOGY. I sort of had it working in the "main" branch, not the default "kate_svn" branch. But for me, I have wetting and drying and it wasn't doing the right thing in the dry parts, nor did it even know if it was wet or dry that day.
by kate
Sun Apr 25, 2021 7:03 pm
Forum: ROMS Discussion
Topic: Blow up
Replies: 2
Views: 125

Re: Blow up

This is when you need to become a detective. Which field is blowing up. You say it is in the land? Is that the only place going bad?
by kate
Sun Apr 25, 2021 7:01 pm
Forum: ROMS Discussion
Topic: A strange error while open the biological options
Replies: 5
Views: 538

Re: A strange error while open the biological options

I haven't used BIO_UMAINE since Hernan changed the advection scheme choices. Good catch, John. As for my github code, I have stopped updating it because updates can be disruptive and I don't want to spend the time sorting them out. Mark Hadfield is keeping his github code up to date. If you find any...
by kate
Sat Apr 24, 2021 7:58 pm
Forum: ROMS Discussion
Topic: A strange error while open the biological options
Replies: 5
Views: 538

Re: A strange error while open the biological options

The check you are running into is in diag.F: IF (maxbio(ibt,ng).gt.max_bio(ibt)) THEN It looks like you have max_bio set to 10000. Does that make sense for the units you are using? You can change it if 10000 is a valid number for your BGC model. I probably added the check for COSINE, but then kept i...
by kate
Fri Apr 23, 2021 10:19 pm
Forum: ROMS Discussion
Topic: problem in open boundary
Replies: 5
Views: 332

Re: problem in open boundary

That's great that the model is running without any problems. You're all set, then, right? You need to look at your solution to see if you are happy with it. Are you happy with it? If not, there are steps you can take to possibly improve your boundary solution, but nothing is guaranteed when it comes...
by kate
Fri Apr 23, 2021 5:41 pm
Forum: ROMS Discussion
Topic: problem in open boundary
Replies: 5
Views: 332

Re: problem in open boundary

How are you getting your ROMS information? There is both a manual and a wiki. The wiki page might help you. The manual has much the same information and also examples containing open boundaries.
by kate
Fri Apr 23, 2021 5:39 pm
Forum: ROMS Discussion
Topic: add passive tracer
Replies: 9
Views: 544

Re: add passive tracer

I see the wiki page on FLOATS is still empty. There is probably something about floats in the manual, though.
by kate
Fri Apr 23, 2021 1:29 am
Forum: ROMS Discussion
Topic: add passive tracer
Replies: 9
Views: 544

Re: add passive tracer

The FLOATS option is a whole different can of worms. Is that what you want? Online or offline?
by kate
Thu Apr 22, 2021 11:33 pm
Forum: ROMS Discussion
Topic: problem in open boundary
Replies: 5
Views: 332

Re: problem in open boundary

What sort of open boundary do you want? In the UPWELLING case, the southern edge is in shallow water and typically has a wall. Are you changing the bathymetry? Are you providing fields on the open boundary? I never use a sponge in the ROMS sense of the word. I do use nudging to a climatology, which ...
by kate
Thu Apr 22, 2021 3:29 am
Forum: ROMS Discussion
Topic: add passive tracer
Replies: 9
Views: 544

Re: add passive tracer

You need these as well:

Code: Select all

# define ANA_BPFLUX        /* analytical bottom passive tracers fluxes */
# define ANA_SPFLUX
unless you want the top and bottom tracer fluxes to be non-zero.
by kate
Wed Apr 21, 2021 5:27 pm
Forum: ROMS Discussion
Topic: add passive tracer
Replies: 9
Views: 544

Re: add passive tracer

First, you turn on T_PASSIVE and set NPT in roms.in. For the initial conditions, ANA_PASSIVE can be used, otherwise you'll have to put the initial conditions into a netcdf file. I've usually initialized with zero via ANA_PASSIVE, then used sources to introduce the passive dye locally. What is your s...
by kate
Mon Apr 19, 2021 6:07 pm
Forum: Job Opportunities
Topic: Position in Colorado
Replies: 0
Views: 355

Position in Colorado

The University of Colorado Cooperative Institute for Environmental Sciences (CIRES) jointly with the NOAA/Physical Sciences Lab is hiring an associate scientist in the area of ocean data analysis/visualization. Here is a brief position description: The chosen candidate will assist in the development...
by kate
Mon Apr 19, 2021 5:50 pm
Forum: Job Opportunities
Topic: Postdoc at Princeton/GFDL
Replies: 1
Views: 559

Re: Postdoc at Princeton/GFDL

This position is closing soon. If you intend to apply, get your application submitted now.
by kate
Fri Apr 16, 2021 4:53 pm
Forum: ROMS Discussion
Topic: Timestamp of ROMS diagnostic files
Replies: 1
Views: 117

Re: Timestamp of ROMS diagnostic files

Yes, the diagnostics are time-averaged.
by kate
Mon Apr 12, 2021 5:28 pm
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 912

Re: segmentation fault when running the COAWST model

Make sure NBCFILES is defined. Some roms.in files have it, some do not. If using COAWST, you need to have it.
by kate
Sat Apr 10, 2021 7:39 pm
Forum: ROMS Discussion
Topic: Direction of Uwind and Vwind
Replies: 2
Views: 162

Re: Direction of Uwind and Vwind

For the incoming 2-D forcing, it depends. ROMS will check the dimensions of each field you provide. If they match that of your ROMS grid, it assumes Uwind and Vwind are aligned in the directions of the grid. If not, it assumes it needs to do the interpolation and the rotation to the ROMS grid.
by kate
Mon Apr 05, 2021 4:14 pm
Forum: ROMS Problems
Topic: How to solve the speed trouble on the open boundaries?
Replies: 3
Views: 613

Re: How to solve the speed trouble on the open boundaries?

I use SODA as my initial condition, no spin-up period. You can try a spin-up and see how much drift you get.
by kate
Sun Apr 04, 2021 5:18 pm
Forum: ROMS Problems
Topic: How to solve the speed trouble on the open boundaries?
Replies: 3
Views: 613

Re: How to solve the speed trouble on the open boundaries?

I typically add a nudging to a monthly climatology in a band around the edge of the domain, ramping to zero in the middle. Note that this will overwrite TNUDG and friends in the default ROMS code.
by kate
Tue Mar 30, 2021 5:07 pm
Forum: ROMS Discussion
Topic: Era5 data is not read by the ROMS
Replies: 5
Views: 548

Re: Era5 data is not read by the ROMS

GET_2DFLD - surface u-wind component, 2013-01-01 00:00:00.00 (Grid=01, Rec=1, Index=1, File: era5_u10_year_2013_roms.nc) (Tmin= 0.0000 Tmax= 30.9583) t = 0.0000 (Min = 1.00000000E+35 Max = -1.00000000E+35) regrid = T GET_2DFLD - surface v-wind component, 2013-01-01 00:00:00.00 (Grid=01, Rec=1, Inde...
by kate
Mon Mar 29, 2021 8:11 pm
Forum: ROMS Discussion
Topic: Era5 data is not read by the ROMS
Replies: 5
Views: 548

Re: Era5 data is not read by the ROMS

NETCDF_GET_TIME_1D - error while inquiring ID for variable: srf_time in input file: era5_msdwswrf_year_2013_roms.nc What does "ncdump -h era5_msdwswrf_year_2013_roms.nc" look like? Is there a time variable? What is it called? ROMS is looking for srf_time and it could have gotten this from (a) the v...
by kate
Sat Mar 27, 2021 2:45 am
Forum: ROMS Discussion
Topic: ROMS blowing up problem
Replies: 15
Views: 3790

Re: ROMS blowing up problem

Same answer as above - use a shorter timestep and look more closely at your model output files.
by kate
Tue Mar 23, 2021 5:00 pm
Forum: ROMS Problems
Topic: Speeding up ROMS Model simulation
Replies: 6
Views: 610

Re: Speeding up ROMS Model simulation

You know your DT and you should be able to figure out how many timesteps per day from that. Say DT=120 leads to 720 steps/day. Three days of steps is therefore 2160 steps.
by kate
Sun Mar 21, 2021 6:25 pm
Forum: ROMS Problems
Topic: Speeding up ROMS Model simulation
Replies: 6
Views: 610

Re: Speeding up ROMS Model simulation

When going coarser, you can take longer timesteps. When going coarser, are you reducing the number of grid points?

Why not use MPI?
by kate
Sun Mar 21, 2021 6:22 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 2728

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

The debugger I use is called totalview. There is also ddt. Both are quite expensive, but the supercomputer center I use pays for the license. For free, there is gdb, which works if you run on just one processor. There may be a gui (ddd?) that can make it easier to use, but using these things is stil...
by kate
Sat Mar 20, 2021 5:27 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 2728

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

We know it's dying for some unknown reason. So maybe it didn't actually finish ana_initial. I would look in there. Did you change anything there?

Personally, I would run this in the debugger to watch it go bad.
by kate
Thu Mar 18, 2021 4:49 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 2728

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

You know it got to ANA_INITIAL because it printed stuff from there. I think the next calls it would make from initial.F are these: CALL set_depth0 (ng, tile, iNLM) CALL set_depth (ng, tile, iNLM) CALL set_massflux (ng, tile, iNLM) I would add print statements before these, then after these, then aft...
by kate
Tue Mar 16, 2021 4:25 pm
Forum: ROMS Problems
Topic: Fatal erorr after update to ubuntu 20.04
Replies: 2
Views: 146

Re: Fatal erorr after update to ubuntu 20.04

What I would do is recompile netcdf-fortran. Did you compile it yourself before or use the one from "apt install"?
by kate
Tue Mar 16, 2021 6:41 am
Forum: ROMS Discussion
Topic: ROMS blowing up problem
Replies: 15
Views: 3790

Re: ROMS blowing up problem

I'll say what I always say: try a shorter timestep and look more closely at your solution to see where it is blowing up.
by kate
Sun Mar 14, 2021 3:04 pm
Forum: ROMS Installation
Topic: upwelling - linux-ifort.mk compilation error
Replies: 5
Views: 377

Re: upwelling - linux-ifort.mk compilation error

Do you have USE_NETCDF4? Do you have netcdf installed properly?
by kate
Sun Mar 14, 2021 7:37 am
Forum: ROMS Installation
Topic: upwelling - linux-ifort.mk compilation error
Replies: 5
Views: 377

Re: upwelling - linux-ifort.mk compilation error

I think you want to have USE_MPIF90 turned on.
by kate
Mon Mar 08, 2021 7:24 pm
Forum: ROMS Discussion
Topic: The defination of Hz
Replies: 2
Views: 210

Re: The defination of Hz

Hz is the height of the vertical grid box, which you could call Hz=∂z/∂σ.
For consistency between the baroclinic and barotropic solutions, you also need D=ΣHz.
by kate
Mon Mar 08, 2021 7:19 pm
Forum: Job Opportunities
Topic: Postdoc at Princeton/GFDL
Replies: 1
Views: 559

Postdoc at Princeton/GFDL

I wanted to let everyone know that Princeton is now accepting applications for a Princeton/GFDL postdoc to further develop the regional capabilities of the MOM6 ocean model: https://puwebp.princeton.edu/AcadHire/apply/application.xhtml?listingId=19521 The position would include algorithmic developm...
by kate
Sun Mar 07, 2021 5:25 pm
Forum: ROMS Discussion
Topic: Modify UPWELLING case to a flat bottom, etc.
Replies: 3
Views: 329

Re: Modify UPWELLING case to a flat bottom, etc.

I would not edit UPWELLING, but rather create a new name for your case. Then search the whole code for UPWELLING so you know what files to look in. For the flat bathymetry, the "else" case on ana_grid.h should take care of that for you. For the wind stress, you should be able to find that in the sam...
by kate
Fri Mar 05, 2021 5:44 pm
Forum: ROMS Discussion
Topic: various time indices in ROMS
Replies: 2
Views: 206

Re: various time indices in ROMS

There is some discussion of the timestepping in the roms manual . Most of them point to the old/new time which alternates between 1 and 2, while 3 is always a working time level for tracers. Note: on restart, the model doesn't save these things and always assumes whatever is true for the first step....
by kate
Tue Mar 02, 2021 7:19 pm
Forum: ROMS Discussion
Topic: MPI_ERR_TRUNCATE
Replies: 4
Views: 302

Re: MPI_ERR_TRUNCATE

initial_ 219 initial.f90
What is on line 219 of your initial.f90? This would be in your scratch directory.
by kate
Tue Mar 02, 2021 6:03 pm
Forum: ROMS Discussion
Topic: nested grid with pyroms: tutorial and python routines
Replies: 15
Views: 1906

Re: nested grid with pyroms: tutorial and python routines

Yes, you got it.

I don't use Matlab if I can help it. You should get the script from the ROMS Matlab repository.
by kate
Tue Mar 02, 2021 8:20 am
Forum: ROMS Discussion
Topic: nested grid with pyroms: tutorial and python routines
Replies: 15
Views: 1906

Re: nested grid with pyroms: tutorial and python routines

Online nesting requires a special 1:3 or 1:5 relationship between the grids, as in Lake Jersey. Both models are run at the same time and the computer spends a lot of time communicating between the domains for you. That's what the contact stuff is for. I did this only once and it was pretty painfully...
by kate
Tue Mar 02, 2021 8:01 am
Forum: ROMS Discussion
Topic: nested grid with pyroms: tutorial and python routines
Replies: 15
Views: 1906

Re: nested grid with pyroms: tutorial and python routines

You are talking about doing online nesting in that case. I have a lot less experience with that than with offline nesting.
by kate
Tue Mar 02, 2021 7:19 am
Forum: ROMS Discussion
Topic: nested grid with pyroms: tutorial and python routines
Replies: 15
Views: 1906

Re: nested grid with pyroms: tutorial and python routines

When doing offline nesting, each domain is a fully independent domain. You run the larger, coarser one first, then use fields from that domain to create the initial/boundary files for the smaller domain. This creation of boundary files for the child grid happens offline, for instance with pyroms, bu...
by kate
Sun Feb 28, 2021 6:34 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 2728

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

I always have USE_MPIF90 on so that it knows where to find mpif.h.
Are you using COAWST? It is recommended to use that code when coupling to SWAN. Not that I've done it with any of the versions.
by kate
Wed Feb 17, 2021 5:53 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 2728

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: 445

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: 17
Views: 1313

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: 37
Views: 2728

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: 37
Views: 2728

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: 668

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: 668

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: 37
Views: 2728

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: 37
Views: 2728

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: 10
Views: 1045

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: 37
Views: 2728

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: 37
Views: 2728

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: 37
Views: 2728

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: 37
Views: 2728

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: 4
Views: 514

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: 452

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: 452

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: 788

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: 788

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: 788

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: 788

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: 788

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: 788

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: 788

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: 216

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: 11
Views: 1355

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: 875

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: 875

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: 875

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: 747

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: 11
Views: 1355

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: 11
Views: 1355

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: 7260

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: 5942

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: 5942

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: 5942

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: 5942

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: 820

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: 820

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: 820

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: 1472

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: 1472

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: 1472

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: 1472

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: 1472

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: 1472

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: 1472

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: 1367

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: 825

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: 825

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: 1367

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: 365

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: 1367

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: 1367

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: 1367

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: 1426

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: 1018

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: 1018

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: 361

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: 5010

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: 200

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: 418

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: 418

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: 550

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: 550

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: 1111

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: 991

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: 254

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: 267

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: 435

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: 435

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: 391

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: 391

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: 466

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: 427

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: 427

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: 427

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: 630

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: 630

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: 630

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: 459

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: 570

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: 643

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: 458

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: 575

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: 421

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: 593

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: 581

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: 501

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: 1149

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: 494

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: 769

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: 769

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: 424

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: 1361

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: 1004

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: 1004

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: 1537

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: 745

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: 745

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: 745

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: 1361

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: 6162

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: 1537

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: 985

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: 1537

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: 1537

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: 585

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: 985

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: 786

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: 2479

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: 2479

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: 896

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: 783

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: 753

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: 965

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: 1846

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: 443

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: 1846

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: 1846

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: 600

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: 15
Views: 1906

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: 1146

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: 1146

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: 844

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: 844

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: 1925

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 1458 times
by kate
Tue Mar 17, 2020 5:24 pm
Forum: ROMS Problems
Topic: Blowing up problem
Replies: 9
Views: 1925

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: 1925

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: 697

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: 629

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: 2118

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: 2118

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: 616

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: 2118

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: 2118

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: 2118

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: 2201

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: 1179

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: 1914

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: 3527

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: 1179

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: 1914

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: 1019

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: 1914

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: 1914

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: 1914

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: 1676

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: 1676

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: 1676

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: 1676

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: 1676

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: 590

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: 1077

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: 1077

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: 1414

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: 1914

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: 1414

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: 1914

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: 1414

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: 12
Views: 2483

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: 12
Views: 2483

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: 397

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: 1739

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: 1739

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: 1739

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: 2861

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: 1812

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: 1812

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: 2659

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: 1812

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: 1812

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: 1812

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: 4542

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: 15
Views: 1906

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: 2035

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: 2035

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: 15
Views: 1906

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: 2035

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: 784

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: 1489

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: 1998

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: 1998

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: 760

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: 760

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: 637

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: 2311

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: 2311

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: 1042

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: 1042

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: 2311

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: 698

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: 2311

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: 2311

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: 1629

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: 1629

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: 1496

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: 579

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: 1214

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: 582

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: 2995

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: 3522

Re: problems with the addition of tide forcing

There are errors and warnings coming from SWAN. Is this from COAWST or not?