Search found 30 matches

by hpftcb
Tue Jul 13, 2021 9:58 pm
Forum: ROMS Discussion
Topic: Shortwave and Longwave radiation.
Replies: 6
Views: 4802

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...
by hpftcb
Sat May 29, 2021 1:50 pm
Forum: ROMS Discussion
Topic: Atmospheric forcing for child grid
Replies: 3
Views: 4338

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
by hpftcb
Mon Oct 26, 2020 12:57 pm
Forum: ROMS Discussion
Topic: 4DVAR codes with Octave
Replies: 0
Views: 978

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...
by hpftcb
Fri Jun 26, 2020 3:26 pm
Forum: ROMS Information
Topic: Curvilinear Coordinates
Replies: 13
Views: 38060

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
by hpftcb
Thu Apr 16, 2020 5:35 pm
Forum: ROMS Discussion
Topic: Roms not processing DSTART
Replies: 0
Views: 1430

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...
by hpftcb
Mon Apr 30, 2018 2:43 pm
Forum: ROMS Discussion
Topic: Sponge Layer visc_factor
Replies: 1
Views: 1716

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...
by hpftcb
Fri Sep 22, 2017 3:09 pm
Forum: ROMS Discussion
Topic: Climatology as open boundary
Replies: 6
Views: 3842

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
by hpftcb
Wed Sep 20, 2017 2:45 pm
Forum: ROMS Discussion
Topic: Climatology as open boundary
Replies: 6
Views: 3842

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...
by hpftcb
Wed Jun 07, 2017 3:33 pm
Forum: ROMS Discussion
Topic: Using Gradient as boundary condition
Replies: 0
Views: 1420

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...
by hpftcb
Wed May 24, 2017 1:58 pm
Forum: ROMS Bugs
Topic: Upwelling Channel with 4 Open Boundaries
Replies: 6
Views: 4191

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...
by hpftcb
Tue May 23, 2017 3:28 pm
Forum: ROMS Bugs
Topic: Upwelling Channel with 4 Open Boundaries
Replies: 6
Views: 4191

Re: Upwelling Channel with 4 Open Boundaries

Thank you for the replies.
I'll also try Gradient BC (Gradient = 0)

regards
by hpftcb
Mon May 22, 2017 5:15 pm
Forum: ROMS Bugs
Topic: Upwelling Channel with 4 Open Boundaries
Replies: 6
Views: 4191

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. ...
by hpftcb
Mon May 22, 2017 1:35 pm
Forum: ROMS Bugs
Topic: Upwelling Channel with 4 Open Boundaries
Replies: 6
Views: 4191

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, ...
by hpftcb
Sun Apr 30, 2017 3:08 pm
Forum: ROMS Discussion
Topic: why kinetic energy raise up?
Replies: 2
Views: 3121

Re: why kinetic energy raise up?

Hi,

did you find the answer for your question?

regards,

Fernando
by hpftcb
Tue Apr 11, 2017 4:00 pm
Forum: ROMS Discussion
Topic: Anonalous currents in shelf regions
Replies: 3
Views: 2314

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...
by hpftcb
Tue Apr 11, 2017 1:20 am
Forum: ROMS Discussion
Topic: Anonalous currents in shelf regions
Replies: 3
Views: 2314

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...
by hpftcb
Fri Mar 31, 2017 1:59 pm
Forum: ROMS Tools and Techniques
Topic: Pyroms remap problem
Replies: 5
Views: 4977

Re: Pyroms remap problem

my email is fernandotcbarreto@gmail.com. Contact me
by hpftcb
Thu Mar 30, 2017 7:30 pm
Forum: ROMS Tools and Techniques
Topic: Pyroms remap problem
Replies: 5
Views: 4977

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
by hpftcb
Mon Mar 13, 2017 2:59 pm
Forum: ROMS Tools and Techniques
Topic: one-way nesting
Replies: 4
Views: 4372

Re: one-way nesting

Thank you, Kate
by hpftcb
Sun Mar 12, 2017 2:14 am
Forum: ROMS Tools and Techniques
Topic: one-way nesting
Replies: 4
Views: 4372

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
by hpftcb
Fri Mar 10, 2017 10:31 pm
Forum: ROMS Problems
Topic: model blow-up due to abnormal velocity
Replies: 8
Views: 8365

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.
by hpftcb
Fri Mar 10, 2017 1:40 pm
Forum: ROMS Problems
Topic: model blow-up due to abnormal velocity
Replies: 8
Views: 8365

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
by hpftcb
Wed Nov 23, 2016 3:41 pm
Forum: ROMS Tools and Techniques
Topic: Pyroms remap problem
Replies: 5
Views: 4977

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...
by hpftcb
Tue Nov 22, 2016 9:24 pm
Forum: ROMS Discussion
Topic: Weights file
Replies: 8
Views: 7552

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.
by hpftcb
Mon Jun 06, 2016 11:25 pm
Forum: ROMS Tools and Techniques
Topic: PYROMS installation
Replies: 45
Views: 72889

Re: PYROMS installation

Hi maxblu33,
I am also having the same error in grid generation. Did you find out the error's source?
Regards
by hpftcb
Sat Jun 04, 2016 6:32 pm
Forum: ROMS Discussion
Topic: Hycom-Ncoda in PYROMS
Replies: 2
Views: 3522

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...
by hpftcb
Fri Jun 03, 2016 11:25 pm
Forum: ROMS Discussion
Topic: Hycom-Ncoda in PYROMS
Replies: 2
Views: 3522

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...
by hpftcb
Wed May 25, 2016 2:00 pm
Forum: ROMS Discussion
Topic: Weights file
Replies: 8
Views: 7552

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...
by hpftcb
Sun May 22, 2016 10:03 pm
Forum: ROMS Discussion
Topic: Weights file
Replies: 8
Views: 7552

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...
by hpftcb
Sun May 22, 2016 7:43 pm
Forum: ROMS Discussion
Topic: Weights file
Replies: 8
Views: 7552

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...