Search found 175 matches

by stef
Fri Mar 15, 2024 4:05 pm
Forum: ROMS Discussion
Topic: 4dvar run hangs using internal lapack module
Replies: 4
Views: 134

Re: 4dvar run hangs using internal lapack module

Oh, I just saw the comment in lapack_mod.F. That explains it, so it seems more a question of code style and easier future maintenance. ! This module includes modernized versions of the selected routines ! ! from thw Linear Algebra Package (LAPACK) library, which are used ! ! in ROMS 4D-Var algorithm...
by stef
Fri Mar 15, 2024 3:48 pm
Forum: ROMS Discussion
Topic: 4dvar run hangs using internal lapack module
Replies: 4
Views: 134

Re: 4dvar run hangs using internal lapack module

Good to know, it was not clear from the trac ticket (at least for me) that NCEP was complaining about the GOTO's. Do you mean their code doesn't compile due to the GOTOs?
by stef
Fri Mar 15, 2024 10:04 am
Forum: ROMS Discussion
Topic: 4dvar run hangs using internal lapack module
Replies: 4
Views: 134

Re: 4dvar run hangs using internal lapack module

I should add that in the successful run, the next lines are (001,007): New Ritz eigenvalues and their accuracy, RitzMaxErr = 1.00000E-01 001 3.6293220E+00 3.6446995E-03 converged (Good=001) 002 8.5684099E+00 9.1391967E-03 converged (Good=002) 003 2.1050775E+01 2.1416478E-03 converged (Good=003) 004 ...
by stef
Fri Mar 15, 2024 9:58 am
Forum: ROMS Discussion
Topic: 4dvar run hangs using internal lapack module
Replies: 4
Views: 134

4dvar run hangs using internal lapack module

I have a toy test case using RBL4DVAR & RPCG with 1 outer loop and 7 inner loops. The run hangs at inner loop 6 at the computation of the Lanczos vectors. This occurs when using the 'develop' branch at git commit b7a47408ba22 ([1]), which uses the internal lapack routines in lapack_mod.F ([2]). ...
by stef
Wed Feb 14, 2024 3:36 pm
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Re: Adjoint blow-up if misfit forcing is too far after initial time

Oh, I did not see the roms_test repository on github! On the 4dvar tutorial page still mentions svn: https://www.myroms.org/wiki/4DVar_Tutorial_Introduction (accessed 2024/02/14) I might have missed some announcements... As long as you keep both repositories (roms and roms_test) updated the tutorial...
by stef
Wed Feb 14, 2024 7:27 am
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Re: Adjoint blow-up if misfit forcing is too far after initial time

Yes, use the ROMS version on GitHub. We will no longer update all the SVN versions after January 1, 2025. I fixed all the indices in the code. Hmm, not sure if I understand. Which ROMS version on github, the 'develop' branch or the 'main' branch? Which revision should we use with the currently avai...
by stef
Tue Feb 13, 2024 5:09 pm
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Re: Adjoint blow-up if misfit forcing is too far after initial time

Sorry may I ask a question about the Tutorial in this thread (it may be related to my problem discussed here): I noticed that the tutorial does not seem to work with git branch 'main'. Currently 'main' is at 21d8ab47ad4b (from June 2023). The debugger complains about what looks like related to the b...
by stef
Tue Feb 13, 2024 8:01 am
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Re: Adjoint blow-up if misfit forcing is too far after initial time

But I suggest staying in the strong-constraint regime since modeling the model error, Q, is not trivial. Yes, I'm trying to look at the simplest case possible. NADJ should be the same size as ntimes. That's what I thought after I first tried to change it. Because if the input is lower than that, it...
by stef
Mon Feb 12, 2024 5:32 pm
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Re: Adjoint blow-up if misfit forcing is too far after initial time

Thanks for the reply! I used RBL4DVAR with RPCG, as in Exercise 3 of the tutorial. It worked in some cases, but then I got blow-ups and didn't now how to investigate the causes. That's when I switched to I4DVAR and modified the NADJ to write out the adjoint solutions more frequently (I'm aware that ...
by stef
Mon Feb 12, 2024 3:41 pm
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Re: Adjoint blow-up if misfit forcing is too far after initial time

Hmm, regarding what I wrote above about preconditioning: In i4dvar.f90 (as noted in the first post I use I4DVAR for debugging), there is a segment ! Convert observation cost function gradient, GRADx(Jo), from model ! space (x-space) to minimization space (v-space) which is executed after the first a...
by stef
Mon Feb 12, 2024 10:27 am
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Re: Adjoint blow-up if misfit forcing is too far after initial time

Hi Ivica, Thanks for the suggestions! Yes I haven't provided enough details about the run, I will try to write it up in a blog post. My present set up is more of a test case for learning 4dvar than a realistic application, i.e. an addition to the tutorial excercises. Alongside I am reading the liter...
by stef
Thu Feb 08, 2024 9:45 am
Forum: ROMS Discussion
Topic: Adjoint blow-up if misfit forcing is too far after initial time
Replies: 14
Views: 1204

Adjoint blow-up if misfit forcing is too far after initial time

I am making simple experiments with 4dvar and sometimes I get blow-ups in the first iteration of the adjoint model (it doesn't actually terminate during the adjoint, but the output is clearly nonsensical). In one experiment I use: * ARGO profiles from a single float, with one profile every 3 days * ...
by stef
Wed Dec 06, 2023 3:50 pm
Forum: ROMS Bugs
Topic: Possible bug in obs_k2z.F
Replies: 0
Views: 13931

Possible bug in obs_k2z.F

I think there is a bug at the following line in obs_k2z.F: https://github.com/myroms/roms/blob/5ca8b9d0559233eaa40e7fc1d04607b5aa8d579e/ROMS/Utility/obs_k2z.F#L187 (src:trac:960, git branch develop at 5ca8b9d05592) Should this r2=REAL(k2-k1,r8)*(Zobs(iobs)-REAL(k1,r8)) be instead r2=REAL(k2-k1,r8)*(...
by stef
Thu Nov 23, 2023 7:46 am
Forum: ROMS Discussion
Topic: AMSR2 collated or uncollated?
Replies: 6
Views: 15532

Re: AMSR2 collated or uncollated?

Thanks for clarifying!

Makes sense to me, because you are probably resolving some diurnal signal in the top layer.

Thanks again for all the help, it is all much clearer now.
by stef
Tue Nov 21, 2023 3:24 pm
Forum: ROMS Discussion
Topic: AMSR2 collated or uncollated?
Replies: 6
Views: 15532

Re: AMSR2 collated or uncollated?

Thanks, great resource!! Hmm, I've skimmed through the paper and there is a sentence: PWR SST is consistent with in situ data during both day and night and thus may be viewed as an estimate of ‘depth’ in situ SST. The PWR SST is expected to be a useful input into L4 SST analyses, especially for foun...
by stef
Wed Nov 15, 2023 5:46 pm
Forum: ROMS Discussion
Topic: AMSR2 collated or uncollated?
Replies: 6
Views: 15532

Re: AMSR2 collated or uncollated?

Thanks for the information! Wow, that's cool! So it probably means you assimilate the ascending and descending passes individually, and maybe the subskin temperature directly as opposed to converting it to a foundation temperature.
by stef
Tue Nov 14, 2023 8:41 am
Forum: ROMS Discussion
Topic: AMSR2 collated or uncollated?
Replies: 6
Views: 15532

AMSR2 collated or uncollated?

Am a beginner in data assimilation and a question came up relating assimilation of microwave SST data, specifically from AMSR2. In Lecture 6 of the 4D-VAR tutorial, it says that the sampling frequency for AMSR2 and WindSat is daily (Table 21.1). According to [1], AMSR2 is near-polar orbiting and sam...
by stef
Tue Jul 11, 2023 5:06 pm
Forum: ROMS Adjoint
Topic: Steric signals and along-track altimetry
Replies: 3
Views: 4327

Re: Steric signals and along-track altimetry

After reading the manuals I guess this is incorrect, with a balance operator one would see both and adjustment of the surface, and an adjustment of dynamic/steric height via tracers. I couldn't test it so far using a balance operator (see below), but looking at the definition of the matrix it seems...
by stef
Tue Jul 11, 2023 1:30 pm
Forum: ROMS Adjoint
Topic: Steric signals and along-track altimetry
Replies: 3
Views: 4327

Re: Steric signals and along-track altimetry

Do I assume correctly that for gridded altimetry, the domain-averaged steric signals should have no effect on the tracer increments (even if not removed), because the model is Boussinesq? After reading the manuals I guess this is incorrect, with a balance operator one would see both and adjustment ...
by stef
Wed Jun 28, 2023 9:06 am
Forum: ROMS Adjoint
Topic: Steric signals and along-track altimetry
Replies: 3
Views: 4327

Re: Steric signals and along-track altimetry

I'm asking because I think in the BoB,the steric seasonal mean sea-level fluctuation can be like 5 cm amplitude. In [1] it seems to be much less, their Fig. 3 suggests roughly 4 cm range (2 cm amplitude) of the seasonal steric component. That's domain-averaged over the northern Bay of Bengal, compu...
by stef
Tue Jun 27, 2023 8:35 am
Forum: ROMS Adjoint
Topic: Steric signals and along-track altimetry
Replies: 3
Views: 4327

Steric signals and along-track altimetry

Hi, Is there a common practice of how to remove seasonal thermo/halosteric signals from along-track altimetry in data assimilation? Do these have to be removed from gridded altimetry? Do I assume correctly that for gridded altimetry, the domain-averaged steric signals should have no effect on the tr...
by stef
Fri May 12, 2023 8:30 am
Forum: ROMS Discussion
Topic: Bottom Friction and How to Slow Down Currents in Shallow Regions
Replies: 11
Views: 7602

Re: Bottom Friction and How to Slow Down Currents in Shallow Regions

Have you tried logarithmic drag? In [1], it says Because the vertical elevation of the velocity in the bottom computational cell will vary spatially and temporally, the inclusion of the elevation provides a more consistent formulation. For testing purposes, it may be possible to disable LIMIT_BSTRES...
by stef
Wed May 03, 2023 1:45 pm
Forum: ROMS Adjoint
Topic: Vertical coordinate for assimilation of CTDs
Replies: 1
Views: 2864

Re: Vertical coordinate for assimilation of CTDs

Looking at the code in extract_obs.F, what I wrote above "obs_depth" in ROMS is "depth of observations below the sea surface" is not true, is it? Do I understand correctly that the varinfo.yaml is out of date, because it says there: variable: obs_depth standard_name: depth_of_obs...
by stef
Thu Apr 27, 2023 2:27 pm
Forum: ROMS Adjoint
Topic: Vertical coordinate for assimilation of CTDs
Replies: 1
Views: 2864

Vertical coordinate for assimilation of CTDs

I'm trying out data assimilation and a question regarding in-situ profiles came up. Do we have to consider dynamic height anomalies when assimilating CTD profiles into ROMS, i.e. do you adjust the "depth" axis with each individual salinity and temperature profile from the cast? In the past...
by stef
Mon Mar 27, 2023 4:56 pm
Forum: ROMS Adjoint
Topic: Possible truncation bug in the binning of observations to baroclinic time steps
Replies: 0
Views: 3433

Possible truncation bug in the binning of observations to baroclinic time steps

I'm trying to learn 4dvar and noticed that some observations got rejected, even though the observation locations and time looked fine. Am I doing the time conversion for observations with the wrong numerical precision, or is there a minor bug in the binning of observations to baroclinic time steps c...
by stef
Thu Jan 26, 2023 7:52 am
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

It's been about a week since I asked at the rda.ucar helpdesk, no answer yet. They seem to re-organize their website, maybe the ticket got lost. It may also be possible that the helpdesk tickets I lodge have lower priority than others (due to my non-institutional email address, no affiliation with l...
by stef
Tue Jan 17, 2023 8:24 am
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

The individual files, e.g., https://rda.ucar.edu/thredds/dodsC/file ... grib2.html (one every 3 hours) have variable Downward_Short-Wave_Radiation_Flux_surface_3_Hour_Average Well, I guess not all of them. The +09 file you linked above https://rda.ucar.edu/thredds/dodsC/files/g/ds083.3/2021/202101/...
by stef
Mon Jan 16, 2023 4:58 pm
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

Or, since the folks at NCAR are usually super helpful, perhaps log a ticket with their helpdesk asking for clarification.
Ok, I'll log a ticket! I'll update this post once I get a response.
by stef
Mon Jan 16, 2023 4:54 pm
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

You are plotting values extracted from individual files Well I'm plotting the files that are accessed by the scripts you mention above, where you wrote roms_{get,write}_gfs_NCARds083_bulkflux.m for accessing the NCAR archive of GFS data to my tools https://github.com/johnwilkin/roms_wilkin. Specifi...
by stef
Mon Jan 16, 2023 9:47 am
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

and if you plot them together you get a nice diurnal cycle resolved at 3-hour intervals? So I'm comparing this with ERA5 data, and also what I would get from applying a (hypothetically necessary) conversion from the 6-hour averages to 3-hour averages of the RDA aggregation/Best fluxes. Below figure...
by stef
Fri Jan 13, 2023 5:13 pm
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

It's my understanding that variables denoted "average" are just that ... they are averages over the 3-hour interval as opposed to snapshots at a time instant. So the question is whether they are 3-hour averages, or mixed 3-hour averages and 6-hour averages, right? Earlier you wrote The co...
by stef
Thu Jan 12, 2023 6:19 pm
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

Wow, great! The codes return everything at 3-hour intervals because NCAR have unpacked those 6-hour accumulated variables before adding to the aggregation. That's interesting, may I ask how you found this? I didn't see it in the manuals. Do you ask at their helpdesk for this type of information? Hmm...
by stef
Wed Jan 11, 2023 9:16 am
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

Hmm, I guess one could easily compute the 3-hourly resolution from the available data: F_a: average flux from +00 to +03 (available from the resources) F_b: average flux from +03 to +06 (missing in the resources) F_c: average flux from +00 to +06 (available from the resources) Can't one compute F_b ...
by stef
Tue Jan 10, 2023 3:48 pm
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

Re: variable frequency in atmospheric forcing fields

Thanks for your reply! I have seen the ERA5 scripts, thanks! However, my question was specific to GDAS/FNL and GFS. In that repository you linked to, do you have scripts that extract hourly data from those as well? The question was whether 1) the radiative fluxes for GDAS/FNL and GFS really are (as ...
by stef
Mon Jan 09, 2023 6:10 pm
Forum: ROMS Discussion
Topic: variable frequency in atmospheric forcing fields
Replies: 15
Views: 20478

variable frequency in atmospheric forcing fields

Hi, I am comparing a simulations forced with ERA5, GDAS/FNL and GFS fields using bulk fluxes (I know it can be problematic [1]). ERA5 has radiation fluxes at the same time steps as wind, albeit temporally staggered because it's an averaged quantity, but effectively all fields have the same resolutio...
by stef
Mon Dec 05, 2022 4:35 pm
Forum: ROMS Discussion
Topic: Manning-Strickler for 2d bottom friction
Replies: 5
Views: 5793

Re: Manning-Strickler for 2d bottom friction

regarding the 'h' variable, in set_vbc.F there is # ifdef UV_QDRAG real(r8), intent(in) :: rdrag2(LBi:UBi,LBj:UBj) real(r8), intent(in) :: zeta(LBi:UBi,LBj:UBj,3) real(r8), intent(in) :: h(LBi:UBi,LBj:UBj) # endif and you could write # ifdef UV_QDRAG real(r8), intent(in) :: rdrag2(LBi:UBi,LBj:UBj) r...
by stef
Mon Dec 05, 2022 4:11 pm
Forum: ROMS Discussion
Topic: Manning-Strickler for 2d bottom friction
Replies: 5
Views: 5793

Re: Manning-Strickler for 2d bottom friction

I never actually used the code snippet above. I'm lazy, what I actually did is "abuse" the existing field "rdrag2" (which is actually the quadratic bottom drag coefficient) in set_vbc.F and initialized it with manning_n values. This is dangerous and can go very wrong, because the...
by stef
Wed Nov 30, 2022 5:44 pm
Forum: ROMS Discussion
Topic: Question about Flather 2d boundary condition for ubar
Replies: 10
Views: 14604

Re: Question about Flather 2d boundary condition for ubar

Ok, thanks for taking the time to explain it! I will test the various options soon. I just started the thread because I was worried that 1/0.5=2 was mistakenly coded as 0.5, but in hindsight I learned a lot.
by stef
Wed Nov 30, 2022 5:00 pm
Forum: ROMS Discussion
Topic: Question about Flather 2d boundary condition for ubar
Replies: 10
Views: 14604

Re: Question about Flather 2d boundary condition for ubar

The distance between zeta(west,j) and zeta(Istr,j) is always a full pm. Now I'm confused. Isn't this inconsistent with set_tides.F? Isn't zeta(west,j) already set to be centered at u-points: 322 ! If appropriate, load tidal forcing into boundary arrays. The "zeta" 323 ! boundary arrays ar...
by stef
Wed Nov 30, 2022 4:16 pm
Forum: ROMS Discussion
Topic: Question about Flather 2d boundary condition for ubar
Replies: 10
Views: 14604

Re: Question about Flather 2d boundary condition for ubar

Thanks for the tip, yes, I have not looked into the time stepping. I will write into my note that I don't fully understand it yet. Maybe I can figure it out in spring or summer when it's less busy. But for the record: I didn't see your first code block in the Flather condition, only in the reduced-p...
by stef
Wed Nov 30, 2022 10:21 am
Forum: ROMS Discussion
Topic: Question about Flather 2d boundary condition for ubar
Replies: 10
Views: 14604

Re: Question about Flather 2d boundary condition for ubar

Thanks for the reply! Ok, I'll test it. Along with testing this, I should test prescribing the barotropic velocity. First I wanted to document carefully what I have been doing up to now (reduced physics). I attached an excerpt from a note I'm writing, in case anybody is interested in this topic. The...
by stef
Mon Nov 28, 2022 8:09 am
Forum: ROMS Discussion
Topic: Question about Flather 2d boundary condition for ubar
Replies: 10
Views: 14604

Re: Question about Flather 2d boundary condition for ubar

I had the same thought about the staggering of zeta_west, i.e. that it's centered on u-points. This is supported by the comments/code in set_tides.F, see below. zeta_w+=0.5(zeta_{i-1}+zeta_i). But assuming that's right, a forward difference would be (zeta_i-zeta_w)/(0.5*dx) = 2*pm*(zeta_i-zeta_w). !...
by stef
Fri Nov 25, 2022 5:15 pm
Forum: ROMS Discussion
Topic: Question about Flather 2d boundary condition for ubar
Replies: 10
Views: 14604

Question about Flather 2d boundary condition for ubar

I'm trying to write up the boundary condition I use and I wanted to confirm whether the following is correct: ! ! Western edge, Flather boundary condition. ! ELSE IF (LBC(iwest,isUbar,ng)%Flather) THEN DO j=Jstr,Jend IF (LBC_apply(ng)%west(j)) THEN #if defined SSH_TIDES && !defined UV_TIDES ...
by stef
Wed Aug 31, 2022 1:47 pm
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

Just a quick callout to a past post of mine viewtopic.php?t=4665 on performance versus memory bandwidth. Thanks, I haven't seen this! ROMS scales well with increasing number of MPI processes (eg. NtileI*NtileJ) until the number of processes exceeds twice or three times the number of memory channels...
by stef
Mon Aug 29, 2022 7:46 am
Forum: ROMS Benchmarks
Topic: reference x86-64 linux images for benchmarking?
Replies: 0
Views: 2223

reference x86-64 linux images for benchmarking?

After reading the thread [1] I am wondering if it wouldn't be nice to have a reference operating system for benchmarking. I don't know much about computers, so I'm not sure if this makes sense. Argument for using clearlinux [2] to produce reproducible OS images for benchmarking x86-64 hardware: Hern...
by stef
Mon Aug 29, 2022 7:12 am
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

Ok, thanks for the tip! Here is a tomshardware article about the AMD threadripper strategy [1], if I understand correctly the prices for the workstation CPUs with quad- and octacore memory do not differ much from server CPUs nowadays. There used to be an "enthusiast" line which they don't ...
by stef
Fri Aug 26, 2022 6:54 pm
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

Wow, the information on the thread you link to is all new to me! I looked at your figure in this thread and wondered why there are so many tiles. Now I see it's explained in the thread. My hopes are up that I can get a significant speedup by configuring everything properly. Hernan is right in the po...
by stef
Fri Aug 26, 2022 4:00 pm
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

Good to know, thanks for all the great advice!
by stef
Fri Aug 26, 2022 2:37 pm
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

Get a high-level gaming mother board. Quad-channel memory system in a MUST. Do not settle with dual-channel memory. Ok! In AMD land that would imply a Ryzen Threadripper I think. They have 4 or 8 channels. In terms of price difference to a server: Doing a quick search, Lenovo offers a ThinkStation ...
by stef
Thu Aug 25, 2022 3:19 pm
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

Good point, and even with a flat internet rate, I guess the compute nodes should be close to wherever the large data archive is. I guess the largest data set is the output data, assuming one wants to keep it all. There is a blog post saying it's ok to run 24/7, except the fans deteriorate easily: ht...
by stef
Tue Aug 23, 2022 9:32 am
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

Thanks for sharing your experience! Yes, SWAN takes up lots of resources. I guess I'll try to run everything on my existing workstation for now.
by stef
Mon Aug 22, 2022 1:40 pm
Forum: ROMS Benchmarks
Topic: Consumer-grade desktop for toy-"operational" model running 24/7?
Replies: 19
Views: 21597

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

I am running a toy "operational" model and consider buying a dedicated consumer-grade desktop for it. The plan is to pick some number of cores (arbitrary, say 6 or 12 cores) and adjust the model configuration and data processing pipeline such that the desktop runs continuously under full l...
by stef
Wed Jul 06, 2022 6:10 am
Forum: ROMS Bugs
Topic: Possible bug in ROMS/Utility/convolve.F
Replies: 2
Views: 4250

Re: Possible bug in ROMS/Utility/convolve.F

Indeed, thanks! Sorry, I thought I had the current release before reporting.

Edit: I take it back, apparently it was fixed after I reported it :D :
https://www.myroms.org/projects/src/ticket/919

Anyway, thanks for the update!!
by stef
Tue Jun 28, 2022 12:23 pm
Forum: ROMS Bugs
Topic: Possible (minor?) bug in wrt_ini.F
Replies: 2
Views: 4200

Re: Possible (minor?) bug in wrt_ini.F

So this is intended behaviour: The units of ocean_time in the ini files can be 'days' when doing forward runs, but *must* be 'seconds' when doing 4dvar runs, otherwise the program will error at runtime (after multiple timesteps have already been performed).

Okay, thanks for the clarification!
by stef
Mon Jun 27, 2022 2:03 pm
Forum: ROMS Bugs
Topic: Possible (minor?) bug in wrt_ini.F
Replies: 2
Views: 4200

Possible (minor?) bug in wrt_ini.F

I have been using ocean_time units of 'days since 1970...' for some time now in the initial conditions netcdf. This worked (even though in varinfo.yaml the units are s). However it does not seem to work when using 4dvar, because the updated initial conditions are always written in seconds (I think a...
by stef
Fri Jun 24, 2022 3:25 pm
Forum: ROMS Bugs
Topic: Possible bug in /ROMS/Utility/time_corr.F
Replies: 0
Views: 2169

Possible bug in /ROMS/Utility/time_corr.F

Following the 4dvar tutorial there was a runtime error: AD_WRT_HIS_NF90 - wrote adjoint fields (Index=1,1) in record = 5 Found Error: ** Line: 525 Source: ROMS/Utility/time_corr.F, time_corr_nf90 TIME_CORR_NF90 - error while writing variable: ubar at time record = 5 into NetCDF file: wc13_tlf.nc The...
by stef
Thu Jun 23, 2022 2:45 pm
Forum: ROMS Adjoint
Topic: Are the observation error covariances in wc13 actually instrument errors?
Replies: 2
Views: 4231

Re: Are the observation error covariances in wc13 actually instrument errors?

Thanks for the information! Searching for references to Desroziers et al., I found Fig. 12a) of [1] where I think the 2 cm error for the 30km experiment are justified. [1] The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems. Part II – Performance and applica...
by stef
Wed Jun 22, 2022 2:32 pm
Forum: ROMS Adjoint
Topic: Are the observation error covariances in wc13 actually instrument errors?
Replies: 2
Views: 4231

Are the observation error covariances in wc13 actually instrument errors?

I was wondering if the observation error covariances in the tutorial represent instrument errors or the actual spread of the data around a valid model. In the tutorial excerise 1, WC13 is described as: While 30 km resolution is inadequate for capturing much of the energetic meso-scale circulation as...
by stef
Wed Jun 15, 2022 3:14 pm
Forum: ROMS Bugs
Topic: Possible bug in ROMS/Utility/convolve.F
Replies: 2
Views: 4250

Possible bug in ROMS/Utility/convolve.F

For tutorial exercise 3 I tried to compile a RBL4DVAR application without RPCG. It fails with cd xxx; /usr/bin/mpif90 -c -frepack-arrays -fallow-argument-mismatch -O3 -ffast-math -fallow-argument-mismatch convolve.f90 convolve.f90:299:37: 299 | CALL wrt_ini (ng, Rnew(ng)) | 1 Error: Missing actual a...
by stef
Sat Jun 11, 2022 8:39 am
Forum: ROMS Bugs
Topic: Undefined/uninitialized local variable eps for WET_DRY
Replies: 2
Views: 4334

Re: Undefined/uninitialized local variable eps for WET_DRY

Thanks for taking the time to provide so much background! I will switch it off, work through more of the tutorials/lectures you provide and at some point revisit your post when I understand more about assimilation (I'll admit I don't yet understand much of what you are saying about the math). I get ...
by stef
Wed Jun 08, 2022 5:05 pm
Forum: ROMS Bugs
Topic: Undefined/uninitialized local variable eps for WET_DRY
Replies: 2
Views: 4334

Undefined/uninitialized local variable eps for WET_DRY

Sorry Hernan I think it's a WET_DRY bug (you mentioned in a previous post that maintaining the option is a burden). No urgency, I can do it myself, but still wanted to report it. Using WET_DRY in a 4dvar results in: cd xxx; /usr/bin/mpif90 -c -frepack-arrays -O3 -ffast-math -fallow-argument-mismatch...
by stef
Tue Jun 07, 2022 10:01 am
Forum: ROMS Bugs
Topic: Typo in matlab/4dvar/d_ssh_obs.m
Replies: 1
Views: 3592

Typo in matlab/4dvar/d_ssh_obs.m

I am tinkering with the assimilation tutorial, and this looks like a typo: provenance.Targo = 6; % ARGO floats temperature from Met Office provenance.Targo = 7; % ARGO floats salinity from Met Office The same lines are in d_sst_obs.m. The .Targo/.Sargo properties are also used in d_ts_metoffice.m, b...
by stef
Fri Jun 03, 2022 7:52 am
Forum: ROMS Discussion
Topic: GFS forcing: forecast at analysis time
Replies: 0
Views: 1559

GFS forcing: forecast at analysis time

Last week I asked a question to ncep.webmaster@noaa.gov about the GFS forcing, but haven't received an answer yet. I was wondering if someone on this forum could help, it might be of interest to many beginners here. I do not understand the difference between the analysis fields and forecast data wit...
by stef
Mon May 30, 2022 10:24 am
Forum: User Applications
Topic: How to harness ECMWF "Open Data"?
Replies: 1
Views: 2303

Re: How to harness ECMWF "Open Data"?

Apparently they may open enough data for WRF downscaling in 2023: For open data, the next plan is to open a 'CAMS Support set' in 2023 which contains more parameters - this would be with a view that users could run WRF. This is their entire answer: Hi Stefan, Apologies that the answer was not suffic...
by stef
Fri May 27, 2022 7:54 am
Forum: User Applications
Topic: How to harness ECMWF "Open Data"?
Replies: 1
Views: 2303

How to harness ECMWF "Open Data"?

Recently ([1]), ECMWF opened access to some of their data. Forecasts are available globally at 0.4 degrees resolution (i.e reduced spatial resolution), in steps of 3h for at least 90 hours. The available fields are listed in [2]. Unfortunately, the list is incomplete for a bulk flux parametrization....
by stef
Tue Apr 26, 2022 7:43 am
Forum: ROMS Discussion
Topic: Flather for u/vbar with reduced physics without tide elevations
Replies: 0
Views: 2000

Flather for u/vbar with reduced physics without tide elevations

I use data from a global model without tides for my u/v/T/S boundary, and (as noted above) use ADD_FSOBC, so the SSH of the global model is used too. I want to make a sensitivity study quantifying the effect of tidal forcing on the circulation with ROMS, and I don't use UV_TIDES. I have Chapman-Flat...
by stef
Tue Apr 12, 2022 3:17 pm
Forum: ROMS Discussion
Topic: Manning-Strickler for 2d bottom friction
Replies: 5
Views: 5793

Re: Manning-Strickler for 2d bottom friction

Great, thanks for your help!
by stef
Tue Apr 12, 2022 9:16 am
Forum: ROMS Discussion
Topic: Manning-Strickler for 2d bottom friction
Replies: 5
Views: 5793

Manning-Strickler for 2d bottom friction

I'm comparing 2d vs 3d tide simulations with ROMS, and many studies with 2d tide models use a Manning-Strickler formula. If I wanted to use a Manning-Strickler type formulation for bottom friction into ROMS in a quick-and-dirty fashion, would I insert something like cff1=0.25_r8*(vbar(i ,j ,krhs)+ &...
by stef
Mon Mar 28, 2022 3:41 pm
Forum: ROMS Documentation
Topic: Missing link to Broquet et al. 2009 journal article
Replies: 2
Views: 5000

Re: Missing link to Broquet et al. 2009 journal article

Thanks for your reply, I went to researchgate and made a request.
by stef
Mon Mar 28, 2022 8:36 am
Forum: ROMS Documentation
Topic: Missing link to Broquet et al. 2009 journal article
Replies: 2
Views: 5000

Missing link to Broquet et al. 2009 journal article

I was trying to download one of the earlier papers on data assimilation [1] from the myroms.org page: https://www.myroms.org/index.php?page=papers Could it be that there is a "bug" on that website: list items 3. and 5. point to the same pdf, which is the excerpt on the Mercator Ocean Quart...
by stef
Thu Feb 17, 2022 2:58 pm
Forum: ROMS Adjoint
Topic: Which ROMS source version to use for the 4DVar tutorial
Replies: 3
Views: 5047

Re: Which ROMS source version to use for the 4DVar tutorial

I confirm that the first excercise ran without error, with the wc13.h and varinfo.dat provided in the test repository, and the latest src version commit 3b74341007c1b2312f091fed3758655d8f2e485f (HEAD -> master, origin/master, origin/HEAD) Author: Hernan G. Arango <arango@marine.rutgers.edu> Date: Su...
by stef
Thu Feb 17, 2022 2:21 pm
Forum: ROMS Adjoint
Topic: Which ROMS source version to use for the 4DVar tutorial
Replies: 3
Views: 5047

Re: Which ROMS source version to use for the 4DVar tutorial

Your best strategy is to start with the latest version of the code.
I see! Thanks!
by stef
Thu Feb 17, 2022 9:59 am
Forum: ROMS Adjoint
Topic: Which ROMS source version to use for the 4DVar tutorial
Replies: 3
Views: 5047

Which ROMS source version to use for the 4DVar tutorial

I had trouble running the first exercise in the 4DVar tutorial. I checked out svn checkout --username joe_roms https://www.myroms.org/svn/src/test my_test as described in [1]. I checked out revision 974 as described in [2]. But they seem to be incompatible, the files in the project directories do no...
by stef
Thu Jan 06, 2022 8:56 am
Forum: ROMS Tools and Techniques
Topic: Extraction of ERA 5 data for custom grid
Replies: 34
Views: 68932

Re: Extraction of ERA 5 data for custom grid

Fantastic! Thanks! I'm sure a lot of us appreciate the work you put into this. It seems to be a very widely used data set, and having the options for direct momenum/heat/salt fluxes and BULK_FLUXES is great. However, something we have noticed in tests is that near the coast, in regions where the ERA...
by stef
Wed Jan 05, 2022 9:03 pm
Forum: ROMS Tools and Techniques
Topic: Extraction of ERA 5 data for custom grid
Replies: 34
Views: 68932

Re: Extraction of ERA 5 data for custom grid

Thanks for the WMO paper, I should read it! These are great. Yes we should look at the rda.ucar documentation. Unfortunately I did not realize that rda.ucar offers the ERA5 data set until it was too late. I like to download the netcdfs via thredds/opendap on rda.ucar. But I had already invested some...
by stef
Wed Jan 05, 2022 5:41 pm
Forum: ROMS Tools and Techniques
Topic: Extraction of ERA 5 data for custom grid
Replies: 34
Views: 68932

Re: Extraction of ERA 5 data for custom grid

I just saw in [1] it says In HRES atmospheric (hourly steps), mean rate precipitation at time=06:00, step=3 represents the average precipitation rate in the 1-hour period from 08:00 to 09:00, so that's clear now. I'm downloading from "cds.climate.copernicus.eu". Not sure which time axis &q...
by stef
Wed Jan 05, 2022 3:37 pm
Forum: ROMS Tools and Techniques
Topic: Extraction of ERA 5 data for custom grid
Replies: 34
Views: 68932

Re: Extraction of ERA 5 data for custom grid

I'm trying to learn how to use atmospheric forcing. Thanks for the matlab scripts on github! They are great. I saw in [1] that e.g. the mean evaporation rate in the era5 reanalysis is This parameter is a mean over a particular time period (the processing period) which depends on the data extracted. ...
by stef
Thu Dec 16, 2021 2:14 pm
Forum: ROMS Discussion
Topic: Unrealistic M2 tide amplitude in Sulu Sea, Indonesian Sea
Replies: 5
Views: 2189

Re: Unrealistic M2 tide amplitude in Sulu Sea, Indonesian Sea

You integrate for 5 days 2400 ntimes Number of timesteps for 3-D equations. 180.000 dt Timestep size (s) for 3-D equations. I would check if the cotide/corange plots are sensitive to the integration period. Depending on how many constituents you are solving for in the LSQR routine, 5 days may be too...
by stef
Tue Dec 07, 2021 8:20 am
Forum: ROMS Discussion
Topic: Differences 2d vs. 3d barotropic simulation with very large AKV_BAK
Replies: 3
Views: 2549

Re: Differences 2d vs. 3d barotropic simulation with very large AKV_BAK

Turns out I probably took the wrong turn here. Most of the difference does not appear to be due to the splitting itself, but due to the fact that that bottom stress was artificially limited in the 3d application. In set_vbc.F, there is: # ifdef LIMIT_BSTRESS cff3=cff*0.5_r8*(Hz(i-1,j,1)+Hz(i,j,1)) b...
by stef
Tue Nov 30, 2021 2:36 pm
Forum: ROMS Discussion
Topic: Differences 2d vs. 3d barotropic simulation with very large AKV_BAK
Replies: 3
Views: 2549

Re: Differences 2d vs. 3d barotropic simulation with very large AKV_BAK

I think the difference is mostly due to the differences in bottom stress application between 2d vs. 3d. From main3d.f90 it seems that the bottom stress is computed at the beginning of the baroclinic time step (set_vbc) and held constant throughout the barotropic stepping, whereas in 2d applications ...
by stef
Mon Nov 29, 2021 5:10 pm
Forum: ROMS Discussion
Topic: Differences 2d vs. 3d barotropic simulation with very large AKV_BAK
Replies: 3
Views: 2549

Re: Differences 2d vs. 3d barotropic simulation with very large AKV_BAK

A run with ndtfast=4 is much similar to the 2d solution. Could it be that the "slow forcing" terms in the 2d momentum equations are the reason? In my case, these are * vertically integrated 3d "rhs" terms (rhs3d.f90) * bottom/surface stress plus (rhs3d.f90) * 3d horizontal viscos...
by stef
Wed Nov 24, 2021 5:09 pm
Forum: ROMS Discussion
Topic: Differences 2d vs. 3d barotropic simulation with very large AKV_BAK
Replies: 3
Views: 2549

Differences 2d vs. 3d barotropic simulation with very large AKV_BAK

I'm running tide simulations and want to compare the results of 2d simulations to 3d barotropic simulations (no tracers, homegeneous density). Vertical viscosity is constant and prescribed by AKV_BAK. Should I expect the 3d solution to converge to the 2d solution for very large AKV, e.g. AKV=1e5? I ...
by stef
Mon Nov 15, 2021 2:46 pm
Forum: ROMS Bugs
Topic: Possible bug in def_tides.F
Replies: 1
Views: 3857

Possible bug in def_tides.F

I have the following error: Found Error: ** Line: 311 Source: ROMS/Utility/wrt_tides.F, wrt_tides_nf90 WRT_TIDES_NF90 - error while writing variable: v_tide into detide harmonics NetCDF file: ... I tried to fix it as follows, but don't have a run long enough to validate the detided output, so I'm no...
by stef
Thu Sep 09, 2021 6:46 am
Forum: ROMS Discussion
Topic: How to define tidal constituent SA?
Replies: 5
Views: 3239

Re: How to define tidal constituent SA?

Again, as I said before, for ocean tides the seasonal effects (tropical year) dominate over the gravitational effects (anomalistic year) so it makes sense to use the tropical year, whereas for earth tides, presumably the opposite is true. Thanks for pointing this out again, I had not fully understo...
by stef
Tue Sep 07, 2021 3:59 pm
Forum: ROMS Discussion
Topic: How to define tidal constituent SA?
Replies: 5
Views: 3239

Re: How to define tidal constituent SA?

Wow, I just saw that appendix A of [1] explains it very well, and it's open access!

[1] Ray, R.D., Loomis, B.D. & Zlotnicki, V. The mean seasonal cycle in relative sea level from satellite altimetry and gravimetry.
https://link.springer.com/article/10.10 ... 21-01529-1
by stef
Fri Aug 20, 2021 2:27 pm
Forum: ROMS Discussion
Topic: Underpredicted tidal amp: improvements with 3d instead of 2d?
Replies: 0
Views: 1713

Underpredicted tidal amp: improvements with 3d instead of 2d?

I'm modelling tides with 2d configuration in the north-eastern Bay of Bengal, as a precursor to a storm-surge modelling exercise. A figure of the northern part of the domain (Ganges-Brahmaputra-Meghna Delta): https://sriha.net/articles/assets/harmAnaBoBPart1/images/map_closeup_overview.png The white...
by stef
Thu Jun 10, 2021 1:44 pm
Forum: ROMS Discussion
Topic: How to define tidal constituent SA?
Replies: 5
Views: 3239

Re: How to define tidal constituent SA?

Very good, thanks for the advice. The question came up indirectly because [1] don't seem to distinguish between "H1" (aka "alpha_2")in the Foreman software, with Cartwright number (2 0 -1 0 0 1), and the "MA_2" constituent proposed by [2], which IMHO should have Cartwri...
by stef
Wed Jun 09, 2021 1:29 pm
Forum: ROMS Discussion
Topic: How to define tidal constituent SA?
Replies: 5
Views: 3239

How to define tidal constituent SA?

I have a general oceanography question about the tidal constituent SA. Does it have the period of a tropical or anomalistic year? I know it's not directly related to ROMS. Sorry, but currently I don't know where else to ask. In the harmonic analysis packages [1] and [2], SA is defined as inverse of ...
by stef
Fri May 14, 2021 11:06 am
Forum: ROMS Bugs
Topic: possible bug in set_tides.F
Replies: 1
Views: 4361

possible bug in set_tides.F

I'm trying to do harmonic analysis and detiding, and get weird results for a child grid when nesting. The results for the parent grid seem ok. Searching for a possible bug in ROMS, I noticed in set_tides.F ! ! If nested grids and detiding, load period and harmonics to other ! grids. ! IF (LprocessTi...
by stef
Wed May 05, 2021 3:40 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

But I do not think that is the cause because that issue seems to affect the input of information from coarse to fine grid and thus the artifacts seen in my results would be present in both one-way and two-way nesting. Please let me know if you agree. Yes, agreed! I have yet to try out one-way nesti...
by stef
Thu Apr 22, 2021 2:32 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

If you aren't using identical bathymetry in the contact regions, the following may be interesting too: viewtopic.php?f=14&t=5438&p=21072#p21072 I haven't updated my ROMS code in a while, so I'm not sure if this is still relevant. I thinkt the problem was that the barotropic volume flux throu...
by stef
Wed Apr 21, 2021 3:14 pm
Forum: ROMS Discussion
Topic: Filtering altimetry data
Replies: 1
Views: 1604

Re: Filtering altimetry data

Ok today I searched for 'Dynamic Atmospheric Correction' and there are a lot of resources! Thanks.
by stef
Mon Apr 19, 2021 3:26 pm
Forum: ROMS Discussion
Topic: Filtering altimetry data
Replies: 1
Views: 1604

Filtering altimetry data

I'm interested in storm surges and I'm trying to understand how to use satellite altimetry. The ROMS 4DVar tutorial [1] mentions 'satellite SSH in the form of a gridded product from Aviso'. I'm wondering if one has to do any preprocessing (filtering?) on that data? I read about how datasets like [2]...
by stef
Mon Mar 22, 2021 9:40 am
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

But if there is a problem with the grid topology, how come the experiment with the same timestepping does run fine. Hmm, it's surprising, isn't it?
by stef
Mon Mar 22, 2021 9:32 am
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

[edit: deleted]
by stef
Thu Mar 11, 2021 5:25 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

There is approximately 5km between the stripes and the barotropic wave speed is 166m/s, so there is a period of 30s, which is in sync with the baroclinic time step? So assuming the signal is generated at the nesting boundary and propagating into the fine domain, maybe check the effect of the cpp opt...
by stef
Wed Feb 10, 2021 10:50 am
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

I was thinking that the masking could be done in the preprocessing step that generates the contact zone input file (e.g., 'roms_ngc.nc') so that the model can be tricked to think that there is an island where the coarse grid is supplanted by the fine grid so that the fine2coarse interpolation ignor...
by stef
Wed Feb 10, 2021 9:03 am
Forum: ROMS Installation
Topic: Compiling Error for COAWST
Replies: 6
Views: 10870

Re: Compiling Error for COAWST

(In case this is directed to me, not the original poster...)

For me everything seems to work using the flag mentioned above! Let me know if you still want the output. Thanks!
by stef
Tue Feb 09, 2021 9:02 am
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

If the parent grid discrete stencil for advection and/or diffusion spans, say, 3 parent grid cells then do we need to match bathymetry precisely over this range, which would be refine_factor times the 3 parent cell stencil width. Yes, any insight into this would be great! If I recall correctly, my ...
by stef
Sat Feb 06, 2021 3:51 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 219606

Re: problem with two-way nesting

[Edit: just saw you replied while I was posting...the following is intended as reply to wilkin's comment..] To test this, what do you think about the following: Mask all coarse grid points in the interior of the fine2coarse contact region, except the ones adjacent to the nest boundary, i.e. the regi...
by stef
Wed Feb 03, 2021 9:18 am
Forum: ROMS Installation
Topic: Compiling Error for COAWST
Replies: 6
Views: 10870

Re: Compiling Error for COAWST

I don't understand much of this, but for the record I'm posting error messages when compiling the following version of Roms Rutgers on a recently updated Arch Linux distribution: >>git log commit 697c77dc4836919704033f970a1b051f24c595d0 (grafted, HEAD -> master, origin/master, origin/HEAD) Author: H...
by stef
Mon Feb 01, 2021 5:03 pm
Forum: ROMS Installation
Topic: Compiling Error for COAWST
Replies: 6
Views: 10870

Re: Compiling Error for COAWST

I had a similar error, using gcc 10. Adding the '-fallow-argument-mismatch' to FFLAGS solved it for me. See also here: https://www.myroms.org/forum/viewtopic.php?f=17&t=5696&p=22037&hilit=fallow+argument+mismatch#p22037 https://github.com/Unidata/netcdf-fortran/issues/212 But then I had ...
by stef
Wed Oct 21, 2020 4:27 pm
Forum: ROMS Discussion
Topic: Bathymetry differences in nested simulations
Replies: 5
Views: 3289

Re: Bathymetry differences in nested simulations

Probably not relevant for your 3 vs 5 refinement, but perhaps interesting to look into if you have very shallow regions or very steep bathymetry gradients along the grid: https://www.myroms.org/forum/viewtopic.php?f=14&t=5438&p=21072#p21072 I haven't updated my ROMS code in a while, so I'm n...
by stef
Wed Oct 21, 2020 4:05 pm
Forum: ROMS Tools and Techniques
Topic: wind vector direction in NAM and NARR data
Replies: 7
Views: 9629

Re: wind vector direction in NAM and NARR data

I noticed this thread just a couple of days ago (when agpc posted), and after reading wilkin's first post from 4 years ago, I'm not sure if I'm doing it right. I'm using the source https://www.ncei.noaa.gov/data/north-american-mesoscale-model/access/historical/analysis/201210/20121027/namanl_218_201...
by stef
Fri May 15, 2020 3:10 pm
Forum: ROMS Discussion
Topic: Instabilities in LF-AM3 caused by wetting/drying?
Replies: 4
Views: 3454

Re: Instabilities in LF-AM3 caused by wetting/drying?

think of the beach, and the tide (or surge) coming in and out. what is the horizontal distance from the high tide line to the low tide line - that is the what you should think about.
Good point, thanks again for your help!
by stef
Fri May 15, 2020 2:44 pm
Forum: ROMS Discussion
Topic: Instabilities in LF-AM3 caused by wetting/drying?
Replies: 4
Views: 3454

Re: Instabilities in LF-AM3 caused by wetting/drying?

I think at 12km your spacial scales are probably not resolving those wet-dry processes. Thanks for the reply, that makes sense. I have a nested toy setup and tried to keep the processing of the grids identical as far as possible. So I have not used a minimum depth at all, even in the coarsest grids...
by stef
Fri May 15, 2020 11:03 am
Forum: ROMS Discussion
Topic: Instabilities in LF-AM3 caused by wetting/drying?
Replies: 4
Views: 3454

Instabilities in LF-AM3 caused by wetting/drying?

I have a grid with 12km resolution and dcrit=0.1. Is this reasonable? I have been looking at a specific blowup and I'm wondering if there is an instability caused by the combination of wet-dry flux-inhibitor and the multiple time levels used in the predictor-corrector scheme. The "three-time&qu...
by stef
Tue May 12, 2020 9:19 am
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

Update: note that there is a problem with the patch in that the wet/dry flux limiter (i.e. suppress outflow from a dry cell) does not work at the boundary. An already dry cell adjacent to the boundary may still export water out of the domain, depending on the boundary condition. It could be fixed si...
by stef
Mon Apr 27, 2020 4:47 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

I should say that the "mod_stef.F" is very messy. Sorry. :roll: The only variables needed are the ones used in step2d.
by stef
Mon Apr 27, 2020 4:43 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

Hi, Thanks for your email, I'll respond tomorrow, meanwhile I attached the patch. Please note: *) Some (much?) of this may be incorrect. I'm a beginner. My approach was "if I get a reasonable result, it's good enough for now". *) I have not really understood MPI yet, in particular, in the ...
by stef
Wed Apr 22, 2020 12:42 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

I wrote a blog post [1] about my experience using ROMS in 2d setups with grid nesting and wetting/drying. Very early stage, but if anybody is interested in this, I appreciate comments/suggestions, in particular if you are working on new code that is being pushed to the repository in the near future....
by stef
Wed Apr 15, 2020 1:13 pm
Forum: ROMS Information
Topic: Curvilinear Coordinates
Replies: 13
Views: 38054

Re: Curvilinear Coordinates

In the first post of this thread it says: u(XI,ETA)=u(LON,LAT)*cos(angle(i,j))+v(LON,LAT)*sin(angle(i,j)) v(XI,ETA)=u(LON,LAT)*sin(angle(i,j))-v(LON,LAT)*cos(angle(i,j)) shouldn't this instead be u(XI,ETA)=u(LON,LAT)*cos(angle(i,j))+v(LON,LAT)*sin(angle(i,j)) v(XI,ETA)=-u(LON,LAT)*sin(angle(i,j))+v(...
by stef
Mon Feb 17, 2020 7:13 pm
Forum: ROMS Discussion
Topic: u/bar boundary condition for nested grid in 3d setup
Replies: 0
Views: 1731

u/bar boundary condition for nested grid in 3d setup

Another question about the boundary condition for barotropic flow for nested grids. This applies to 3d setups. In 'u/v2dbc_im.F', there are lines like: cff=0.5_r8*GRID(ng)%on_u(Istr,j)* & & (GRID(ng)%h(Istr-1,j)+zeta(Istr-1,j,kout)+ & & GRID(ng)%h(Istr ,j)+zeta(Istr ,j,kout)) cff1=GR...
by stef
Fri Feb 14, 2020 11:01 am
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

So maybe in a 2d setup, one could store the volume flux of the donor after each barotropic corrector step, and impose that during the 3 (or 5 etc...) time steps made in the receiver grid? Hmm, I guess one should store two time levels of the donor volume flux along the boundary, and then do a tempor...
by stef
Mon Feb 10, 2020 6:31 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

I guess that we never accounted for a User needing nesting in a 2D application. The motivation would be to do storm surge modelling, in a way that is roughly consistent with other storm surge models such as SLOSH or ADCIRC (both mostly 2d, at least until a couple of years ago...I'm not sure about t...
by stef
Fri Feb 07, 2020 8:14 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

I see, thanks. That DU_Avg field is actually appearing in u2dbc_im.F, for the "nested" type boundary conditions (in a 2D application). I substituted it there by a dummy value to test, and then the dummy value appears in my history file. Unfortunately I can only continue testing on Monday. ...
by stef
Fri Feb 07, 2020 6:20 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Re: Possible bug in wetdry.F for nested apps

Thanks for your replies! I work on a 2d application, and another potentially relevant thing I just noticed is that the nesting boundary conditions for the refined grid use fields called e.g. REFINED(cr)%DU_avg2(1,m,tnew) I think this is the case even in 2d applications. However, I think that this fi...
by stef
Fri Feb 07, 2020 1:26 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 19357

Possible bug in wetdry.F for nested apps

Hi,

I read wetdry.F, and was wondering if the index boundaries should account for the contact regions in nested applications? Currently, there are ranges like "JstrR,JendR" etc., but shouldn't they be "JstrT,JendT" instead?

Thanks!

Regards, Stefan
by stef
Wed Jan 29, 2020 8:44 am
Forum: ROMS Bugs
Topic: Possible bug in nesting.F for 2D apps with wetting/drying
Replies: 1
Views: 1953

Re: Possible bug in nesting.F for 2D apps with wetting/drying

I wonder if in 2d applications, the multi-valued mask can just be converted to a binary mask like so
I mean only in fine2coarse2d(), not elsewhere....
by stef
Wed Jan 29, 2020 8:42 am
Forum: ROMS Bugs
Topic: Possible bug in nesting.F for 2D apps with wetting/drying
Replies: 1
Views: 1953

Possible bug in nesting.F for 2D apps with wetting/drying

Hi, This is unrelated to the incorrect netcdf header reproduced in the wiki (as reported earlier in https://www.myroms.org/forum/viewtopic.php?f=17&t=5425 ), but it's related to nesting in general. In 2-d applications with wetting/drying, the u/vmasks supplied to fine2coarse2d() are in their mul...
by stef
Mon Jan 27, 2020 9:18 am
Forum: ROMS Problems
Topic: Three nested refined grids - Ngrids vs. Ncontact
Replies: 2
Views: 4564

Re: Three nested refined grids - Ngrids vs. Ncontact

I think that we, as ROMS developers, have freely provided a very sophisticated modeling framework to the ocean community without asking anything in return. I highly appreciate your effort and that of your team, and the appreciation and admiration grows every day I work with the software you provide...
by stef
Fri Jan 24, 2020 6:03 pm
Forum: ROMS Problems
Topic: Three nested refined grids - Ngrids vs. Ncontact
Replies: 2
Views: 4564

Three nested refined grids - Ngrids vs. Ncontact

Hi, I was able to do a refinement using 2 grids, using a self-written python script to generate the contact file. It's a 2d application and although I hit a bug with the wet masks using wetting/drying, I was able to patch it and proceed. However, now I want to use 3 refined nests (i.e. one telescopi...
by stef
Wed Dec 04, 2019 2:57 pm
Forum: ROMS Tools and Techniques
Topic: netCDF browser for remote work
Replies: 7
Views: 6072

Re: netCDF browser for remote work

Thanks, that sounds perfect! I have never heard about this!

This way I don't have to set up the server, plus it occurred to me I would have needed to set up a VPN too. Much simpler!
by stef
Tue Dec 03, 2019 7:06 pm
Forum: ROMS Tools and Techniques
Topic: netCDF browser for remote work
Replies: 7
Views: 6072

Re: netCDF browser for remote work

Wow, thanks for this! Looks great! I like the zooming with the mouse wheel and the panning.

I'm not sure if I need a time axis with calendar units, but maybe I can tweak it. I'll read the docs! Thanks!
by stef
Tue Dec 03, 2019 6:39 pm
Forum: ROMS Tools and Techniques
Topic: netCDF browser for remote work
Replies: 7
Views: 6072

Re: netCDF browser for remote work

Also, I just found a paper on "Accessing HDF5 Files via SSHFS".

https://support.hdfgroup.org/pubs/paper ... -final.pdf

They say it's comparable to NFS, and useful for subsetting remote files. So it may be comparable to using the THREDDS server anyways.
by stef
Tue Dec 03, 2019 5:28 pm
Forum: ROMS Tools and Techniques
Topic: netCDF browser for remote work
Replies: 7
Views: 6072

Re: netCDF browser for remote work

Thanks for your response! Yes that would be great, but does THREDDS have an option to deliver a UI to the clients so they can browse data, in the sense that they can specify e.g. which x-y slice at fixed z and t they want to look at? This way only the figure would be sent over the internet, not the ...
by stef
Mon Dec 02, 2019 10:50 am
Forum: ROMS Tools and Techniques
Topic: netCDF browser for remote work
Replies: 7
Views: 6072

netCDF browser for remote work

I was wondering if anybody has experience with visualtization tools for remote work over the internet (WWW, not LAN). I used to mount remote directories on my local machine via sshfs, so any netcdf browser will work. However, now I sometimes make larger simulations, and the UI of the browser becomes...
by stef
Thu Jul 04, 2019 7:40 am
Forum: ROMS Bugs
Topic: typos in radiation_stress.F
Replies: 8
Views: 7071

Re: typos in radiation_stress.F

Also, there is a typo in the computation of the radiation stress component derivative for the vbar forcing: >> git diff master mybranch ROMS/Nonlinear/nearshore_mellor08.h ... @@ -1116,9 +1116,9 @@ DO j=JstrV,Jend DO i=Istr,Iend cff=(Syyl(i,j )- & - & Syyl(i,j-1))*on_v(i,j)+ & - & (S...
by stef
Tue Jun 11, 2019 2:00 am
Forum: ROMS Problems
Topic: Stokes velocity fills dry cell uphill
Replies: 3
Views: 2945

Re: Stokes velocity fills dry cell uphill

I tested the spurious wetting of high elevation grid points now with the 3d version of COAWST. No spurious wetting occurred in my test, neither using WEC_VF nor WEC_MELLOR. Browsing through the source code and Warner et al. (2013) I currently think the following (unverified/untested hypothesis): In ...
by stef
Thu May 30, 2019 10:21 pm
Forum: ROMS Problems
Topic: Stokes velocity fills dry cell uphill
Replies: 3
Views: 2945

Re: Stokes velocity fills dry cell uphill

Thanks for your reply! Yes, I should have mentioned it's a 2d application. I'm trying to simulate storm surges. At some point I want to have at least a bottom boundary layer, but having 2d simulations for reference and comparison is interesting. (also I just noticed that I have the WET_DRY_LIMIT dir...
by stef
Thu May 30, 2019 5:25 am
Forum: ROMS Problems
Topic: Stokes velocity fills dry cell uphill
Replies: 3
Views: 2945

Stokes velocity fills dry cell uphill

I have tried ROMS+SWAN coupling with wetting/drying and the model blew up. jcwarner already had told me to use the COAWST version for ROMS+SWAN coupling, but I wasn't able to compile it (see below for the error. Also, I'm not sure where to report errors.). So I kept working with the Rutgers version....
by stef
Thu May 23, 2019 11:53 pm
Forum: ROMS Discussion
Topic: Alternative wave coupling for efficient storm surge modeling
Replies: 0
Views: 1641

Alternative wave coupling for efficient storm surge modeling

I'm interested in modeling storm surges in the context of probabilistic risk assessment, e.g. using the Joint Probability Method as described e.g. in [1]. See also ROMS forum post [3]. I just found a paper [2] describing an efficient alternative to SWAN. The authors of [2] note that, as opposed to S...
by stef
Thu May 23, 2019 5:16 am
Forum: ROMS Bugs
Topic: typos in radiation_stress.F
Replies: 8
Views: 7071

Re: typos in radiation_stress.F

Thanks John, I sent the email! Sooner or later I wanted to learn about COAWST anyways, so that's great. For the near future, I have to keep it simple as I'm still learning nearshore processes, wind waves and small scales in general. I looked for the most simple method to compute wave setup. SWAN has...
by stef
Wed May 22, 2019 11:15 pm
Forum: ROMS Bugs
Topic: typos in radiation_stress.F
Replies: 8
Views: 7071

Re: typos in radiation_stress.F

I'm trying to set up a minimal storm surge model and hit the same bug. It seems to affect any 2d application using both NEARSHORE_MELLOR0*. I don't know much about wind waves yet, still reading literature. But this is still a bug, isn't it?
by stef
Thu Apr 18, 2019 12:53 am
Forum: ROMS Discussion
Topic: Coastal flood risk assessment with ROMS?
Replies: 4
Views: 13783

Re: Coastal flood risk assessment with ROMS?

I should add that an example for the 1D transects can e.g. be found on page 25 in: Federal Emergency Management Agency, 2013: Flood Insurance Study (preliminary). City of New York, New York. https://hazards.fema.gov/femaportal/prelimdownload/prelim/ProductsDownLoadServlet?pfiProdId=16353 This public...
by stef
Thu Apr 18, 2019 12:41 am
Forum: ROMS Discussion
Topic: Coastal flood risk assessment with ROMS?
Replies: 4
Views: 13783

Coastal flood risk assessment with ROMS?

Hi, May I kindly ask for your opinions on using ROMS for storm surge modelling? I'm particularly interested in risk assessment, similar to what has been done by FEMA for coastal flood hazard studies. Here is a link to the New York - New Jersey region: http://www.region2coastal.com/ They are using AD...
by stef
Thu Jun 08, 2017 1:59 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

Just a quick warning in case anybody (like me) considers making oblique Mercator projections with Proj version (4.9.3): be careful :wink:
by stef
Tue May 30, 2017 5:44 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

It turns out that "oblique Mercator projection" is exactly what I need, I think. It lets you define two points in lat,lon coordinates through which the central line for the cylindrical projection goes (that's the line at which the projection cylinder touches the sphereoid), and this can be...
by stef
Mon May 29, 2017 3:12 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

Wow, Lambert conformal conic! On Wikipedia it says: By scaling the resulting map, two parallels can be assigned unit scale, with scale decreasing between the two parallels and increasing outside them. This gives the map two standard parallels. In this way, deviation from unit scale can be minimized ...
by stef
Thu May 25, 2017 7:21 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

figure_grid_rotation_tutorial_mercator.png I'd like to show off a figure about one way to correctly (?) rotate the grid, where I'm trying to follow the steps from above post: For pm, pn, I have an initial value based on the flat projection space. Then I later remap to the sphere so that I can assig...
by stef
Tue May 23, 2017 7:41 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

Very good, thanks for your patience and the information, it helped a lot. I will try to process all of this during the next couple of weeks :)
by stef
Tue May 23, 2017 6:01 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

Ok, but ROMS doesn't care whether it's orthogonal on the sphere, right? What it wants is orthogonality in the plane. And this cannot be fulfilled if we measure pm and pn on the sphere and then use it on the plane. Is this right? Or is there a conformal projection that can do that? The way I see it, ...
by stef
Tue May 23, 2017 4:38 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

Thanks for the help! Now I also found a presentation , where you explain this on page 4 (may be helpful for other readers). However, strictly speaking, when we measure pm, pn on the sphere and use it in the horizontal plane, we don't get an orthogonal grid, right? So if we computed the ds quantities...
by stef
Mon May 22, 2017 3:45 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

Ok thanks a lot for your time Kate. Let me see if I understand correctly: The real question to ask is what to use for pm and pn in the first place. This reflects how the projection is made from the sphere to the plane. It already involves a trade off: If one chooses Mercator, one doesn't conserve vo...
by stef
Fri May 19, 2017 2:37 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

After spending some time today browsing through the SCRIP manual, I think I should better explain the motivation of my initial question. If I understand correctly :? none of the bilinear and bicubic methods take into account the sphericity, neither in SCRIP nor in ROMS. What they do is interpolate i...
by stef
Thu May 18, 2017 6:38 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Re: Roms grid rotation and interpolation

Very good, thanks! Meanwhile I also found some Python code, specific to ROMS: https://github.com/ESMG/pyroms/blob/master/pyroms/pyroms/remapping/remap.py with "conservative" and "bicubic" remapping. It seems they use the "SCRIP" library you linked to, so that starts to ...
by stef
Mon May 15, 2017 3:22 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 14550

Roms grid rotation and interpolation

I'm trying to configure a model to practice data assimilation, but stumbled across a very basic question about interpolation (maybe embarrassing for a geophysicist :oops: , but I never had to do with that problem before.). When I interpolate data on an equidistant lat/lon grid, do I need to do some ...
by stef
Mon Jan 23, 2017 2:00 am
Forum: ROMS Benchmarks
Topic: Benchmarks
Replies: 31
Views: 61044

Re: Benchmarks

met_c48xlarge.png Hi, I'm tinkering with benchmarks on Amazon Web Services (AWS). The attached figure shows time (in seconds) spent per process for the ROMS ”large” benchmark test (benchmark3.in), as function of the number of processes. Axes are logarithmic with base 2. Computations are performed o...
by stef
Thu Jan 12, 2017 11:05 pm
Forum: ROMS Installation
Topic: Anyone tried running ROMS on HPC Cloud (e.g. Amazon)?
Replies: 8
Views: 11037

Re: Anyone tried running ROMS on HPC Cloud (e.g. Amazon)?

Yes I'm trying AWS. At the moment I'm doing some benchmark testing, see https://github.com/poidl/awsroms and https://poidl.github.io/awsroms/ . I still have few results, and they are not interpreted properly. I'd like to become an independent (sub-) contractor, so I'm also very interested in estimat...
by stef
Thu Sep 29, 2016 4:10 pm
Forum: ROMS Bugs
Topic: make caught in infinite loop
Replies: 2
Views: 4062

Re: make caught in infinite loop

That's it, thanks. Setting the right path in the platform dependent file solves it!

On a side note, this still leaves me wondering why the line after the prerequisites is executed (sfmakedepend), if the prerequisites are not fulfilled? Not urgent, though.
by stef
Tue Sep 27, 2016 8:23 am
Forum: ROMS Bugs
Topic: make caught in infinite loop
Replies: 2
Views: 4062

make caught in infinite loop

I'm trying to run the upwelling example, and make gets caught in an infinite loop. To reproduce: I change the compiler and netcdf paths in the makefile: >> diff makefile path/to/original/makefile 142c142 < FORT ?= gfortran --- > FORT ?= pgi 312,313c312,313 < NETCDF_MODFILE := ../../include/netcdf.mo...
by stef
Mon Sep 24, 2007 5:11 pm
Forum: ROMS Tools and Techniques
Topic: analytical coriolis with netcdf bathymetry
Replies: 2
Views: 3517

Thanks for the guidance! I'll try that.
--stefan
by stef
Mon Sep 24, 2007 1:33 pm
Forum: ROMS Tools and Techniques
Topic: analytical coriolis with netcdf bathymetry
Replies: 2
Views: 3517

analytical coriolis with netcdf bathymetry

Hello! We use a grid generated with 'Roms Tools', and we changed the code in the m-files to have constant, user defined coriolis acceleration. Everytime we make a run with different f we have to generate a new grid, although the rest is the same as before. Is there an easy way to analytically define...
by stef
Thu Aug 30, 2007 5:51 pm
Forum: ROMS Discussion
Topic: Collaboration, releases, SVN, Trac, flavors of ROMS...
Replies: 11
Views: 19845

Relating to this discussion I want to point to a different post. The following is *very* instructive:
viewtopic.php?t=622
Thanks!
by stef
Fri Aug 24, 2007 8:49 am
Forum: ROMS Discussion
Topic: Collaboration, releases, SVN, Trac, flavors of ROMS...
Replies: 11
Views: 19845

Thank you for the replies! Ok, we will try to isolate the user-modified parts of the code and avoid using svn for anything else than checking out a working copy from the roms/toms-group repository for every user. This seems to be easy and gives every user a chance to learn how to use svn. When we ge...
by stef
Thu Aug 23, 2007 8:59 am
Forum: ROMS Discussion
Topic: Collaboration, releases, SVN, Trac, flavors of ROMS...
Replies: 11
Views: 19845

We are in the process of switching to the new ROMS version and I have some questions concerning administration: Imagine that on a cluster of our university, there are three people using ROMS and making setups, all working on the same project, but all going into a slightly different direction (one is...
by stef
Mon Jul 16, 2007 11:02 am
Forum: ROMS Installation
Topic: An urgent problem about downloading ROMS by using SVN
Replies: 5
Views: 5665

Hi Zhou Wei!
If it's because of the blank, mybe you can try type your username 'Zhou\ Wei'. the backslash is use in linux (and I guess unix) to mark a blank.
good luck
--Stefan
by stef
Wed Jul 11, 2007 7:38 am
Forum: ROMS Installation
Topic: running roms with debugger
Replies: 12
Views: 11408

Hi Mark! Yes, I will read the manual for ifort to see what kind of options it has. I also will ask our sysadmin to install g95 and idb. SALINITY is undefined. I think the error was caused by my altering the files over and over again and making some copy-paste error with emacs, because the error show...
by stef
Tue Jul 10, 2007 5:31 pm
Forum: ROMS Installation
Topic: running roms with debugger
Replies: 12
Views: 11408

Hi!
The post I was referring to is the one from Mon May 28, 2007 10:34 pm in:
viewtopic.php?t=538

There seem to be a lot of other posts on this topic, but I haven't read them yet...
by stef
Tue Jul 10, 2007 8:34 am
Forum: ROMS Installation
Topic: running roms with debugger
Replies: 12
Views: 11408

Hi! I recently ran into troubles with segmentation faults and it was the first time that I tried the debug option. I'm using ifort on Linux. As m.hadfield suggested in another post, I also added a '-traceback' option the the '-g -check bounds' in the Linux-ifort makefile. Now when I start oceanG, th...
by stef
Fri Jun 15, 2007 2:52 pm
Forum: ROMS Problems
Topic: barotropic wave in the channel (no forcing, no dissipation)
Replies: 9
Views: 10003

Concerning the pressure gradient: I'm looking at prsgrd31.h (probably not good, but maybe for flat ocean bottom it's ok) and there is a term called fac2 that is multiplied by a gradient of z_w(i,j,N). It looks like the shallow-water-part of the prsgrd. z_w(i,j,N) is set in set_depth (called at the l...
by stef
Thu Jun 14, 2007 8:16 am
Forum: ROMS Problems
Topic: barotropic wave in the channel (no forcing, no dissipation)
Replies: 9
Views: 10003

I see. So the baroclinic mode is effectively doing nothing. I'm not sure, but my guess is that is has to do with the time-filtering of the barotropic mode in the coupled system. I don't know if time-filtering is used in the standalone 2d-code, but I guess not, since it serves to remove instabilities...
by stef
Wed Jun 13, 2007 7:44 pm
Forum: ROMS Problems
Topic: barotropic wave in the channel (no forcing, no dissipation)
Replies: 9
Views: 10003

Hi!
Sorry, I don't know the answer, but I wonder how the vertical velocity-profile looks like? Is v z-independent, like it's supposed to be?
by stef
Tue Jun 05, 2007 8:26 pm
Forum: ROMS Information
Topic: Grid cell volume
Replies: 7
Views: 10126

I think that what I wrote above only makes sense for 2nd order centered advection, both for velocity and tracers.
by stef
Tue Jun 05, 2007 7:36 pm
Forum: ROMS Information
Topic: Grid cell volume
Replies: 7
Views: 10126

Hi! I still use version 2.2, and don't know if anything has changed in the new version. In mod_grid.F there are some variables defined, three of them are: Hz: Thickness of vertical rho-Points Huon: Total u-Momentum Flux term: interpolated(Hz)*u*dy (dy is the same as 1/pn, at least for a rectilinear ...
by stef
Mon May 28, 2007 8:36 am
Forum: ROMS Documentation
Topic: Numerical Issues
Replies: 7
Views: 13016

I just reread the scrum manual and there it just says "finite diffences", not finite volume. I looked at the 2d-momentum equation and if I did everything right, then for a rectilinear grid the finite-volume and centered finite difference formulation is the same thing. But for a curvilinear...
by stef
Sun May 27, 2007 12:51 pm
Forum: ROMS Documentation
Topic: Numerical Issues
Replies: 7
Views: 13016

Hello! Thanks for the reply! I am also searching for information on ROMS finite volume formulations on the c-grid. Specifically, I have troubles understanding how the variables on the boundary of the gridcells are reconstructed to do the integration of the fluxes. It seems like there are many ways t...
by stef
Mon May 07, 2007 7:52 pm
Forum: ROMS Documentation
Topic: Numerical Issues
Replies: 7
Views: 13016

Hi! I was wondering if someone could give me hints about where to find documentation of the 3d momentum time-stepping (Adams-Bashforth 3) used in version 2.2. I tried to understand the Shchepetkin-McWilliams paper, and could extract some info about the barotropic LF-AM3, but I think they don't cover...
by stef
Tue Apr 10, 2007 3:53 pm
Forum: ROMS Usage
Topic: Coupling of m3-m2 initial conditions
Replies: 8
Views: 8230

I'm going to try that as soon as I can. My most stable setup so far is: Everything clamped in north, m3, m2 and tracer clamped in south, with free surface radiation.
by stef
Wed Apr 04, 2007 8:57 am
Forum: ROMS Usage
Topic: Coupling of m3-m2 initial conditions
Replies: 8
Views: 8230

Hi! Hmm, that's interesting. When I define m2 to be zero and define a m3 profile, e.g. u=w=0 and v=v(z), with all values negative from zeta down to the bottom (so vbar should be negative), the entry in the history files for the first calculated baroclinic timestep (NHIS=2) containes a vertical v pro...
by stef
Mon Apr 02, 2007 9:08 pm
Forum: ROMS Usage
Topic: Coupling of m3-m2 initial conditions
Replies: 8
Views: 8230

Thanks Kate, In my version (2.2) there is only one possibility to define 3d velocity initial conditions analytically: ANA_INITIAL. If I define it, I have to insert analytical initial fields for BOTH m2 and m3 velocities (and tracers). That's why I don't think that there is a call to any routine that...
by stef
Sat Mar 31, 2007 3:52 pm
Forum: ROMS Usage
Topic: Coupling of m3-m2 initial conditions
Replies: 8
Views: 8230

Coupling of m3-m2 initial conditions

My configuration is a flow-channel with open boundaries in north/south and closed walles in the east and west. I want to have velocity boundary conditions identical to velocity initial conditions. All of these conditions are analytically defined. If I want inflow from the north with u=w=0 and a vert...
by stef
Sat Mar 24, 2007 7:05 pm
Forum: ROMS Discussion
Topic: Simplest vertical mixing
Replies: 2
Views: 3526

I see. Thanks. Another question: In my mathematical model, I will have water flowing in a channel from north to south, with u=w=0 and a vertical shear in v (just like in the paper above). This implies a horizontal density gradient wich I want to impose by using a linear equation of state. A tilted f...
by stef
Fri Mar 23, 2007 5:04 pm
Forum: ROMS Documentation
Topic: Typos
Replies: 19
Views: 19235

In the vertical boundary conditions for z=zeta there should be the total time derivative of zeta for the w condition. If I understand correctly, the advective terms are missing.
by stef
Mon Mar 19, 2007 8:56 am
Forum: ROMS Discussion
Topic: Simplest vertical mixing
Replies: 2
Views: 3526

Simplest vertical mixing

I would like to set up a thermal wind balanced channel flow. The channel will have flat bottom for now, later I will put a cylinder or an island into the domain. I basically follow the setup in the article "Island Wakes in Deep Water" by Changming Dong. http://www.atmos.ucla.edu/~cdong/pap...