Ocean Modeling Discussion

langlangxie

i think it works without adding the river. give it a try. Thanks. Sorry, I am still confused. I think if I add the river as the point source, it needs to switch on the SAND_Ltsrc in the sediment.in. And then a river netcdf file should be also set to SSFNAME (i.e. SSFNAME == my_rivers.nc) in the ...

langlangxie

I have not used this in a while, but there is a LtracerSrc option. Can you look in the ocean.in for this? There is a description at the bottom of the file. -j Thanks. I looked in the description of LtracerSrc in the ocean.in. It seems that it needs to add river runoff as a point source. However ...

langlangxie

Hi, I am running ROMS model with a sediment component. I would like to implement one point source of sediment at different level in the ocean. There is a post similar to my question. The post address is “https://www.myroms.org/forum/viewtopic.php?f=20&t=3787”. The author of this post considered the p ...

langlangxie

Hi,everyone. I have had some trouble in my ROMS. Look :NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00 (Grid 01, File: roms_ini_Y2001M1.nc, Rec=0001, Index=1) - free-surface (Min = 3.77269772E-01 Max = 9.62686659E-01) - vertically integrated u-momentum component (Min = -2.90209194E-01 ...

langlangxie

rsoutelino wrote:Did you try to increase the stack limit in the machine you're running? Just a guess, but since you decreased the time step, there is more memory being used in between netcdf time records.
Code: Select all
ulimit -s unlimited

Thanks,can you tell me how to find ulimit.I change my NRST,NAVG and so on.
thanks again!

langlangxie

hello,everyone! I have a trouble in my roms.in. the first version: NTIMES == 1036800 ！50 years DT == 600.0d0 NDTFAST == 180 NRREC == 0 LcycleRST == T NRST == 51840 NSTA == 1 NFLT == 1 NINFO == 51840 LDEFOUT == T NHIS == 0 NDEFHIS == 0 NTSAVG == 1 NAVG == 4320 NDEFAVG == 51840 NTSDIA == 1 NDIA == 0 N ...

langlangxie

Did you just update to a new version of ROMS? If that's the situation, you can read this thread[https://www.myroms.org/forum/viewtopic.php?f=14&t=2511]. This happened to me before, but my simulation broken at the very beginning. So before updating your compiler, you can turn off the average output ...

langlangxie

after some seconds,the ROMS is broken.the problem is: Lateral Boundary Conditions: NLM ============================ Variable Grid West Edge South Edge East Edge North Edge --------- ---- ---------- ---------- ---------- ---------- zeta 1 Chapman Imp Chapman Imp Chapman Imp Chapman Imp ubar 1 Flather ...

langlangxie

So you are saying: * With QCORRECTION, SST looks good, SSH looks bad. * Without QCORRECTION, SST looks bad, SSH looks good. Is that right? I don't know how that could be. kate,I don't mean 'Without QCORRECTION, SST looks bad, SSH looks good'.without QCORRECTION,the SST and SSH are all bad,but if ...

langlangxie

Is this a surface heating problem or a boundary condition problem? How are you providing the surface heat fluxes? You've commented out BULK_FLUX, so you must be providing some heat flux from outside. If so, you should also provide the dQdSST term. Hello,kate. I made a forced field data by agr ...

langlangxie

Dear all I have ran into a problem.My result of Roms is abnormal,the temp is more high in the northern of my areas. My test_south.in is : ! Keyword Lateral Boundary Condition Type ! ! Cha Chapman ! Cla Clamped ! Clo Closed ! Fla Flather _____N_____ j=Mm ! Gra Gradient | 4 | ! Nes Nested | | ! Nud ...

langlangxie

I recently had some zetas get too negative because my air pressure units were off. ROMS wants air pressures in millibars, so around a value of 1000. Or you could post your ROMS output to see if anyone sees anything obvious in there - but I'm going on vacation for a few days. thanks kate,I will try ...

langlangxie

It is very hard to understand what you mean. Although SSH is more relevant to your netcdf driving file than to .in file, you might want to check your driving force. I'm sorry to reply to you too lately .I mean that I think my ssh which is wrong,so I want to find a way to solve my error.I check my ...

langlangxie

This my ssh picture,in the picture ,ssh is -8.32 - -7.32.The phenomenon shows sea water is only ouput without into in the area.And I find the distribution is wrong ,the high ssh should be around Luzon Strait,but not in the northern of south china sea. I can't find the error in my roms.in,so I want ...

langlangxie

This is worrisome: GET_NGFLD - 2D u-momentum western boundary condition, t = 15 00:00:00 (Rec=0000001, Index=1, File: scs_bry.nc) (Tmin= 15.0000 Tmax= 345.0000) (Min = NaN Max = NaN) GET_NGFLD - 2D v-momentum western boundary condition, t = 15 00:00:00 (Rec=0000001, Index=1, File: scs_bry.nc) (Tmin ...

langlangxie

when I are running ROMS,I run into a problem.follows: Basin information for Grid 01: Maximum grid stiffness ratios: rx0 = 2.071374E-01 (Beckmann and Haidvogel) rx1 = 5.425647E+00 (Haney) Initial basin volumes: TotVolume = 6.3267250293E+15 m3 MinVolume = 8.3358286833E+08 m3 MaxVolume = 3.2792063556E ...

langlangxie

Are you sure it's blowing up on the boundary? Have you looked at the boundary conditions you are providing for it? Have you tried without the volume conservation? Were you running something similar before with more success? What did you change this time? Sorry,I don't say totally,before I change ...

langlangxie

I want to make my area are all open,this is VolCons(west) == T VolCons(east) == T VolCons(south) == T VolCons(north) == T Then I should change lateral boundary condition type: ! ! Cha Chapman ! Cla Clamped ! Clo Closed ! Fla Flather _____N_____ j=Mm ! Gra Gradient | 4 | ! Nes Nested | | ! Nud ...

langlangxie

kate wrote:I wrote about segmentation faults here. ROMS allocates big arrays right after printing the cppdefs.

The run with T=100 must be some other run? Unless that's your initial condition?

Thanks a lot,kate!

langlangxie

Hi,can you help me? I make my boundary all is T,but my roms can't running.the error is South China Sea model Operating system : Linux CPU/hardware : x86_64 Compiler system : pgi Compiler command : /dcfs2/program/mpi/mvapich2/1.8/pgcc.pgf90/bin/mpif90 Compiler flags : -O3 -Mfree Input Script : SVN ...

langlangxie

Can you help me to determine the "Open boundaries switches" Following obc = [1 1 1 1]; % open boundaries (1=open , [S E N W])(using agrifroms) My ROMS area is 104°E-120°E，2°N-24°N I make obc = [1 1 0 1],that is wrong. ??? Error using ==> geost_currents_bry at 261 Replace your South boundary Error in ...

langlangxie

I recently heard a quote that goes something like "I used to worry about that, but now I don't". I have also tried ndtfast=30, but when I need to shorten my main dt for stability, sometimes I reduce ndtfast. Very rarely (never?) is my model blowing up because my dtfast is too long - it finds many ...

langlangxie

You're asking how to pick ndtfast? I usually have it at 20 since my baroclinic timesteps are short enough for this to be stable. I could probably get away with smaller values, but haven't really pushed it. Yes,I am asking help about how to pick Ndtfast.Kate,you choose ndtfast=20 in your every case ...

langlangxie

Hi
I know dt=Ndtfast*dtfast,and dafast depends on the maximum depth and grid resolution.what's more,dtfast is limited by the barotropic CFL criterion.
So ,I want to know how to determine dafast.
Can you help me !

langlangxie

Hi, You can take a look at this: https://www.myroms.org/forum/viewtopic.php?f=19&t=2512&hilit=dtfast to find out about the relationship between ntimes, dt and ndtfast. ntimes is the total number of timesteps so if you have for example: NTIMES=1440 DT=300 NRST,NHIS,...=288 then your model will run ...

langlangxie

I told you to try NINFO=1. With dt=60 and that many steps, you are asking it to run for 36 days, saving an averages file every 3 days. You shouldn't get 20 averages files, just 12. Right? I've never heard of anyone using ndtfast as large as 600. You must have told it that your initial file is a ...

langlangxie

Are you changing NAVG as well as NINFO? You should read up on these parameters - they are described briefly at the bottom of the ocean.in file. As for ROMS blowing up, there are a few tips on the wiki. Dear kate Firstly,thanks again! My NAVG is not equal to NINFO.NAVG=4320,and NINFO=51840. When I ...

langlangxie

Dears,
I doubt there is a big relationship between some variables,such as ntimes,ndtfast,nhis,navg,ndia and nrst.
Can you give me some explainings?
Thank you very much!

langlangxie

The best value for NINFO is 1 (as above). That way, ROMS will check for NaN and other such values every time-step. Your model is going bad and you need to figure out where and why. Thanks! I am finding my errors in my ROMS.can you give me some directions of going bad? Today,I get a result ,and I ...

langlangxie

Zeta is the free surface elevation and it evolves during a model run. It would be in your ini file and also in the restart file, plus whatever other outputs you've asked for. I suggest ncview to look at the netcdf input and output files. Thanks for your help! Firstly ,I got a error is "WRT_AVG ...

langlangxie

Have you set NINFO to 1? If you do, you will get a diagnostics report each timestep. If the model is blowing up, you would then see it there first. If NINFO is large, you can get a NaN and ROMS will keep running. I believe you have NaN or some such number in your zeta field, which is why it won't ...

langlangxie

kate wrote:NINFO is set in your ocean_xx.in file. Try setting it to 1. If you can get things to run stably, then you can set it to a large value again if you like. I always leave it at 1 and put up with the extra output.

Thanks a lot. Thanks a million!

langlangxie

Have you set NINFO to 1? If you do, you will get a diagnostics report each timestep. If the model is blowing up, you would then see it there first. If NINFO is large, you can get a NaN and ROMS will keep running. I believe you have NaN or some such number in your zeta field, which is why it won't ...

langlangxie

Hi,I have a problem in computing ROMS. This is : WRT_AVG - error while writing variable: zeta into averages NetCDF file for time record: 1 Elapsed CPU time (seconds): Node # 0 CPU: 1009.411 Node # 4 CPU: 1009.387 ROMS/TOMS - Output NetCDF summary for Grid 01: Node # 5 CPU: 1009.386 Node # 10 CPU ...

langlangxie

kapitzah wrote:The argument to option -np seems to have a misplaced '-character. It is

-np 'cat${PBS_NODEFILE}|' wc -l

while it should be something like

-np `cat ${PBS_NODEFILE}| wc -l`

Please note the backticks ` instead of '

Thank you very much!

langlangxie

Still, how far does it get? What is the last thing ROMS outputs? No time. + NSLOTS='cat${PBS_NODEFILE}| wc -l' + echo 'This jobs is 4681.node200@general' This jobs is 4681.node200@general + cd /dcfs2/users/langlangxie/Projects/upwelling + /dcfs2/program/mpi/mvapich2/1.8/pgcc.pgf90/bin/mpirun_rsh ...

langlangxie

kate wrote:You're not giving us much to go on. How does it not compute? Can you provide the error output from ROMS or did you not get any?

I don't get any error.The ROMS just stop computing.
Thanks for your help!

langlangxie

Now,I want to use ROMS to make a model of South China Sea.First,I product some files in using agrif_roms.such as *_bry.nc,*_ini.nc,*_grd.nc,*_frc.nc.*_clm.nc,then I take these files into my project,and compile ROMS,lastly,I put these into cluster computing.In processing,I'm sure compiling and Job ...

langlangxie

kate wrote:Hmm, I don't know that one. Here's what I wrote about ROMS blowing up.

Thanks your help!

langlangxie

Could you please be more specific about what you are trying to do and what is going wrong? You don't need all those files for the upwelling case. Have you looked through the documentation on the wiki? There is also a manual. Thanks! I am now doing a model of the South China Sea in this case on the ...

langlangxie

Could you please be more specific about what you are trying to do and what is going wrong? You don't need all those files for the upwelling case. Have you looked through the documentation on the wiki? There is also a manual. Thanks! I am now doing a model of the South China Sea in this case on the ...

langlangxie

Kindhearted persons:
I already have some files,such as .ini,.grd,.clm,.frc and so on.Now I want to calculate my ROMS(It is just upwelling example),but it can't be calculated.I don't know How to modify the file(such upwelling.h and ocean_upwelling.in).
Thanks for everyone help!

Search found 42 matches