Ocean Modeling Discussion

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

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

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

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

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

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

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

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

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

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

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

regards

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

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

Hi,

did you find the answer for your question?

regards,

Fernando

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

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

my email is fernandotcbarreto@gmail.com. Contact me

Hi Luis,
I didn't, sorry. If you are interested I can send you the code I'm using for interpolation hycom2roms.

Regards

Thank you, Kate

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

Hi Kate,

no, I'm not using WetDry and my shelf region has very shallow parts. I am gonna try that.

Thank you.

Hi, I have the same problem with high velocity near the coast (resulting in model blow up). Did you solve yours?

Regards,

Fernando

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

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.

Hi maxblu33,
I am also having the same error in grid generation. Did you find out the error's source?
Regards

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

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

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

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

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