Search found 3799 matches

by kate
Wed Sep 23, 2020 4:25 pm
Forum: ROMS Bugs
Topic: Error with wclock_on subroutine
Replies: 1
Views: 219

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: 1
Views: 74

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Re: An error for two-way nesting with restart

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Re: Segmentation fault

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

Re: d_ecmwf2roms.m

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

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

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

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

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

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

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

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

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

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

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

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

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

Re: Very early blowing up problem

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

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

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

Re: Parallel IO not scaling over Serial IO for ROMS

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Re: problems with the addition of tide forcing

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

Re: problems with the addition of tide forcing

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

Re: problems with the addition of tide forcing

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

Re: problems with the addition of tide forcing

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

Re: Problem with setting up inlet_test

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

Re: Error while compiling ROMS

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

Re: Error while compiling ROMS

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

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

Re: Huon_eastward and Hvom_northward

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

Re: Tracers don't radiate out at boundary

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

Re: Huon_eastward and Hvom_northward

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

Re: Shchepetkin (Shc) boundary condition reference?

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

Re: Tracers don't radiate out at boundary

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

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

Re: Finished compiling COAWST Installation - Ubuntu 14.04

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

Re: how to improve vertical profile velocity

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

Re: Compile error

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

Re: how to improve vertical profile velocity

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

Re: Asking for cpp and control files

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

Re: Amplification in zeta after adding river forcing

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

Re: PYROMS error

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

Re: Some confusion about boundary condition

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

Re: Some confusion about boundary condition

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

Re: Some questions about river input

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

Re: Some questions about river input

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

Re: WET_DRY option results in too much land_mask?

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

Re: erro about MPDADA

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

Re: open boundary for future projection of CMIP6 downscaling

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

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

Re: open boundary for future projection of CMIP6 downscaling

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

Re: open boundary for future projection of CMIP6 downscaling

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

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

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

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

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

Re: DIAG speed trouble / gridbuilder

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

Re: Error during compilation of build file

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

Re: ROMS processing water sluice in upstream

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

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

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

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

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

Re: The xi_position and eta_position of rivers.nc input

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

Re: surface boundary layer depth

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

Re: Interannual simulation

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

Re: Interannual simulation

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

Re: problem with pyroms.grid.Gridgen

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

Re: i have a question

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

Re: Interannual simulation

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

Re: Interannual simulation

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

Re: nudging/relaxation time scale field

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

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

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

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

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

Re: blow up the model

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

Re: Shortwave and Longwave radiation.

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

For the shortwave, it might be worth trying what you say. I use an albedo option because the ice model wants to have its own albedo.
by kate
Thu Aug 22, 2019 4:01 pm
Forum: ROMS Discussion
Topic: blow up the model
Replies: 3
Views: 442

Re: blow up the model

You need to investigate some more. Have you tried a shorter timestep? Have you looked at where the thing died? Is it surface? Boundary? Middle or what?
by kate
Fri Aug 16, 2019 3:53 pm
Forum: ROMS Discussion
Topic: Abnormal Cooling at Shallow Water Area
Replies: 7
Views: 1025

Re: Abnormal Cooling at Shallow Water Area

That's better, but still, what are your BODYFORCE parameters? I haven't used that option since DAMEE, which was a long time ago. It applies the surface stress not just at the surface, but over some number of points in the vertical. It was for a time when we couldn't resolve the surface boundary laye...
by kate
Fri Aug 16, 2019 5:04 am
Forum: ROMS Discussion
Topic: Abnormal Cooling at Shallow Water Area
Replies: 7
Views: 1025

Re: Abnormal Cooling at Shallow Water Area

Code: Select all

#undef  TS_C4VADVECTION
#define BODYFORCE
These and others of your choices put you into a realm where I have no experience at all. Have you done any bathymetry smoothing?
by kate
Mon Aug 12, 2019 6:42 pm
Forum: ROMS Sediment
Topic: Calculating sediment flux across transect
Replies: 4
Views: 1065

Re: Calculating sediment flux across transect

Are you aware of how the whole grid looks? See for example the figure here. To get something in the centers, try something like: Huon_sand_01 = 0.5 * (Huon_sand_01[:,:,1:-1,1:] + Huon_sand_01[:,:,1:-1,:-1]) # u ==> rho Hvom_sand_01 = 0.5 * (Hvom_sand_01[:,:,1:,1:-1] + Hvom_sand_01[:,:,:-1,1:-1]) # v...
by kate
Sun Aug 04, 2019 4:23 pm
Forum: ROMS Problems
Topic: Error climatology file: not include initialization time
Replies: 6
Views: 1342

Re: Error climatology file: not include initialization time

Are you starting from a restart file? What is the time in that file? I have a Python script for changing the time(s) in a file because sometimes you just need to.
by kate
Sat Aug 03, 2019 3:15 pm
Forum: ROMS Installation
Topic: upwelling: - Configuration error ..... exit_flag: 5
Replies: 3
Views: 645

Re: upwelling: - Configuration error ..... exit_flag: 5

Here is the real error message:
CHECKDEFS - no bottom stress formulation is selected.
The upwelling case should have one selected in the roms_upwelling.in file.
by kate
Fri Aug 02, 2019 3:09 pm
Forum: ROMS Sediment
Topic: ana_sediment.h:27:0: fatal error: tile.h: No such file
Replies: 8
Views: 775

Re: ana_sediment.h:27:0: fatal error: tile.h: No such file

I don't know if you need a wave model, but if you do, you should contact John Warner to get his version of the wave-coupling code. You might well need a wave model in that region to get the fresh water influx to mix down in a physical way. You will need to think about how to get the river inputs int...
by kate
Fri Aug 02, 2019 5:50 am
Forum: ROMS Sediment
Topic: ana_sediment.h:27:0: fatal error: tile.h: No such file
Replies: 8
Views: 775

Re: ana_sediment.h:27:0: fatal error: tile.h: No such file

Have you tried running a really simple case? Turning on sediments does not explain why you are failing to compile ROMS. What exactly are you doing? Operating system? ROMS version? Build script choices?
by kate
Thu Aug 01, 2019 6:33 pm
Forum: ROMS Sediment
Topic: ana_sediment.h:27:0: fatal error: tile.h: No such file
Replies: 8
Views: 775

Re: ana_sediment.h:27:0: fatal error: tile.h: No such file

You say it's running - was it running before you turned on the sediments? What you show is a failure to compile. I have a ROMS/Include/tile.h and so should you. During the build, it will look for include files in ROMS/Include, thanks to: includes := ROMS/Include vpath %.h $(includes) in the makefile...
by kate
Mon Jul 08, 2019 4:56 pm
Forum: ROMS Installation
Topic: Problem with the file makefile
Replies: 11
Views: 1546

Re: Problem with the file makefile

There was a recent update from Hernan in which the makefile was changed. It used to make "oceanM" and so forth, now it makes "romsM" and family. The final letter is "S" for serial, "M" for parallel, and "G" for debug mode (serial or parallel).
by kate
Sun Jul 07, 2019 6:57 am
Forum: ROMS Source
Topic: passive tracers - segmentation error ......help? :(
Replies: 6
Views: 864

Re: passive tracers - segmentation error ......help? :(

It is possible for problems to show up in debug mode that are not fatal when optimized. I still like to get to the bottom of all such issues, if I can.
by kate
Fri Jul 05, 2019 9:02 pm
Forum: ROMS Source
Topic: passive tracers - segmentation error ......help? :(
Replies: 6
Views: 864

Re: passive tracers - segmentation error ......help? :(

Did you change anything in your build except USE_DEBUG? I assume you cleaned out the old build first. Anyway, it's a link error, not a compile error. The netCDF library wants to link to the curl library. Sometimes, "nc-config --flibs" will want to bring in the curl library. You can add it explicitly...
by kate
Fri Jul 05, 2019 7:45 pm
Forum: ROMS Source
Topic: passive tracers - segmentation error ......help? :(
Replies: 6
Views: 864

Re: passive tracers - segmentation error ......help? :(

If you recompile with USE_DEBUG=on, you should get filenames and linenumbers instead of those annoying UNKNOWNS: oceanS 000000000059C85A Unknown Unknown Unknown oceanS 000000000059BF55 Unknown Unknown Unknown Without line numbers, we're pretty blind!
by kate
Wed Jul 03, 2019 4:31 am
Forum: ROMS Discussion
Topic: Marchesiello's RSUP3 in ROMS?
Replies: 9
Views: 2421

Re: Marchesiello's RSUP3 in ROMS?

You might try it first in the UPWELLING problem to see if it works when there is no land mask. If so, then try one of the simpler, smaller domains with a land mask. Build your own small tests, with and without land. If it is some trouble with the land mask, go through the equations in the paper and ...
by kate
Tue Jul 02, 2019 7:09 pm
Forum: ROMS Discussion
Topic: Marchesiello's RSUP3 in ROMS?
Replies: 9
Views: 2421

Re: Marchesiello's RSUP3 in ROMS?

Have you tried the split advection scheme on TS only? For your ellipse a, you may have boundary water too close to the land mask. Terms in the OBC schemes can include the first and second points in from the boundary and if the latter is land, it is not well-behaved. Your ellipses b and d have one-gr...
by kate
Tue Jul 02, 2019 6:38 pm
Forum: ROMS Discussion
Topic: SST less than -10 ℃
Replies: 1
Views: 391

Re: SST less than -10 ℃

Have you looked at the boundary file? Is that where you have inflow going on? It does look like a boundary problem to me.
by kate
Tue Jul 02, 2019 2:19 am
Forum: ROMS Discussion
Topic: Marchesiello's RSUP3 in ROMS?
Replies: 9
Views: 2421

Re: Marchesiello's RSUP3 in ROMS?

Have you looked to see where and how it is blowing up? Is it going for three steps no matter what the timestep?
by kate
Mon Jul 01, 2019 5:43 pm
Forum: ROMS Discussion
Topic: Marchesiello's RSUP3 in ROMS?
Replies: 9
Views: 2421

Re: Marchesiello's RSUP3 in ROMS?

The easiest thing to try is a shorter timestep. Have you tried that?
by kate
Sun Jun 30, 2019 1:12 am
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

I like to use ncview to get a quick and dirty view of what's in a netCDF file. Looking at 'ncdump -h' on a file can also be useful. If you didn't get temp and salt, then something didn't go right in the boundary extraction. If you like the idea of using these Python scripts, perhaps it's time to lea...
by kate
Sun Jun 30, 2019 12:25 am
Forum: ROMS Installation
Topic: Problem with the file makefile
Replies: 11
Views: 1546

Re: Problem with the file makefile

Did you create the out directory for it?
by kate
Fri Jun 28, 2019 5:52 pm
Forum: ROMS Discussion
Topic: Alongshore currents with periodic alongshore boundaries
Replies: 9
Views: 1363

Re: Alongshore currents with periodic alongshore boundaries

To answer here or there? You don't need to ask things twice on the ROMS forum. With a rigid-lid model, you can impose the longshore transport via the streamfunction boundary conditions. We no longer have a rigid-lid model, though. If you have too much flow over your shallow areas, perhaps you can ma...
by kate
Fri Jun 28, 2019 2:50 pm
Forum: ROMS Installation
Topic: Problem with the file makefile
Replies: 11
Views: 1546

Re: Problem with the file makefile

If you have a romsS, it is a serial executable. Run it with:

Code: Select all

romsS < roms_upwelling.in
The way you had it, it was waiting for you to type in the whole of roms_upwelling to standard-in (stdin).
by kate
Fri Jun 28, 2019 2:47 pm
Forum: ROMS Bugs
Topic: I/O syntax error
Replies: 2
Views: 619

Re: I/O syntax error

If you recompile with USE_DEBUG, then rerun this, it should give you file names and line numbers rather than all those "Unknown"s. Note that you don't want to do a long run with USE_DEBUG but it can be so helpful for solving these problems.
by kate
Thu Jun 27, 2019 11:26 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

What you want is to call do_file for each file in the list. Did it print the list of files to operate on? Is it what you expected? Then you can comment out the chunk with processes = 4 to the end. Add a loop going through the list of files, calling do_file(file, src_grd, dst_grd), something like: fo...
by kate
Thu Jun 27, 2019 5:46 pm
Forum: ROMS Problems
Topic: ROMS different src version leading to unalike results
Replies: 7
Views: 1277

Re: ROMS different src version leading to unalike results

You can hunt this down by using the equivalent of "git bisection", where you check out the version in the middle and see which camp it falls into. Pick a simulation of only a day or so to compare to. You can iteratively try intermediate versions until you find out which commit changed your results, ...
by kate
Thu Jun 27, 2019 5:41 pm
Forum: ROMS Discussion
Topic: Nudging time scales working on bry or clim
Replies: 7
Views: 1606

Re: Nudging time scales working on bry or clim

First, what is the land mask doing right there? You might want to change:

______
.........|_____

to

________________

filling in the leftmost water points with dirt. That's what I would try first. Next, look to see what the velocities are doing in that corner.
by kate
Wed Jun 26, 2019 9:11 pm
Forum: ROMS Installation
Topic: Problem with the file makefile
Replies: 11
Views: 1546

Re: Problem with the file makefile

Now you have to tell it where to find the netCDF libraries at link time, both via the makefile and also via LD_LIBRARY_PATH or whatever the equivalent is on your OS (I live in a Linux world now).
by kate
Wed Jun 26, 2019 4:39 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Maybe this is biting you:

Code: Select all

memorylocked    16384 kbytes
by kate
Wed Jun 26, 2019 12:36 am
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

You can read up on multiprocessing online. Something wasn't right about your last attempt. I would check your system limits. With a bash shell: (snowdrifts) [kshedstrom@chinook02 Tidal_bay]$ ulimit -a core file size (blocks, -c) 0 data seg size (kbytes, -d) unlimited scheduling priority (-e) 0 file ...
by kate
Tue Jun 25, 2019 11:38 pm
Forum: ROMS Installation
Topic: Problem with the file makefile
Replies: 11
Views: 1546

Re: Problem with the file makefile

cp -f /usr/include/netcdf.mod /home/al/Projects/upwelling/Build_romsG cp -f /usr/include/typesizes.mod /home/al/Projects/upwelling/Build_romsG It is looking for the netcdf.mod file which got copied into the Build directory. You know those files got copied, but what compiler was used to make them? I...
by kate
Tue Jun 25, 2019 4:41 pm
Forum: ROMS Discussion
Topic: ROMS 1D for Mac needed
Replies: 2
Views: 492

Re: ROMS 1D for Mac needed

The way to run in 1D is to set up a very small horizontal grid of say 3x3 or 4x4 with doubly periodic boundary conditions. The grid can be analytic and the boundary conditions don't need any inputs. Apply the same surface forcing to all points and you should be good to go. It should run on a Mac, no...
by kate
Mon Jun 24, 2019 6:07 pm
Forum: ROMS Problems
Topic: one question about boundary condition in an idealized test
Replies: 2
Views: 578

Re: one question about boundary condition in an idealized te

Dynamically, the zeta outside is not playing a direct role, though it can be used in Flather type open boundary condition. You don't say what open boundary options you have chosen. In the model zeta is a second-order field in that it is computed from lateral convergences and divergences - the barotr...
by kate
Mon Jun 24, 2019 4:17 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Can you try on a system with more memory?
by kate
Mon Jun 24, 2019 6:34 am
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Does it object to processes = 1 the same way as processes = 4, or was there a different error? What if you don't let it kill the job?
by kate
Mon Jun 24, 2019 5:18 am
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Maybe you should try a pool size of 1, not running in parallel. I'm not familiar with that exact display, but the swap seems high, like it doesn't all fit in memory and none of them are making any progress.
by kate
Sun Jun 23, 2019 7:18 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Each system is different. I run on one with front-end nodes containing 8 cores, while the compute nodes contain 24-28 cores. I run these scripts on the front-end nodes because they have access to the filesystems containing these files, but I only ask for 4 cores in my "pool" so other people can get ...
by kate
Sun Jun 23, 2019 4:53 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

year = int(sys.argv[1]) Did you give it a year as a command-line argument? It is expecting one. There's also a commented out line: #lst_year = sys.argv[1:] that can be enabled, allowing you to give it a list of years. Or you can edit the script to provide the year(s) internally. For the other, I th...
by kate
Sat Jun 22, 2019 3:40 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

I can give you the PALAU1 chunk of gridid.txt, but then it points to a grid file which you don't have either. Why don't you put in a chunk for your grid, since that is what you want anyway? id = PALAU1 name = PALAU1 grdfile = /import/AKWATERS/kshedstrom/gridpak/Palau/grid_Palau_1.nc N = 50 grdtype =...
by kate
Sat Jun 22, 2019 4:03 am
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Sometimes Python evolves so that something that used to run no longer runs. Sometimes when that happens I can figure out how to change it and get it pushed to github. You need to update your code through either "git pull" or just editing it in place.
by kate
Sat Jun 22, 2019 3:35 am
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

My latest version of that troublesome line is: h[j,i] = depth_bnds[int(bottom[j,i])+1] Sometimes it pays to "git pull". Step 5 happens with this script: #!/bin/ksh print -n "Enter year to process: "; read year; print "" leap=`echo $(($year % 4))` if [ $leap == 0 ] ; then nday=366 else nday=365 fi #n...
by kate
Wed Jun 19, 2019 5:10 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Are you sure you are downloading on the GLBa0.08 HYCOM grid? They have moved that to the "inactive" status. It was on a tripole grid and needed lat,lon to be functions if both i and j. The newer HYCOM products are on a simpler grid, where lat is a function of j only and lon is a function of i only. ...
by kate
Tue Jun 18, 2019 7:13 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

No, I don't think you have a version problem, unless it's a HYCOM version problem.

Code: Select all

lon_t = lon[1:-1,1:-1]
Longitude is often a 1-D array, not a 2-D array. What do you have?
by kate
Tue Jun 18, 2019 5:20 am
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

Looks like it should be 'depth', not 'z'.
by kate
Mon Jun 17, 2019 6:25 pm
Forum: ROMS Tools and Techniques
Topic: Create boundary and initial file from seagrid and GLBa0.08
Replies: 32
Views: 3023

Re: Create boundary and initial file from seagrid and GLBa0.

I have successfully made ROMS IC and BC files from HYCOM using the scripts in pyroms/examples/Arctic_HYCOM/. First you run make_remap_weights_file, then make_ic_file/make_bdry_file. A more equatorial example is in Palau_HYCOM. For the vertical, set the number of vertical points in the gridid.txt fil...
by kate
Sat Jun 15, 2019 9:59 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS error
Replies: 12
Views: 1775

Re: PYROMS error

The only use I know of for lpsolve is to smooth the bathymetry. I usually use pyroms for other tasks, not grid generation. What do you need pyroms to do for you?
by kate
Sat Jun 15, 2019 4:22 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS error
Replies: 12
Views: 1775

Re: PYROMS error

Spectral has a capital "S". See for example this Stackoverflow thread. I don't know how you ended up with a lower case spectral.

As for the other, my solution is that I don't use lpsolve. I think I once copied the thing to my current directory and that worked.
by kate
Fri Jun 14, 2019 5:44 pm
Forum: ROMS Discussion
Topic: Floats in nested configuration
Replies: 3
Views: 578

Re: Floats in nested configuration

I assume that would work. Why don't you try it and let us know?
by kate
Fri Jun 14, 2019 5:29 pm
Forum: ROMS Installation
Topic: Tutorial: Installing ROMS in Windows XP / Vista (via Cygwin)
Replies: 20
Views: 13305

Re: Tutorial: Installing ROMS in Windows XP / Vista (via Cyg

I have tried to provide all you need to know (and more!) for Unix in the manual.
by kate
Fri Jun 14, 2019 5:26 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS error
Replies: 12
Views: 1775

Re: PYROMS error

This troublesome bit: -Inf-config not yet implemented for cmake builds/include is coming from: INCDIR = $(shell $(NC_CONFIG) --prefix)/include which as you point out is from nf-config. It is executing "nf-config --prefix" which on my machine gives: chinook04.rcs.alaska.edu 543% nf-config --prefix /h...
by kate
Thu Jun 13, 2019 4:49 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS error
Replies: 12
Views: 1775

Re: PYROMS error

No, don't go backwards! We always had problems with python finding those libraries, so going backwards is not the answer. Stick with Python 3 and read the README where I tell you what I have to do for python to find the libraries - I just have copies in three different places, with hopes that at lea...
by kate
Thu Jun 13, 2019 4:46 pm
Forum: ROMS Discussion
Topic: Floats in nested configuration
Replies: 3
Views: 578

Re: Floats in nested configuration

No, it has not. I think you'd be better off doing floats offline with dedicated float-tracking software instead. Adding multi-grid features would still be tricky there too.
by kate
Thu Jun 13, 2019 4:35 am
Forum: ROMS Tools and Techniques
Topic: PYROMS error
Replies: 12
Views: 1775

Re: PYROMS error

Type "git branch" to see which branch you have. The "python3" branch is the current default, so that is probably what you have already.

As for those libraries, you need to make sure your LD_LIBRARY_PATH environment variable includes the directory they live in.
by kate
Wed Jun 12, 2019 6:45 am
Forum: ROMS Discussion
Topic: pan-Arctic simulations
Replies: 7
Views: 2125

Re: pan-Arctic simulations

You want to be working in a polar stereographic projection so that there is no pole problem. The domain is roughly a rectangle over the pole. The curvilinear orthogonal grids are made in a flat geometry, so you need software that can do map projections and inverse map projections to get back to the ...
by kate
Wed Jun 12, 2019 6:19 am
Forum: ROMS Discussion
Topic: pan-Arctic simulations
Replies: 7
Views: 2125

Re: pan-Arctic simulations

What do you want to know? I'm running a pan-Arctic domain now, which is described in the ROMS manual.
by kate
Tue Jun 11, 2019 6:38 pm
Forum: ROMS Discussion
Topic: Time type in nc file
Replies: 5
Views: 754

Re: Time type in nc file

It doesn't matter. Just be consistent across all your inputs.
by kate
Tue Jun 11, 2019 6:35 pm
Forum: ROMS Discussion
Topic: unrealistic discontinuity of surface salinity
Replies: 2
Views: 495

Re: unrealistic discontinuity of surface salinity

I would try LMD_SHAPIRO as a first step.
by kate
Mon Jun 10, 2019 3:34 pm
Forum: ROMS Discussion
Topic: Slow running of nesting simulation
Replies: 5
Views: 638

Re: Slow running of nesting simulation

Hernan once said that he thought the way to do it is fine-to-coarse, with a fine grid the size of your larger domain. However, as you find, one must use coarse-to-fine because that's the matlab script you have. I made all the grids and did the bathymetry the way I wanted to. What I can help you with...
by kate
Sat Jun 08, 2019 5:48 am
Forum: ROMS Bugs
Topic: WC13/R4DVAR Compiling Problem: rp_set_vbc.f90
Replies: 1
Views: 421

Re: WC13/R4DVAR Compiling Problem: rp_set_vbc.f90

You need to fix rp_set_vbc.F, not the .f90 file.
by kate
Fri Jun 07, 2019 9:35 pm
Forum: ROMS Discussion
Topic: Slow running of nesting simulation
Replies: 5
Views: 638

Re: Slow running of nesting simulation

Yes, the online nesting is extremely slow. I have not seen any ROMS development in the myroms code to lead me to expect that this has changed since those old posts. I am instead doing offline nesting because it gives more flexibility in the child grids.
by kate
Fri Jun 07, 2019 9:29 pm
Forum: ROMS Discussion
Topic: bulk flux advice from clivar community
Replies: 3
Views: 1117

Re: bulk flux advice from clivar community

The latest ROMS has the WIND_MINUS_CURRENT, which is what I was talking about. Actually, the CLIVAR people now say that one should not use WIND_MINUS_CURRENT with the new JRA55-do forcings because of how they used the satellite-derived winds in the data assimilation. To get realistic variance, I wou...
by kate
Thu Jun 06, 2019 10:36 pm
Forum: ROMS Bugs
Topic: Error when i try to run the upwelling.in
Replies: 2
Views: 508

Re: Error when i try to run the upwelling.in

We need more information. Are you asking ROMS for USE_NETCDF4? When you say you downloaded netcdf, what exactly do you mean by that? You can download the source code and compile it or you can find a way to get it prebuilt for some systems. So what operating system is this anyway? This is a problem w...
by kate
Wed Jun 05, 2019 5:34 pm
Forum: ROMS Problems
Topic: Segmentation fault after 10 model years of running
Replies: 5
Views: 878

Re: Segmentation fault after 10 model years of running

I don't know how to fix this. I have always used restarts because I operate on a supercomputer with queue length restrictions. You might want to be using PERFECT_RESTART.
by kate
Tue Jun 04, 2019 10:00 pm
Forum: ROMS Problems
Topic: Segmentation fault after 10 model years of running
Replies: 5
Views: 878

Re: Segmentation fault after 10 model years of running

It's possible that there's a memory leak, so that the program grows over time until it hits some limit. I just ran something in debug mode and got a *ton* of these messages: ==6434==ERROR: LeakSanitizer: detected memory leaks Direct leak of 34816 byte(s) in 8 object(s) allocated from: #0 0x7f1d00486...
by kate
Tue Jun 04, 2019 6:01 pm
Forum: ROMS Tools and Techniques
Topic: History/Averages file output for a single sigma layer
Replies: 2
Views: 511

Re: History/Averages file output for a single sigma layer

You should check the Qout values in roms.in and turn on the Quick output feature. It already has the Usur, Vsur, and surT for surface fields.
by kate
Tue Jun 04, 2019 5:41 am
Forum: ROMS Problems
Topic: NESTING RUNTIME ERROR code 59
Replies: 5
Views: 808

Re: NESTING RUNTIME ERROR code 59

The oceanG with USE_DEBUG can certainly be the one using MPI. You can't tell from the name of the executable, just from how you built it. It should be possible to run it in an MPI environment.
by kate
Mon Jun 03, 2019 10:45 pm
Forum: ROMS Problems
Topic: NESTING RUNTIME ERROR code 59
Replies: 5
Views: 808

Re: NESTING RUNTIME ERROR code 59

The generic advice to figure out what's going on is: "recompile with USE_DEBUG, rerun to get the line numbers, look at the lines referred to in the output". This is the preferred method. The other thing people do is post their roms.in file to this list and hope that someone will look at it and tell ...
by kate
Mon Jun 03, 2019 10:41 pm
Forum: ROMS Problems
Topic: How to adjust roms temperature with observation (Bias)
Replies: 2
Views: 577

Re: How to adjust roms temperature with observation (Bias)

You've told us about your vertical and horizontal diffusion, but what is the vertical boundary condition? In other words, the surface heat flux? Can you output it and look at it? Is it computed internally (bulk fluxes) or provided externally?
by kate
Wed May 29, 2019 3:27 pm
Forum: ROMS Problems
Topic: Noisy SST
Replies: 9
Views: 1134

Re: Noisy SST

As for MIX_S_TS, it's a really bad idea over a sloping bottom and can lead to unphysical cross-isopycnal mixing. I hear you! But there was no way to solve the issue of punctual alternating high/low temperatures at the slope when mixing along z-planes which would even cause the model to blow up. The...
by kate
Tue May 28, 2019 4:10 pm
Forum: ROMS Discussion
Topic: Extracting temperature
Replies: 9
Views: 906

Re: Extracting temperature

The ROMS default is to write out the grid file into every output netcdf file. I make many, many output netcdf files and have told ROMS not to put the grid into any of them. I know the grid because it is one of the ROMS input files. If you are using someone else's model output, you have to ask them f...
by kate
Tue May 28, 2019 4:07 pm
Forum: ROMS Problems
Topic: Noisy SST
Replies: 9
Views: 1134

Re: Noisy SST

The equivalent horizontal smoothing for LMD is LMD_SHAPIRO, to horizontally filter the mixed layer depth. I've heard from some other colleagues that they too need to smooth the mixed layer depth to avoid that sort of noise. As for MIX_S_TS, it's a really bad idea over a sloping bottom and can lead t...
by kate
Tue May 28, 2019 4:03 pm
Forum: ROMS Problems
Topic: How to adjust temperature mix, salinity advection and overall vertical mixing in roms
Replies: 3
Views: 723

Re: How to enhance vertical mixing in roms MY2.5

I've heard it's a difficult problem. We also have too little vertical mixing from all the fresh water input into the Gulf of Alaska. One idea is that we need to couple to a wave model and enable the wind-wave surface mixing in COAWST. Easier said than done; we haven't tried it yet.
by kate
Tue May 28, 2019 4:00 pm
Forum: ROMS Problems
Topic: ROMS restart time steps
Replies: 4
Views: 809

Re: ROMS restart time steps

That's weird, that's not what I get. Is there something about your forcing that ends before the end of the year? Do you see the NaNs in the ROMS output for the last few steps?
by kate
Tue May 28, 2019 3:56 pm
Forum: ROMS Discussion
Topic: Extracting temperature
Replies: 9
Views: 906

Re: Extracting temperature

I think you are saying that the index is xi=91, eta=28. Yes, you can extract the temperature at that location to get a vertical column of numbers for each time.
by kate
Mon May 27, 2019 4:13 am
Forum: ROMS Information
Topic: Error compiling COAWST
Replies: 4
Views: 732

Re: Error compiling COAWST

You need to change yours to look like mine for it to compile.
by kate
Sun May 26, 2019 4:52 am
Forum: ROMS Discussion
Topic: ad passive tracer
Replies: 10
Views: 1337

Re: ad passive tracer

Are you by chance asking for Aout(idTTav)? I have never done so and can't guarantee that it would work.
by kate
Sat May 25, 2019 7:09 am
Forum: ROMS Discussion
Topic: Extracting temperature
Replies: 9
Views: 906

Re: Extracting temperature

There should be a grid file somewhere with lat and lon as a function of xi and eta. The algorithm for finding the xi, eta (or i,j) of a given lat, lon is implemented inside ROMS in the routine hindices inside interpolate.F.
by kate
Sat May 25, 2019 7:03 am
Forum: ROMS Discussion
Topic: ad passive tracer
Replies: 10
Views: 1337

Re: ad passive tracer

Which mod_ncparam.f90 error are you getting and what is on the line number in question? We want you to look at the file. Maybe you should also explain exactly what it is you want to do with your passive tracer.
by kate
Fri May 24, 2019 4:17 pm
Forum: ROMS Discussion
Topic: Extracting temperature
Replies: 9
Views: 906

Re: Extracting temperature

If you only have a few such location from which you want to extract fields, the easiest thing is to ask ROMS for a STATIONS file. You provide the lat,lon in the stations input file. If you don't have a chance to rerun the model, you will have to build tools for extracting from the ROMS output, using...
by kate
Fri May 24, 2019 4:13 pm
Forum: ROMS Information
Topic: Error compiling COAWST
Replies: 4
Views: 732

Re: Error compiling COAWST

Could you look at your diag.F to see how Nreduce is set and used? My code has: #ifdef DISTRIBUTE # ifdef SOLVE3D # ifdef BIOLOGY integer :: Nreduce # else integer, parameter :: Nreduce = 6 # endif integer, parameter :: Ncourant = 7 # else integer, parameter :: Nreduce = 4 integer, parameter :: Ncour...
by kate
Thu May 23, 2019 6:40 pm
Forum: ROMS Problems
Topic: Blowing up problem
Replies: 9
Views: 1083

Re: Blowing up problem

Minimum Z-grid spacing, DZmin = -2.17046894E-01 m This would worry me. It means you need to be using WET_DRY or make the minimum depth deeper. 1)shorten my DT from 10 seconds to 1 second. The standard output is shown below. The maximum barotropic Courant Number is smaller than 1. Minimum barotropic...
by kate
Thu May 23, 2019 4:52 pm
Forum: ROMS Discussion
Topic: ad passive tracer
Replies: 10
Views: 1337

Re: ad passive tracer

You've got one clue here:

Code: Select all

At line 2008 of file mod_ncparam.f90
Can you tell us what is at line 2008 of your mod_ncparam.f90 file?

What is NRREC in your roms.in file? The restart file is not necessarily the same as a file needed for initialization, especially when you go changing settings.
by kate
Thu May 23, 2019 4:03 am
Forum: ROMS Discussion
Topic: low velocity in deep water
Replies: 1
Views: 352

Re: low velocity in deep water

We need a good deal more information about what you are doing! By the way, you should be looking at the T and S water mass properties as well. There's a known problem on the basin scale with at least most advection schemes being too dissipative, especially over multi-decadal timescales.
by kate
Wed May 22, 2019 9:09 pm
Forum: ROMS Discussion
Topic: Nudging time scales working on bry or clim
Replies: 7
Views: 1606

Re: Nudging time scales working on bry or clim

Your values look like you should use 'second-1' as the units. 'day-1' is also a valid option. See get_nudgcoef.F to check that. You say you want 1 day at the edge and one year in the middle. I think if it rather as strong (1 dayish) at the edge to infinite time (zero coefficient) over some width, sa...
by kate
Mon May 20, 2019 3:55 pm
Forum: ROMS Discussion
Topic: Lateral nudging and nudging to climatology
Replies: 2
Views: 598

Re: Lateral nudging and nudging to climatology

This was discussed just days ago right here.