Search found 12 matches

by lfparejar
Mon Jan 29, 2018 12:45 am
Forum: ROMS Discussion
Topic: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles
Replies: 8
Views: 1493

Re: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles

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...
by lfparejar
Sun Jan 28, 2018 11:31 pm
Forum: ROMS Discussion
Topic: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles
Replies: 8
Views: 1493

Re: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles

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.
by lfparejar
Sun Jan 28, 2018 5:09 pm
Forum: ROMS Discussion
Topic: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles
Replies: 8
Views: 1493

Re: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles

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 ...
by lfparejar
Sat Jan 27, 2018 10:54 pm
Forum: ROMS Discussion
Topic: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles
Replies: 8
Views: 1493

Re: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles

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 ...
by lfparejar
Fri Jan 26, 2018 7:53 pm
Forum: ROMS Discussion
Topic: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles
Replies: 8
Views: 1493

Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles

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...
by lfparejar
Fri Jan 26, 2018 2:02 am
Forum: ROMS Problems
Topic: the same wind but the wave height is different
Replies: 1
Views: 1001

Re: the same wind but the wave height is different

I would double check that the direction in wind_3h.dat is in the same coordinate system as when you prescribe WIND 5 45
by lfparejar
Tue Sep 05, 2017 1:51 pm
Forum: ROMS Discussion
Topic: About non-breaking wave enhanced vertical mixing
Replies: 2
Views: 978

Re: About non-breaking wave enhanced vertical mixing

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...
by lfparejar
Mon May 15, 2017 4:07 pm
Forum: ROMS Discussion
Topic: mellor instead of VF
Replies: 2
Views: 1204

Re: mellor instead of VF

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 ...
by lfparejar
Wed Apr 05, 2017 1:11 am
Forum: ROMS Information
Topic: stability criterion
Replies: 2
Views: 2306

Re: stability criterion

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,...
by lfparejar
Wed Apr 05, 2017 1:02 am
Forum: ROMS Installation
Topic: tide velocity problem
Replies: 3
Views: 2292

Re: tide velocity problem

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...
by lfparejar
Wed Apr 05, 2017 12:24 am
Forum: ROMS Installation
Topic: an error about coupled SWAN and ROMS model
Replies: 1
Views: 1615

Re: an error about coupled SWAN and ROMS model

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...
by lfparejar
Mon Apr 03, 2017 2:20 pm
Forum: ROMS Sediment
Topic: Wave-current BBL drag coefficient as output
Replies: 4
Views: 2712

Wave-current BBL drag coefficient as output

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