Ocean Modeling Discussion

Wave-current interaction and wave effect on currents are the same concept. You are trying to use WEC_VF with wave forcing variables, which is essentially what SWAN will provide in a streamlined way. Coupling with SWAN may sound daunting but it's actually more straightforward for wind generated waves ...

In theory you can calculate Lwave based on the dispersion relation for surface gravity waves, but you would be missing a term of of k (dot) u for the Doppler effect because you have no currents. In SWAN, you may get bathymetric refraction (PROPTHETA) but not accurate frequency shift (PROPSIGMA), as ...

Just to confirm: you are using COAWST (e.g., https://github.com/DOI-USGS/COAWST) or ROMS (https://github.com/myroms)? It's probably easier to make a SWAN domain based on your ROMS grid, figure out a SWAN input file, and run the models coupled - that way you would have all the wave variables you need ...

One way to look at this is through your flags. Depending on your flags, forcing files may or may not be used. You have: #define BULK_FLUXES #ifdef BULK_FLUXES # define ANA_WINDS # define ANA_CLOUD # define ANA_TAIR # define ANA_HUMIDITY # define ANA_PAIR # define ATM_PRESS # define PRESS_COMPENSATE ...

I made a netdcf file from a .cdl template like this below. They're on rho points. Not sure how to treat 'time' here. What should I do so that I can use the forcing file akin to the analytical option (without having to change the forcing file if I change NTIMES in ocean.in) --- NetCDF-3 Classic atmos ...

Ok, I see. The pcolor is of Umag wind speed magnitude, though.

I didn't make this grid, don't know how the angle was computed. I got the grid .nc file and configured everything else from scratch.

thanks I defined CURVGRID and tried rotating the winds in ana_winds: I took a look at set_data.f and found how the clockwise rotation matrix is applied. Looks something like this: cff1=Uwind*CosAngler + Vwind*SinAngler cff2=Vwind*CosAngler - Uwind*SinAngler Uwind(i,j)=cff1 Vwind(i,j)=cff2 ---- So I ...

I have no curvgrid on the project header file. In ana_winds.h I have this (no rotation): Wmag=13.0_r8 DO j=JstrT,JendT DO i=IstrT,IendT Uwind(i,j)=0.0_r8 Vwind(i,j)=10.0_r8 END DO END DO In ana_smflux.h I have: DO j=JstrT,JendT DO i=IstrP,IendT sustr(i,j)=0.0_r8 svstr(i,j)=0.1_r8/rho0 END DO END DO ...

I'm trying to force a shallow bay with uniform, northward wind. When I use ana_smflux and prescribe svstr =0.1/rho0 and no eastward wind, I get this (ok) pattern: https://i.imgur.com/uwz0iVV.jpg?1 But when I use ana_winds and prescribe Vwind=10.0_r8 and no eastward wind, I get vectors that seem to ...

I would double check that the direction in wind_3h.dat is in the same coordinate system as when you prescribe WIND 5 45

I don't think ROMS/COAWST has a special treatment of breaking and non-breaking waves when it comes to mixing and TKE injection in the water column. Deep-water breaking (whitecapping) is a poorly understood phenomenon that remains a challenge to parameterize. I know of two options to account for wave ...

Sorry to necropost on this issue, but was this solved? I get the same compilation error when I #define WEC_MELLOR and try to compile. I'm using COAWST, though, so I'm not sure if this is the right place to ask. I got the code in August 2016 (last time I did svn checkout). I must say that the WEC_VF ...

As far as I know, the CFL criterion (or C) is calculated for u, v and w in ROMS. This is a somewhat old page on that topic but I think it's the same now? https://www.myroms.org/projects/src/ticket/486 ------------- C = Cu + Cv + Cw where Cu = ABS(u(i,j,k) + u(i+1,j,k)) * dt * pm(i,j) Cv = ABS(v(i,j ...

One way could be to decrease the value of Zob in ocean.in. Zos is the surface roughness length, and you wouldn't need to tweak that here. Another would be to decrease rdrg2, also in the ocean.in you show. And another way could be to prescribe a constant, lower, bottom drag coefficient to increase ...

Ideally, you might want to have winds in both ROMS and SWAN for the coupled runs. If you provide SWAN with wind forcing, ROMS will not know unless you create a ROMS wind forcing file too. That, or you can set a surface stress in ROMS using ANA_SMFLUX. Regarding the SWAN error, SWAN outputs PRINT and ...

I'm looking at wave-current stresses in an estuarine BBL. I've been using ANA_DRAG, UV_QDRAG, and UV_DRAG_GRID to set a constant drag coefficient everywhere (CD). Now I have these questions: 1 - I want to get the time and space-dependent drag coefficient (CD) as output in the history file (for ...