Ocean Modeling Discussion

Thank you Parker! :-) I will try your idea using some history files. Just like you thought, I used <temp> and I calculated <v> (using <zeta>) both from the averages files. That is because I have a big domain, tidally forced and I don't have enough space storage to save a year of hourly history files ...

Hello everyone, I'm also interested in testing my tracer time rate of change. In particular, I'm trying to contrast the temperature time rate of change with the time derivative of temperature, after setting the time averaging window for both the diagnostic (NDIA) and average (NAVG) files to be the ...

Hi Rodrigo,

I think the main problem is not the proportionality (aspect ratio) between the coordinate axes. If that were the case, the appearance of the zeta contours would also be affected so that they could never be aligned with the velocity vectors.

Thank you for your comments and work! :-) I was doing some tests with vector plots over a map projection and I think could be useful to share it. I attach two Figures with the first areas that I showed (at the beginning of this post) be part of whole West Antarctic Peninsula domain. The left panels ...

OK Hernan, I made a new plot with barotropic velocities at every grid point. The vectors are plotted using "quiver.m" at the same fashion of " roms_quivergrd.m ", I included a vector's scale over the Antarctic Peninsula and I think vector 's tails are located at [lon,lat] grid cells.

I attached here a new plot that shows an example of model output's bottom velocities (subsampled each 2 grid points), including bathymetry contours every 100 m up to 500 m depth (dashed lines) and every 500 m up to 5000 m depth (solid lines). As Mike said, I would expect that velocities would follow ...

Thank you John, I appreciate your effort. I downloaded your code but the results are almost the same, the only differences are velocity magnitudes due to used precision at gravity of Earth. On the other hand, the complete domain of my application include West Antarctic Peninsula. In this sense, I ...

I'm agree the expression "rho2u_2d(1./pm)" is not exactly the same as "1./rho2u_2d(pm)" [which, according to your comments, seems to me to be equal to "om_u(i,j)=2.0_r8/(pm(i-1,j)+pm(i,j))" at metrics.F]. But that difference is quite small and could not explain the presented issue. I attached here ...

Thank you John by your feedback! I think we may do the overall calculus at different order, but the results will not be very different. By the way, I made de calculus as you said according the following MATLAB code lines: % Spatial derivates (adimensional) dzetadxi = diff(zeta,1,2) ./ rho2u_2d(1./pm ...

Hello guys, I'm trying to do some post processing of a curvilinear application, such as computing the surface geostrophic velocity. To do that, I'm derivating zeta field along it's cols and rows (which correspond to xi and eta coordinates), dividing by it's metric resolution (1/pm, 1/pn) and then ...

Oh, that sound completely new to me. If you access ERA5 datasets at the Climate Data Store of Copernicus, you will not see the concept of 'analysis' and 'forecast'. You need to request all variables, including radiation variables (fluxes), like 'analysis' fields. However, that doesn't mean all ...

Thank you John! I'm not sure to understand that answer pretty well. But, in my interpretation, the background reason to request different data sources (analysis vs forecast) would be related to store physical limitation at the NCAR. That is correct or have a sense for you? Cheers, David. Responses ...

Hello Mr. John, I hava a conceptual question about roms_get_era5_NCARds633_bulkflux.m script. In that file, why you use forecast fields to 'msdwlwrf','msnlwrf','msnswrf','mtpr' and analysis fields to 'msl','t2','d2','u10','v10'? Wouldn't it be better to just use analysis fields? Thanks in advance ...