Search found 30 matches
- Tue Jul 13, 2021 9:58 pm
- Forum: ROMS Discussion
- Topic: Shortwave and Longwave radiation.
- Replies: 6
- Views: 5284
Re: Shortwave and Longwave radiation.
Hi Rafael, posting the cpp flags I usually use in case it helps. One of the times I had problem with a colder surface temperature was due to humidity (wrong unit). /* momentum equations */ #define UV_ADV #define UV_COR #define UV_QDRAG #define UV_U3HADVECTION #define UV_C4VADVECTION #define TS_U3HAD...
- Sat May 29, 2021 1:50 pm
- Forum: ROMS Discussion
- Topic: Atmospheric forcing for child grid
- Replies: 3
- Views: 5330
Re: Atmospheric forcing for child grid
Hi Ash,
From the error it does not seem related to the forcing file, as ROMS performs the interpolation internally. Have you checked if the dimensions in the .in file from the son grid are correct?
Regards
From the error it does not seem related to the forcing file, as ROMS performs the interpolation internally. Have you checked if the dimensions in the .in file from the son grid are correct?
Regards
- Mon Oct 26, 2020 12:57 pm
- Forum: ROMS Discussion
- Topic: 4DVAR codes with Octave
- Replies: 0
- Views: 1145
4DVAR codes with Octave
Hi all, I am trying to run matlab 4dvars routines (https://www.myroms.org/svn/src/matlab/4dvar/) using Octave. I had sucess until line '[status]=c_observations(S)' ( a function from c_observations.m file). The problem is that this function does not have the option to not use MEXNC like the other one...
- Fri Jun 26, 2020 3:26 pm
- Forum: ROMS Information
- Topic: Curvilinear Coordinates
- Replies: 13
- Views: 39598
Re: Curvilinear Coordinates
You are right Stef, you can see in the obc_roms2roms.m code % Receiver Grid Orientation: rotate from TRUE East/North to % (XI,ETA) coordinates. B.(ufield) = u .* cos(angle.(edge)) + ... v .* sin(angle.(edge)); B.(vfield) = v .* cos(angle.(edge)) - ... u .* sin(angle.(edge)); Regards
- Thu Apr 16, 2020 5:35 pm
- Forum: ROMS Discussion
- Topic: Roms not processing DSTART
- Replies: 0
- Views: 1594
Roms not processing DSTART
Hi all, I am using ROMS3.8 and so far I have been setting my time with time_ref only (DSTART=0). After changing the value of DSTART I realized that Roms is not processing this paremeter, using only time_ref for initialization. As you can check in the log file. Has anyone found the same pattern? Rega...
- Mon Apr 30, 2018 2:43 pm
- Forum: ROMS Discussion
- Topic: Sponge Layer visc_factor
- Replies: 1
- Views: 1774
Sponge Layer visc_factor
Hello, I'm trying to set a sponge layer in roms3.6 as described in this post https://www.myroms.org/projects/src/ticket/627. I have set LuvSponge == T ! horizontal momentum LtracerSponge == T T T ! temperature, salinity, inert The problem is Roms isn't reading visc_factor and diff_factor from my gri...
- Fri Sep 22, 2017 3:09 pm
- Forum: ROMS Discussion
- Topic: Climatology as open boundary
- Replies: 6
- Views: 3958
Re: Climatology as open boundary
Thank you for the replies.
I tried to put Rad instead of RadNud in order to use clm as boundary conditions, but it still asked for bry file. In there a way to use only clm file as boundary condition.
Thank you
I tried to put Rad instead of RadNud in order to use clm as boundary conditions, but it still asked for bry file. In there a way to use only clm file as boundary condition.
Thank you
- Wed Sep 20, 2017 2:45 pm
- Forum: ROMS Discussion
- Topic: Climatology as open boundary
- Replies: 6
- Views: 3958
Climatology as open boundary
Hi all, I'd like to know how to use only climatology file as open boundary (not nudging, as boundary conditions). Do I have to create a bry file too? If I have a bry and a clm file, which one is used as open boundary? My BC are: LBC(isFsur) == Clo Che Che Che ! free-surface LBC(isUbar) == Clo Shc Sh...
- Wed Jun 07, 2017 3:33 pm
- Forum: ROMS Discussion
- Topic: Using Gradient as boundary condition
- Replies: 0
- Views: 1498
Using Gradient as boundary condition
Hi all, I'm trying to use gradient as BC in the following configuration: LBC(isFsur) == Che Che Che Che ! free-surface LBC(isUbar) == Gra Gra Gra Gra ! 2D U-momentum LBC(isVbar) == Gra Gra Gra Gra ! 2D V-momentum LBC(isUvel) == Gra Gra Gra Gra ! 3D U-momentum LBC(isVvel) == Gra Gra Gra Gra ! 3D V-mo...
- Wed May 24, 2017 1:58 pm
- Forum: ROMS Bugs
- Topic: Upwelling Channel with 4 Open Boundaries
- Replies: 6
- Views: 4390
Re: Upwelling Channel with 4 Open Boundaries
I tried with grad = 0 in all boundaries and the domain energy kept rising (velocity rising) until my case blew up. Does anyone know if what I did is right, if I can't use grad to termohaline or something like that. Thank you in advance LBC(isFsur) == Gra Gra Gra Gra ! free-surface LBC(isUbar) == Gra...
- Tue May 23, 2017 3:28 pm
- Forum: ROMS Bugs
- Topic: Upwelling Channel with 4 Open Boundaries
- Replies: 6
- Views: 4390
Re: Upwelling Channel with 4 Open Boundaries
Thank you for the replies.
I'll also try Gradient BC (Gradient = 0)
regards
I'll also try Gradient BC (Gradient = 0)
regards
- Mon May 22, 2017 5:15 pm
- Forum: ROMS Bugs
- Topic: Upwelling Channel with 4 Open Boundaries
- Replies: 6
- Views: 4390
Re: Upwelling Channel with 4 Open Boundaries
Hi, thanks for the reply Yes,It's an idealized experiment (bathymetry is attached). It's a 3D case with 30 sigma layers (vt=2,vs=4), with salinity and temp field build from one WOA point (uniform horizontally), zero current velocity in all boundaries, and 5 m/s eastern wind input with ana_smflux.h. ...
- Mon May 22, 2017 1:35 pm
- Forum: ROMS Bugs
- Topic: Upwelling Channel with 4 Open Boundaries
- Replies: 6
- Views: 4390
Upwelling Channel with 4 Open Boundaries
Hi all, I'D like to know if there is a way to set up an upwelling channel case with four open boundaries (on the contrary of Roms' upwelling test case) and zero inflow and outflow current. I tried the following configuration, yet the wind generated current did not manage to radiate from the domain, ...
- Sun Apr 30, 2017 3:08 pm
- Forum: ROMS Discussion
- Topic: why kinetic energy raise up?
- Replies: 2
- Views: 3582
Re: why kinetic energy raise up?
Hi,
did you find the answer for your question?
regards,
Fernando
did you find the answer for your question?
regards,
Fernando
- Tue Apr 11, 2017 4:00 pm
- Forum: ROMS Discussion
- Topic: Anonalous currents in shelf regions
- Replies: 3
- Views: 2371
Re: Anonalous currents in shelf regions
Hi Kate, thanks for your reply It's a baroclinic case (although I found the same problem when testing the barotropic mode). I aim to build an idealized case to study baroclinic instability problems. In this specific case the BC e IC where built from a parametrization of the Brazil Current with the T...
- Tue Apr 11, 2017 1:20 am
- Forum: ROMS Discussion
- Topic: Anonalous currents in shelf regions
- Replies: 3
- Views: 2371
Anonalous currents in shelf regions
Hi Roms forum, I'm running a case in which my boundary conditions were built from my initial condition (initial_condition.png). It's a gausian current flowing along shelf break with 0 velocity elsewhere. The problem is I'm getting strong current velocities in the shallowest parts of my region (in th...
- Fri Mar 31, 2017 1:59 pm
- Forum: ROMS Tools and Techniques
- Topic: Pyroms remap problem
- Replies: 5
- Views: 5226
Re: Pyroms remap problem
my email is fernandotcbarreto@gmail.com. Contact me
- Thu Mar 30, 2017 7:30 pm
- Forum: ROMS Tools and Techniques
- Topic: Pyroms remap problem
- Replies: 5
- Views: 5226
Re: Pyroms remap problem
Hi Luis,
I didn't, sorry. If you are interested I can send you the code I'm using for interpolation hycom2roms.
Regards
I didn't, sorry. If you are interested I can send you the code I'm using for interpolation hycom2roms.
Regards
- Mon Mar 13, 2017 2:59 pm
- Forum: ROMS Tools and Techniques
- Topic: one-way nesting
- Replies: 4
- Views: 4525
Re: one-way nesting
Thank you, Kate
- Sun Mar 12, 2017 2:14 am
- Forum: ROMS Tools and Techniques
- Topic: one-way nesting
- Replies: 4
- Views: 4525
Re: one-way nesting
Hi,
I would appreciate a lot to know the answer to this question. For a one-way nesting I compile with #define ONE-WAY?
regards,
Fernando Barreto
I would appreciate a lot to know the answer to this question. For a one-way nesting I compile with #define ONE-WAY?
regards,
Fernando Barreto
- Fri Mar 10, 2017 10:31 pm
- Forum: ROMS Problems
- Topic: model blow-up due to abnormal velocity
- Replies: 8
- Views: 8607
Re: model blow-up due to abnormal velocity
Hi Kate,
no, I'm not using WetDry and my shelf region has very shallow parts. I am gonna try that.
Thank you.
no, I'm not using WetDry and my shelf region has very shallow parts. I am gonna try that.
Thank you.
- Fri Mar 10, 2017 1:40 pm
- Forum: ROMS Problems
- Topic: model blow-up due to abnormal velocity
- Replies: 8
- Views: 8607
Re: model blow-up due to abnormal velocity
Hi, I have the same problem with high velocity near the coast (resulting in model blow up). Did you solve yours?
Regards,
Fernando
Regards,
Fernando
- Wed Nov 23, 2016 3:41 pm
- Forum: ROMS Tools and Techniques
- Topic: Pyroms remap problem
- Replies: 5
- Views: 5226
Pyroms remap problem
Hi all, I'm using pyroms/Palau_hycom to generate .bry .ini files to my grid (Hycom data were downloaded through get_hycom_GLBa0.08..py ). The problem is, during remap (dst_varz = pyroms.remapping.remap(src_varz, wts_file, spval=spval) line in remap_bry.py), in which hycom data is horizontally interp...
- Tue Nov 22, 2016 9:24 pm
- Forum: ROMS Discussion
- Topic: Weights file
- Replies: 8
- Views: 8161
Re: Weights file
Hi João,
thanks a lot for your help. Sorry for the late reply, I was working in other projects and completely forgot the thread.
You are completely right,
Thank you
Regards,
Fernando.
thanks a lot for your help. Sorry for the late reply, I was working in other projects and completely forgot the thread.
You are completely right,
Thank you
Regards,
Fernando.
- Mon Jun 06, 2016 11:25 pm
- Forum: ROMS Tools and Techniques
- Topic: PYROMS installation
- Replies: 45
- Views: 86230
Re: PYROMS installation
Hi maxblu33,
I am also having the same error in grid generation. Did you find out the error's source?
Regards
I am also having the same error in grid generation. Did you find out the error's source?
Regards
- Sat Jun 04, 2016 6:32 pm
- Forum: ROMS Discussion
- Topic: Hycom-Ncoda in PYROMS
- Replies: 2
- Views: 3586
Re: Hycom-Ncoda in PYROMS
Thank you for your reply. I've downloaded hycom-Ncoda data in .nc format for an area bigger than my roms area. Pyroms is reading and interpolating, the problem is that the output data is wrong, it is like the program is not interpolating for the right points. At least until after flooding, everythin...
- Fri Jun 03, 2016 11:25 pm
- Forum: ROMS Discussion
- Topic: Hycom-Ncoda in PYROMS
- Replies: 2
- Views: 3586
Hycom-Ncoda in PYROMS
Hi all, I am trying to generate boundary and initial conditions to ROMS from Hycom-Ncoda Analysis-Reanalysis through PYROMS. I am using as background examples/Arctic-HYCOM. I have installed everything correctly (I think :D),so I reckon my problem is not related with Pyroms installation. The thing is...
- Wed May 25, 2016 2:00 pm
- Forum: ROMS Discussion
- Topic: Weights file
- Replies: 8
- Views: 8161
Re: Weights file
Hi Matheus, thanks a lot for your reply. After checking the flood_fast.py file In realized that 4 2d arrays are generated in the section, MASK, C1, C2, C3. The shape of C2 and C3 is 119-126 (equal Hycom's), although the shape of MASK and C1 is 117-124, which causes the error. I set MASK and C1 to C2...
- Sun May 22, 2016 10:03 pm
- Forum: ROMS Discussion
- Topic: Weights file
- Replies: 8
- Views: 8161
Re: Weights file
Thank you, I managed to generate the weight files. During make_bdry I got the following error in the temperature part: ''ValueError: operands could not be broadcast together with shapes (117,124) (119,126)'' I tracked down the error as belonging to flood_fast.py in the following part: if kk == 0: c3...
- Sun May 22, 2016 7:43 pm
- Forum: ROMS Discussion
- Topic: Weights file
- Replies: 8
- Views: 8161
Weights file
Hi all, I am trying to generate boundary file from HYCOM-NCODA data using as background pyroms/examples. My grid was generate through Seagrid. The problem is I get an error while this process, which I think is probably related with the weight file ("wts_file"). The question: how do I get t...