Ocean Modeling Discussion

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

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

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