Ocean Modeling Discussion

Hi all,


I need to configure the roms model for forecast run . The initial, clim and bdry files are from HYCOM data(real-time data) using pyroms.I need to add the forcing files from GFS data.Can you share the script or idea how to create make_frc.py for GFS .
Any help will be highly appreciated

Hi Kate,
Thanks for the help provided during our previous discussion .Happy journey!!

Hi Kate, Actually once I run the make_bdry.py ,this error came. zeta = remap_bdry(file, 'ssh', src_grd, dst_grd, dst_dir=dst_dir) Can you help me to understand what is causing this kind of error src grid=hycom_glby grid dimensions: ocean_time = UNLIMITED ; // (1 currently) lon = 438 ; lat = 751 ; z ...

Hi Kate, Thanks for the reply. My input hycom netcdf file is having negative fill values,it created problem in masking .Now I am facing another error while running another script remap_bdry.py ValueError: operands could not be broadcast together with remapped shapes [original->remapped]: (103,) and ...

Hi Kate,

Thanks for reply,actually i havnt changed the code, ~is present in the new version of pyroms.
"~" returns true,false,1,0
Can you help with the 'mask' variable,based on what condition ,masking happens
mask_t = np.array(~var[:].mask, dtype='int')

Hi Kate, The variables list along with error >>> print(lon.shape) (438,) >>> print(depth.shape) (40,) >>> print(lat.shape) (751,) >>> print(ssh.shape) (751, 438) >>> print(var.shape) (40, 751, 438) >>> mask_t = np.array(~var[:].mask, dtype='int') >>> print(mask_t) 1 >>> z_t = np.tile(depth,(mask_t.s...

Thanks Kate, for the reply I think some issues with my src_gridfile!! The src file is as follows:- netcdf HYCOM_GLBy0.08_2022_065 { dimensions: lon = 438 ; lat = 751 ; ocean_time = UNLIMITED ; // (1 currently) z = 40 ; variables: float lon(lon) ; lon:long_name = "longitude" ; lon:units = &...

Hi Kate,
I am getting following error while running the get_nc_Grid_HYCOM2.py


z_t = np.tile(depth,(mask_t.shape[2],mask_t.shape[1],1)).T
IndexError: tuple index out of range

How to overcome this problem

Hi Kate,
Thanks for the reply,I am running Palau_HYCOM example for another study domain.Using gridbuilder ,I created grid for the study region,once I run the make_remap_weights_file.py ,I am facing above error

Hi everyone, I am facing the error of segmentation fault while executing the make_remap_weights_file.py. Assuming spherical is integer b'T' <class 'numpy.ma.core.MaskedArray'> Load geographical grid from file grid shape 338 640 grid shape 338 639 grid shape 337 640 Segmentation fault Any help with t...

Hi all, I am facing same error and Iam using roms version 3.9. Found Error: 01 Line: 361 Source: ROMS/Utility/close_io.F, close_out ROMS/TOMS - Blows up ................ exit_flag: 1 Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FL...