Ocean Modeling Discussion

Hi, All, I was wondering if anyone could give me some suggestions on experiencing segmentation fault when reads in bry.nc. The boundary conditions I used is for applying zeta and ubar/vbar at South,North,East boundaries. My idealized rectangular domain is closed at west. I'm using the Chapman/Flathe...

Hi, I was wondering if anyone has the experience of imposing Holland wind field in ROMS. I have an idealized rectangular domain (250 km * 350 km, resolution = 2km) with tidal forcing propagating in from the east boundary and west boundary closed. I created a forcing netcdf file using Holland model a...

Hi, I was wondering if the Ocean Modeler position is still open? I didn't find any NOAA/ERT related jobs through the provided linkage. Has that already been filled?

Thanks and please let me know.

Good to know all of these. Thanks a lot, Kate. I masked some cells that blowing up to be dry and now it works fine. I am also trying to make a new grid just as you said making a transition from the shallow to deep. But what I have found that is the cell during this transition portion is very weird, ...

Hi, I have tried with BODYFORCE, but there is no much improvement. The model still blows up either at western boundary or at the junction between shallow and deep water, due to the speed trouble. Now I am testing svstr = 0.5N/m^2. And the current velocity along y_rho direction often gets into a very...

I got you. Thank you, Kate. Let me play with it and see how thing's going.

And now I am trying BODYFORCE. In order to do that, I was wondering if the following steps are correct: 1) define BODYFORCE in cpp option 2) Set LEVSFRC = 8 in ocean.in Besides these two steps, I was wondering if I need to generate a new netcdf file for the wind? Also, I was confused about the index...

zeta along the west boundary is no more than 1.6e-04 m. What happened is the V-component is too high (over 35 m/s) along the west boundary. If Coriolis effect is not considered. Everything went fine. However, if f is not zero, the cell gets blow up.

Hi, Kate, Here I attached my ideal study domain. The point(1,1,8) is located at boundary of the shallow water. But I was wondering if grid cell (1,1) is ghost cell or not? Basically when I create the domain, I didn't mask any cells dry because it is a simple test for an assumption in our real case s...

Thank you, Kate. I adjusted svstr(i,j)=0.1_r8, however, it only lasted 2 steps more. STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME C => (i,j,k) Cu Cv Cw Max Speed 0 0 00:00:00 0.000000E+00 4.861661E+03 4.861661E+03 6.219716E+13 (000,000,0) 0.000000E+00 0.000000E+00 0.000000E+00 0.0...

Hi, Dear All, During my simulation, I encountered a new problem when trying to add constant winds through analytical function in ROMS. I ran my simulations for five days with 1 min time step, however, after I define svstr=10 in the enana_smflux.h, I encountered such problem: NL ROMS/TOMS: started ti...

Thanks, Dr.Kate.

Thank you, Dr.John.Now it works. I have two additional question related to this ideal simulation: I found a header file named as ana_m2obc.h, which looks to me that it also deals with O.B.C. Since in my test case, the M2 component is imposed on the EAST boundary, I thought it is necessary to call th...

Hi,All, I am now using the COAWST v3.1. I tried to input an analytical M2 tide through the east open boundary conditions for an ideal test case. My domain is rectangular. What I did is add the following sentence into ana_fsobc.h and ana_m2obc.h: # elif defined IDEAL fac=1.0_r8 ! zeta0 omega=2.0_r8*p...

Thanks, John. Since I am not specifying OUTPUT in swan.in, only results from last time step is expected. If SWAN blows up in the middle, unlike in ROMS where the last time step can be saved in rst.nc, nothing can be printed out ?

I see. I actually have another simulation with ROMS and SWAN coupled. I used the following sentence, but nothing been generated from SWAN. BLOCK 'COMPGRID' NOHEADER 'rtp.mat' LAY 4 RTP 1. BLOCK 'COMPGRID' NOHEADER 'depth.mat' LAY 4 DEPTH 1. BLOCK 'COMPGRID' NOHEADER 'vel.mat' LAY 4 VEL 1. BLOCK 'COM...

Hi, did anyone encounter issue with opening the .mat files generated by SWAN during the ROMS-SWAN coupling process? I used the following script in swan.in BLOCK 'COMPGRID' NOHEADER 'rtp.mat' LAY 4 RTP 1. OUTPUT 20110822.000000 30 SEC BLOCK 'COMPGRID' NOHEADER 'depth.mat' LAY 4 DEPTH 1. OUTPUT 201108...

Thank you. Let me check the wave length and wave peak from ROMS first.

Thanks, John. I looked into the code, and found the basic differences did exist in whether Pwave_top are needed or not. However, the Pwave_top is provided by the SWAN. In the ROMS code, there are two places when Pwave shows up besides the Bulk_flux.f90 and they are either import the value from SWAN ...

Dear all, I tried to use ROMS coupled with SWAN to simulate a real hurricane. My problem is by using Taylor and Yelland option to calculate the sea surface roughness, my simulation went all well. However, by using the OOST option, the simulation quickly blew up after running about 1 day. My time ste...