Ocean Modeling Discussion

After interpolation, all variables no longer have NaN values in all input files, even for land points. The interpolated values of variable is shown in the figure. (takign temp as the example)

Metrics information for Grid 01: =============================== Minimum X-grid spacing, DXmin = 4.54195740E+00 km Water points = 4.54195740E+00 km Maximum X-grid spacing, DXmax = 5.18572491E+00 km Water points = 5.18572491E+00 km Minimum Y-grid spacing, DYmin = 4.44440448E+00 km Water points = 4.4...

Hello, everyone. Recently, I have prepared for a several-years simulation, after making all needed input files, I tried to run it firstly. But model blows up at beginning, I found there are some NaNs at step 0. TIME-STEP YYYY-MM-DD hh:mm:ss.ss KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME C => (i,j,...

thank you, professor, wish you a happy happy life work!

Dear Prof.Jamie Thank you for your prompt reply and contribution to all ROMS learners, it's really helpful. And Professor, I have reviewed your posts https://www.myroms.org/forum/viewtopic.php?f=29&t=4665 , one of these results is that dual CPU performs better. So, I think a dual EPYC9654 is goo...

Hello,everyone! Recently, I plan to build my own server in the office. And when I was selecting CPUs，I was torn between choosing AMD's Epyc or Intel's Xeon. EPYC 9754 has more cores and higher cost performance than Xeon(Dual EPYC 9754 has 256 cores and 512 threads and I think this already meets my n...

Hello, everyone, wishing you all the best. Recently, I have been learning how to incorporate sediment models into my case. I have tried to find some information about it from wiki roms and forums, but I am still a little confused about adding it. My steps are as follows: 1. modify .in file(set path ...

Kate, thanks for your prompt reply, it's really helpful.

Hello, everyone! Wish you all the best! I added codes to decay the passive tracer in step3d_t.F, as shown in the figure. 1.png And I set a cpp option in .h file. 1698063626366.png When I ran the case, the screen displayed as follow: 1698064990729.png There is no cpp option of TRACER DECAY in the log...

Hello everyone, when I run the TPXO9v5_ROMS_script , it reported some errors. The functions or variables 'ellipses' and 'cj_map' are not recognized. I think there are some libraries I missed, could you please help me point them out? Thank you advanced! subplot(2,1,1) z2a(maskR==0)=nan; contourf(lonR...

Hello, zduvims, sorry for my late reply. It's Panoply, the soft can quick read your .nc file directly.

Thank you, Kate, more information I have got. My simulation is still running which indicates that there are no obvious issues with bathymetry. So I'd like to keep stationary for now.

Dear Kate, Thank you for your reply, besh wishes to you! :D My grid has 15 layers in vertical direction. Actually, I don't want to simulate the situation during a typhoon, it's accidental. So I changed my start time to a calm day which is helpful, and now, with a DT of 20s, the mode has run successf...

Hello, all guys, I checked my forcing files, I found the high values in Uwind and Vwind. So I guess the time I started simulation was during a typhoon, and my simulation domain is quete small which may leads the blowup. Now I would change the start time and try it again, I'm not sure if this help, o...

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000051840) GET_2DFLD_NF90 - surface u-wind component, 2022-07-01 03:00:00.00 (Grid=01, Rec=2, Index=2, File: MERRA_Uwind_3hours_2022_07.nc) (Tmin= 738336.0000 Tmax= 738366.8750) t = 738336.1250 (Min = -1.33804505E+01 Max = ...

The log as shown below: (base) enzo@enzo:~/roms/project/yangjiang_0.3km$ ./romsO < roms_riverplume1_uv.in -------------------------------------------------------------------------------- Model Input Parameters: ROMS/TOMS version 4.1 Tuesday - June 27, 2023 - 9:28:47 PM ------------------------------...

Hello everyone! I am using ROMS to simulate the diffusion of tracers over several days. But I encountered the BLOW UP REASON: MaxSpeed = 2.1E+01 at early steps(71). The maximum water depth of research area is 142m, the resolution of grid is 0.3km. I checked the rst.nc ,I found the max speed in V(-22...

Dear, abiola and Kate, thank you for your kindly and helpfule suggestions! I am so sorry about my late reply. :cry: I rerun the scripts, it works after I added the import words as you said!!!! :D :D This error confused me for several days, and abiola and Kate, you have all helped me more than once. ...

Hello, everyone! Thanks for your attention. Due to my hard drive breaking, I had to rebuild the ROMS environment. I have already successfully installed ROMS and PYROMS, but when I run the make_remap_weights_file.py , it reported an error: (base) enzo@enzo:~/pyroms/inputfiles/yangjiang$ python make_r...

These are remapping .py scripts. remap_weights_file.py: import matplotlib matplotlib.use('Agg') import pyroms import pyroms_toolbox # load the grid srcgrd = pyroms_toolbox.Grid_HYCOM.get_nc_Grid_HYCOM('/home/zeng/soft/pyroms_soft/pyroms/zhuhai_HYCOM/HYCOM_GLBy0.08_ZHUHAI_grid_0.18km.nc') dstgrd = py...

Is that the salt initial condition you want? I doubt it. I think you need to sort out the basic solution before you worry about dyes. And why all the little lakes in the land? Thank you, doctor. I found the reason why some salt concentrations in avg.nc file are abnormally 0. Because there are some ...

You could simplify your test case by choosing Rad instead of RadNud for the dye open boundary condition ... then ROMS should not need use (or even read) the boundary dye values. You have OBCFAC = 0.0, which I think disables the nudging on inflow. Maybe start with OBCFAC = 1.0 so you have boundary n...

Hellow, everyone! I want to simulate the diffusion of some kind of liquid substance(dye, sewage or liquid metal salts etc..) in the ocean. After brief learning on the wiki website and forum, I think adding passive tracer into the run is a good way(probably?) . But I encountered some problems in this...

Hi Zeng, I am not sure if the timeout you added is necessary, though, or probably it is necessary for python3, but I don't think so. You might remove the timeout and try again. Sometimes, it shows some error, especially if you go in bulk. You can try each file once. If it doesn't work still, you mi...

Hi Zeng, You first need to register and activate your account on the data site. Go to https://urs.earthdata.nasa.gov/profile, click on applications, and check your authorized apps. You need to activate or authorize if you are yet to be authorized to use GESDISC and GESDISC Test Data Archive. I hope...

Sorry, it seems like that the URL cannot be displayed completely.. In fact, the URL created by the code is exactly the same as the URL displayed on the data website (url https:// goldsmr4.gesdisc.eosdis.nasa.gov/opendap/MERRA2/M2T1NXRAD.5.12.4/2022/01/MERRA2_400.tavg1_2d_rad_Nx.20220101.nc4). :cry: ...

Hello,everyone! I am a beginer of ROMS, and when I followed the videos to learn how to creat forcing files with Pyroms, I encountered one problem. (base) zeng@zeng-virtual-machine:~/soft/pyroms_soft/pyroms/MERRA-2$ python3 get_MERRA_albedo_from_nasa_opendap_daily.py 2022 url https://goldsmr4.gesdisc...