Search found 1171 matches

by jcwarner
Mon Mar 18, 2024 1:14 pm
Forum: ROMS Discussion
Topic: some problems about coawst_ini.nc,spinup period and wave input
Replies: 3
Views: 175

Re: some problems about coawst_ini.nc,spinup period and wave input

you can modify those m files as needed. and you can certainly search around for other tools. for the wave forcing, as i said already, we used to run swan, save the data, create a roms forc file, and then run roms. in theory it used to work, and it might now. but i dont run the modeling system that w...
by jcwarner
Mon Mar 18, 2024 1:06 pm
Forum: ROMS Discussion
Topic: Some problems about adding wave data in ROMS case
Replies: 1
Views: 116

Re: Some problems about adding wave data in ROMS case

these are settings for surface flux of tke due to wave breaking, see
https://www.sciencedirect.com/science/a ... 0312000042
by jcwarner
Mon Mar 18, 2024 12:28 pm
Forum: ROMS Problems
Topic: Error in inp_decode.f90 while running the roms_upwelling.in
Replies: 1
Views: 22

Re: Error in inp_decode.f90 while running the roms_upwelling.in

look at
line 307 of file inp_decode.f90
and see what the code is trying to read.
by jcwarner
Fri Mar 15, 2024 12:15 pm
Forum: ROMS Discussion
Topic: some problems about coawst_ini.nc,spinup period and wave input
Replies: 3
Views: 175

Re: some problems about coawst_ini.nc,spinup period and wave input

1- i am not sure what is causing the blips at the coast. not sure i have seen this before. is it also in u and v? i would think so. are you nudging to ubar vbar? if not, then it is not an issue. 3 - we used to run roms forced with the output from swan. to do that, we would run swan, make a roms forc...
by jcwarner
Tue Mar 05, 2024 1:01 pm
Forum: ROMS Installation
Topic: COAWST JOE_TC case compile error
Replies: 9
Views: 14276

Re: COAWST JOE_TC case compile error

it would be best to post this on the coawst issue site:
https://github.com/DOI-USGS/COAWST/issues
and include the full output, not just those last few lines.
by jcwarner
Tue Feb 06, 2024 4:06 pm
Forum: ROMS Problems
Topic: Problem with Blowing-up, REASON: MaxSpeed = 2.5E+01
Replies: 5
Views: 830

Re: Problem with Blowing-up, REASON: MaxSpeed = 2.5E+01

you need to do something there.
not sure what the bath looks like. maybe a little smoothing there. maybe increase visc2 a little.
maybe decrease time step. try a few things.
by jcwarner
Tue Feb 06, 2024 1:07 pm
Forum: ROMS Problems
Topic: Problem with Blowing-up, REASON: MaxSpeed = 2.5E+01
Replies: 5
Views: 830

Re: Problem with Blowing-up, REASON: MaxSpeed = 2.5E+01

"WRT_RST_NF90 - wrote re-start fields (Index=2,2) in record = 1"

have you viewed the rst file?
use ncview or similar tool to quickly look in the file, see where the problem exists.
by jcwarner
Sat Jan 27, 2024 12:16 am
Forum: ROMS Tools and Techniques
Topic: different pm and pn using COAWST and Rutgers Toolkits
Replies: 2
Views: 985

Re: different pm and pn using COAWST and Rutgers Toolkits

yes sorry.
i was changing the way to compute the metrics to be more in line with gridbuilder and i messed up the x y values used to compute dx dy and angle.
sorry.
i reverted my method back to what i had before. if you do a git pull it should be better now.

-j
by jcwarner
Sun Jan 21, 2024 5:41 pm
Forum: ROMS Discussion
Topic: how to add sediment model into my simulation
Replies: 1
Views: 12568

Re: how to add sediment model into my simulation

sorry for the late reply on this one.
are you still trying to incorporate the sediment.?
You can use the coawst model for now, we have more updated methods there. we are working to include all the sediment into the rutgers version now.
-j
by jcwarner
Mon Jan 08, 2024 2:29 pm
Forum: ROMS Discussion
Topic: Tidal dissipation and topographic wave drag in ROMS
Replies: 7
Views: 15715

Re: Tidal dissipation and topographic wave drag in ROMS

ok. here is some 'no responsibility' guidance. Look in ROMS/Nonliear/set_vbc.F and there is a section for the bottom stress near line 637 or so: # elif defined UV_QDRAG ! ! Set quadratic bottom stress. ! DO j=Jstr,Jend DO i=IstrU,Iend cff1=0.25_r8*(v(i ,j ,1,nrhs)+ & & v(i ,j+1,1,nrhs)+ &...
by jcwarner
Mon Jan 08, 2024 1:31 pm
Forum: ROMS Discussion
Topic: Tidal dissipation and topographic wave drag in ROMS
Replies: 7
Views: 15715

Re: Tidal dissipation and topographic wave drag in ROMS

ok. So does this term need to be added throughout the water column, or just as an extra term that could be added to the bottom stress?
by jcwarner
Thu Jan 04, 2024 8:09 pm
Forum: ROMS Discussion
Topic: Tidal dissipation and topographic wave drag in ROMS
Replies: 7
Views: 15715

Re: Tidal dissipation and topographic wave drag in ROMS

This could be added as a body force, or it could be the same way we added vegetation drag throughout the water column.
What is the T?
by jcwarner
Thu Dec 21, 2023 2:53 pm
Forum: ROMS Discussion
Topic: Abnormal termination: NetCDF INPUT
Replies: 13
Views: 22334

Re: Abnormal termination: NetCDF INPUT

ok. now you know what is happening. next is to figure out why it is happening. it blew up in 12 time steps, so save every his step, and watch what happens in that corner. what are the sal and temp in your bc files? I personally do not use A4 advection, i use U3. maybe try different advection or diff...
by jcwarner
Tue Dec 19, 2023 1:44 pm
Forum: ROMS Discussion
Topic: Abnormal termination: NetCDF INPUT
Replies: 13
Views: 22334

Re: Abnormal termination: NetCDF INPUT

12 2022-07-01 12:18:00.00 2.693400E-02 5.341908E+03 5.341935E+03 4.910017E+13 (588,005,03) 4.670190E-01 2.073180E-01 2.452327E+00 2.479567E+01 DIAG speed trouble 24.795668439173749 DIAG speed ijk 20.680908753216158 589 5 1 1 Blowing-up: Saving latest model state into RESTART file REASON: MaxSpeed = ...
by jcwarner
Sun Dec 17, 2023 3:18 pm
Forum: ROMS Discussion
Topic: Abnormal termination: NetCDF INPUT
Replies: 13
Views: 22334

Re: Abnormal termination: NetCDF INPUT

i think you can start out by first adjusting the time step. and see if you can get the configuration you have to run.
by jcwarner
Fri Dec 15, 2023 1:07 pm
Forum: ROMS Discussion
Topic: Abnormal termination: NetCDF INPUT
Replies: 13
Views: 22334

Re: Abnormal termination: NetCDF INPUT

192 ntimes Number of timesteps for 3-D equations. 900.000 dt Timestep size (s) for 3-D equations. 30 ndtfast Number of timesteps for 2-D equations between each 3D timestep. roms time step should be a lot smaller. i see you have ~1km spacing. try 90 sec or something like that to start.
by jcwarner
Thu Dec 14, 2023 8:40 pm
Forum: ROMS Discussion
Topic: Abnormal termination: NetCDF INPUT
Replies: 13
Views: 22334

Re: Abnormal termination: NetCDF INPUT

here is some help.
viewtopic.php?t=2930
by jcwarner
Thu Dec 14, 2023 1:27 pm
Forum: ROMS Discussion
Topic: Abnormal termination: NetCDF INPUT
Replies: 13
Views: 22334

Re: Abnormal termination: NetCDF INPUT

keep looking in your file, it says

"INQUIRY - unable to find requested variable: EminusP
in file:
Projects/UPWELLING/ROMS_forcing.nc"
by jcwarner
Wed Dec 13, 2023 2:08 pm
Forum: ROMS Discussion
Topic: Blow up before starting
Replies: 5
Views: 14307

Re: Blow up before starting

GET_NGFLD - free-surface western boundary condition, 1947-12-12 00:00:00.00 (Grid= 01, Rec=12, Index=2, File: croco_bry.nc) (Tmin= 15.0000 Tmax= 345.0000) t = 345.0000 (Min = NaN Max = NaN) GET_NGFLD - free-surface northern boundary condition, 1947-12-12 00:00:00.00 (Grid= 01, Rec=12, Index=2, File:...
by jcwarner
Wed Dec 06, 2023 1:18 pm
Forum: ROMS Discussion
Topic: WRF quilting support in the WRF-ROMS setup in COAWST
Replies: 3
Views: 10653

Re: WRF quilting support in the WRF-ROMS setup in COAWST

this is best to post a message on the coawst forum:
https://github.com/DOI-USGS/COAWST/issues
and the answer is yes, it should work.
by jcwarner
Fri Nov 24, 2023 2:15 pm
Forum: ROMS Installation
Topic: COAWST JOE_TC case compile error
Replies: 9
Views: 14276

Re: COAWST JOE_TC case compile error

for coawst issues, specific to compiling or running multiple models together, it is best to post an issue here: https://github.com/DOI-USGS/COAWST/issues your issue is that you selected a compiler for openmp and mpi. (dm+sm) For cowast, you can only use mpi. if you have further issues please post a ...
by jcwarner
Thu Nov 02, 2023 11:30 pm
Forum: ROMS Usage
Topic: "thin dams" in ROMS
Replies: 4
Views: 18220

Re: "thin dams" in ROMS

can you just set the u-mask to be zero at a point and that would prevent the flow across that interface (similar for v). ? .
by jcwarner
Fri Sep 29, 2023 10:35 pm
Forum: ROMS Discussion
Topic: Boundary Zeta Instability
Replies: 2
Views: 2623

Re: Boundary Zeta Instability

look at the temp and salt, they might be contributing.
what did you set for the bc's in the ocean.in file.
by jcwarner
Fri Aug 25, 2023 1:22 pm
Forum: ROMS Discussion
Topic: How to download HYCOM GLBa0.08 expt60.5 data?
Replies: 14
Views: 14773

Re: How to download HYCOM GLBa0.08 expt60.5 data?

a couple of things here. 1 - there is a depth in your grid exactly = 0 (h(52,4) ==0) and that caused some troubles in the routine to comute ubar and vbar cc=roms_zint_mw(u,gn); ubar=rho2u_2d_mw(u2rho_2d_mw(cc)./gn.h); cc=roms_zint_mw(v,gn); vbar=rho2v_2d_mw(v2rho_2d_mw(cc)./gn.h); %== Rotate the vel...
by jcwarner
Mon Aug 21, 2023 3:54 pm
Forum: Job Opportunities
Topic: Coastal Ocean Modeler USGS Woods Hole
Replies: 0
Views: 9119

Coastal Ocean Modeler USGS Woods Hole

The USGS - Woods Hole Coastal and Marine Science Center (WHCMSC) is seeking a coastal ocean modeler. The USGS conducts a wide range of interdisciplinary research activities in marine, estuarine, lacustrine, and terrestrial regions of the United States, including the Atlantic Ocean, Pacific Ocean, Gu...
by jcwarner
Fri Aug 18, 2023 1:05 pm
Forum: ROMS Discussion
Topic: Linear Bottom Drag & Horizontal Viscosity
Replies: 2
Views: 2913

Re: Linear Bottom Drag & Horizontal Viscosity

1. Look in ROMS/Nonlinear/set_vbc.F The bottom stress for a 2D application is # if defined UV_LDRAG ! ! Set linear bottom stress. ! DO j=Jstr,Jend DO i=IstrU,Iend bustr(i,j)=0.5_r8*(rdrag(i-1,j)+rdrag(i,j))* & & ubar(i,j,krhs) ... etc and for 3D it is # elif defined UV_LDRAG ! ! Set linear b...
by jcwarner
Wed Aug 09, 2023 12:52 pm
Forum: ROMS Problems
Topic: Problem with dimensions.
Replies: 4
Views: 12855

Re: Problem with dimensions.

not sure, but it looks like you put ! in front of a lot of lines. Even if you are not using that info, the inp script looks for that data. for example ! NestLayers = 1 < this line needs to be read ! Number of grids in each nesting layer [1:NestLayers]. ! GridsInLayer = 1 < this line needs to be read...
by jcwarner
Sun Aug 06, 2023 1:06 pm
Forum: ROMS Bugs
Topic: cannot work parallel nodes
Replies: 1
Views: 2534

Re: cannot work parallel nodes

ROMS/TOMS: Wrong choice of grid 01 partition or number of parallel nodes. NtileI * NtileJ = 2 must be equal to the number of parallel processes = 1 Change -np value to mpirun or change domain partition in input script. if you set to run with -np X, where X is an integer, then you need to change ocea...
by jcwarner
Tue Jul 18, 2023 12:59 pm
Forum: ROMS Bugs
Topic: Nsrc init
Replies: 2
Views: 5109

Re: Nsrc init

I had a mix of options that were not correct and it caused an issue. There is a catch in mod_arrays to prevent my problem: IF (LuvSrc(ng).or.LwSrc(ng).or.ANY(LtracerSrc(:,ng))) THEN CALL allocate_sources (ng) END IF so this is all good. sorry to raise the flag. happy modeling. -j
by jcwarner
Fri Jul 14, 2023 2:12 pm
Forum: ROMS Bugs
Topic: Nsrc init
Replies: 2
Views: 5109

Nsrc init

I am having an issue with a memory allocation, and it came back to the init for the var Nsrc. In ROMS/Modules/mod_sources.F we have #ifndef ANA_PSOURCE .... Nsrc(mg)=var_Dsize(1) ! first dimension .... #else ... Msrc(ng)=200 Nsrc(ng)=Msrc(ng) #endif ... ! Allocate point Sources/Sinks variables. ! al...
by jcwarner
Wed Jul 12, 2023 12:47 pm
Forum: ROMS Discussion
Topic: Random inconsistencies on the temperature and salinity variables
Replies: 5
Views: 4624

Re: Random inconsistencies on the temperature and salinity variables

ok do me a favor.
cd ROMS/Nonlinear
mv gls_corstep.F gls_corstep.F_orig
go here:
https://myroms.org/forum/viewtopic.php?p=24570#p24570
and download their alt_gls_corstep.F,
mv alt_gls_corstep.F gls_corstep.F
then recompile and see if that helps.
-j
by jcwarner
Mon Jul 10, 2023 11:43 pm
Forum: ROMS Discussion
Topic: Random inconsistencies on the temperature and salinity variables
Replies: 5
Views: 4624

Re: Random inconsistencies on the temperature and salinity variables

if you set the coefficients to use k-e or k-gen, does the spots still appear?
by jcwarner
Mon Jul 10, 2023 12:18 am
Forum: ROMS Discussion
Topic: Random inconsistencies on the temperature and salinity variables
Replies: 5
Views: 4624

Re: Random inconsistencies on the temperature and salinity variables

what kind of vertical mixing do you have on?
do you have any wave-ocn coupling?
by jcwarner
Wed Jun 28, 2023 9:32 pm
Forum: ROMS Discussion
Topic: The magnitude of vertical mixing and advection terms
Replies: 4
Views: 4081

Re: The magnitude of vertical mixing and advection terms

those are about the same dynamics that we simulated for H Ivan.
Zambon, J.B., He, R., and Warner, J.C. (2014). Investigation of Hurricane Ivan using the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) Model, Ocean Dynamics, 64, 1535-1554. doi:10.1007/s10236-014-0777-7.
by jcwarner
Tue Jun 13, 2023 4:13 pm
Forum: ROMS Discussion
Topic: DIAGNOSTICS_TS in step3d_t
Replies: 3
Views: 3157

Re: DIAGNOSTICS_TS in step3d_t

i think it has to deal with the time stepping and the Hz grid cell thickness. in prestep3d we have cff1=Hz(i,j,k)*t(i,j,k,nstp,itrc) ... # ifdef DIAGNOSTICS_TS DiaTwrk(i,j,k,itrc,iTrate)=cff1 so the iTrate is t(nstp)*Hz ~= told*Hzold in step3dt we do 2 things: DO idiag=1,NDT DiaTwrk(i,j,k,itrc,idiag...
by jcwarner
Tue Jun 13, 2023 12:34 pm
Forum: ROMS Discussion
Topic: Model blows up before it starts
Replies: 1
Views: 2580

Re: Model blows up before it starts

the last line of your file says "Could not read valid record axe from file" i had to search for this one. it is in SWAN, call agnc_set_error( 'Could not read valid record axe from file') looks like you are trying to read a netcdf file. I suggest you use SWAN by itself, try to run that appl...
by jcwarner
Thu Jun 08, 2023 12:41 pm
Forum: ROMS Installation
Topic: COAWST Compiling cowastM failed
Replies: 2
Views: 2916

Re: COAWST Compiling cowastM failed

For compile errors with coawst, please submit them to https://github.com/jcwarner-usgs/COAWST/issues i took a quick look at your compile log. there are lots of issues. first, it shows "f951: Warning: Nonexistent include directory ‘/main’ [-Wmissing-include-dirs] f951: Warning: Nonexistent inclu...
by jcwarner
Fri Jun 02, 2023 12:34 pm
Forum: ROMS Source
Topic: Using rain as LwSrc
Replies: 0
Views: 3229

Using rain as LwSrc

We are using ROMS to simulate rainfall and urban flooding during hurricanes. Currently the 'rain' input from the surface forcing is used for several physics such as a BC for salinity. This will allow the water to get fresher. But the volume of that rain is not accounted for. The LwSrc point source a...
by jcwarner
Mon May 15, 2023 12:58 pm
Forum: ROMS Discussion
Topic: about Bottom Friction
Replies: 1
Views: 2313

Re: about Bottom Friction

if you have a uniform current on a non-uniform grid, the stress would be uniform because the current is uniform.
the areas are taken into account when the time stepping of the full equations.
you can dig into the details here:
https://www.myroms.org/wiki/Documentation_Portal
by jcwarner
Thu May 11, 2023 11:55 am
Forum: ROMS Discussion
Topic: Bottom Friction and How to Slow Down Currents in Shallow Regions
Replies: 11
Views: 7310

Re: Bottom Friction and How to Slow Down Currents in Shallow Regions

Do you think it is due to the vegetation, or just that the water is shallow? you could use the VEGETATION option that we have in coawst. it extracts momentum from the water column up to the veg height, and this has been shown to reduce flow in those areas. Beudin, A., Kalra, T. S., Ganju, N. K., and...
by jcwarner
Thu May 04, 2023 1:15 am
Forum: ROMS Discussion
Topic: Interpolation of bed stress for step_floats
Replies: 8
Views: 11084

Re: Interpolation of bed stress for step_floats

i think 0 is at the lower left psi point. maybe Hernan is better to answer. if so,then int(track(ixgrd)) would give the left side of the box at a u-point, and int(track(ixgrd))+1 would be the right side u-point. but you have i1 = MIN(INT (track(ixgrd,nfp1,1)+0.5_r8) ,Lm(ng)+1) i2 = MAX(i1 + 1,0) so ...
by jcwarner
Thu Apr 27, 2023 8:16 pm
Forum: ROMS Bugs
Topic: An error in ROMS-SWAN coupling and INLET_TEST problem
Replies: 4
Views: 4267

Re: An error in ROMS-SWAN coupling and INLET_TEST problem

yes, use coawst for roms+swan coupling, and that system also has the option to also use wrf. Rutgers has WRF coupling and that is with ESMF. i have not used that yet, but it will be happening soon. Also, Rutgers is looking to couple WAVEWATCH III to ROMS via ESMF, but that is just starting. So for n...
by jcwarner
Thu Apr 27, 2023 9:14 am
Forum: ROMS Bugs
Topic: An error in ROMS-SWAN coupling and INLET_TEST problem
Replies: 4
Views: 4267

Re: An error in ROMS-SWAN coupling and INLET_TEST problem

if you are using coawst, please post a message here:
https://github.com/jcwarner-usgs/COAWST/issues

you would need to use coawst.bash
we include a user manual.doc that has lots of info.
thanks
j
by jcwarner
Sun Apr 16, 2023 12:47 am
Forum: ROMS Problems
Topic: Tidal flood-ebb asymmetry
Replies: 3
Views: 3810

Re: Tidal flood-ebb asymmetry

looks like model is taking longer to ebb, that could be due to friction too high on the landward side of this site.
you need to look at an animation and see how the flow looks.
by jcwarner
Fri Apr 07, 2023 5:17 pm
Forum: ROMS Discussion
Topic: Error while using contact.m : Unrecognized field name "contact".
Replies: 2
Views: 3547

Re: Error while using contact.m : Unrecognized field name "contact".

when you made the child grid, did you use coarse2fine:

Istr=22; Iend=60; Jstr=26; Jend=54;
ref_ratio=3;
roms_child_grid='Sandy_roms_grid_ref3.nc';
F=coarse2fine('Sandy_roms_grid.nc','Sandy_roms_grid_ref3.nc', ...
ref_ratio,Istr,Iend,Jstr,Jend);
by jcwarner
Tue Apr 04, 2023 12:47 pm
Forum: ROMS Discussion
Topic: Issue with HISTORY file logical switch for Hbbl
Replies: 6
Views: 5472

Re: Issue with HISTORY file logical switch for Hbbl

in ROMS/Utility/def_his.F I see # ifdef LMD_SKPP ! ! Define depth of surface boundary layer. ! IF (Hout(idHsbl,ng)) THEN ......... status=def_var(ng, iNLM, HIS(ng)%ncid, HIS(ng)%Vid(idHsbl), & & NF_FOUT, nvd3, t2dgrd, Aval, Vinfo, ncname) IF (FoundError(exit_flag, NoError, __LINE__, MyFile))...
by jcwarner
Mon Apr 03, 2023 12:46 pm
Forum: ROMS Discussion
Topic: wrf-roms coupling for tropical cyclone simulation duration
Replies: 1
Views: 2394

Re: wrf-roms coupling for tropical cyclone simulation duration

yes i saw your post on the other site and just responded the same thing there. The Sandy test case is a test case. I do that for purposes to demonstrate a procedure for setting up and running the coupled system. The grids are very coarse, and the simulation is short, just to show how to do things. E...
by jcwarner
Fri Mar 24, 2023 3:50 pm
Forum: ROMS Discussion
Topic: How to download HYCOM GLBa0.08 expt60.5 data?
Replies: 14
Views: 14773

Re: How to download HYCOM GLBa0.08 expt60.5 data?

i updated the script (get_ijrg.m) a few months ago (?) to use: % % Read HYCOM lon lat depth % display(['getting HYCOM grid data from ', url]) lonlat_fullgrid=1; try numX=ncread(url,'X'); numY=ncread(url,'Y'); hycom_lon=ncread(url,'Longitude'); hycom_lat=ncread(url,'Latitude'); hycom_depth=ncread(url...
by jcwarner
Wed Mar 15, 2023 1:37 pm
Forum: ROMS Bugs
Topic: Possible Coding Error in gls_corstep.F: application of gls_Kmin and gls_Pmin to control advection
Replies: 4
Views: 12017

Re: Possible Coding Error in gls_corstep.F: application of gls_Kmin and gls_Pmin to control advection

ok thanks for digging in that deep. it helps a lot.
let me look into it. I do remember those min values were only suggestions, from the literature. But if there is an issue with the min being applied to tke*Hz and not to tke (and gls), then that needs to be resolved.
-j
by jcwarner
Fri Mar 03, 2023 12:12 pm
Forum: ROMS Installation
Topic: undefined reference to `crosskss_ While Compiling COAWST
Replies: 2
Views: 3789

Re: undefined reference to `crosskss_ While Compiling COAWST

this is a coawst issue that has been resolved.
please do a git pull.
in the future please post these issues on the coawst site here:
https://github.com/jcwarner-usgs/COAWST/issues

thanks
j
by jcwarner
Sun Feb 19, 2023 3:03 am
Forum: ROMS Discussion
Topic: Model could not recognize variables
Replies: 8
Views: 8068

Re: Model could not recognize variables

i see the mud_time is the same as salt_time. but i dont think that is the problem. you should be able to list any name, as long as the dimension information is correct. you also need to look in the varinfo.dat. the mud vars should be listed as using mud_time. maybe try to change varinfo from mud_tim...
by jcwarner
Sat Feb 18, 2023 2:56 pm
Forum: ROMS Discussion
Topic: Model could not recognize variables
Replies: 8
Views: 8068

Re: Model could not recognize variables

what does this give:
ncdisp('D:\Data\NC\coase_1010\out_4m\coase_21_bry_nonorth_sed05.nc')
by jcwarner
Fri Feb 17, 2023 1:46 pm
Forum: ROMS Discussion
Topic: Model could not recognize variables
Replies: 8
Views: 8068

Re: Model could not recognize variables

are these 2 files the same:

D:\Data\NC\coase_1010\out_4m\coase_21_bry_nonorth_sed05.nc
and
/public1/home/sc90163/project_hzhy/coase_21/4m_in/coase_21_bry_nonorth_sed05.nc
by jcwarner
Fri Feb 03, 2023 1:32 pm
Forum: ROMS Discussion
Topic: ROMS crashing when using LwSrc
Replies: 18
Views: 31087

Re: ROMS crashing when using LwSrc

"The run seems to stop halfway when ninfo is too small, but I'll try to make it as small as possible." When you set Ninfo == 1, roms will make a call to Nonlinear/diag.F every time step IF (MOD(iic(ng)-1,ninfo(ng)).eq.0) THEN in here , roms computes max values for vel and dens, compares th...
by jcwarner
Wed Feb 01, 2023 2:58 pm
Forum: ROMS Discussion
Topic: problem in using Hvom and Huon
Replies: 8
Views: 10398

Re: problem in using Hvom and Huon

yes i figured that. run the job for a short duration and look at the diagnositcs.nc file. you will need to spend some time, start with the depth avg terms. they are easier to figure out. try to add up the terms at one of the time steps: ubar_accel = ubar_sustr+ubar_bustr+ubar_prsgrd+ubar_cor+ubar_ha...
by jcwarner
Wed Feb 01, 2023 2:02 pm
Forum: ROMS Discussion
Topic: problem in using Hvom and Huon
Replies: 8
Views: 10398

Re: problem in using Hvom and Huon

the wiki does not have a full explanation, but there are other sources to get the meaning of the diagnostics. https://www.myroms.org/wiki/Model_Diagnostics here is a list of the equations https://www.myroms.org/wiki/Numerical_Solution_Technique Kates manual has a nice more full expansion of the term...
by jcwarner
Wed Feb 01, 2023 1:17 pm
Forum: ROMS Discussion
Topic: problem in using Hvom and Huon
Replies: 8
Views: 10398

Re: problem in using Hvom and Huon

this can get complicated in the code. Suggest you use the diagnostics, they are set up to compute all this for you.
by jcwarner
Tue Jan 31, 2023 1:46 pm
Forum: ROMS Discussion
Topic: ROMS crashing when using LwSrc
Replies: 18
Views: 31087

Re: ROMS crashing when using LwSrc

ok but can you still try my other suggestions: 1- set Ninof to be 1, so we can see if there is more information. 3 what it the temp coming in for this point source? 4- what is the time series of Qin (river flow) at the problematic point? can you plot that? 5- use KANTHA_CLAYSON, not CANUTO. my exper...
by jcwarner
Mon Jan 30, 2023 1:41 pm
Forum: ROMS Discussion
Topic: ROMS crashing when using LwSrc
Replies: 18
Views: 31087

Re: ROMS crashing when using LwSrc

so i can offer a few suggestions again 1- set Ninof to be 1, so we can see if there is more information. 2- use HSIMT for all the tracers including salt + temp 3 you say "Another note: forcing for vertical sources varies on a monthly basis. Thus, the input transport, temperature, biogeochem, et...
by jcwarner
Fri Jan 20, 2023 5:51 pm
Forum: ROMS Installation
Topic: Problem compiling version 3.7 of COAWST
Replies: 5
Views: 4874

Re: Problem compiling version 3.7 of COAWST

at the top of your coawst.bash you have
parallel=0
clean=1
dprint=0
cleanwrf=0
cleanwrfhydro=0
cleanww3=0

switch that to
parallel=0
clean=1
dprint=0
cleanwrf=1
cleanwrfhydro=1
cleanww3=1

this should fix it.
by jcwarner
Fri Jan 20, 2023 1:55 pm
Forum: ROMS Installation
Topic: Problem compiling version 3.7 of COAWST
Replies: 5
Views: 4874

Re: Problem compiling version 3.7 of COAWST

ok now i see your log file.
but that log file is for "william_onlyroms.h"
And that built correctly.
So that is not the correct log file for me to look at.
by jcwarner
Fri Jan 20, 2023 1:50 pm
Forum: ROMS Installation
Topic: Problem compiling version 3.7 of COAWST
Replies: 5
Views: 4874

Re: Problem compiling version 3.7 of COAWST

the problem probably occurred way near the beginning of the build. and then it got to the end and all the modules were not there. When WRF builds, if it gets an error, it keeps going. So you need to scroll way up to the top and start from there to look for the 'error'. You can output the build using...
by jcwarner
Mon Jan 16, 2023 2:16 pm
Forum: ROMS Discussion
Topic: Sandy COAWST - Bathymetry
Replies: 1
Views: 2750

Re: Sandy COAWST - Bathymetry

the bathy in that test case is very crude. it is only intended as a guide.
Strongly suggest you use a different source for the bathy.
here is one source: https://download.gebco.net/
here is another: https://www.ncei.noaa.gov/products/etop ... lief-model

-j
by jcwarner
Fri Jan 13, 2023 2:30 pm
Forum: ROMS Discussion
Topic: ROMS crashing when using LwSrc
Replies: 18
Views: 31087

Re: ROMS crashing when using LwSrc

here are a few thoughts / suggestions. 1 what are the magnitudes of mass transport for the LwSRC? i see: GET_NGFLD_NF90 - river runoff mass transport, 2021-02-20 12:00:00.00 (Grid= 01, Rec=5, Index=2, File: rivers.nc) (Tmin= 18674.5000 Tmax= 18681.5000) t = 18678.5000 (Min = -9.00360189E+03 Max = 1....
by jcwarner
Fri Dec 09, 2022 3:24 pm
Forum: ROMS Discussion
Topic: About the CPP option for the WRF-ROMS-SWAN coupling model
Replies: 1
Views: 2388

Re: About the CPP option for the WRF-ROMS-SWAN coupling model

you really should post this on the coawst site here: https://github.com/jcwarner-usgs/COAWST/issues because Rutgers ROMS does not have all the wave physics yet. We are porting them over but it will take some time. Suggest you look at the COAWST user manual, test cases of Inlet_test, Rip_current, Duc...
by jcwarner
Wed Dec 07, 2022 1:34 pm
Forum: ROMS Problems
Topic: [SIGSEGV: Segmentation fault] when running SWAN only
Replies: 1
Views: 3047

Re: [SIGSEGV: Segmentation fault] when running SWAN only

this post is about SWAN, and it would best placed on the coawst ticket site
https://github.com/jcwarner-usgs/COAWST/issues
as a short reply, you need to look at any Err* or PRINT* files for more info.
if you are going to follow up, post a new instance on the coawst issues.
-j
by jcwarner
Mon Dec 05, 2022 2:21 pm
Forum: ROMS Problems
Topic: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic
Replies: 7
Views: 6175

Re: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic

Can you look at the different components of the heat budget: lwrad, swrad, latent, sens, and compare the coupled to uncoupled of those 4 terms?
by jcwarner
Fri Dec 02, 2022 3:15 pm
Forum: ROMS Problems
Topic: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic
Replies: 7
Views: 6175

Re: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic

ok, that figure looks better and i can see what you are talking about. Suggest you look at the different components of the heat budget: lwrad, swrad, latent, sens, and compare the coupled to uncoupled of those 4 terms. Try undef ANA_SRFLUX, although i dont think that is the problem. i am not familia...
by jcwarner
Fri Dec 02, 2022 1:18 pm
Forum: ROMS Problems
Topic: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic
Replies: 7
Views: 6175

Re: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic

yes i understand what T2 is.
I see an image, with lots of colors, i dont know geographically where it is.
I dont know where the 'anamoly' is, etc. Do the temp issue appear right away, is it on land, on the ocean,
I do not have enough info to understand the issue.
by jcwarner
Thu Dec 01, 2022 7:07 pm
Forum: ROMS Problems
Topic: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic
Replies: 7
Views: 6175

Re: COAWST - ROMS+WRF Coupled Strange BIAS Atlantic

not sure what i am looking at here.
what does the SST in roms look like?
by jcwarner
Wed Nov 23, 2022 6:18 pm
Forum: ROMS Discussion
Topic: weakly and partially stratified estuaries using ROMS
Replies: 2
Views: 2304

Re: weakly and partially stratified estuaries using ROMS

i remember setting this up. there was some trial-and-error to keep the salt where i wanted it. there is a balance between bathy, Qriver, and mixing (basically the tidal forcing). So as you adjust that balance, the salt field will evolve. For that particular very idealized case (constant dh/dx) the b...
by jcwarner
Sat Nov 19, 2022 2:47 pm
Forum: Meetings/Workshops
Topic: EGU2023 - OS1.2/AS2: Air-sea interactions through modeling and observation on multiple temporal and spatial scales.
Replies: 0
Views: 3451

EGU2023 - OS1.2/AS2: Air-sea interactions through modeling and observation on multiple temporal and spatial scales.

EGU2023 - OS1.2/AS2: Air-sea interactions through modeling and observation on multiple temporal and spatial scales. We would like to invite submissions to session OS1.2/AS2: Air-sea interactions through modeling and observation on multiple temporal and spatial scales. This session aims to gather res...
by jcwarner
Mon Oct 31, 2022 1:29 pm
Forum: ROMS Discussion
Topic: How to download HYCOM GLBa0.08 expt60.5 data?
Replies: 14
Views: 14773

Re: How to download HYCOM GLBa0.08 expt60.5 data?

so i just tried the second one you have, and it worked for me: disp('getting roms grid, hycom grid, and overlapping indices') [gn, clm]=get_ijrg(url, modelgrid, theta_s, theta_b, Tcline, N, Vtransform, Vstretching); getting roms grid, hycom grid, and overlapping indices getting roms grid dimensions ...
by jcwarner
Thu Sep 01, 2022 2:34 pm
Forum: ROMS Installation
Topic: Error compiling ROMS in COAWST
Replies: 5
Views: 3988

Re: Error compiling ROMS in COAWST

there is a lot of info in the COAWST_User manual.doc. if you only select the roms model, then coawstM is essentially romsM. if you select additional other models, only 1 executable is created, coawstM. That would run roms + wrf, or whatever. this is a different approach than what Rutgers does, so if...
by jcwarner
Thu Sep 01, 2022 12:01 pm
Forum: ROMS Installation
Topic: Error compiling ROMS in COAWST
Replies: 5
Views: 3988

Re: Error compiling ROMS in COAWST

use the coawst.bash file in the main dir. I have modified the build to be able to compile all the different models.
if you have more issues please post here:
https://github.com/jcwarner-usgs/COAWST/issues
thanks, j
by jcwarner
Tue Aug 23, 2022 12:40 pm
Forum: ROMS Discussion
Topic: Issue with the interpolation between time steps in forcing file
Replies: 3
Views: 2038

Re: Issue with the interpolation between time steps in forcing file

seems odd.
when you did the test with the extra time step in the middle, what did roms report out when it read the file?
it always reports out the min, max, times, etc. check that.
by jcwarner
Sat Jul 30, 2022 11:34 pm
Forum: ROMS Usage
Topic: An error occurred when running ROMS in COAWST
Replies: 10
Views: 12739

Re: An error occurred when running ROMS in COAWST

DIAG speed trouble 93.106400122773167 DIAG speed ijk 131.67231528779021 96 64 1 1 Found Error: 01 Line: 341 Source: ROMS/Nonlinear/main3d.F Found Error: 01 Line: 331 Source: ROMS/Drivers/nl_ocean.h, ROMS_run Blowing-up: Saving latest model state into RESTART file REASON: MaxSpeed = 9.3E+01 WRT_RST -...
by jcwarner
Sat Jul 30, 2022 1:52 pm
Forum: ROMS Usage
Topic: An error occurred when running ROMS in COAWST
Replies: 10
Views: 12739

Re: An error occurred when running ROMS in COAWST

I see still: GET_GRID - bathymetry at RHO-points: h (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc) (Min = 5.95202990E-03 Max = 3.35037820E+03) that is a min of 0.005 m . I dont think you need a min of 5 mm. try again. in the grid file to set h to have a minimum depth. maybe something like ...
by jcwarner
Fri Jul 29, 2022 8:51 pm
Forum: ROMS Discussion
Topic: Issue with HISTORY file logical switch for Hbbl
Replies: 6
Views: 5472

Re: Issue with HISTORY file logical switch for Hbbl

for gls i typically use

# define GLS_MIXING
# ifdef GLS_MIXING
# define RI_SPLINES
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# endif

it can get more options, for surface mixing etc, not sure what type of applications you are using.

for LMD i will let someone else answer.
by jcwarner
Fri Jul 29, 2022 7:49 pm
Forum: ROMS Discussion
Topic: Issue with HISTORY file logical switch for Hbbl
Replies: 6
Views: 5472

Re: Issue with HISTORY file logical switch for Hbbl

gls and tke will only be output if you activate GLS_MIXING or maybe also MY25_MIXING.

for the hsbl, those go with the LMD options.

-j
by jcwarner
Fri Jul 29, 2022 2:58 pm
Forum: ROMS Usage
Topic: An error occurred when running ROMS in COAWST
Replies: 10
Views: 12739

Re: An error occurred when running ROMS in COAWST

GET_GRID - bathymetry at RHO-points: h (Grid = 01, File: Projects/Katrina/roms_grid_Katrina.nc) (Min = -2.15744419E+02 Max = 3.33282292E+03) you need to activate wet dry if you have negative h values. Suggest you take a look at your grid and decide if you are really resolving wet/dry processes. if n...
by jcwarner
Wed Jul 27, 2022 12:36 am
Forum: ROMS Discussion
Topic: Didn't create average files
Replies: 3
Views: 2076

Re: Didn't create average files

how many time steps did you run it for?
by jcwarner
Mon Jul 25, 2022 11:45 am
Forum: ROMS Problems
Topic: MaxDensity error where the h changes rapidly
Replies: 4
Views: 5473

Re: MaxDensity error where the h changes rapidly

where are the locations for the : LwSrc == T
by jcwarner
Mon Jul 25, 2022 11:37 am
Forum: ROMS Problems
Topic: ROMS<kshedstrom version> blow at the fist timestep
Replies: 9
Views: 7102

Re: ROMS<kshedstrom version> blow at the fist timestep

we have MAIN: Abnormal termination: BLOWUP. REASON: MaxDensity = 3.4E+02 can you look in the ocean_rst.nc and see where that density is large? is it at the open boundary? at a river location? maybe try to turn off the rivers, and turn off clm nudging. see if it runs without them to help diagnose the...
by jcwarner
Sun Jul 24, 2022 1:10 pm
Forum: ROMS Problems
Topic: ROMS<kshedstrom version> blow at the fist timestep
Replies: 9
Views: 7102

Re: ROMS<kshedstrom version> blow at the fist timestep

ok so this is a different user on the same post. to zhaoyun- all you are telling us is that it blows up on first time step. You need to provide more information. Look at the std out and see if there is any information like h<0 and no wetdry, or nans in input forc fields, etc. the warnings are ok, th...
by jcwarner
Sat Jul 23, 2022 1:42 pm
Forum: ROMS Problems
Topic: ROMS<kshedstrom version> blow at the fist timestep
Replies: 9
Views: 7102

Re: ROMS<kshedstrom version> blow at the fist timestep

did you also update your ocean.in file to include the BC's and advection choices?
-j
by jcwarner
Sat Jul 23, 2022 1:41 pm
Forum: ROMS Bugs
Topic: Model running error: REASON: KEchar = NaN, PEchar = NaN
Replies: 3
Views: 4162

Re: Model running error: REASON: KEchar = NaN, PEchar = NaN

there can be many reasons why the model is not starting. i remember seeing this GET_GRID - bathymetry at RHO-points: h (Grid = 02, File: Projects/Mangkhut/Mangkhut_roms_grid_ref3.nc) (Min = 0.00000000E+00 Max = 3.49541846E+03) if your min depth is 0, you need to define wet_dry. -j
by jcwarner
Fri Jul 22, 2022 1:33 pm
Forum: ROMS Discussion
Topic: The roms simulates a 6 hour difference in tide levels.
Replies: 2
Views: 2216

Re: The roms simulates a 6 hour difference in tide levels.

looks like the simulation starts with zeta = -0.5 but the obs are at +0.5. so it will take a while for that barotropic mode to slosh out of the model. after 2 days the water levels look like they are starting to get into phase. if that storm at time 40 is your main interest, i suggest you start the ...
by jcwarner
Sun Jul 17, 2022 2:56 pm
Forum: ROMS Discussion
Topic: How to download HYCOM GLBa0.08 expt60.5 data?
Replies: 14
Views: 14773

Re: How to download HYCOM GLBa0.08 expt60.5 data?

oops, not enough coffee yet. those ecmwf is for bulk fluxes. there are other posts on getting hycom data, such as: https://www.myroms.org/forum/viewtopic.php?p=23944&hilit=hycom#p23944 https://www.myroms.org/forum/viewtopic.php?p=23938&hilit=hycom#p23938 https://www.myroms.org/forum/viewtopi...
by jcwarner
Sun Jul 17, 2022 2:45 pm
Forum: ROMS Discussion
Topic: How to download HYCOM GLBa0.08 expt60.5 data?
Replies: 14
Views: 14773

Re: How to download HYCOM GLBa0.08 expt60.5 data?

doing griddata u for HYCOM level 1 ... doing griddata u for HYCOM level 2 catch u Unable to download HYCOM u data at17-Jul-2022 23:08:30 if it worked for level 1, but not for level 2, that is an internet or Thredds issue. this is common. you can try again, or try some other way. Also, if you want to...
by jcwarner
Sun Jul 17, 2022 1:36 pm
Forum: ROMS Discussion
Topic: How to download HYCOM GLBa0.08 expt60.5 data?
Replies: 14
Views: 14773

Re: How to download HYCOM GLBa0.08 expt60.5 data?

the tools we have there are for aggregated data

https://tds.hycom.org/thredds/catalog.html
i dont see the 60.5 on that list.

-j
by jcwarner
Sat Jul 09, 2022 1:30 am
Forum: ROMS Information
Topic: LIMIT_VVISC
Replies: 5
Views: 6955

Re: LIMIT_VVISC

i dont think GLS has changed much in several years, so it should be the same.
-j
by jcwarner
Fri Jul 08, 2022 10:01 pm
Forum: ROMS Information
Topic: LIMIT_VVISC
Replies: 5
Views: 6955

Re: LIMIT_VVISC

ok. First of all let's tackle the AKV_limit issue. This option is intended to prevent the values of Akt and Akv from getting very large. I do a lot of estuary work with lots of mixing, so I recommend setting these limiter values be at something like AKV_limit=1.0d0 and Akt_limit=1.0d0. IF you are un...
by jcwarner
Wed Jul 06, 2022 12:36 am
Forum: ROMS Tools and Techniques
Topic: obc_roms2roms.m
Replies: 1
Views: 4258

obc_roms2roms.m

Is anyone else using obc_roms2roms.m? I am getting a size mismatch at lines 720-723 B.(ufield) = u .* cos(angle.(edge)) - ... v .* sin(angle.(edge)); the ufield = 'u_west' size(u) = 2 327 40 size(cos(angle.(edge))) = 2 327 40 size(v) = 2 326 40 size(sin(angle.(edge))) = 2 327 40 obviously u and v ar...
by jcwarner
Sat Jul 02, 2022 3:17 pm
Forum: ROMS Discussion
Topic: The problem of creating initial boundary conditions
Replies: 1
Views: 2058

Re: The problem of creating initial boundary conditions

the variable names change over the years. I have no control in that. I just pushed out an update to the m files. please do an update and try again. thanks. if you have more troubles, please post a message on the coawst git site, as this issue for those tools. https://github.com/jcwarner-usgs/COAWST/...
by jcwarner
Tue Jun 14, 2022 6:55 pm
Forum: ROMS Bugs
Topic: LwSrc i,j bug
Replies: 1
Views: 4339

LwSrc i,j bug

For the LwSrc, in step3d_t.F near lines 1146 we need to change cff=dt(ng)*pm(i,j)*pn(i,j) to cff=dt(ng)*pm(Isrc,Jsrc)*pn(Isrc,Jsrc) and the same near lines 1343 from cff=dt(ng)*pm(i,j)*pn(i,j) to cff=dt(ng)*pm(Isrc,Jsrc)*pn(Isrc,Jsrc) -j
by jcwarner
Mon Jun 13, 2022 11:51 pm
Forum: ROMS Discussion
Topic: How to improve modelled salinity and salinity initial condition
Replies: 14
Views: 15406

Re: How to improve modelled salinity and salinity initial condition

what bottom stress are you using? UV_LOGDRAG?
what is your Zob?
by jcwarner
Mon Jun 13, 2022 2:51 pm
Forum: ROMS Discussion
Topic: How to improve modelled salinity and salinity initial condition
Replies: 14
Views: 15406

Re: How to improve modelled salinity and salinity initial condition

i think you need more Volume to be filled on the flood, and then this will allow increased flows on the ebb. I was hoping to see more overtides in that water level time series (ie quicker rise to peak water level and longer tail to lower water level). do you see that further into the estuary? i see ...
by jcwarner
Mon Jun 13, 2022 1:04 pm
Forum: ROMS Discussion
Topic: How to improve modelled salinity and salinity initial condition
Replies: 14
Views: 15406

Re: How to improve modelled salinity and salinity initial condition

some quick looks: -how long did you run this for? looks like the salt is getting better as time goes on. -in the attached ppt, the vel in Charcoal bay does not get as strong (peaky) as the observations. So the salt will not get advected as far at those peaks, and that is exactly what you see in the ...
by jcwarner
Sun Jun 12, 2022 3:19 pm
Forum: ROMS Discussion
Topic: too high bottom temperature result in estuary
Replies: 6
Views: 9700

Re: too high bottom temperature result in estuary

you need to look at how the temperature is evolving over time. does it get warmer from the surface net fluxes and then mix downwards? does it get warmer from the offshore boundary and then propagate in? i see: GET_3DFLD - potential temperature climatology, t = 45 15:45:04 (Rec=0000002, Index=2, File...
by jcwarner
Mon Jun 06, 2022 11:48 am
Forum: ROMS Discussion
Topic: Questions about wrf-roms coupling
Replies: 20
Views: 27968

Re: Questions about wrf-roms coupling

glad it is working for you!
-j
by jcwarner
Fri Jun 03, 2022 12:30 pm
Forum: ROMS Sediment
Topic: unable to find variable: sand in input NetCDF
Replies: 8
Views: 6606

Re: unable to find variable: sand in input NetCDF

not exactly sure what is going on from that information. but it looks like you are doing a restart, and the model is looking for fields in that restart file. did you do a run with sediment and now restarting that? or did the first run not have sediment? if the first run did not have sediment, then y...
by jcwarner
Thu Jun 02, 2022 8:09 pm
Forum: ROMS Discussion
Topic: Questions about wrf-roms coupling
Replies: 20
Views: 27968

Re: Questions about wrf-roms coupling

ok i spent some time on this. I think the issue is a bathy mismatch between the child and the parent, especially near the boundary. At the location of that vortex, the child grid has a bathymetric rise, but the parent grid did not see it. So i recommend that you use average child bathy to replace th...
by jcwarner
Wed Jun 01, 2022 12:48 pm
Forum: ROMS Discussion
Topic: Questions about wrf-roms coupling
Replies: 20
Views: 27968

Re: Questions about wrf-roms coupling

I looked at the grids and they seem ok. the bathy matches well. the dmde looks large in the parent grid. are you running with curvgrid on? also the parent has spherical ='T' and the child has spherical = 1. this is probably not an issue. how did you set the init conditions? are the T/S fields consis...
by jcwarner
Tue May 31, 2022 12:01 pm
Forum: ROMS Discussion
Topic: How to calculate turbulent shear rate / energy dissipation
Replies: 3
Views: 4185

Re: How to calculate turbulent shear rate / energy dissipation

the turbulence and the sediment should be separate. Both MY25 and GLS will output tke and gls. IF you use MY25, then tke is 0.5q^2 (i think, it has been a while) and gls will be the q^2l. IF you use GLS , then tke is tke and gls will be the set depending on the p m n parameters (gls can be dissipati...
by jcwarner
Mon May 30, 2022 10:54 pm
Forum: ROMS Discussion
Topic: Questions about wrf-roms coupling
Replies: 20
Views: 27968

Re: Questions about wrf-roms coupling

looks like the only differences for roms is along the southern boundary.
with nesting, does the vortex start and grow, is it formed during the experiment at the boundary and advect in?
does the child bathy match the parent?
i assume this probably gets deep along that child southern boundary.
-j
by jcwarner
Fri May 27, 2022 12:20 am
Forum: Job Opportunities
Topic: USGS Woods Hole Technician V Coastal Modeler Support
Replies: 0
Views: 2853

USGS Woods Hole Technician V Coastal Modeler Support

Akima is hiring a Technician V position to support the USGS Woods Hole Coastal and Marine Science Center (WHCMSC), Woods Hole, Massachusetts. The Technician V will assist the WHCMSC in developing, implementing, testing, and applying a coupled ocean-atmosphere-wave-sediment transport modeling system ...
by jcwarner
Mon May 23, 2022 5:36 pm
Forum: ROMS Discussion
Topic: How to run COAWST model with real-time data
Replies: 6
Views: 8211

Re: How to run COAWST model with real-time data

if you are having trouble getting hycom data, suggest that you use ECMWF data. John Wilkin has a nice set of tools to do this:
roms_get_era5_NCARds633_bulkflux.m
and
roms_write_era5_NCARds633_frcfile.m
from here:
https://github.com/johnwilkin/roms_wilkin
by jcwarner
Mon May 23, 2022 12:55 pm
Forum: ROMS Discussion
Topic: Questions about wrf-roms coupling
Replies: 20
Views: 27968

Re: Questions about wrf-roms coupling

i am not worried about the free surface, if that is all dry anyway. do you have wet dry on? did you set hmin = 5 in both the parent and child? is this vortex created at the boundary and advected in? can you do me a favor, i made a change to my version of nesting.F, Are you using COAWST? if so please...
by jcwarner
Fri May 20, 2022 3:47 pm
Forum: ROMS Installation
Topic: Memory leaks - Upwelling
Replies: 8
Views: 9618

Re: Memory leaks - Upwelling

all those 'mem leaks' are for calls to allocate state variables such as allocate ( OCEAN(ng) % rubar(LBi:UBi,LBj:UBj,2) ) allocate ( OCEAN(ng) % rvbar(LBi:UBi,LBj:UBj,2) ) allocate ( OCEAN(ng) % rzeta(LBi:UBi,LBj:UBj,2) ) allocate ( OCEAN(ng) % ubar(LBi:UBi,LBj:UBj,3) ) allocate ( OCEAN(ng) % vbar(L...
by jcwarner
Fri May 20, 2022 1:22 am
Forum: ROMS Installation
Topic: Memory leaks - Upwelling
Replies: 8
Views: 9618

Re: Memory leaks - Upwelling

ok i need another one. Build/mod_ocean.f90
-j
by jcwarner
Thu May 19, 2022 7:15 pm
Forum: ROMS Installation
Topic: Memory leaks - Upwelling
Replies: 8
Views: 9618

Re: Memory leaks - Upwelling

not sure what these are, but if you got a netcdf output file, then the model ran fine.
can you post
Build/mod_arrays.f90
that seems to come up a lot.
thanks
j
by jcwarner
Thu May 19, 2022 3:04 pm
Forum: ROMS Discussion
Topic: Questions about wrf-roms coupling
Replies: 20
Views: 27968

Re: Questions about wrf-roms coupling

I think this person is using COAWST. And yes they should add SST to cover non-Roms areas, and they probably are. We have had this option since the beginning of the coupling. What we do is use the SST from the WRFlow file, and replace that with ROMS data wherever the ocean is. We do that for each rom...
by jcwarner
Tue May 17, 2022 11:44 am
Forum: ROMS Discussion
Topic: Questions about wrf-roms coupling
Replies: 20
Views: 27968

Re: Questions about wrf-roms coupling

is this the nest or the parent grid?
if this is the nest, are the vortices advected into the nest from the parent?
what do the temp and salt fields look like? you are using radnudge on the south, so what are you nudging to?
by jcwarner
Tue May 10, 2022 12:08 pm
Forum: ROMS Discussion
Topic: How to create the hycom_time with the real-time data
Replies: 4
Views: 3847

Re: How to create the hycom_time with the real-time data

it is not possible for me to tell what is going on. The issue is with the array sizes. You will need to put in break points, see what the arrays are, the sizes, the values, and follow the arrays as they are passed to each routine. it is not that hard, just put in break points. when matlab stops at t...
by jcwarner
Fri May 06, 2022 12:03 pm
Forum: ROMS Discussion
Topic: How to create the hycom_time with the real-time data
Replies: 4
Views: 3847

Re: How to create the hycom_time with the real-time data

for example here https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/FMRC/runs/GLBy0.08_930_FMRC_RUN_2022-05-05T12:00:00Z.html the variable list will show for time: units: hours since 2022-04-30 12:00:00.000 UTC so the 120 is hours since that time. You need to look for these attributes on the dat...
by jcwarner
Wed May 04, 2022 11:50 am
Forum: ROMS Discussion
Topic: An ocean_coupler problem when compiled
Replies: 2
Views: 2100

Re: An ocean_coupler problem when compiled

sorry. i think you posted a message on the coawst issues and i did not get back to you.
i think you have a misspelled cpp def

define ATM2OCE_FLUXES
should be
define ATM2OCN_FLUXES

then do a ./coawst.bash and see if that compiles better.
-john
by jcwarner
Wed Apr 27, 2022 7:42 pm
Forum: ROMS Discussion
Topic: River runoff problem
Replies: 9
Views: 13514

Re: River runoff problem

https://www.myroms.org/wiki/River_Runof ... _Eposition
is a good description of the river locations.
it needs to be coming out of the land masked cell.
by jcwarner
Wed Apr 27, 2022 4:47 pm
Forum: ROMS Discussion
Topic: River runoff problem
Replies: 9
Views: 13514

Re: River runoff problem

it could be i, j, or transport.
make sure the Q is positive for flow to the right.
look at the u values and ensure the rivers are coming from the correct locations.
those locations should be masked points. the river flow will overwrite the zero values at the masked locations.

-j
by jcwarner
Wed Apr 27, 2022 12:56 pm
Forum: ROMS Discussion
Topic: converting roms grid to swan grid
Replies: 7
Views: 4907

Re: converting roms grid to swan grid

if you are using coawst, you can post a message here:
https://github.com/jcwarner-usgs/COAWST/issues

also i suggest you look at the Inlet_test case and the User manual for descriptions on how to set these up.
-j
by jcwarner
Wed Apr 27, 2022 12:54 pm
Forum: ROMS Discussion
Topic: How do I solve the simulation speed anomaly?
Replies: 3
Views: 3101

Re: How do I solve the simulation speed anomaly?

ok. if you min depth = 5 meters, and then wet/dry helped, i suggest you look at the zeta values to make sure the tides are ok. what is the tide range in this area? for the next problem, that is probably a density issue. closely look at what is being put along that boundary for temp salt u v etc. you...
by jcwarner
Tue Apr 26, 2022 1:26 pm
Forum: ROMS Problems
Topic: Runtime error line 322 of inp_decode.f90
Replies: 1
Views: 2531

Re: Runtime error line 322 of inp_decode.f90

you could look at
At line 322 of file inp_decode.f90
and see what is trying to be read. looks like it wants a real but getting an integer or something like that.
by jcwarner
Fri Apr 22, 2022 11:54 am
Forum: ROMS Discussion
Topic: How do I solve the simulation speed anomaly?
Replies: 3
Views: 3101

Re: How do I solve the simulation speed anomaly?

what is the bathy in that corner? i see you have 2 Vtransform S-coordinate transformation equation. 4 Vstretching S-coordinate stretching function. 5.0000E+00 theta_s S-coordinate surface control parameter. 4.0000E-01 theta_b S-coordinate bottom control parameter. 50.000 Tcline S-coordinate surface/...
by jcwarner
Thu Apr 21, 2022 2:31 pm
Forum: ROMS Problems
Topic: At line 1895 of file mod_ncparam.f90 Fortran runtime error
Replies: 6
Views: 5298

Re: At line 1895 of file mod_ncparam.f90 Fortran runtime error

not sure if this is the correct answer, but you need to be consistent with salinity and NT value settings.

If you use
#define SALINITY
then make sure ocean.in has NAT = 2

If you use
#undef SALINITY or not even mention salinity,
then make sure ocean.in has NAT = 1

-j
by jcwarner
Tue Apr 19, 2022 4:48 pm
Forum: ROMS Discussion
Topic: tracers(temp and salt )distributed noisy
Replies: 2
Views: 2238

Re: passive temp and salt distributed unphysically

not sure what you mean by the visc values. Are you using horiz viscosity: #define UV_VIS2 #define MIX_S_UV then you can set VISC2 == 2.0d0 ! m2/s ??? or something like that Perhaps the tracer horiz diffusivity is more what you need: # define TS_DIF2 # define MIX_GEO_TS TNU2 == 0.5d0 ! m2/s ??? or so...
by jcwarner
Wed Apr 13, 2022 12:15 am
Forum: ROMS Discussion
Topic: "gls" in the output file, and calculation of TKE dissipation
Replies: 7
Views: 5915

Re: "gls" in the output file, and calculation of TKE dissipation

oh ok. let me look at that glossary, it might not be worded clearly.
can you paste that part of the glossary here?
if you choose k-e coeffs, then gls will be e (dissipation).
by jcwarner
Tue Apr 12, 2022 8:18 pm
Forum: ROMS Discussion
Topic: "gls" in the output file, and calculation of TKE dissipation
Replies: 7
Views: 5915

Re: "gls" in the output file, and calculation of TKE dissipation

here is a good paper with lots of info: Warner, J.C., Sherwood, C., Arango, H., and Signell, R. (2005). “Performance of Four Turbulence Closure Models Implemented Using a Generic Length Scale Method.” Ocean Modelling, v. 8/1-2, p. 81-113. the GLS scheme is a method that transports tke and 'gls' = ge...
by jcwarner
Fri Apr 08, 2022 6:53 pm
Forum: ROMS Installation
Topic: Finished COAWST installation ubuntu 18.04
Replies: 2
Views: 4558

Re: Finished COAWST installation ubuntu 18.04

that is only a warning.
For WRF you need to
export NETCDF=/location/of/netcdf
did you set that?
by jcwarner
Wed Apr 06, 2022 12:27 pm
Forum: ROMS Discussion
Topic: Couple WRF and ROMS:Could not find input parameter: NestLayers
Replies: 1
Views: 1957

Re: Couple WRF and ROMS:Could not find input parameter: NestLayers

this is more of a coawst configuration issue, and should be posted here:
https://github.com/jcwarner-usgs/COAWST/issues
can you post this and the coupling.in file over there.
-j
by jcwarner
Tue Apr 05, 2022 12:09 pm
Forum: ROMS Discussion
Topic: ROMS not bit-reproducible: northern boundary?
Replies: 36
Views: 28578

Re: ROMS not bit-reproducible: northern boundary?

i had a little trouble following your logic, but it is not really clear what is causing the differences in the tke. Can you look at fields of temp, u, v, AKv, etc to see if there is something noticeably different in the runs? like something happening along the open boundary. Are there rivers (with f...
by jcwarner
Mon Apr 04, 2022 2:22 pm
Forum: ROMS Installation
Topic: COAWST installation through Cygwin on Windows 10
Replies: 4
Views: 4464

Re: COAWST installation through Cygwin on Windows 10

this is a common issue with gfortran
you need to edit the Makefile.conf (i think) and add a fflag
-fallow-argument-mismatch

-j
by jcwarner
Sat Apr 02, 2022 2:44 pm
Forum: ROMS Installation
Topic: Error when installing SCRIP_COAWST
Replies: 1
Views: 2503

Re: Error when installing SCRIP_COAWST

your executable is created here: -o ./scrip_coawst so it does not have a .exe on the end. do a ls -ltr and do you see ..... -rw-rw-r-- 1 jwarner domain users 103960 Jan 10 2020 remap_write.o -rw-rw-r-- 1 jwarner domain users 19413 Jan 10 2020 scrip.fpp -rw-rw-r-- 1 jwarner domain users 20302 Jan 10 ...
by jcwarner
Fri Apr 01, 2022 12:53 pm
Forum: ROMS Installation
Topic: COAWST installation through Cygwin on Windows 10
Replies: 4
Views: 4464

Re: COAWST installation through Cygwin on Windows 10

make sure you have
c:/cygwin64/usr/include/mpif.h
if not then you missed something during the setup, like libopenmpifh40 or something like that
cygwinlib.png
by jcwarner
Mon Mar 28, 2022 2:37 pm
Forum: ROMS Discussion
Topic: EW_periodic,NS_periodic
Replies: 2
Views: 2180

Re: EW_periodic,NS_periodic

All the boundary condition options have been moved to the ocean.in file. You dont need the boundary cpp options like EW_PERIODIC or FSCLAMPED etc. it is all in the ocean.in For Example: ! Set lateral boundary conditions keyword. Notice that a value is expected ! for each boundary segment per nested ...
by jcwarner
Sun Mar 27, 2022 6:19 pm
Forum: ROMS Installation
Topic: Fortran runtime error: Assignment of scalar to unallocated array
Replies: 7
Views: 10115

Re: Fortran runtime error: Assignment of scalar to unallocated array

Sunny up here but looks like rain coming your way Hernan. Enjoy the day! -j
by jcwarner
Sun Mar 27, 2022 3:16 pm
Forum: ROMS Installation
Topic: Fortran runtime error: Assignment of scalar to unallocated array
Replies: 7
Views: 10115

Re: Fortran runtime error: Assignment of scalar to unallocated array

maybe you could help a bit and track that message, what is on

/mnt/c/Users/lyras/OneDrive/install_roms/folder_tes/Build_romsG/mod_scalars.f90:2595
by jcwarner
Thu Mar 24, 2022 4:56 pm
Forum: ROMS Problems
Topic: Signed Integer Overflow
Replies: 10
Views: 14226

Re: Signed Integer Overflow

i am not sure the difference between those 2 output files. the server.log has - bathymetry at RHO-points: h (Grid = 01, File: bowl_grid.nc) (Min = 2.00018817E-02 Max = 1.00000000E+03) Minimum Z-grid spacing, DZmin = 6.53658218E-04 m that is a small dz somewhere. you have thetas =10. maybe adjust the...
by jcwarner
Thu Mar 24, 2022 4:49 pm
Forum: ROMS Installation
Topic: compile ROMS error
Replies: 12
Views: 16241

Re: compile ROMS error

both the Netcdf libs and roms need to be built with the same compiler and version of that compiler.

type this:
nf-config --all

and this
mpif90 -v

find out what compiler was used to build the Netcdf libs, or have someone rebuild the netcdf with the same mpif90 that you are using.
by jcwarner
Wed Mar 23, 2022 12:46 pm
Forum: ROMS Problems
Topic: Signed Integer Overflow
Replies: 10
Views: 14226

Re: Signed Integer Overflow

how about you post the whole output, not just the last part.
there might be something wrong with the init conditions.
by jcwarner
Mon Mar 21, 2022 1:08 pm
Forum: ROMS Discussion
Topic: Different boundary scenarios
Replies: 15
Views: 16663

Re: Different boundary scenarios

I am nots really sure what all to suggest here. A few images back , you had zeta on the west side at -8m. Can you figure out why that happened? Even a super strong wind would have trouble to do that. I thought the winds might not be rotated to the grid correctly. Suggest you set roms to write out th...
by jcwarner
Wed Mar 16, 2022 6:43 pm
Forum: ROMS Discussion
Topic: Different boundary scenarios
Replies: 15
Views: 16663

Re: Different boundary scenarios

do you have
#define CURVGRID
that might help with some of the advection.
take a look at the winds in that ncviewer. how strong are they on the west side?

-j
by jcwarner
Mon Mar 14, 2022 2:00 am
Forum: ROMS Ice
Topic: blow-up in first time step, using glorys and ERA5
Replies: 8
Views: 11135

Re: blow-up in first time step, using glorys and ERA5

does it run without the clm nudging?
look at the entire output file, there might have been more info before the nan.
by jcwarner
Mon Mar 14, 2022 1:56 am
Forum: ROMS Bugs
Topic: An error was reported after only a few hours of operation
Replies: 2
Views: 2623

Re: An error was reported after only a few hours of operation

1581120 2008-07-02 00:00:00.00 0.000000E+00 1.453619E+04 1.453619E+04 1.282427E+15 (000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 DEF_HIS - creating history file, Grid 01: Ev08_his_coupled_00184.nc WRT_HIS - wrote history fields (Index=1,1) in record = 0000001 DEF_AVG - creating av...
by jcwarner
Mon Mar 14, 2022 12:17 am
Forum: ROMS Discussion
Topic: Different boundary scenarios
Replies: 15
Views: 16663

Re: Different boundary scenarios

think you are changing too many things here. ERA winds should not affect the tides. and dont rely on ncview to give you the max/min. click at a point and get a time series. try this: You have this: Variable Grid West Edge South Edge East Edge North Edge --------- ---- ---------- ---------- ---------...
by jcwarner
Fri Mar 11, 2022 1:26 pm
Forum: ROMS Discussion
Topic: Different boundary scenarios
Replies: 15
Views: 16663

Re: Different boundary scenarios

1.JPG
on this page, top figure, what is 'wrong' with the zeta?
by jcwarner
Wed Mar 09, 2022 1:44 pm
Forum: ROMS Results
Topic: ROMS output. Best technique to interpolate in the vertical
Replies: 8
Views: 7859

Re: ROMS output. Best technique to interpolate in the vertical

looking better.
there are 'extrapolation' options, which need to be used with caution.
but give that a try.
maybe someone else has a better idea.
by jcwarner
Wed Mar 09, 2022 12:14 am
Forum: ROMS Results
Topic: ROMS output. Best technique to interpolate in the vertical
Replies: 8
Views: 7859

Re: ROMS output. Best technique to interpolate in the vertical

it just wants you to enter the data as a column. lets say this is a 2D field in x and z. lets also say that you have 10 vertical layers and we are looking at the y index of 35 (just picking any number here). N=10; yidx=35; X=repmat(lon_rho(:,yidx),N); Z=z_rho(:,yidx,:); S=squeeze(salt(:,yidx,:,tidx)...
by jcwarner
Tue Mar 08, 2022 8:03 pm
Forum: ROMS Results
Topic: ROMS output. Best technique to interpolate in the vertical
Replies: 8
Views: 7859

Re: ROMS output. Best technique to interpolate in the vertical

" linear interpolation over one column at the time" - that is probably not a good idea.
try a 2D interpolation scheme. Are you using matlab?
help scatteredInterpolant
F = scatteredInterpolant(x_rho,z_rho,salt);
saltnew= F(xnew,znew);
by jcwarner
Mon Mar 07, 2022 1:43 pm
Forum: ROMS Discussion
Topic: SIGSEGV, segmentation fault occurred (with Empty log file)
Replies: 6
Views: 8214

Re: SIGSEGV, segmentation fault occurred (with Empty log file)

take a look at the file Build/mod_ncparam.f90 at line 1039 and see what it is looking for.

romsS 00000000005696B9 mod_ncparam_mp_in 1039 mod_ncparam.f90
by jcwarner
Sat Mar 05, 2022 6:48 pm
Forum: ROMS Discussion
Topic: forrtl: severe (174): SIGSEGV, segmentation fault occurred
Replies: 9
Views: 11265

Re: forrtl: severe (174): SIGSEGV, segmentation fault occurred

"forrtl: severe (408): fort: (2): Subscript #2 of the array FRC has value 1 which is greater than the upper bound of -1" i looked in your ocean.in and you have: !NFFILES = 1 ! number of unique forcing files !FRCNAME == /disk3/chenzhixing/ROMS/cases/laizhou/laizhou-frc.nc ! forcing file 1, ...
by jcwarner
Thu Mar 03, 2022 1:14 pm
Forum: ROMS Discussion
Topic: DIAG speed trouble 20.125553145124616. Blowing-up: Saving latest model state into RESTART file
Replies: 6
Views: 3801

Re: DIAG speed trouble 20.125553145124616. Blowing-up: Saving latest model state into RESTART file

this configuration is for ROMS by itself. the model runs for a while, as you say, and then 13510 2010-01-04 03:03:20.00 5.896200E-02 1.757339E+04 1.757345E+04 1.160390E+16 (127,016,32) 1.571469E-02 1.362909E-02 7.575998E+00 2.012555E+01 DIAG speed trouble 20.125553145124616 ... Blowing-up: Saving la...
by jcwarner
Thu Mar 03, 2022 1:08 pm
Forum: ROMS Discussion
Topic: boundary condition
Replies: 2
Views: 2266

Re: boundary condition

" Cha for free-surface and Rad for momentum in 2D and 3D mode, the boundary file is not read" Yes, those BC's do not need any external data. so a file does not need to be read. " RadNud for free-surface or Fla mode for momentum in 2D mode, the input boundary file is usually read"...
by jcwarner
Sat Feb 26, 2022 3:35 pm
Forum: ROMS Discussion
Topic: shallow water blow up due to large wind speed or too much cooling?
Replies: 3
Views: 4107

Re: shallow water blow up due to large wind speed or too much cooling?

difficult to assess why the model would blow up with that limited information. if you want to run the simulation and have the blow up occur, look at the rst file, his file time steps right before the blow up, etc. that can help to diagnose what may have been going on. look at water levels, currents,...
by jcwarner
Fri Feb 04, 2022 2:52 am
Forum: ROMS Discussion
Topic: Curvilinear grid and velocity at rho-points problem
Replies: 2
Views: 1998

Re: Curvilinear grid and velocity at rho-points problem

what does the grid angle look like?
also, when you make those pcolor plots, are they
pcolorjw(lon_rho,lat_rho,angle)
and use
pcolorjw(lon_u,lat_u,ubar)
etc

You said the u and v are ok, but the u_east and v_north are not? so that might be an angle issue.
by jcwarner
Wed Jan 05, 2022 12:18 am
Forum: ROMS Bugs
Topic: SINGLE_PRECISION bug
Replies: 0
Views: 2841

SINGLE_PRECISION bug

small but tricky bug when #define SINGLE_PRECISION 1) in get_idata we have real(r8) :: time_save = 0.0_r8 but it needs to be real(dp) :: time_save = 0.0_dp 2) And then there are multiple uses of: time_save=time(ng) time(ng)=8640000.0_r8 .... call get2d.. time(ng)=time_save I think this should be : t...
by jcwarner
Tue Jan 04, 2022 2:30 pm
Forum: ROMS Problems
Topic: Question about tracer routine
Replies: 4
Views: 5208

Re: Question about tracer routine

you need to be careful in here, as some routines (like MPDATA) are a little different than others.
but at the bottom of step3d_t is where the units become Tracer only. it can happen in the diffusion step, and here again you need to be careful if splines are on or not.
by jcwarner
Thu Dec 30, 2021 10:49 pm
Forum: ROMS Discussion
Topic: Adding a river runoff in simulation
Replies: 2
Views: 4090

Re: Adding a river runoff in simulation

i just downloaded the code the other day, and the RIVERPLUME1 test case worked fine for me. give that a try. i dont understand if you are using RIVERPLUME1 or trying to use that as a guide. if it is a guide, then you need to define your own application in ana_psource and copy those lines of code to ...
by jcwarner
Wed Dec 29, 2021 3:23 pm
Forum: ROMS Discussion
Topic: qeustion about ROMS+SWAN(COAWST)
Replies: 1
Views: 2520

Re: qeustion about ROMS+SWAN(COAWST)

sorry for the very long delay on this response. i had lost this thread. when coupled to swan, the wec breaking is the wave dissipation from swan. it comes from swan as a variable DISSURF and put into roms variable Dissip_break. That Dissip_break is used: - #define WEC_VF then it is used to drive cur...
by jcwarner
Tue Dec 28, 2021 3:31 pm
Forum: ROMS Discussion
Topic: usage of both LuvSrc and LwSrc
Replies: 4
Views: 8217

Re: usage of both LuvSrc and LwSrc

starting with the most recent version of Rutgers code, I used the RIVER_PLUME1 test case because it has 1 UV_Psource location for a river entering a N-S Periodic domain with East and West walls. Therefore, the entire volume change within the domain is due to the river inflow. I ran the test case 'as...
by jcwarner
Tue Dec 28, 2021 2:13 pm
Forum: ROMS Installation
Topic: why my roms src has no WAVE?
Replies: 1
Views: 2504

Re: why my roms src has no WAVE?

if you want to use a coupled roms+swan, then you need to use COAWST
https://github.com/jcwarner-usgs/COAWST
-j
by jcwarner
Mon Dec 27, 2021 2:29 pm
Forum: ROMS Discussion
Topic: bottom boundary layer height calculation (Hbbl)
Replies: 2
Views: 2144

Re: bottom boundary layer height calculation (Hbbl)

i did a grep -nir hbbl ROMS/* and got a lot of results, but i think the most relevant is in ROMS/Nonlinear/lmd_bkpp.F and it says "This routine determines the depth of bottom oceanic boundary ! ! layer, hbbl, as the deepest depth where the bulk Richardson ! ! number is equal to the critical val...
by jcwarner
Mon Dec 27, 2021 2:27 pm
Forum: ROMS Discussion
Topic: Inlet_test:the Swan output has no result
Replies: 1
Views: 1882

Re: Inlet_test:the Swan output has no result

if you are going to use ROMS + SWAN coupling, then please use the COAWST code here:
https://github.com/jcwarner-usgs/COAWST

that should run 'as is' and give the outputs, this is an example test case.
thanks,
john
by jcwarner
Wed Dec 22, 2021 2:19 pm
Forum: ROMS Discussion
Topic: READ_PHYPAR - Invalid dimension parameter
Replies: 11
Views: 9946

Re: READ_PHYPAR - Invalid dimension parameter

there are probably a few things not correct here, you will need to try and figure them out. for swan Err file01 -001 : this is for grid 1, tile 1 (so probably left side) it has Error : invalid boundary point Error : ( -99 -99) is outside computational grid so most likely one of the Boundspec locatio...
by jcwarner
Tue Dec 21, 2021 1:44 pm
Forum: ROMS Discussion
Topic: READ_PHYPAR - Invalid dimension parameter
Replies: 11
Views: 9946

Re: READ_PHYPAR - Invalid dimension parameter

half right, half not. for ini and clm: you can use the parent ini and clm to make the first child ini and clm. Then use those first child ini and clm to make the files for the 3rd grid. for the contact: roms only wants 1 contact file, so to create a contact file for 3 grids you need to make all 3 gr...
by jcwarner
Mon Dec 20, 2021 1:36 pm
Forum: ROMS Discussion
Topic: READ_PHYPAR - Invalid dimension parameter
Replies: 11
Views: 9946

Re: READ_PHYPAR - Invalid dimension parameter

this is now a SWAN error. I see in Errfile02-001 (this is for swan grid 2, tile partition 1) Error : invalid boundary point Error : ( -99 -99) is outside computational grid Error : invalid boundary point Error : ( -99 -99) is outside computational grid Error : File is already opened: Projects/Tembin...
by jcwarner
Sat Dec 18, 2021 3:43 pm
Forum: ROMS Discussion
Topic: READ_PHYPAR - Invalid dimension parameter
Replies: 11
Views: 9946

Re: READ_PHYPAR - Invalid dimension parameter

this is for ROMS.
It looks like your ocean.in has

Ngrids = 2
NestLayers = 2

but it must have been compiled without nesting. the error is saying that you can only have NestLayers = 1 for non-nesting.
think you need to add
#define NESTING
to *.h file and recompile.
-j
by jcwarner
Wed Dec 15, 2021 5:46 pm
Forum: ROMS Problems
Topic: Cannot run 2D only mode with ROMS3.9 or higher
Replies: 9
Views: 9137

Re: Cannot run 2D only mode with ROMS3.9 or higher

wow! thanks for diggin in. But i see in mod_ncparam.F

ic=5
DO i=1,MT
ic=ic+1
isTvar(i)=ic
END DO
....
isWvel=ic+1

are you saying ifort can not do "DO i=1,0" ?

there still should be a value for ic =5.
by jcwarner
Mon Dec 13, 2021 8:39 pm
Forum: ROMS Problems
Topic: Cannot run 2D only mode with ROMS3.9 or higher
Replies: 9
Views: 9137

Re: Cannot run 2D only mode with ROMS3.9 or higher

what is on line 904 of mod_ncparam.f90
(the f90, not the F)
mod_ncparam_mp_in 904 mod_ncparam.f90
by jcwarner
Tue Nov 30, 2021 1:25 pm
Forum: ROMS Installation
Topic: COAWST compilation error -- wec_vf.f90
Replies: 4
Views: 4527

Re: COAWST compilation error -- wec_vf.f90

can you post the issue here:
https://github.com/jcwarner-usgs/COAWST/issues
and include the output during the run or during the build?
thanks,
j
by jcwarner
Mon Nov 22, 2021 5:37 pm
Forum: ROMS Discussion
Topic: usage of both LuvSrc and LwSrc
Replies: 4
Views: 8217

Re: usage of both LuvSrc and LwSrc

thanks Hernan. yes of course. the user would create a single file for all the points (LuvSrc and Lwsrc). Then for each point, the user would define a Dsrc value. If Dsrc is 0, then it is for a u point (LuvSrc). If Dsrc is 1, then it is for a v point (LuvSrc). If Dsrc is -1, then it is for rho point ...
by jcwarner
Mon Nov 22, 2021 3:42 pm
Forum: ROMS Discussion
Topic: usage of both LuvSrc and LwSrc
Replies: 4
Views: 8217

usage of both LuvSrc and LwSrc

Dear ROMS- I am trying to use both LuvSrc and LwSrc at the same time. However, the code currently is not allowing that configuration, but i need to make it happen. Can we make some simple modifications, such as defining a new direction of -1 to be used for the LwSrc? This will need some minor code c...
by jcwarner
Mon Nov 08, 2021 1:04 pm
Forum: ROMS Discussion
Topic: Vertical and horizontal mixing
Replies: 1
Views: 2164

Re: Vertical and horizontal mixing

the horizontal viscosity is usually a function of grid size and horizontal velocities. the horiz mixing also has several choices, https://www.myroms.org/wiki/Horizontal_Mixing the vertical viscosity is usually a function of the vertical stratification (temp, salt, sediment sometimes) and the vertica...
by jcwarner
Mon Nov 01, 2021 11:41 pm
Forum: ROMS Bugs
Topic: tiling and divide by zero error
Replies: 10
Views: 17334

Re: tiling and divide by zero error

is ​there any more relevant info? i assume you change the NtileI NtileJ values accordingly.
also, looks like the error might be in output. did you change anything in wrt_his or any other place to write things out for your setup?
-j
by jcwarner
Mon Nov 01, 2021 5:21 pm
Forum: ROMS Discussion
Topic: The intensity of tropical cyclone
Replies: 3
Views: 3723

Re: The intensity of tropical cyclone

i guess i was not really clear in my suggestion, sorry. The plots you just posted show the initial conditions for the TC centered at 12Z02Nov2017. What do the initial conditions look like west of that region, where the TC will make landfall? The WRF CTL SST looks warmer than the WRF-sstupdate and th...
by jcwarner
Sun Oct 31, 2021 2:00 am
Forum: ROMS Discussion
Topic: The intensity of tropical cyclone
Replies: 3
Views: 3723

Re: The intensity of tropical cyclone

What did the SST for ROMS look like before the TC crossed into the domain, and compare that to the SST from the WRF control.
How did you set the initial vertical stratification for ROMS?
-j
by jcwarner
Tue Oct 26, 2021 12:35 pm
Forum: ROMS Installation
Topic: COAWST compilation error when coupling ROMS with other models
Replies: 2
Views: 2894

Re: COAWST compilation error when coupling ROMS with other models

can you post the entire output of the build? there might have been something earlier.
and can you create a new post for this issue here:
https://github.com/jcwarner-usgs/COAWST/issues

thanks,
john
by jcwarner
Tue Sep 28, 2021 1:03 pm
Forum: ROMS Installation
Topic: output files of SWAN model
Replies: 7
Views: 8057

Re: output files of SWAN model

i think for some versions of gfortran you need to add

FFLAGS += -fconvert=big-endian

see if that works.
what do you get for
gfortran --version
by jcwarner
Thu Sep 23, 2021 1:44 pm
Forum: ROMS Discussion
Topic: Question about the typhoon-induced sea surface cooling in ROMS-WRF coupled simulation
Replies: 3
Views: 4617

Re: Question about the typhoon-induced sea surface cooling in ROMS-WRF coupled simulation

it looks like the ocean SST is higher in roms all around than in the GHRSST. so maybe your initial conditions are not matching the SST data even before the storm. How did you set the initial stratification? perhaps that could be too thin of a surface mixed layer allowing the cooler temp to mix upwar...
by jcwarner
Thu Sep 23, 2021 1:39 pm
Forum: ROMS Problems
Topic: Overmixing in stripes
Replies: 5
Views: 5739

Re: Overmixing in stripes

all this limit does is put a maximum value on the AKt and AKv. if you set it to 1e-5 (a very low example) then the mixing will never get above those values (not recommended). If you set it to 10, then the mixing will never go above 10. It doesnt mean that the mixing is 10 everywhere, is just means t...
by jcwarner
Wed Sep 22, 2021 5:26 pm
Forum: ROMS Problems
Topic: Overmixing in stripes
Replies: 5
Views: 5739

Re: Overmixing in stripes

you stated that the stripes go away if you undef the AKV limit. i see your ocean.in has AKT_LIMIT == 1.0d-3 1.0d-3 ! m2/s AKV_LIMIT == 1.0d-3 1.0d-3 ! m2/s if you undef the limit, what are the values of AKV? if they are above 1e-3, then your limiters are preventing the mixing. that could be the issu...
by jcwarner
Fri Sep 10, 2021 12:23 am
Forum: ROMS Discussion
Topic: u and u_stokes & zeta and zeta_w clarification
Replies: 2
Views: 2140

Re: u and u_stokes & zeta and zeta_w clarification

this was a discussion point, and we decided to always have u- v- ubar- vbar be Eulerian. then we also can write out u_stokes, v_stokes, ubar_stokes, vbar_stokes. the lagrangian would then be eulerian + stokes. "In the history file, does u = u_e + u_stokes, v = v_e + v_stokes, and zeta_c = zeta ...
by jcwarner
Fri Sep 10, 2021 12:17 am
Forum: ROMS Discussion
Topic: Horizontal surface pressure gradient forcing in 2D and 3D
Replies: 18
Views: 23752

Re: Horizontal surface pressure gradient forcing in 2D and 3D

GRho0 * z_w = 1000.0_r8 * GRho * z_w = 1000.0_r8* g/rho0 * z_w so this is g*z_w * 1000/rho0 and you want it to be g*z_w well, the difference is 1000/1025 = 0.97, as you said. But it does not have the wrong units. The units are correct. The 1000/rho0 is density over density and those units cancel. I ...
by jcwarner
Thu Sep 09, 2021 9:42 pm
Forum: ROMS Discussion
Topic: Horizontal surface pressure gradient forcing in 2D and 3D
Replies: 18
Views: 23752

Re: Horizontal surface pressure gradient forcing in 2D and 3D

I think the original code was correct. P is only an internal variable for this routine. Maybe you want to re-define it for the tides, but the way it was originally seemed correct. It was: P(i,j,N(ng))=GRho0*z_w(i,j,N(ng))+ & & GRho*(rho(i,j,N(ng))+cff2)* & & (z_w(i,j,N(ng))-z_r(i,j,N...
by jcwarner
Wed Aug 25, 2021 6:47 pm
Forum: ROMS Discussion
Topic: Thin Wall Between Grid Cells
Replies: 3
Views: 3181

Re: Thin Wall Between Grid Cells

just try the adjustment to just mask _u _v for now and see if that works. it may look strange, as there could be dzeta/dx gradients across the rho cells, but no flow. but give it a try.
-j
by jcwarner
Wed Aug 25, 2021 4:59 pm
Forum: ROMS Discussion
Topic: Thin Wall Between Grid Cells
Replies: 3
Views: 3181

Re: Thin Wall Between Grid Cells

Typically to provide a permanent wall, then use the mask _u _v _rho variables. ROMS is solved on a staggered grid, so when you set the mask_rho, then edit_mask (or similar) can determine the appropriate settings for mask_u and mask_v and mask_psi. However, it seems you want a special case, where 2 a...
by jcwarner
Mon Aug 09, 2021 11:30 pm
Forum: ROMS Discussion
Topic: grid cell switches from dry to wet while (u,v) are NaN
Replies: 1
Views: 1946

Re: grid cell switches from dry to wet while (u,v) are NaN

the cells will become wet when ever the total water depth is > Dcrit. The water will move from all the processes in the model, such as pressure gradients from gravity or density, wind, horiz mixing, etc. At the model init, how did you set the water level (zeta) on the grid? did you use an init file,...
by jcwarner
Wed Aug 04, 2021 1:56 pm
Forum: ROMS Discussion
Topic: How to erase sediment-variables in ROMS?
Replies: 9
Views: 4492

Re: How to erase sediment-variables in ROMS?

if go back to your original setup, and activate #define ANA_SEDIMENT this will create the vars needed. to do this, i suggest you - copy the ROMS/Functionals/ana_sediment.h to your project folder -edit build.bash or whatever to set the path to the ANALYTICAL files -edit ANA_SEDIMENT and copy or chang...
by jcwarner
Mon Aug 02, 2021 6:44 pm
Forum: ROMS Information
Topic: conversion of sigma levels in to depth
Replies: 4
Views: 8340

Re: conversion of sigma levels in to depth

wow, i originally answered this in 2006. there are some different options now. i am not sure if that m file is still around, but you can also use get_roms_grid.m and you can enter the grid name, or this ocean_his name, and the time index from a history file. that m file has example input commands at...
by jcwarner
Fri Jul 16, 2021 12:20 pm
Forum: ROMS Problems
Topic: SWAN encountered 'Unexpected end of file while reading DUMMY', produces constant SWH
Replies: 3
Views: 3902

Re: SWAN encountered 'Unexpected end of file while reading DUMMY', produces constant SWH

not sure what you mean '2 more lines for the date." please read the swan manual for the proper format.
you need to have data for each cell
424*24*227*181*2=836202624
or something close to taht. you have about half of that, so maybe that is why it stops short.
by jcwarner
Thu Jul 15, 2021 12:45 pm
Forum: ROMS Problems
Topic: SWAN encountered 'Unexpected end of file while reading DUMMY', produces constant SWH
Replies: 3
Views: 3902

Re: SWAN encountered 'Unexpected end of file while reading DUMMY', produces constant SWH

i can give this a try. suggest you look in the swan manual for the syntax on the wind lines. There are more lines for the wind, something like: && KEYWORD TO CREATE WIND GRID && INPGRID WIND REGULAR -105 10 0 220 160 0.25 0.25 & NONSTATIONARY 20121028.000000 6 HR 20121031.000000 ...
by jcwarner
Wed Jul 14, 2021 1:53 pm
Forum: ROMS Discussion
Topic: Plotting vector velocity u and w and bedload_coeff
Replies: 1
Views: 2111

Re: Plotting vector velocity u and w and bedload_coeff

sorry for late reply. 1. i think u-east and v-north have been averaged to rho points (?). but normally the u- and v- points are on a staggered grid, see here: https://www.myroms.org/wiki/Fractional_Coordinate_System_(%CE%BE_-_%CE%B7_space) and similar pages. you can interpolate the velocities to the...
by jcwarner
Sat Jul 10, 2021 12:12 pm
Forum: ROMS Bugs
Topic: Error to run roms_clm_roms_master_climatology_coawst_mw
Replies: 6
Views: 5465

Re: Error to run roms_clm_roms_master_climatology_coawst_mw

if it works sometimes and not others, then it is your internet connection or the thredds server on the other end.
suggest you try from a hard wired connection (not wifi), and try to acquire data for a short period of time.

-j
by jcwarner
Thu Jun 24, 2021 12:29 pm
Forum: ROMS Sediment
Topic: Bed armoring
Replies: 4
Views: 4407

Re: Bed armoring

yes the armoring effect would not work well if the top layer is very thick. if you are interested in having multiple layers you may need to run some test simulations to see what would work best for you. there is a flag 'new bed layer thick' set in the sediment.in that tells the model when to combine...
by jcwarner
Wed Jun 23, 2021 1:37 pm
Forum: ROMS Sediment
Topic: Bed armoring
Replies: 4
Views: 4407

Re: Bed armoring

i am not sure what NONCOHESIVE_BED1 is, but for bed armoring- there is not an explicit armoring routine. Instead, as the finer grain material is eroded, the surface will become composed of coarser grain sizes. this coarsening will increase the bottom(i,j,itauc), and thus create a higher critical str...
by jcwarner
Mon Jun 14, 2021 3:28 pm
Forum: ROMS Usage
Topic: Queries regarding possible use of modified test_chan case in ROMS
Replies: 4
Views: 6660

Re: Queries regarding possible use of modified test_chan case in ROMS

maybe it needs a leading /
change
home/sbasak08/roms/Projects/test_chan/External/sediment_test_chan.in
to
/home/sbasak08/roms/Projects/test_chan/External/sediment_test_chan.in

??
by jcwarner
Mon Jun 14, 2021 2:39 pm
Forum: ROMS Usage
Topic: Queries regarding possible use of modified test_chan case in ROMS
Replies: 4
Views: 6660

Re: Queries regarding possible use of modified test_chan case in ROMS

I think you are modifying the ana_sediment in ROMS/Functionals ??? what did you set for the path to the analytical files, in the build.bash or whatever you are using? export MY_HEADER_DIR=${MY_PROJECT_DIR}/?????? export MY_ANALYTICAL_DIR=${MY_PROJECT_DIR}/???? - what did you set here??? I suggest yo...
by jcwarner
Thu Jun 03, 2021 2:25 am
Forum: ROMS Bugs
Topic: dogbone composite nesting: contact()%Kdg in nesting.f90 not assigned?
Replies: 2
Views: 5617

Re: dogbone composite nesting: contact()%Kdg in nesting.f90 not assigned?

i think there is a issue with the updated calls for z_weights. There is a check in set_contact that has IF (.not.ANY(Lcoincident).and.ANY(Lcomposite)) THEN get_Vweights=.TRUE. ELSE get_Vweights=.FALSE. END IF i think we need this to allow get_vweights to be true. then at top of nesting we have IF ((...
by jcwarner
Mon May 31, 2021 11:05 pm
Forum: ROMS Problems
Topic: rectangles delimited by anomalous values
Replies: 4
Views: 5110

Re: rectangles delimited by anomalous values

i have seen this problem if the ocean.in does not have the updated section to select the tracer advection schemes. Check to make sure your input file is current.
by jcwarner
Thu May 20, 2021 12:50 pm
Forum: ROMS Installation
Topic: Error to compile WRF+SWAN
Replies: 1
Views: 3285

Re: Error to compile WRF+SWAN

ok. so you have set #undef ROMS_MODEL #undef NESTING #define WRF_MODEL #define SWAN_MODEL #undef WW3_MODEL #define MCT_LIB #define MCT_INTERP_OC2AT #define MCT_INTERP_WV2AT #define MCT_INTERP_OC2WV if you only use swan and wrf, then you do not need INTERP_OC2AT or INTERP_OC2WV so set this as #undef ...
by jcwarner
Mon May 17, 2021 5:26 pm
Forum: ROMS Discussion
Topic: Fatal error in MPI_Send: Other MPI error, error stack:
Replies: 1
Views: 3357

Re: Fatal error in MPI_Send: Other MPI error, error stack:

i see "MCT::m_SparseMatrixPlus::initFromRoot_:: FATAL--length of vector y different from row count of sMat.Length of y = 128250 Number of rows in sMat = 127516 01E.MCT(MPEU)::die.: from MCT::m_SparseMatrixPlus::initFromRoot_()" so this is probably due to the scrip weights not being compute...
by jcwarner
Mon May 10, 2021 12:53 pm
Forum: ROMS Discussion
Topic: Question about SST in WRF-ROMS coupling
Replies: 6
Views: 6412

Re: Question about SST in WRF-ROMS coupling

something is making roms get 2-3 degrees cooler in 6 hour is quite extreme. I do not think it is the waves.
what if you ran roms by itself, with surface forcing from wrf or other source?

-j
by jcwarner
Fri May 07, 2021 1:54 pm
Forum: ROMS Discussion
Topic: Question about SST in WRF-ROMS coupling
Replies: 6
Views: 6412

Re: Question about SST in WRF-ROMS coupling

Looking back at the original post, i see that the SST from NCEP and the SST from ROMS are very different. So the 'coupling' is not the issue. It seems that the issue is "how did you set up the initial conditions for ROMS?" If ROMS is starting out 2-3 degrees cooler, then you cant expect si...
by jcwarner
Wed Apr 28, 2021 12:28 pm
Forum: ROMS Discussion
Topic: How to run COAWST model with real-time data
Replies: 6
Views: 8211

Re: How to run COAWST model with real-time data

you may have to do a little work for yourself, modify existing tools, or create new tools for your own particular application. For example, the HYCOM variable 'u' has now been replaced in some locations with the variable 'water_u', etc. If you look here: http://tds.hycom.org/thredds/catalog.html thi...
by jcwarner
Wed Apr 28, 2021 1:03 am
Forum: ROMS Discussion
Topic: COAWST zeta output in VF mode
Replies: 1
Views: 2222

Re: COAWST zeta output in VF mode

from my understanding, the output 'zeta' does not include zetaw, bh, or qsp. Those terms are separate, but i have not looked at Kumar's paper in a long time. Suggest you look here for more reference. I can spend some time to look deeper at some point in the future. https://www.sciencedirect.com/scie...
by jcwarner
Wed Apr 28, 2021 12:55 am
Forum: ROMS Discussion
Topic: Using wind data to model wind induced current when coupling ROMS and SWAN
Replies: 1
Views: 1966

Re: Using wind data to model wind induced current when coupling ROMS and SWAN

Leah- i am not completely clear of what you are saying. If you are looking to couple ROMS and SWAN, that has been accomplished. It has taken lots of work and many years. I would recommend you use the COAWST system here: https://github.com/jcwarner-usgs/COAWST and read Kumar's paper: https://www.scie...
by jcwarner
Thu Apr 22, 2021 12:43 pm
Forum: ROMS Problems
Topic: Error aobut COAWST parallel mode
Replies: 2
Views: 2700

Re: Error aobut COAWST parallel mode

I think you are barking up the wrong tree. it is most likely not an mpi issue. Now that you changed to mpich you need to make sure that the mpirun call points to mpich, and that all the libraries (hdf5, netcdf ,etc ...) are all with that flavor of mpi. Suggest you go back to openmpi and fix the mess...
by jcwarner
Tue Apr 20, 2021 3:33 pm
Forum: ROMS Discussion
Topic: Error when run restart file in cygwin
Replies: 2
Views: 2523

Re: Error when run restart file in cygwin

I have been running roms in cygwin for a long time, mostly for debugging small test cases when we add new physics. You do need to be careful when an ascii character file is convert to windows, the ascii line character endings can be problematic, sometimes tabs appear (tabs are bad), etc. there shoul...
by jcwarner
Mon Apr 19, 2021 2:51 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 205091

Re: problem with two-way nesting

it might be on the call to set_depth before the call to fine2coarse of the r3dvars in nesting. just a long shot but if i could get a copy of this setup, i can dig into this.
-j
by jcwarner
Fri Apr 16, 2021 1:23 pm
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 19154

Re: segmentation fault when running the COAWST model

it is not that big of an effort. Just do a diff of your ocean.in to lets say the Sandy ocean.in. that is all. It will probably not be the issue.

Suggest you compile in debug mode, and run that.
by jcwarner
Thu Apr 15, 2021 12:18 pm
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 19154

Re: segmentation fault when running the COAWST model

not sure. if you are using COAWST, make sure you start with a COAWST ocean input file. Compare yours to the one in Inlet Test or Sandy or something like that. We have a few differences to Rutgers, such as NBCfiles and NCLmfiles, etc. Another thing to try is to compile in debug=on and that gives a &q...
by jcwarner
Wed Apr 14, 2021 1:20 pm
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 19154

Re: segmentation fault when running the COAWST model

good that you are making progress. but at this point, there is not much information there to help. look through that entire log file, are there any NaN's? does the data look correct? use ncview or ncdump or matlab or anything to look at the fields in your "mycase-ini.nc" file. maybe start ...
by jcwarner
Tue Apr 13, 2021 4:55 pm
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 19154

Re: segmentation fault when running the COAWST model

so the ana_sediment and sediment.in are different files. If you need sediment, then you can provide initial information in a netcdf file, or by using ana_sediment. - Suggest you copy ana_sediment from the Sandy or inlet_test case, and move that file to the location of your project.h file. -Edit the ...
by jcwarner
Tue Apr 13, 2021 11:51 am
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 19154

Re: segmentation fault when running the COAWST model

the sediment is required if you set SSW_BBL or any other BBL.
you can use ana_sediment and set some basic properties in there.
-j
by jcwarner
Mon Apr 12, 2021 5:08 pm
Forum: ROMS Discussion
Topic: Increase Bottom Vertical Grid Resolution
Replies: 4
Views: 4819

Re: Increase Bottom Vertical Grid Resolution

here is a good start:
https://www.myroms.org/wiki/Vertical_S-coordinate

there is also an m file "roms_getnplot_grid_rutgers" that you can use to see what an adjustment to the various parameters will do to the vertical stretching.
-j
by jcwarner
Mon Apr 12, 2021 12:33 pm
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 19154

Re: segmentation fault when running the COAWST model

it probably has to do with reading the names or numbers of boundary files.
make sure you do not have NBCFILES listed twice or something like that.
by jcwarner
Mon Apr 12, 2021 1:21 am
Forum: ROMS Problems
Topic: segmentation fault when running the COAWST model
Replies: 16
Views: 19154

Re: segmentation fault when running the COAWST model

couple of things.
the error are asking about the Build/*.f90 files, not the .F.

can you look at the correct lines in the .f90 's ?

1629 inp_decode.f90
1557 read_phypar.f90
by jcwarner
Fri Mar 26, 2021 7:50 pm
Forum: ROMS Problems
Topic: Blowing up problem when using sediment
Replies: 12
Views: 11732

Re: Blowing up problem when using sediment

this gets into a world of research. Suggest you spend a good bit of time looking thru literature, and doing some model sensitivity tests. this all depends on fractions of non-cohesive, cohesive, is it mixed, etc. There are other sediment bed options in the code now for mixed beds. Here is an example...
by jcwarner
Tue Mar 23, 2021 1:17 pm
Forum: ROMS Discussion
Topic: Error to create the ww3partition_2TPAR
Replies: 1
Views: 2779

Re: Error to create the ww3partition_2TPAR

can you please post this message on the coawst forum:
https://github.com/jcwarner-usgs/COAWST/issues
thanks, -j
by jcwarner
Thu Mar 18, 2021 8:08 pm
Forum: ROMS Discussion
Topic: Fatal error in MPI_Recv
Replies: 13
Views: 14492

Re: Fatal error in MPI_Recv

not sure if you figured this out yet, but maybe i misled you . For the scrip file you need to list the full size of the SWAN grid in your scrip file use SWAN_COORD(1)='../../Projects/Mangkhut/swan_coord.grd', SWAN_BATH(1)='../../Projects/Mangkhut/swan_bathy.bot', SWAN_NUMX(1)=525, SWAN_NUMY(1)=351, ...
by jcwarner
Sun Mar 07, 2021 3:42 pm
Forum: ROMS Problems
Topic: Blowing up problem when using sediment
Replies: 12
Views: 11732

Re: Blowing up problem when using sediment

time to put on your investigator hat. keep digging.
what else is weird there? what is forcing this? winds? tides? what did bottom stress look like? ubar? zeta?
the answer is in there.
by jcwarner
Fri Mar 05, 2021 1:54 pm
Forum: ROMS Discussion
Topic: ** Warnin: Quadruplets will be updated per iteration
Replies: 5
Views: 5320

Re: ** Warnin: Quadruplets will be updated per iteration

this is becoming non-roms related, so can you post this message here:
https://github.com/jcwarner-usgs/COAWST/issues
thanks,
j
by jcwarner
Wed Mar 03, 2021 6:17 pm
Forum: ROMS Discussion
Topic: ** Warnin: Quadruplets will be updated per iteration
Replies: 5
Views: 5320

Re: ** Warnin: Quadruplets will be updated per iteration

for the quadruplets, swanmain.F says "! ! Here the various problems with quadruplets are checked ! The combination of quadruplets and sectors is an error ! in the calculation of quadruplets when the SECTOR option is ! used in the CGRID command. This error should be corrected ! in the future ! I...
by jcwarner
Tue Mar 02, 2021 7:05 pm
Forum: ROMS Problems
Topic: Blowing up problem when using sediment
Replies: 12
Views: 11732

Re: Blowing up problem when using sediment

with any blowup, you need to look at the rst file, his file, and stdout to see what information it is telling you. if you look at the rst file, where is the problem? at the boundary? at the wet/dry masking (for example). try to restart it before the solution occurred. you can typically restart from ...
by jcwarner
Sun Feb 28, 2021 10:05 pm
Forum: ROMS Discussion
Topic: Problem in getting started with SWAN+ROMS
Replies: 1
Views: 3218

Re: Problem in getting started with SWAN+ROMS

Neel- It seems like you could use our COAWST model, it has the same ROMS version and we have added many updates to SWAN Coupling and other capabilities. If you want to get that code here: https://github.com/jcwarner-usgs/COAWST we have a user manual that describes how to setup some basic cases like ...
by jcwarner
Sun Feb 28, 2021 2:50 pm
Forum: ROMS Problems
Topic: Problem when I try to compile ROMS with "make" command
Replies: 1
Views: 2135

Re: Problem when I try to compile ROMS with "make" command

i am not sure about that tutorial, but to build roms you need to use one of the build scripts.

go here
https://www.myroms.org/wiki/Getting_Started
and read that page.

also make sure you see the section about
https://www.myroms.org/wiki/build_Script

that should help a bit.
by jcwarner
Fri Feb 26, 2021 12:03 am
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 205091

Re: problem with two-way nesting

either way please let us know. would like to figure this one out.
by jcwarner
Thu Feb 25, 2021 2:49 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 205091

Re: problem with two-way nesting

we have a 3 nest of USeast - Carolinas - Outer Banks and i am not seeing this behavior of the currents along the boundaries. can you put the grids and the contact files somewhere accessible?
by jcwarner
Wed Feb 24, 2021 6:22 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 205091

Re: problem with two-way nesting

the exterior points of the child grid are replaced with information from the contact file during initialization, in get_grid. At those points, data is interpolated from the parent and the nesting needs characteristics of the parent at those locations. I have found some issues with this for the maski...
by jcwarner
Mon Feb 22, 2021 3:16 am
Forum: ROMS Discussion
Topic: Fatal error in MPI_Recv
Replies: 13
Views: 14492

Re: Fatal error in MPI_Recv

in your scrip file use (1) SWAN_COORD='../../Projects/Mangkhut/swan_coord.grd', SWAN_BATH='../../Projects/Mangkhut/swan_bathy.bot', SWAN_NUMX=524, SWAN_NUMY=350, CARTESIAN=0, SWAN_COORD(1)='../../Projects/Mangkhut/swan_coord.grd', SWAN_BATH(1)='../../Projects/Mangkhut/swan_bathy.bot', SWAN_NUMX(1)=5...
by jcwarner
Mon Feb 22, 2021 3:12 am
Forum: ROMS Discussion
Topic: Error to run script_coawst
Replies: 5
Views: 5673

Re: Error to run script_coawst

i downloaded all of your files, compiled SCRIP, and I can duplicate your error. but there is more information: " Error in netWarning! ***HDF5 library version mismatched error*** The HDF5 header files used to compile this application do not match the version used by the HDF5 library to which thi...
by jcwarner
Sat Feb 20, 2021 3:23 pm
Forum: ROMS Discussion
Topic: Fatal error in MPI_Recv
Replies: 13
Views: 14492

Re: Fatal error in MPI_Recv

the sparse matrix calls are reading the scrip interpolation weights.
What do you have in the scrip.in to compute the interpolation weights?
-j
by jcwarner
Tue Feb 09, 2021 10:21 pm
Forum: ROMS Discussion
Topic: Error to run script_coawst
Replies: 5
Views: 5673

Re: Error to run script_coawst

sorry for the late reply.
did you figure this out?
can you attach the scrip.in file?
by jcwarner
Tue Feb 09, 2021 9:48 pm
Forum: ROMS Installation
Topic: Compiling Error for COAWST
Replies: 6
Views: 10506

Re: Compiling Error for COAWST

are you still having this problem?
can you send the full output of the std out.
-j
by jcwarner
Tue Feb 09, 2021 1:22 pm
Forum: ROMS Discussion
Topic: problem with two-way nesting
Replies: 58
Views: 205091

Re: problem with two-way nesting

I have been working on the nesting for a while now. There are some simple steps that i can take to check on a few things. Is your setup easy to share out?
what ts_advection are you using? time steps?
-j
by jcwarner
Mon Feb 08, 2021 2:51 pm
Forum: ROMS Discussion
Topic: Fatal error in MPI_Recv
Replies: 13
Views: 14492

Re: Fatal error in MPI_Recv

if your roms grid size(h) is 525 351

then, for the SWAN INPUT file, you need to have
INPGRID BOTTOM CURVILINEAR 0 0 524 350 EXC 9.999000e+003
by jcwarner
Thu Feb 04, 2021 9:14 pm
Forum: ROMS Installation
Topic: COAWST compilation error -- wec_vf.f90
Replies: 4
Views: 4527

Re: COAWST compilation error -- wec_vf.f90

looks like a couple of issues.
a year ago or so i changed some parts of the makefile and coawst.bash files and the Compilers files.
are you using most recent versions of them?
which Compilers/ ? file are you using.
-j
by jcwarner
Wed Feb 03, 2021 1:36 pm
Forum: ROMS Discussion
Topic: Fatal error in MPI_Recv
Replies: 13
Views: 14492

Re: Fatal error in MPI_Recv

if you look in the output, you will see that you basically have the same error with different values MCT::m_SparseMatrixPlus::initFromRoot_:: FATAL--length of vector y different from row count of sMat.Length of y = 185152 Number of rows in sMat = 184275 031.MCT(MPEU)::die.: from MCT::m_SparseMatrixP...
by jcwarner
Tue Feb 02, 2021 7:39 pm
Forum: ROMS Discussion
Topic: Fatal error in MPI_Recv
Replies: 13
Views: 14492

Re: Fatal error in MPI_Recv

please read all of the Err* files. This is half way through: "MCT::m_SparseMatrixPlus::initFromRoot_:: FATAL--length of vector y different from row count of sMat.Length of y = 184275 Number of rows in sMat = 183400 010.MCT(MPEU)::die.: from MCT::m_SparseMatrixPlus::initFromRoot_()" so the ...
by jcwarner
Tue Feb 02, 2021 3:16 pm
Forum: ROMS Discussion
Topic: Fatal error in MPI_Recv
Replies: 13
Views: 14492

Re: Fatal error in MPI_Recv

there is not enough information here to figure out what the issue is.
can you dig in a bit deeper?
by jcwarner
Tue Feb 02, 2021 3:13 pm
Forum: ROMS Installation
Topic: Compiling Error for COAWST
Replies: 6
Views: 10506

Re: Compiling Error for COAWST

yes, the compiler flags might be set too strict, and reducing that can help. But fundamentally there is a potential issue that an integer is being broadcasted as a real. I did a search in distribute.F for "CALL mpi_bcast (A, Asize" and the only place that shows up is once near the bottom i...
by jcwarner
Sat Jan 23, 2021 3:42 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 110842

Re: Error while executing ROMS for a self created grid region

those other errors were initiated by the first one, so get the init file method figured out then try again.
by jcwarner
Fri Jan 22, 2021 5:43 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 110842

Re: Error while executing ROMS for a self created grid region

looks closely at the output text, and you will see INP_PAR:FIND_FILE - empty file name string for standard input script KeyWord: ININAME so the code is looking for a file listed in the ININAME field in the ocean.in If you use #define ANA_INITIAL then it will not look for the ininame and you can use ...
by jcwarner
Thu Jan 21, 2021 3:30 pm
Forum: ROMS Problems
Topic: Blowing up problem when using sediment
Replies: 12
Views: 11732

Re: Blowing up problem when using sediment

I am not really clear what you are trying to do but here are some thoughts. The estuary test case has grid cells on the order of 100m (or so, i dont remember exactly). So therefore it can not generate ripples that are on the order of 10's of cm. Instead, we have coded different methods to 'predict' ...
by jcwarner
Wed Jan 20, 2021 12:54 pm
Forum: ROMS Problems
Topic: netcdf output from SWAN
Replies: 3
Views: 5464

Re: netcdf output from SWAN

that is a known issue for me. Currently you can only have netcdf output for one of the grids when using my swan nesting. it is on my list to correct soon. sorry. also, for future reference, can you and others who read this please post non-roms issues on this other github site: https://github.com/jcw...