Ocean Modeling Discussion

hbzong

Recently, I've been tuning an expensive model. To save time, I compiled the model with the SINGLE_PRECISION option. The initialization time of my simulation is 2/1/2021 12:00:00. After turning on SINGLE_PRECISION option, the time of the first step became 2/1/2021 12:01:04, which is 64 seconds ahead ...

hbzong

Did you try to let the model calculate upward longwave radiation (define LONG_OUT)?

In Bohai Sea, WTYPE is also a sensitive parameter for temperature simulation in the summer season.

hbzong

You could read lon_rho, lat_rho, and h from the _roms.nc then use mat2roms_rps.m to generate a *.nc file.

hbzong

Usually, I use the following steps to convert Delft3d grid to ROMS 1. Read the Delft3D grid (see example codes) %Read delft3D grid %Get PSI-points location try fid = fopen(d3dgrd); catch disp('ERROR: Cannot Open file') end % for m = 1:5 tline = fgetl(fid); end % eta_num = fscanf(fid,'%d',1); xi_num ...

hbzong

During running a longterm simulation, I met an error several times. This error occurred in December of model time when restarting the model from the same year. The error says the model couldn't read zeta_time in boundary file. By checking the input, the zeta_time has no problem. If I run the model ...

hbzong

! Vertical mixing coefficients for biological tracers for nonlinear ! model and basic state scale factor in adjoint-based algorithms: ! [1:NBT,Ngrids]. AKT_BAK == 12*1.0d-6 ! m2/s AKT_BAK at here is invalid because model doesn't calculate AKT for passive tracers and uses AKT of salinity to calculate ...

hbzong

I think that this error was caused by NaN values in South/North boundary.

hbzong

It seems that you use Matlab scripts in the directory roms_clm of COAWST to generate boundary file. In their case, west and north boundary are close boundaries. You need to edit create_roms_netcdf_bndry_mw.m and updatbdry_coawst_mw.m to write west boundary informations in netCDF file.

hbzong

When 0 is selected, analyzed data are available. When 3,6,9 and 12 is selected, forecast data are available. If you follow the parameters in d_ecmwf2roms.m, both analyzed and forecast fields are available. And the time array in the netCDF file is look like this: analyzed data time + forecast data ...

hbzong

Thanks, Kate.

hbzong

I'm running a case with ROMS(svn 666). At beginning, I used Flather boundary conditions for both Ubar and Vbar. The model worked. Then I changed the both boundary conditions to Clamped. The model crashed. And I found that the model could run if set Clamped condition for Ubar and Flather boundary for ...

hbzong

The model could not find a input file.

hbzong

There is a possible bug in routines fennel_mod.h. ! Set number of diagnostics terms. ! NDbio3d=2 NDbio2d=0 # ifdef DENITRIFICATION NDbio2d=NDbio2d+1 # endif # ifdef CARBON NDbio2d=NDbio2d+2 # endif # ifdef OXYGEN NDbio2d=NDbio2d+1 # endif IF DIAGNOSTICS_BIO is defined and DENITRIFICATION, CARBON and ...

hbzong

I'm trying to setup a tidal model with wetting-drying scheme. The minimum bathymetry in my grid is 0.1m, and the minimum depth for wetting and drying was set to 0.5m. Blowing-up occurs at the step 16979. Only tide forcing was added to the model. DT == 60.0 NDTFAST == 10 Minimum barotropic Courant ...

hbzong

I updated the bathymetry of a coupling case. Before I updated the bathymetry, the model could run successfully. But now, it can’t run. Firstly, I updated the bathymetry and the mask of ROMS grid file. Then, generated SWAN grid file, bathymetry file and wind files with new ROMS grid file. When only 1 ...

hbzong

Did you consider the time lag between Beijing Time and GMT when you compared your model result with observation?

hbzong

John, I'm using ROMS-SWAN coupling model to study wind wave induced sediment resuspension in the Yellow River Estuary of China.

hbzong

feroda wrote: By Matlab,
>> nc{'river'}(:)

ans =

1
2

[/quote]
Set the river ID starting from 0.

nc{'river'}(:) = [0,1];

hbzong

There is "READ_WATER" option in the ROMS Cpp file. Does it mean that ROMS can read the gird without the longitude and latitude information of dry points? I tried it.The longitude and latitude of dry points were defined as zero. The model can't run because TotValume = NaN, However, the MinVolume and ...

hbzong

Hi,Lyon.Could you leave your email?

hbzong

Same error during compiling.
I created a new directory and used svn checkout to get the totally new codes. And same problem ...

hbzong

I added #define MCT_LIB to my application header file:

#define SWAN_COUPLING
#ifdef SWAN_COUPLING
# define MCT_LIB
# define SSW_BBL
#endif

The error still occurred during compiling.

hbzong

I see. Thank you very much.

hbzong

New update was applied this morning. However, I couldn't compile it either. fortcom: Error: inp_par.f90, line 494: The ARRAY argument of the ALLOCATED inquiry intrinsic function shall be an allocatable array. [ALLOCATED] IF (.not.allocated(Export(Iocean)%code)) THEN ...

hbzong

I updated ROMS this morning and then I couldn't compile it successfully with swan coupling. A array named 'OutField' in roms_export.f90 (line 91)was not be declared.

hbzong

t seems that the error occurred when extrapolating the values to the empty grid points. I also did a four points case, but i didn't find the same error. You could comment the following lines when you do four points case: % [x,y] = size(GRD_LON); % nn = 1; % for m = 1:x % for n = 1:y % if ~isnan(GRD ...

hbzong

Thanks everyone.
Rob, how did you extrapolate the values to the empty grid points. I read your pyroms code. Because I'm not familiar with python, I did found how you do it. I used the MATLAB griddata function, but the result doesn't seem good.

hbzong

In order to use gridgen conveniently, I wrote a MATLAB-based GUI tool for creating input file for gridgen, reading output file of gridgen and displaying grid. It is easy to create a polygon for gridgen by using it. See http://www.coastal.udel.edu/~hbzong https://www.myroms.org/links/t_762/step1.jpg

hbzong

The output of gridgen contains 'nan' values when using gridgen to generate grid for the domain shown below. How to deal with these 'nan' value?

hbzong

Hi, everyone I'm trying to use Roms_tools to create initial fields file, forcing fields file and climatology fields file as input files of ROMS3.0. As you known, some variables' name in Roms_tool are not the same in ROMS3.0. I have checked the codes in Preprocessing_tools of Roms_tools. I just found ...

hbzong

Maybe you should install OpenMP into your system and turn on USE_OpenMP while you compile ROMS.

hbzong

There are two output files (*.mat-001,*mat-002,PRINT-001,PRINT-002) when NnodesWAV = 2.
I removed all *mat-001,*mat-002, PRINT*, swaninit and etc files, then rerun it. There were same errors.

hbzong

I'm running it on a Linux cluster. The ROMS is the latest version.
As the error messages shown, it seems that a MPI node of ocean model(p0, ocean model use p0,p1,p2,p3) died first, and then caused the communication problem. But i don't know what make p0 node death.

hbzong

X=6 , NnodesOCN = 4 NnodesWAV = 2 The error message are shown below: ------------------------------------------------------------------------------------------------------------------------------- NL ROMS/TOMS: started time-stepping:( TimeSteps: 00000001 - 00001440) == SWAN sent wave fields and ...

hbzong

I have done what you said. My X = NnodesOCN + NnodesWAV I can set X=3 , NnodesOCN = 2, NnodesWAV = 1 or X=5 , NnodesOCN = 4, NnodesWAV = 1 or any others while NnodesWAV = 1 But I can't set X=4 , NnodesOCN = 2, NnodesWAV = 2 or X=6 , NnodesOCN = 4, NnodesWAV = 2 or any others while NnodesWAV > 1 ...

hbzong

When i run the TESTHEAD example on a cluster with the mpi, I can't set more than one CPU for SWAN. ROMS can use more than one CPU, but SWAN can't. If I set more than one CPU for SWAN, there was a error during running. Could you give me some suggestions?

Ocean Modeling Discussion

Search found 36 matches

Model time issue with SINGLE_PRECISION option

Re: High Temperature in Summer

Re: Delft3D grids to NetCDF

Re: Delft3D grids to NetCDF

One possible bug with restart

Invalid parameter in bio_Fennel.in

Re: Model blowing up using elevation at the boundary

Re: unable to find requested variable: u_west

Re: Selection with regard to ECMWF ERA datasets

Re: Possible problem of clamped condition for Vbar

Possible problem of clamped condition for Vbar

Re: I/O error in running

Possible bug in writing out biological diagnostics file

Blowing-up problem

Problem on coupling

Re: urgent problem on time ref and tide_start !!!!!!!

Re: ROMS-SWAN Coupled System Users

Re: problems after updating to latest revision(321)

Question about "READ_WATER" option

Re: Delft3D grids to NetCDF

Compiling problem with swan coupling

A MATLAB-based GUI for gridgen

A question about the gridgen

Use Roms_tools to creat input files for ROMS

SWAN can't use more than one CPU