Ocean Modeling Discussion

I'm searching for a good conceptual description of the RadNud boundary condition, specifically for tracers and 3-D momentum. I understand the nudging inflow timescales set in the *.in file and how the boundary is updated. I'm having trouble with the concept of radiating out of the boundary and the ...

I'm using 3D data from RTOFS for my initial conditions file for my 1km nearshore Florida ROMS grid. Some of my ROMS grid is land masked. Zeta is set to 0 over the land mask as are 2D momentums. My confusion is on the 3D tracers and momentums. I have more resolution in the nearshore than is supplied ...

Hi Kate, This topic has always been one I'm curious about and still feel unsure with diagnosing problems. The descriptions you gave the person that started the post are clear, but what does one do with the information they find out with the tests you suggested? Is there a way to say definitively ...

I'm using RTOFS as my LBCs for a central east Florida shelf simulation and I'm curious about bias correcting RTOFS. RTOFS already assimilates some in situ data into it. I've seen a couple powerpoints on the web of users bias correcting HYCOM. How are users doing this and are there references for a ...

I've been running ROMS successfully with a 1km rectangular grid along the central east Florida shelf. I'm focusing on summer time cold water upwelling so I'm starting on January 1st and running through October 31st; a ~150 day spin up period before comparing to in situ data. I'm using NAM for ...

Hello Kate,

My rxo is 0.69 and my rx1 is 10.4. I tried to smooth with the smooth_bath Matlab file but it the values wouldn't go any lower. I have 78 in the I direction and 248 in the J and 20 vertical layers.

I did get the mex files compiled for seagrid.

I have my grid shown below of the central east Florida shelf. I'm trying to run 214 days and it runs 154 days and then becomes unstable and blows up in a few time steps. My 3d timestep is 150s and my 2d is 10s. When I check the restart file, the temperature and zeta values, at the Ci,Cj,Cw index ...

While running my simulation, two days into my run I receive the following error message. I have it set to run 183 days. SET_2DFLD - current model time exceeds ending value for variable: Pair TDAYS = 2.0023 Data Tmin = 0.0000 Data Tmax = 183.7500 Data Tstr = 0.0000 Data Tend = 0.2500 TINTRP1 = 2.0000 ...

I did originally have NaNs in my u and v velocities. I used the native netcdf tools in Matlab. I've since corrected my mistake and no longer have NaNs in my u and v in my initial conditions file. I looking into diag.F to determine where the problem could be on my end. Does the u and v used in diag.F ...

I solved most of my problems that I was having in this thread. Thank you for grep and print statements. I had a error in my zeta variable, so I no longer have a NaN for my POTEN_ENRG and NET_VOLUME. I still have NaNs for my KINETIC and TOTAL_ENRG and I traced it to the calculation of u2v2 in diag.F ...

I redesigned my grid to be rectangular, and I was able to compile all the mex files in seagrid successfully. My setup is the same as the posted .in file at the start of this thread. I did check my bathymetry for 0s and my land mask for NaNs. When i ran ROMS using oceanM it blew up at the same point ...

Thank you all for the input. Kate, I don't appear to have compiled the mex files properly for seagrid. I'm trying again but MATLAB continues to crash when it tries to call the mexrect command in the test_mexrect.m file. Mexinside and mexsepeli are compiling fine. Given my area of study, the Florida ...

The holes in the time step graphic are the locations of my mask_rho points. I'm not sure why the hslope m-file is configured to only apply the mask to that graphic.

I just used SeaGrid out of the box from the svn repository.

I will try a rectangular grid.

Thank you,

Chris

Here is the pm and pn figures. I don't know what to look for in these variables.





Chris

Thanks for the reply Kate, I've checked what I thought could be the reason for the NaNs. My initial conditions file doesn't have a mask variable, but there are no zeros or NaNs in the h variable. Same for h and hraw in my grid file. I checked the masks in my grid file an they were 0s and 1s where ...

Thank you for the input so far. The standing questions, where I'm stuck, are: why would the energies and net volume of the first output time step be NaNs and would the B&H stiffness of 0.85 be the reason for it blowing up at the first time step?

Below is my ROMS output log for my first attempts at running ROMS for a coarse grid of the central east Florida Shelf for my dissertation. I've smoothed my bathymetry using the smooth_bath MATLAB script and I made sure to not have any 0 values for my bathymetry in the grid file. Nothing jumps out at ...

I'm using RTOFS for my initial conditions and boundary conditions for a small grid along the central east Florida shelf. I'm looking for information on how others have handled the interpolation of RTOFS, or similarly HYCOM, mixed vertical coordinate to ROMS sigma coordinate. Thank you for your ...

I have successfully run the test cases in ROMS on my MacBook Air. Upon looking at the Activity Monitor for the OceanS file during the run, it didn't show up as a 64 bit process. My question is, is there a different setup required than that shown on the WikiROMS page to have ROMS run in 64 bit? I ...

Thank you for the information, the espresso was a nice example to see what they have used. My biggest problem is that there is a lot of software available and it's difficult to know what to use. I'm unfortunately the only ROMS user at my school. I will start with a few direct questions and that ...

I'm working on my dissertation to model the cold water upwelling that occurs most every summer along the central east Florida shelf. I have a long ADCP record in the nearshore and a few summers of bottom mounted ADCP data across the shelf to the shelf break, which I will use as a diagnostic ...

I'm using a Macbook Air with OS X 10.7.2, gfortran 4.6.2, GNU make version 3.8.1, and NetCDF version 3.6.3. I was able to successfully run the ROMS upwelling example by following the directions from the WikiROMS site under the tutorial of "Installing ROMS on a Macintosh". I haven't attempted to run ...

I'm using a Macbook Air with OS X 10.7.2, gfortran 4.6.2, GNU make version 3.8.1, and NetCDF version 3.6.3. I was able to successfully run the ROMS upwelling example by following the directions from the WikiROMS site under the tutorial of "Installing ROMS on a Macintosh". I haven't attempted to run ...