Search found 71 matches

by stef
Fri May 15, 2020 3:10 pm
Forum: ROMS Discussion
Topic: Instabilities in LF-AM3 caused by wetting/drying?
Replies: 4
Views: 227

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

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

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" Adams-...
by stef
Tue May 12, 2020 9:19 am
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 976

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

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

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 end of ste...
by stef
Wed Apr 22, 2020 12:42 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 976

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

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

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=GRID(ng)%on_u(Istr...
by stef
Fri Feb 14, 2020 11:01 am
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 976

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

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

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. Anyways, t...
by stef
Fri Feb 07, 2020 6:20 pm
Forum: ROMS Bugs
Topic: Possible bug in wetdry.F for nested apps
Replies: 14
Views: 976

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

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

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

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 multi-v...
by stef
Mon Jan 27, 2020 9:18 am
Forum: ROMS Problems
Topic: Three nested refined grids - Ngrids vs. Ncontact
Replies: 2
Views: 397

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

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

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

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

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

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

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

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)+ & - & (Sxy_psi(i+1,j)- &...
by stef
Tue Jun 11, 2019 2:00 am
Forum: ROMS Problems
Topic: Stokes velocity fills dry cell uphill
Replies: 3
Views: 515

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

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

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

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

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

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: 2
Views: 674

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: 2
Views: 674

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

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

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 tilted wi...
by stef
Mon May 29, 2017 3:12 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 4151

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

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

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

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

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

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

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

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 make sense to me :) Cheers, St...
by stef
Mon May 15, 2017 3:22 am
Forum: ROMS Discussion
Topic: Roms grid rotation and interpolation
Replies: 16
Views: 4151

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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 grid you ...
by stef
Sun May 27, 2007 12:51 pm
Forum: ROMS Documentation
Topic: Numerical Issues
Replies: 7
Views: 7426

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

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

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

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

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

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

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

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

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/papers/iswake...