Ocean Modeling Discussion

LiuZHQ

The Southern University of Science and Technology (SUSTECH) in ShenZhen (next to Hong Kong) China has two immediate postdoc openings. We seek highly motivated postdoctoral researchers to join our team in the Department of Ocean Science and Engineering. The selected candidates are expected to 1 ...

LiuZHQ

Postdoc positions in Department of Ocean Science and Engineering, Southern University of Science and Technology Job Description The Southern University of Science and Technology has several immediate postdoc openings in coastal and shelf circulation modeling. We seek highly motivated postdoctoral ...

LiuZHQ

Job Description The Southern University of Science and Technology (SUSTC) has an immediate postdoc opening in coastal and shelf circulation modeling. We seek highly motivated postdoctoral researchers to join our ocean modeling team in the Department of Ocean Science and Engineering in Shenzhen (next ...

LiuZHQ

Wondering what version of ROMS you're using. You defined the OBC for three-dimensional velocity as Red, which supposedly represented a reduced-physics OBC. I checked the ROMS code (Version 645 2013-01-22 23:21:54Z), and didn't find that option. If the actual OBC is a closed one, due to, possibly ...

LiuZHQ

looks like you're adapting Flahter and Chapman to simulate the tidal current. To me, it is important that we provide tidal current (suggested by Kate) to this sort of OBC in Flather form (an "active" OBC in fact) to include both the tidal current and the tidal wave correctly. Or, if you insist, use ...

LiuZHQ

I mean, you may want to increase the length of your domain to the south by adapting more grids. So that, the southern boundary will be far away from the river mouth to prevent that contaminantion from Periodic OBCs.

regards,
ZQ

LiuZHQ

Since the plume propogates very fast along the coastline, and it can easily reach the southern boundary in several days. Perodic OBC brings the plume into your domain through northern boundary, which contaminant your solution. That's why I suggest you extend the domain to the south.

Regards,
ZQ

LiuZHQ

I think, that periodic OBC for southern and northern boundaries should be better than the others. You may want to extend your grid southward.

Regards,
ZQ

LiuZHQ

Hi liuy, Looks like you want to hide those corners with extremely large orthogonality error in the land points, but even though, the orthogonality error in the areas you interested in is still too large (~5 degree) to be acceptable. I suggest you start from a linear grid with gradually relaxiation ...

LiuZHQ

Hi there, Yes, everyone gets the liberty to organize their own study, and I'm not trying to offense you personally. In fact, although it may take a lot of effort to generate this grid, I'm not quite sure about the orthogonality, which is critically essential as one of the basic law of ROMS. You ...

LiuZHQ

It is not related to the grid resolution. To be frank, this grid againsted the basic law of ROMS, and it was wrong! Before you run the model, I suggest you find a book (e.g. Dale's book) or read the user manual provided by the wikiROMS, instead of taking the model as a black box and throw all those ...

LiuZHQ

Hi there,
kind of strange in your OBC configuration, why did you adapt chapman to zeta, Flather to the barotrophical part, closed for baroclinic part and clamped for T and S? this is totally controdiction.

Is there any reason? Or what kind of fluid you want to simulate?

regrads,
ZQ

LiuZHQ

Hi often, How about you compare the solutions of zeta_east, zeta_south,zeta_north in your boundary file with SODA? it there big difference? I think you included the tidal forcing in your case. Did the model you want to nest with included tide also? Or, according to your description, you nested with ...

LiuZHQ

Branching of the KC to Tsushima Strait is, nowadays, one of the major scientific questions related to interactions between the ECS shelves and the KC. It is not easy to answer that question: why Tsushima Warm Current and Taiwan Warm Current getting smaller during your modelling time. The way I will ...

LiuZHQ

Hi Zhao,
You can not expect that vertical temperature profile to be similar with the one illustrated in the WikiROMS, since you didn't plot it in the right way.

You should use h.*s_rho in MATLAB to find out the yaxis in the figure you demonstrated.

Good Luck

LiuZHQ

Hi there,
How about you initate the model run from rest?

LiuZHQ

Well, it's tricky to get the proper boundary solution. You can not use WOA as boundary forcing since it can definitely not resolve the find structure of the plume. Even worse, it is affected by many factors, not only the forcing you're considering here. Looks like you want to do some generalized ...

LiuZHQ

It looks like an OBC-related problem. Two suggestions: 1 (most likely). I noticed that you included the plume from PRE, and it looks like during the winter time. Frehser and colder plume arrives the open boundary and since you adapt M3RADIATION without nudging. Disturbances are reflected back into ...

LiuZHQ

That *_VOLCON options do nothing to your problem here, because it only conserves the inward and outward volume transport at the open boundaries. That non-conserved estuary circulation looks strange, and I never see that kind of problem in my idealized estuary cases with buoyancy forcing only. I ...

LiuZHQ

Hi Qianfang,
Please refer to this:
https://www.myroms.org/wiki/index.php/Tidal_Forcing
for more information.

In fact, if you download the data that contains the Harmonic Constants and the software from OTIS, you will find those files and know how to generate the tidal forcing.

happy modelling

LiuZHQ

Hi Shou, According to the figure and question "xiaozhu557" posted here. The so called "over-heated" coastal waters are his/her major concern. Yes, the waters are warmer here. But, the temperature is still in reasonable range, and this implies that the wrong thing is not that strongly related to the ...

LiuZHQ

Hi there, I don't think that over-heated temperature in Bohai Sea during the summertime is because of any possible problem in ROMS itself, and the wrong thing is your configuration of the model and forcing you adapted in this domain. Please re-review the related papers and you'll find the solution ...

LiuZHQ

Hi vviocean,

you may want to try Delft3D to generate the grid. Someone posted that topic.

LiuZHQ

please also post your configure file to help..

LiuZHQ

Did you smooth your bathymetry?

LiuZHQ

Was these figures illustrating your model results or just the forcing?

Did you also include the other forcings? wind, river, or remote forcing from the major currents? Kuroshio, NEC etc...

what is your horizontal mixing algorithm? harmonic or biharmonic, what is the coefficient?

ZQ

LiuZHQ

How about just calculate the barotropic PGF from the standard zeta output? I think it is better than changing the ROMS codes.
you can get baroclinic PGF by subtracting barotropic PGF from total PGF.

LiuZHQ

there is a free software call "VisualSVN", which can be utilized to download the source code under windows OS.

LiuZHQ

looks like the OBC didn't handle the coastal trapped waves well.
Since you seleted Flather, did you also include tidal forcing in your application?

LiuZHQ

leonjld, thanks for your suggestion, and I'll check it.

LiuZHQ

leonjld and kate, thanks for your reply. kate was right, and to me, it looks better if we can find some method to output it directly from the model. His file can not give us the information related to time, which means we can not obtain the terms expressed as parcial something parcial time. I'd ...

LiuZHQ

Well, the diagnostic terms involve a time-integral like the averages files. I usually recommend to have NAVG=NDIA to facilitate computing budgets and fluxes. The history fields are instantaneous :!: Therefore, the time assigned to the time-averaged fields is always at the middle of the time ...

LiuZHQ

Dear All, I noticed that the standard momentum terms output file or the "dia" file is following the type of "avg", which is averaged during the time we defined. My question is, can we also obtain the terms at a specific time, say the time of "his" for example? Or is there any toolbox which can ...

LiuZHQ

I wrote one program by myself to calculate the energy flux from velocity and elevation fields.
anyway, thanks for your attention all the same.
ZQ

LiuZHQ

Dear all,
Is it possible to also output the energy flux crossing the OBCs into the standard "his" files? I didn't find any related options in the "in" files, and the version I'm using is ROMS 3.4.

Thanks,
ZQ

LiuZHQ

I follow the question of Hetland, "what is your OBC?". EW_PERIODIC?
How about also try f plane, instead of beta plane?

LiuZHQ

Hi Bian,
I want to remind you that the grid num. of ROMS should be started from 0, instead of 1 you read in the matlab maxtrix (Array Editor), which means you may already put the river correctly this time, even it looks like on the "land".
regards,
ZQ

LiuZHQ

Hi Bian, From my point of view, the abnormal "salinity" problem is also largely related to the grid or the land-sea mask at the river points. Plz also try to modify your "mask" and make the river more "smooth" or you may want to further dig the river land-ward straightly, and artificially give more ...

LiuZHQ

Dear All, I'm doing some direct simulations of the circulations in western north Pacific Ocean using ROMS. After QuickSCAT retired, I can not find a very good wind speed/stress data to force my domain except some meteorologic model results, and some of my meteorology colleagues suggested the remote ...

LiuZHQ

Bian, plz also chang the "sustr:scale_factor = 0.0005 ;" and "svstr:scale_factor = 0.0005 ;" to "sustr:scale_factor = 1 ;" and "svstr:scale_factor = 1 ;" respectively, "svstr:units = "Newton meter-2 * 1000" ;" to "svstr:units = "Newton meter-2" ;" and then generate the new wind forcing files. good ...

LiuZHQ

Hi Bian, I remind you to recheck the "frc_uvstress.cdl" file before generating the forcing netcdf files. I noticed that there is some defination say "svstr:units = "Newton meter-2 * 1000;" in this file, which indicates that you should time the svstr with 1000 after converting from the wind speed V ...

LiuZHQ

Hi kate,
you are right, I forgot to difine it in the grid file
now, the "lon_rho" appeared.
thanks very much

LiuZHQ

Did you tell it you have a "spherical" grid? If spherical is false, it assumes you are on a flat idealized x,y plane. How did you provide the grid in the first place? Hi kate, I tried the including/excluding "spherical" from ".h" file cases and there are still not "lon_rho" in the out files. I also ...

LiuZHQ

They will be there in your output netcdf files, unless you define 'NO_WRITE_GRID' in your CPP configuration file. Hi feroda, thanks for your reply, I re-checked the "avg" and "his" files and found only the "x_rho" and "y_rho" instead of "lon_rho" and "lat_rho" which I needed, and I didn't define ...

LiuZHQ

Hi all,
Is there any swith cotrolling outputing the variables "lon_rho", "lat_rho" to every standard output "nc" files, or the "his" and "avg" files?
Thanks and any suggestion will be welcomed.

LiuZHQ

I got the same problem say "*** glibc detected *** ./oceanM: free(): invalid next size (normal): 0x00000000352ae2f0 ***
"
and I skip the new version of ROMS and back to the ROMS148, and the problem disappeared. This is a possible solution before anyone find out the real reason.

Search found 46 matches