Search found 813 matches

by jcwarner
Wed Nov 20, 2019 7:32 pm
Forum: ROMS Benchmarks
Topic: Speed-up advice?
Replies: 1
Views: 95

Re: Speed-up advice?

we have been using TS_HSIMT in COAWST. it is a pos def scheme, faster than mpdata. would be nice to know how the tracer advection compares between the 2 schemes. Have you tried different tiling? how many cores / node does it have?
by jcwarner
Tue Nov 12, 2019 3:10 am
Forum: ROMS Discussion
Topic: COAWST grid script help
Replies: 2
Views: 257

Re: COAWST grid script help

in create_roms_xygrid it should have projection='mercator'; and then it saves a file (i just called it temp_jcw33.mat) then mat2roms_mw loads that file into a structure s. as s = load(theMatFile); so in the mat2roms_mw there should be a s.projection, and it should = mercator. do you know how to chec...
by jcwarner
Thu Nov 07, 2019 1:56 pm
Forum: ROMS Sediment
Topic: unable to find variable: sand in input NetCDF
Replies: 5
Views: 429

Re: unable to find variable: sand in input NetCDF

not sure. if you ncview the ini file, does bed thick look ok? do you have bed_thick defined with the correct dimensions
bed_thickness(1:xi_rho,1:eta_rho,1:Nbed,1:length(init_time))
-j
by jcwarner
Tue Nov 05, 2019 12:59 pm
Forum: ROMS Installation
Topic: ROMS-WRF compiling problem
Replies: 5
Views: 2633

Re: ROMS-WRF compiling problem

"***************************************************************************** No environment variable NETCDF set. Stopping ***************************************************************************** " Please go through the log file closely. WRF needs to have the locations of the netcdf libs. Of c...
by jcwarner
Sun Nov 03, 2019 1:44 am
Forum: ROMS Discussion
Topic: error in comile a file
Replies: 2
Views: 145

Re: error in comile a file

probably sounds like you did not set the correct compiler in the build.bash script.
need to set
export FORT=ifort
or
export FORT=gfortran
or
export FORT=pgi
by jcwarner
Thu Oct 31, 2019 6:15 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

this is the Rutgers way of specifying the coupling:
Ocean Export: bath:SSH:Ubar:Vbar:ZO
Waves Export: Wdir:Wamp:Wlen:Wptop:Wpbot:Wubot

i will send you the coawst way in a minute.
-j
by jcwarner
Wed Oct 30, 2019 5:29 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

we have been advancing the coupling to swan a lot more, it will be a more recent version of swan, we also added more wave current interactions, more sediment algorithms, vegetation/marsh, wave watch 3 coupling, wrf coupling, infragravity wave model, new wave-current asymmetry bedload, ... bunch of s...
by jcwarner
Wed Oct 30, 2019 5:07 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

i suggest that you should be using COAWST for roms + swan coupling. We distribute a manual that has explanations of the inlet test and many other test cases. also how to compile and install mct. we dont use arpack. if you need coawst, send me an email to jcwarner@usgs.gov the same swan issue will be...
by jcwarner
Wed Oct 30, 2019 12:45 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

ok i can see the problem now. for some reason, when your code does the cpp preprocessing it changes line 14 from ! equal to REAL2 ) 30.72 to 30.72 the pre-processor tries to clean up the code a bit, and for some reason it is messing this up. It looks like the rest of the Do me a favor and edit SWAN/...
by jcwarner
Wed Oct 30, 2019 3:01 am
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

no i need to see the .f90
i need to see what the compiler is seeing on that line.
by jcwarner
Wed Oct 30, 2019 1:57 am
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

ok the search continues. it has cd /Users/venuscheung/ocean_model/inlet_test_romsG/Build_roms; /usr/local/Cellar/open-mpi/4.0.2/bin/mpif90 -c -frepack-arrays -O3 -I/usr/local/MCT/include -ffixed-form ocpmix.f90 ocpmix.f90:14:9: 14 | 30.72 | 1 Error: Invalid character in name at (1) can you send to m...
by jcwarner
Tue Oct 29, 2019 5:29 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

ok so you need to edit
Compilers/Darwin-gfortran.mk

can you send the whole build.log? so i can see what flags are being issued? Look at the FFLAGS and make sure it is sending the Fixed format options.
by jcwarner
Mon Oct 28, 2019 7:37 pm
Forum: ROMS Discussion
Topic: Readings for different CPP options
Replies: 1
Views: 178

Re: Readings for different CPP options

we have a coupled system call COAWST that has coupled roms + wrf +swan + ww3 and a few other submodules. If you want access to that code, you can send me an email at
jcwarner@usgs.gov
we have a manual with a lot of information and lists about 50 pubs that have used the system.
-john
by jcwarner
Mon Oct 28, 2019 7:32 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

if you type
uname -s
what comes back?
can you send the whole build log?
-j
by jcwarner
Mon Oct 28, 2019 4:08 pm
Forum: ROMS Discussion
Topic: Problem with setting up inlet_test
Replies: 23
Views: 405

Re: Problem with setting up inlet_test

the problem here is that SWAN has characters after column 72 that list when the code was updated. for example, the offending line in ocpmix.F is ! logical function EQREAL ( Checks whether REAL1 is appr. 30.72 and as you can see, the 30.72 is actually swan jargon for the version update. To fix this y...
by jcwarner
Sun Sep 29, 2019 11:14 pm
Forum: ROMS Discussion
Topic: SWAN can't use more than one CPU
Replies: 16
Views: 5096

Re: SWAN can't use more than one CPU

not sure. are you writing to multiple netcdf files? this part is where swan opens files for each processor for writing. comment out swan writing of output just to see if that makes is work.
by jcwarner
Mon Sep 23, 2019 1:30 pm
Forum: ROMS Installation
Topic: Asking for cpp and control files
Replies: 3
Views: 97

Re: Asking for cpp and control files

if you are just trying to run roms, then stick with applications that just use roms. Upwelling is great. Once you get that to work, then start building your own application. There is no need to keep running test cases. If you can find a test case that is close to what you want to do, then run that c...
by jcwarner
Mon Sep 23, 2019 12:38 pm
Forum: ROMS Installation
Topic: Asking for cpp and control files
Replies: 3
Views: 97

Re: Asking for cpp and control files

what is your long term goal? do you want to run a coupled application? or are you just trying to compile the differnt test cases. if you are running a coupled application of roms+swan, then i suggest you use COAWST. you can send me an email at jcwarner@usgs.gov. If you are just running that test to ...
by jcwarner
Mon Sep 16, 2019 5:28 pm
Forum: ROMS Tools and Techniques
Topic: restart error with nested grid
Replies: 4
Views: 111

Re: restart error with nested grid

notice the "\"
by jcwarner
Mon Sep 16, 2019 3:51 pm
Forum: ROMS Tools and Techniques
Topic: restart error with nested grid
Replies: 4
Views: 111

Re: restart error with nested grid

the mud and sand LBC's are in the sediment.in file, such as for Projects/Inlet_test/Refined

LBC(isTvar) == Gra Clo Gra Gra \ ! idsed(:), compact
Nes Nes Nes Nes

Yours would be different on the first line, but the 2nd line should be Nes.
by jcwarner
Mon Sep 16, 2019 3:17 pm
Forum: ROMS Discussion
Topic: NGCname
Replies: 2
Views: 120

Re: NGCname

for non-nesting make sure you have
Ngrids = 1
NestLayers = 1
GridsInLayer = 1
What is the actual error?
by jcwarner
Tue Sep 10, 2019 12:02 pm
Forum: ROMS Bugs
Topic: DIAG speed trouble
Replies: 4
Views: 303

Re: DIAG speed trouble

you have several options activated. suggest you try a more simple first approach and then add complexity. Try it without the river forcings and maybe without the volcons stuff and see if that works to start.
-j
by jcwarner
Mon Sep 09, 2019 2:44 pm
Forum: ROMS Discussion
Topic: Lat_u in ROMs output
Replies: 2
Views: 97

Re: Lat_u in ROMs output

sounds strange. check to see if you can get other variables, it might just be the server is not getting the data to you. if you can get lon_rho and lat_rho, or even better lon_psi and lat_psi, then you can compute the lon_u and lat_u and lon_v and lat_v to a reasonable degree of accuracy as compared...
by jcwarner
Wed Sep 04, 2019 11:44 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : obvious error in the results
Replies: 4
Views: 282

Re: ROMS + SWAN with nesting : obvious error in the results

BLOCK 'COMPGRID' NOHEADER 'NorthSea_wind.mat' LAY 4 WIND 1. OUTPUT 20110101.000000 1 HR
BLOCK 'COMPGRID' NOHEADER 'NorthSea_nest_wind.mat' LAY 4 WIND 1. OUTPUT 20110101.000000 1 HR

how about pcolor or quiver of the winds?
-j
by jcwarner
Wed Sep 04, 2019 12:13 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : obvious error in the results
Replies: 4
Views: 282

Re: ROMS + SWAN with nesting : obvious error in the results

ok. i think i understand what the issue is. If you just run SWAN by itself with NESTING (swan with 2 grids) you are getting an Hsig that is very large compared to SWAN by itself and other similar tests. Here are some suggestions: - do a pcolor of Hsig from the parent, hold on, pcolor hsig from child...
by jcwarner
Thu Aug 01, 2019 12:36 pm
Forum: ROMS Installation
Topic: Error with upwelling case
Replies: 2
Views: 131

Re: Error with upwelling case

did you build this with MPI=on ?? if so, then you need to run it with something like mpirun -np X ./romsG roms_upwelling.in where X is the number of total compute cores X = NtileI * NtileJ , as set in roms_upwelling.in if you did not use MPI=on, then it would be something like ./romsG < roms_upwelli...
by jcwarner
Tue Jul 30, 2019 2:51 am
Forum: ROMS Problems
Topic: variable not properly initialized in ROMS svn978 ?
Replies: 1
Views: 174

Re: variable not properly initialized in ROMS svn978 ?

if you are using roms+swan, you should use our coawst system. We also do not use the Mellor option anymore, we use vortex force.
if you need access to coawst you can send me an email.
jcwarner@usgs.gov

-j
by jcwarner
Mon Jul 01, 2019 1:54 pm
Forum: ROMS Usage
Topic: Could I only couple the parent domains ROMS and WRF?
Replies: 1
Views: 417

Re: Could I only couple the parent domains ROMS and WRF?

When COAWST runs it will look to see that you have 2 WRF grids selected and then it will try to get data from both wrf grids to roms. So i suggest you do the weights. You can try it without the weights but i dont think it will work.
-j
by jcwarner
Wed Jun 19, 2019 2:19 am
Forum: ROMS Sediment
Topic: implement the point source of sediment in the ocean
Replies: 4
Views: 298

Re: implement the point source of sediment in the ocean

i think it works without adding the river. give it a try.
by jcwarner
Tue Jun 18, 2019 8:42 pm
Forum: ROMS Sediment
Topic: implement the point source of sediment in the ocean
Replies: 4
Views: 298

Re: implement the point source of sediment in the ocean

I have not used this in a while, but there is a LtracerSrc option. Can you look in the ocean.in for this? There is a description at the bottom of the file.
-j
by jcwarner
Mon Jun 17, 2019 1:01 pm
Forum: ROMS Discussion
Topic: Unstable with MY25_MIXING, but GLS works
Replies: 2
Views: 447

Re: Unstable with MY25_MIXING, but GLS works

it has been a while since I was in this part of the code, but one quick thought is that we have surface flux bc's for the gls term that is not included in the my25 method. # if defined CHARNOK Zos_eff(i)=MAX(charnok_alpha(ng)/g*0.5_r8* & & SQRT((sustr(i,j)+sustr(i+1,j))**2+ & & (svstr(i,j)+svstr(i,j...
by jcwarner
Fri Jun 07, 2019 5:28 pm
Forum: ROMS Information
Topic: error in compiling scrip_coawst
Replies: 6
Views: 625

Re: error in compiling scrip_coawst

let just see if you can create a file in that directory. try this command

touch help.txt

and then do
ls -ltr
do you see that file?
by jcwarner
Thu May 30, 2019 2:52 pm
Forum: ROMS Problems
Topic: Stokes velocity fills dry cell uphill
Replies: 3
Views: 288

Re: Stokes velocity fills dry cell uphill

Thank you for the post. A couple of things: - looks like you are compiling for just a 2D application? because the error was in main2d. I have not tried to run a 2d application in a long time. i really suggest that you do SOLVE3D. Many of the physics we have added are for 3D applications, like turbul...
by jcwarner
Tue May 28, 2019 2:52 pm
Forum: ROMS Problems
Topic: Noisy SST
Replies: 8
Views: 440

Re: Noisy SST

do the atm forcing fields look smooth ?
you could try mixing from my25 or gls and then activate N2S2_HORAVG.

-j
by jcwarner
Thu May 23, 2019 1:51 am
Forum: ROMS Bugs
Topic: typos in radiation_stress.F
Replies: 8
Views: 2715

Re: typos in radiation_stress.F

we dont really use the Mellor approach any more. suggest that you use the vortex force method (based on wave energy dissipation instead of wave energy gradients). if you want to use the WEC method, we have that in our coawst system. send me an email and i can add you to the distribution. jcwarner@us...
by jcwarner
Fri May 17, 2019 3:09 am
Forum: ROMS Problems
Topic: Performance issues using floats
Replies: 4
Views: 284

Re: Performance issues using floats

as you increase the number of procs for the job, then the tile size decreases. 300x300 tiled 8x14 gives a tile 38 x 20 cells. that is not very big. the job is spending more time exchanging information than performing computations. here are some previous posts in this issue https://www.myroms.org/for...
by jcwarner
Wed May 15, 2019 1:08 pm
Forum: ROMS Sediment
Topic: unable to find variable: sand in input NetCDF
Replies: 5
Views: 429

Re: unable to find variable: sand in input NetCDF

have you been able to work on this? If you activate SEDIMENT, then the model needs some initial fields of bed thickness, porosity, etc. you can provide this information in ANA_SEDIMENT, or you need to include that information in the roms init nc file. we have some m files to create roms init with se...
by jcwarner
Mon May 13, 2019 6:42 pm
Forum: ROMS Problems
Topic: Tide Phase Problem
Replies: 4
Views: 324

Re: Tide Phase Problem

ok. i forgot how the input went. mine does the same: tide_Ephase in degress in netcdf, varinfo converts to rads.
Do you use ramp_tides? What if you adjust tide_start in the ocean.in.
-j
by jcwarner
Mon May 13, 2019 6:21 pm
Forum: ROMS Problems
Topic: Tide Phase Problem
Replies: 4
Views: 324

Re: Tide Phase Problem

basic Q: but are you adding the phase in radians? Ephase would have the most impact. are you adding to all the freqs or just some?
by jcwarner
Tue Apr 23, 2019 1:09 am
Forum: ROMS Usage
Topic: generating nested grid contacts
Replies: 2
Views: 395

Re: generating nested grid contacts

so i know that we tested all this stuff before the workshop in Feb, but maybe we missed something. I dont do the exact procedure you show, i dont use the F=grid_extract.. command. so i am not sure if that is the problem. but here is what worked for our H Sandy test case. after you make the parent gr...
by jcwarner
Fri Apr 19, 2019 5:02 pm
Forum: ROMS Problems
Topic: Segmentation Fault when executing oceanM
Replies: 11
Views: 4449

Re: Segmentation Fault when executing oceanM

- i have not used pgi in awhile, so i did a Google on "Error: no such instruction: `vinserti128 $1,%xmm0,%ymm0,%ymm1'" and found this: https://www.pgroup.com/userforum/viewtopic.php?t=4904&sid=b1ec7d574c2f357163260c5965caa1f6 looks like a pgi issue. - Another thought is that you would not need paral...
by jcwarner
Thu Apr 18, 2019 2:50 pm
Forum: ROMS Discussion
Topic: Coastal flood risk assessment with ROMS?
Replies: 2
Views: 508

Re: Coastal flood risk assessment with ROMS?

ROMS has a wetting/drying capability and we use many roms+swan coupled simulations regularly. I know others have looked at inundations during storms, for example: Beudin, A., Ganju, N. K., Defne, Z., & Aretxabaleta, A. L. (2017). Physical response of a back-barrier estuary to a post-tropical cyclone...
by jcwarner
Wed Apr 17, 2019 5:50 pm
Forum: ROMS Problems
Topic: Segmentation Fault when executing oceanM
Replies: 11
Views: 4449

Re: Segmentation Fault when executing oceanM

i am not sure what you mean by a similar problem. How about we start over, and you explain the issue with as much information as you can. the last few lines of stdout saying "#0 0xFFFFFFFFFFFFFFFF" dont really tell me anything. so if you have something similar i would need to see the whole std out. ...
by jcwarner
Fri Apr 12, 2019 7:43 pm
Forum: ROMS Usage
Topic: tide time corrections
Replies: 1
Views: 264

Re: tide time corrections

set_tides.F cff=2.0_r8*pi*(time(ng)-tide_start*day2sec) DO itide=1,NTC ... omega=cff/Tperiod(itide) DO, .. DO ... Etide(i,j)=Etide(i,j)+ & & ramp*SSH_Tamp(i,j,itide)* & & COS(omega-SSH_Tphase(i,j,itide)) so it takes ocean_time-tide_start to get omega. if you have dstart and tide_start but start at 0...
by jcwarner
Wed Apr 10, 2019 1:10 pm
Forum: ROMS Installation
Topic: compiling upwelling
Replies: 10
Views: 826

Re: compiling upwelling

well that looks worse.
can you attach the build.bash. want to make sure you set the compiler.
by jcwarner
Wed Apr 10, 2019 11:58 am
Forum: ROMS Installation
Topic: compiling upwelling
Replies: 10
Views: 826

Re: compiling upwelling

so when i build i get something like this: [jwarner@poseidon-l1 breach338]$ ./coawst.bash rm -f -r core *.ipo ./Build /vortexfs1/usgs/jwarner/make_macros.mk makefile:248: INCLUDING FILE /vortexfs1/usgs/jwarner/make_macros.mk WHICH CONTAINS APPLICATION-DEPENDENT MAKE DEFINITIONS cp -f /vortexfs1/apps...
by jcwarner
Mon Apr 08, 2019 12:26 pm
Forum: ROMS Discussion
Topic: blow up after first time step
Replies: 1
Views: 425

Re: blow up after first time step

so if you look at the history file using ncview (or similar), what do the fields look like? you are going to have to dig into the output to see what is happening. Is it having trouble at the open boundary, along the coast, etc. i see you have theta_s = 7 and that might be high. not sure. does it run...
by jcwarner
Mon Apr 08, 2019 11:32 am
Forum: ROMS Installation
Topic: compiling upwelling
Replies: 10
Views: 826

Re: compiling upwelling

you might get to a problem with netcdf, and perhaps the -lnetcdff flag is needed, but right now your problem is: "cp -f /home/kzf/softwares/packman.active/include/netcdf.mod /home/kzf/myroms/projects/upwelling/Build_romsG cp -f /home/kzf/softwares/packman.active/include/typesizes.mod /home/kzf/myrom...
by jcwarner
Fri Mar 29, 2019 12:32 pm
Forum: ROMS Adjoint
Topic: ad passive tracer!
Replies: 1
Views: 320

Re: ad passive tracer!

1. you have ! Logical switches (TRUE/FALSE) to activate tracers point Sources/Sinks ! (like river runoff) and to specify which tracer variables to consider: ! [1:NAT+NPT,Ngrids]. See glossary below for details. LtracerSrc == F F F ! temperature, salinity, inert maybe this needs to be F F T 2. also s...
by jcwarner
Wed Mar 27, 2019 2:11 pm
Forum: ROMS Problems
Topic: Blowing up problem when using sediment
Replies: 3
Views: 652

Re: Blowing up problem when using sediment

How about you look at the history file for these 600000 steps and see what is happening. is the model blowing up at the boundary (upstream or downstream), somewhere in the domain, what does the flow look like, etc. you may have to dig in and see how the solution looks to help figure out what is goin...
by jcwarner
Mon Mar 11, 2019 11:56 am
Forum: ROMS Problems
Topic: Run error
Replies: 1
Views: 304

Re: Run error

this looks like a problem with the paths of your system. Your work environment cannot find some shared objects (.so) for netcdf. This might be due to the PATH, LD_LIBRARY_PATH, or other environment variable. You should ask your computer administrator for some help. maybe something like this: https:/...
by jcwarner
Fri Mar 08, 2019 2:23 pm
Forum: ROMS Bugs
Topic: Pair and MPI version has problem with regridding
Replies: 10
Views: 878

Re: Pair and MPI version has problem with regridding

what happens if you do MPI with one processor?
i can look at the press compensate, we want to use that for our hurricane simulations.
-j
by jcwarner
Thu Mar 07, 2019 2:24 pm
Forum: ROMS Problems
Topic: Test Case
Replies: 1
Views: 294

Re: Test Case

whatever test case you are running, the code can not find the grid
"READ_PHYPAR - Grid 01, could not find input file: roms_grd.nc"

i am not sure the application, but the code can not find that file.
-john
by jcwarner
Thu Mar 07, 2019 2:21 pm
Forum: ROMS Usage
Topic: BLOWING UP during nesting
Replies: 7
Views: 1003

Re: BLOWING UP during nesting

are you using wet/dry for the nest? how shallow do the water depths get?
-j
by jcwarner
Wed Feb 13, 2019 1:20 pm
Forum: ROMS Sediment
Topic: Negative suspended load concentrations near river runoff
Replies: 17
Views: 1211

Re: Negative suspended load concentrations near river runoff

think i found it! in sediment_inp.h there is a section of code CASE ('Qout(iMmass)') Npts=load_l(Nval, Cval, NCS, Ngrids, Lmud) DO ng=1,Ngrids DO itrc=1,NCS i=idBmas(itrc) Hout(i,ng)=Lmud(itrc,ng) END DO END DO but that 6th line should be Qout(i,ng)=Lmud(itrc,ng) So the Hout was being overwritten. C...
by jcwarner
Tue Feb 12, 2019 7:45 pm
Forum: ROMS Sediment
Topic: Negative suspended load concentrations near river runoff
Replies: 17
Views: 1211

Re: Negative suspended load concentrations near river runoff

does it report in the std out at the top with all the other zob, tcline, etc etc stuff, then it gets to sediment, it should report the sed classes, then does it say it will Hout().... and write out the sediment mass??
by jcwarner
Tue Feb 12, 2019 7:09 pm
Forum: ROMS Sediment
Topic: Negative suspended load concentrations near river runoff
Replies: 17
Views: 1211

Re: Negative suspended load concentrations near river runoff

are you doing mud or sand?
NNS =
NCS =
by jcwarner
Tue Feb 12, 2019 4:38 pm
Forum: ROMS Sediment
Topic: Negative suspended load concentrations near river runoff
Replies: 17
Views: 1211

Re: Negative suspended load concentrations near river runoff

i get bed mass written out.
make sure you set it for sand or mud, whichever you have activated.
by jcwarner
Tue Feb 12, 2019 2:32 am
Forum: ROMS Discussion
Topic: How to get boundary condition from new data Hycom for COAWST
Replies: 3
Views: 781

Re: How to get boundary condition from new data Hycom for CO

we have updated all of our m files in the latest coawst distribution (feb 11, 2018). so you can give them a try. As ptimko said, look here http://tds.hycom.org/thredds/catalog.html -=-=-=-=-=-=-=-=-=-=-=-ACTIVE-=-=-=-=-=-=-=-=-=-=-=- -- Folder GLBy0.08/latest (Present + FORECASTS)/ -- Folder GLBy0.0...
by jcwarner
Tue Feb 12, 2019 1:32 am
Forum: ROMS Installation
Topic: test case
Replies: 4
Views: 989

Re: test case

sorry i did not get back to you on this earlier. i am using coawst, which i just updated to the latest roms into our system, and this test case works. The i saw what you were doing: ./oceanS < ROMS/External/sediment_test_chan.in Sorry but you need to point to the ocean.in! well, now it is called rom...
by jcwarner
Tue Feb 12, 2019 1:03 am
Forum: ROMS Sediment
Topic: sediment initial conditions for median grain diameter
Replies: 2
Views: 458

Re: sediment initial conditions for median grain diameter

if you are serious about using the sediment capabilities, i suggest you use our coawst system. we have added more capabilities such as a more advanced bed load formulation, a floc model, and a faster positive definite advection scheme. we distribute m files for init conditions. you can send me an em...
by jcwarner
Mon Feb 11, 2019 4:31 pm
Forum: ROMS Sediment
Topic: Negative suspended load concentrations near river runoff
Replies: 17
Views: 1211

Re: Negative suspended load concentrations near river runoff

so a suggestion would be to change define TS_U3HADVECTION to define TS_MPDATA and see how that looks. Also, just a question: are the river sources all putting flow into the domain? some people use tidal river forcing that oscillates +/- to simulate a tidal flow. when you say negative sediment next t...
by jcwarner
Mon Feb 11, 2019 3:38 pm
Forum: ROMS Sediment
Topic: Negative suspended load concentrations near river runoff
Replies: 17
Views: 1211

Re: Negative suspended load concentrations near river runoff

There are many reasons why there can be negative concentration values, but most likely it may be from the advection scheme. Suggest you use TS_MPDATA, that should remain positive. Also, we have added a new scheme called TS_HSIMT into COAWST that is also positive definite, it is based on some flux li...
by jcwarner
Tue Feb 05, 2019 2:30 pm
Forum: ROMS Problems
Topic: RIP currents
Replies: 1
Views: 359

Re: RIP currents

our COAWST system has swan + roms coupled and we have a rip_current test case as shown in: Kumar, N., Voulgaris, G., Warner, J.C., and M., Olabarrieta (2012). Implementation of a vortex force formalism in the coupled ocean-atmosphere-wave-sediment transport (COAWST) modeling system for inner-shelf a...
by jcwarner
Sat Feb 02, 2019 3:38 pm
Forum: ROMS Usage
Topic: units of variable 'Dwave'
Replies: 6
Views: 744

Re: units of variable 'Dwave'

Well, for some reason SWAN uses the word 'Cartesian' for 2 separate things. SET NAUTICAL = winds/waves dir from or SET CARTESIAN = winds/waves dir towards and now for something completely different COORDINATES CARTESIAN = uses x-y in meters or COORDINATES SPHERICAL = uses lat + lon more info in the ...
by jcwarner
Fri Feb 01, 2019 12:49 pm
Forum: ROMS Usage
Topic: units of variable 'Dwave'
Replies: 6
Views: 744

Re: units of variable 'Dwave'

Lawrence- if you are using roms + swan, i recommend that you use COAWST instead, we have been updating swan there, adding new wave physics, etc. You can send me an email at jcwarner@usgs.gov to get that code. In short, Dwave output is in degrees, direction that the waves are from. (because we set NA...
by jcwarner
Thu Jan 31, 2019 1:10 pm
Forum: ROMS Results
Topic: ROMS Granular model result in all variables
Replies: 8
Views: 1397

Re: ROMS Granular model result in all variables

yes and try
#define N2S2_HORAVG
by jcwarner
Wed Jan 30, 2019 9:14 pm
Forum: ROMS Results
Topic: ROMS Granular model result in all variables
Replies: 8
Views: 1397

Re: ROMS Granular model result in all variables

does it start out like that or evolve to that? check the fields in your grid file like dmde dndx. check angle.
are you using N2S2_HORAVG?

-j
by jcwarner
Tue Jan 29, 2019 1:11 pm
Forum: ROMS Discussion
Topic: Correct CPP options to use ANA_WWAVE?
Replies: 4
Views: 667

Re: Correct CPP options to use ANA_WWAVE?

The ana_wwave is for hwave, dwave, pwave_bot _top, etc = from wind waves, ie periods from ~ 2 to 20 seconds. Is this the wave band that you are interested in? when i hear 'gravity waves' i usually think of longer periods, but just wanted to make sure you knew what we intended for the names.
-john
by jcwarner
Wed Dec 19, 2018 1:22 pm
Forum: ROMS Discussion
Topic: High water speed values in hisfile
Replies: 4
Views: 706

Re: High water speed values in hisfile

the error at the end of the output log is 132000 0001-02-25 00:00:00.00 1.203129E-02 7.265249E+03 7.265261E+03 5.839635E+14 (075,001,15) 7.890818E-04 3.869123E-03 2.404095E-02 1.595520E+00 NETCDF_PUT_FVAR_1D - error while writing variable: ocean_time in input file: /media/STORAGE/prueba_forc_build/o...
by jcwarner
Sat Dec 15, 2018 8:30 pm
Forum: ROMS Discussion
Topic: High water speed values in hisfile
Replies: 4
Views: 706

Re: High water speed values in hisfile

how about looking at some forcing variables near those locations. what does sustr svstr zeta rho look like?? something is making the water move.
-j
by jcwarner
Sat Dec 15, 2018 8:27 pm
Forum: ROMS Bugs
Topic: Equation error in GLS model coded into gls_corstep.F
Replies: 5
Views: 983

Re: Equation error in GLS model coded into gls_corstep.F

Hey Hernan. Before you update the code, how about you let me take a look at it and get back to you next week. I spent a good bit of time to code this in, and in a separate email chain Hans is looking at this as well. Give the people who coded it a chance to look it over first.

thanks,
john
by jcwarner
Fri Dec 07, 2018 1:54 pm
Forum: ROMS Usage
Topic: Operational Coupled Forecast Systems?
Replies: 3
Views: 780

Re: Operational Coupled Forecast Systems?

Forecast systems: Hydro and Agro Informatics Institute http://www.thaiwater.net/v3/wrfroms/rain_forecast_pre/tab1/image1 USGS Woods Hole http://woodshole.er.usgs.gov/project-pages/cccp/public/COAWST.htm North Carolina State University http://omgsrv1.meas.ncsu.edu:8080/ocean-circulation-useast2/ Made...
by jcwarner
Tue Dec 04, 2018 1:09 pm
Forum: ROMS Sediment
Topic: Initial sediment layers
Replies: 3
Views: 548

Re: Initial sediment layers

The sediment part in roms can be run with just roms. The swan coupling provides wave information that can be used by the ocean model to: - compute enhanced surface mixing due to tke wave breaking - computed wave driven flow using vortex force approach - compute enhanced bottom stresses due to combin...
by jcwarner
Tue Dec 04, 2018 12:42 am
Forum: ROMS Usage
Topic: Accessing forcing data from NCEP NOMADS
Replies: 1
Views: 527

Re: Accessing forcing data from NCEP NOMADS

Thanks John. Yes we did take some steps earlier to take into account the http to https issue from NCEI etc. We got hit with that earlier when matlab failed, but Christie and Steve & others spent some time to update some m files, here is a post on that https://www.myroms.org/forum/viewtopic.php?f=23&...
by jcwarner
Sun Dec 02, 2018 3:39 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

can send to me the scrip file you used to create the scrip weights?
and are the grids the same size for the parent roms/swan and the child roms/swan?
you need to make sure you enter the swan grid sizes correctly.
-j
by jcwarner
Sat Dec 01, 2018 2:16 am
Forum: ROMS Problems
Topic: Inlet_test and MCT coupler
Replies: 1
Views: 404

Re: Inlet_test and MCT coupler

looks like the compiler can not find your mct libraries. I have some explanation to build MCT in the coawst manual. Then you need to set MCT_INCDIR and MCT_LIBDIR, as described in the manual. give that a try and see if it helps.
-j
by jcwarner
Sat Dec 01, 2018 2:09 am
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

if i had to guess, i would ask what the time steps are being used for roms.
by jcwarner
Fri Nov 30, 2018 3:24 am
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

So i can get this to work. The grid files seem to be ok (2xbathy.bot and 2x coord.grd). They are ok. Here are some changes i did to the INPUT files: 1) I did not have your wind file, so i commented out the wind commands and added WIND 10. 0. this puts a 10 m/s wind from north on both grids. 2) You h...
by jcwarner
Thu Nov 29, 2018 1:11 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

swan with nesting by itself (no roms) does not need the scrip weights, so there is an issue with the swan child grid. Can you send to me the 2 *.bot and the 2 *.grd swan files (4 files)
jcwarner@usgs.gov
thanks
-j
by jcwarner
Wed Nov 28, 2018 1:42 pm
Forum: ROMS Tools and Techniques
Topic: River input problem forrtl: severe (174): SIGSEGV,
Replies: 9
Views: 1201

Re: River input problem forrtl: severe (174): SIGSEGV,

this section of mod_netcdf is newly modified. which one of these is line 852: DO i=1,MIN(Matts,n_gatt) att_id=i status=nf90_inq_attname(ncid, nf90_global, att_id, & & att_name(i)) IF (status.eq.nf90_noerr) THEN status=nf90_inquire_attribute(ncid, nf90_global, & & TRIM(att_name(i)), & & xtype = att_k...
by jcwarner
Wed Nov 28, 2018 1:21 am
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

i cant see the images. can you upload png files?
-j
by jcwarner
Tue Nov 27, 2018 7:33 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

I am not really clear on these results, and perhaps i was not very clear either. 1) if you just run swan by itself, do you get results on both the parent and child that are consistent and look correct? 2) currents continuous = good, roms refinement works. 3) roms+swan coupled: you say that with it c...
by jcwarner
Tue Nov 27, 2018 1:10 pm
Forum: ROMS Tools and Techniques
Topic: River input problem forrtl: severe (174): SIGSEGV,
Replies: 9
Views: 1201

Re: River input problem forrtl: severe (174): SIGSEGV,

a few things: - does it work with just LuvSrc == T ! horizontal momentum transport and the rest == F ? that way we know it is just for river info, and not the tracers. - the mod sources has MyVarName = TRIM(Vname(1,idRxpo)), & so it is trying to read idRxpo, and trying to read attributes of that var...
by jcwarner
Tue Nov 27, 2018 12:54 am
Forum: ROMS Tools and Techniques
Topic: River input problem forrtl: severe (174): SIGSEGV,
Replies: 9
Views: 1201

Re: River input problem forrtl: severe (174): SIGSEGV,

so this is telling you something: oceanG 0000000002469367 mod_netcdf_mp_net 852 mod_netcdf.f90 oceanG 000000000253DC6C mod_sources_mp_al 109 mod_sources.f90 look in Build/mod_sources.f90 line 109 and it is probably calling mod_netcdf with some information. Look in Build/mod_netcdf.f90 line 852 and s...
by jcwarner
Thu Nov 22, 2018 11:56 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

some things to try: - run it with just swan with 2 grids, not coupled. - do the currents in both grids look continuous? - try to run the swan files with INIT in both of them. this will help the child grid to have some init conditions. if you run the swan without any init commands, it can take a long...
by jcwarner
Thu Nov 22, 2018 10:51 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

the winds should be important for the child.
can you send to me both of your swan input files?
jcwarner@usgs.gov
by jcwarner
Thu Nov 22, 2018 9:58 pm
Forum: ROMS Results
Topic: ROMS + SWAN with nesting : error on the child grid
Replies: 21
Views: 2091

Re: ROMS + SWAN with nesting : error on the child grid

how long into the model run is this? it may take some time for information to propagate into the child.
this depends on the initial conditions. Also, do you have the same winds on both swan grids? you can use copy the parent wind grid to put on the child.
-john
by jcwarner
Tue Nov 20, 2018 12:57 pm
Forum: ROMS Discussion
Topic: Dye boundary doesn't work when using passive tracer
Replies: 4
Views: 686

Re: Dye boundary doesn't work when using passive tracer

i have not added a dye BC in a while, but I think you are missing a few things. - i assume you have something like NPT = 2 ! Number of inactive passive tracers - i think you need to add: ! Each tracer variable requires [1:4,1:NAT+NPT,Ngrids] values. Otherwise, ... LBC(isTvar) == Clo Per Gra Per \ ! ...
by jcwarner
Thu Nov 15, 2018 2:50 pm
Forum: ROMS Information
Topic: TPAR wave directions: correct or not on the ROMS grid angle
Replies: 2
Views: 584

Re: TPAR wave directions: correct or not on the ROMS grid an

no.
The TPAR files go to SWAN, and SWAN knows how to deal with them on a curvilinear grid.
Then Dwave comes from SWAN to ROMS, and the coupler knows how to deal with the grid.
So in summary, just put the data in the TPAR files as you would normally.
-john
by jcwarner
Wed Nov 14, 2018 4:23 pm
Forum: ROMS Bugs
Topic: tiny little bug in get_state.F for time evolving bathymetry
Replies: 2
Views: 566

Re: tiny little bug in get_state.F for time evolving bathyme

yes ROMS/Utility/get_state.F line 496 should be
WRITE (stdout,80) string, TRIM(Vname(1,idbath)), &
you can put that as a src ticket for Hernan to update.
thanks,
-john
by jcwarner
Tue Nov 13, 2018 1:41 pm
Forum: ROMS Sediment
Topic: tracer surface flux
Replies: 2
Views: 598

Re: tracer surface flux

try using this in the project.h file

# define ANA_SMFLUX
# define ANA_BPFLUX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_SPFLUX
# define ANA_SRFLUX
# define ANA_SSFLUX
# define ANA_STFLUX
by jcwarner
Mon Nov 12, 2018 12:23 am
Forum: ROMS Discussion
Topic: Temperature significantly near land areas
Replies: 4
Views: 912

Re: Temperature significantly near land areas

a couple of things jump out here. - first it is best to put temp scales the same so we can see the differences. - second it seems a lot of cooling on east side, looks like a BC issue. - for the temp spikes near land: how is wrf data interpolated to roms grid? we had issues with the coupling at first...
by jcwarner
Wed Nov 07, 2018 2:26 pm
Forum: Meetings/Workshops
Topic: COAWST modeling system Training 2019
Replies: 0
Views: 1008

COAWST modeling system Training 2019

We are announcing a COAWST Training workshop to be held Mon – Thurs Feb 25-28, 2019 on the campus of the North Carolina State University, in Raleigh, NC. If you would like to attend the Training, please send an email to jcwarner@usgs.gov with a full list the people that would like to attend (Name, a...
by jcwarner
Wed Nov 07, 2018 1:27 am
Forum: ROMS Discussion
Topic: NetCDF: Start+count exceeds dimension bound
Replies: 2
Views: 526

Re: NetCDF: Start+count exceeds dimension bound

yes the dimensions are strictly enforced.
If you are new, i would recommend to use m files from someone/somewhere that already work. It will save you a lot of time.
-john
by jcwarner
Tue Nov 06, 2018 1:27 pm
Forum: ROMS Problems
Topic: Transport adjustment on the boundary
Replies: 6
Views: 974

Re: Transport adjustment on the boundary

if you use flather boundary conditions (Fla), it does not guarantee that the mass flux imposed by the ubar on the boundary will be imposed at the boundary. For example one of the choices ends up with : ! ! Western edge, Flather boundary condition. ! .... bry_val=BOUNDARY(ng)%ubar_west(j) cff=1.0_r8/...
by jcwarner
Tue Nov 06, 2018 1:19 pm
Forum: ROMS Installation
Topic: compling build.bash
Replies: 2
Views: 618

Re: compling build.bash

the compilation is using some ana* files (looks like ana_grid.h and ana_vmix.h) I suggest you: - copy those 2 files from ROMS/Functionals and place them in the directory with your project.h file - edit those files and provide the necessary information (they have examples in them). - make sure your b...
by jcwarner
Tue Nov 06, 2018 2:55 am
Forum: ROMS Problems
Topic: Transport adjustment on the boundary
Replies: 6
Views: 974

Re: Transport adjustment on the boundary

i am not really clear what the full story is. I dont use obc_volcons, that is to conserve volume in the domain. but it seems you are not using that. other than that, there is no guarantee that the volume flux across all open boundaries will balance. For example, if you had an estuary with only one o...
by jcwarner
Mon Nov 05, 2018 3:14 pm
Forum: Job Opportunities
Topic: USGS Mendenhall postdoc
Replies: 0
Views: 571

USGS Mendenhall postdoc

We are pleased to announce opportunities for the USGS Mendenhall post doctoral positions. The main web site is here: https://geology.usgs.gov/postdoc/ with specific attention to Research Opportunity Number 17-19 "Hydrodynamics-based numerical modeling for predicting coastal change" where we are seek...
by jcwarner
Fri Nov 02, 2018 12:51 pm
Forum: ROMS Problems
Topic: Growing oscillations when using TS_MPDATA
Replies: 4
Views: 898

Re: Growing oscillations when using TS_MPDATA

are the files for this mostly ana*?
can you send them to me? i could try it here quick.
-j
by jcwarner
Wed Oct 31, 2018 7:33 pm
Forum: ROMS Usage
Topic: How to generate contact netcdf?
Replies: 7
Views: 1124

Re: How to generate contact netcdf?

what exactly is the problem? do you get an error message? the commands would be something like this: 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); Gnames={'Sandy_roms_grid.nc','Sandy_roms_grid...
by jcwarner
Mon Oct 29, 2018 11:33 pm
Forum: ROMS Problems
Topic: Growing oscillations when using TS_MPDATA
Replies: 4
Views: 898

Re: Growing oscillations when using TS_MPDATA

i cant open the ps file. how about a png?
how many grid cells is the domain. mpdata has a big footprint. needs at least 5 cells.
-j
by jcwarner
Tue Oct 16, 2018 12:57 pm
Forum: ROMS Problems
Topic: Cannot define variable, h when creating restart file
Replies: 1
Views: 556

Re: Cannot define variable, h when creating restart file

not sure what is going on here. Can you try it without the changes that you added?
-john
by jcwarner
Thu Oct 04, 2018 8:44 pm
Forum: ROMS Usage
Topic: What would happen if LMD_BKPP is undefined in LMD scheme?
Replies: 5
Views: 865

Re: What would happen if LMD_BKPP is undefined in LMD scheme

Here is another paper that i remembered- https://marine.rutgers.edu/po/Papers/Durski_JGR04.pdf See what Scott did in that paper. I think he undefined the bottom bbl. The code looks to me like lmd_bkpp computes the eddy diffusivity and viscosity, and will replace those values where "! Compute boundar...
by jcwarner
Thu Oct 04, 2018 7:52 pm
Forum: ROMS Usage
Topic: What would happen if LMD_BKPP is undefined in LMD scheme?
Replies: 5
Views: 865

Re: What would happen if LMD_BKPP is undefined in LMD scheme

here is a place to start:
https://www.myroms.org/wiki/Vertical_Mi ... erizations
suggest you also dig into the code a bit. you can grep for those options and they are mostly in Nonlinear/lmd_...

-j
by jcwarner
Wed Sep 26, 2018 1:51 pm
Forum: ROMS Discussion
Topic: ROMS history file has suddenly included restart record
Replies: 7
Views: 913

Re: ROMS history file has suddenly included restart record

i am glad you found this out. I was looking around at all the recent updates and could not see any reason for that change in behavior. if you want to stop the addition of _0001, you could look in Nonlinear/output.F for the code: ifile=(iic(ng)-1)/ndefHIS(ng)+1 IF (Master) THEN WRITE (HIS(ng)%name,10...
by jcwarner
Tue Sep 11, 2018 1:04 pm
Forum: ROMS Usage
Topic: Compiling ROMS error
Replies: 2
Views: 666

Re: Compiling ROMS error

if you set export USE_MPI=on # distributed-memory parallelism export USE_MPIF90= # compile with mpif90 script export USE_OpenMP= # shared-memory parallelism or export USE_MPI=on # distributed-memory parallelism export USE_MPIF90=on # compile with mpif90 script export USE_OpenMP= # shared-memory para...
by jcwarner
Tue Sep 04, 2018 2:37 pm
Forum: ROMS Bugs
Topic: MPI Tiling - Dumb Question
Replies: 5
Views: 937

Re: MPI Tiling - Dumb Question

what does your mpirun command line look like? maybe your cluster also needs some slurm info like
#SBATCH --ntasks=108 # Number of MPI ranks
#SBATCH --nodes=3
#SBATCH --ntasks-per-node=36
or something ????
-j
by jcwarner
Wed Aug 22, 2018 2:09 pm
Forum: ROMS Problems
Topic: forrtl: severe (151): allocatable array is already allocated
Replies: 17
Views: 3434

Re: forrtl: severe (151): allocatable array is already alloc

yeah that is what i was going to look for. TABS = bad.
sometimes an errant character can get stuck in the file and cause an issue.
by jcwarner
Wed Aug 22, 2018 12:06 pm
Forum: ROMS Problems
Topic: forrtl: severe (151): allocatable array is already allocated
Replies: 17
Views: 3434

Re: forrtl: severe (151): allocatable array is already alloc

can you post your ocean.in that causes the trouble?
is there a way I can download it. dont just copy it to this widget.
-j
by jcwarner
Fri Jul 20, 2018 12:26 pm
Forum: ROMS Discussion
Topic: Blowing up problem using GLS mixing
Replies: 8
Views: 1492

Re: Blowing up problem using GLS mixing

if the closure is making it unstable, then look at temp, salt, u, v, and AKv and AKt. perhaps the options of #define LIMIT_VDIFF --> use to impose an upper limit on vertical diffusion ** #define LIMIT_VVISC --> use to impose an upper limit on vertical viscosity ** might help ??? Also, when the depth...
by jcwarner
Fri Jul 13, 2018 2:54 pm
Forum: ROMS Sediment
Topic: Wave-current BBL drag coefficient as output
Replies: 4
Views: 2432

Re: Wave-current BBL drag coefficient as output

Parker- the routine set_vbc.F computes the bottom stresses (in case you want to look at the code). For QDRAG, it essentially uses (for u-direction) bustr = rdrag2 * u(z=1) * |V| For this option, as u(z=1) changes elevation, the roughness (rdrag2) stays the same. Another option is to use UV_LOGDRAG b...
by jcwarner
Fri Jul 13, 2018 12:43 pm
Forum: ROMS Tools and Techniques
Topic: wind vector direction in NAM and NARR data
Replies: 5
Views: 3030

Re: wind vector direction in NAM and NARR data

possibly. they all get data from here: https://www.ncdc.noaa.gov/nomads/data-products but there had been some changes at the noaa site. I remember they add the 's' to https and that was a problem in matlab. Matblab had been comparing the first 5 letters of a file, and if it started with http:, then ...
by jcwarner
Fri Jul 13, 2018 12:28 pm
Forum: ROMS Discussion
Topic: Blowing up problem using GLS mixing
Replies: 8
Views: 1492

Re: Blowing up problem using GLS mixing

this is a situation where you will need to dig in a little deeper. You say with n2s2 avg turned on, it blows up near the beginning. Suggest you see how many steps it takes before it blows up, lets say it is a 100. then save the model output in a his file every 5 times (or so). that will give you 20 ...
by jcwarner
Fri Jul 06, 2018 11:59 am
Forum: ROMS Discussion
Topic: on the wave mixing
Replies: 4
Views: 909

Re: on the wave mixing

i believe you are understanding correctly.
-j
by jcwarner
Thu Jul 05, 2018 2:54 pm
Forum: ROMS Discussion
Topic: on the wave mixing
Replies: 4
Views: 909

Re: on the wave mixing

These features are best described in Kumar, N., Voulgaris, G., Warner, J.C., and M., Olabarrieta (2012). Implementation of a vortex force formalism in a coupled modeling system for inner-shelf and surf-zone applications. Ocean Modelling, 47, 65-95. and Carniel, S., Warner, J.C., Chiggiato, J., and S...
by jcwarner
Thu Jul 05, 2018 12:39 pm
Forum: ROMS Discussion
Topic: Problems with SWAN in the COAWST system
Replies: 6
Views: 1439

Re: Problems with SWAN in the COAWST system

do me a favor and lets post this issue here:
https://coawstmodel-trac.sourcerepo.com ... /newticket
by jcwarner
Thu Jun 14, 2018 1:34 pm
Forum: ROMS Discussion
Topic: Problems with SWAN in the COAWST system
Replies: 6
Views: 1439

Re: Problems with SWAN in the COAWST system

ok.. 2 things here. 1) Comment the Boundary lines, and then you get the 'recordaxe' error. Comment out the writing of netcdf files for swan. I have a problem when both the parent and child are writing out netcdf files. I need to correct this but just did not get to it yet. You can write out mat file...
by jcwarner
Wed Jun 13, 2018 1:07 pm
Forum: ROMS Discussion
Topic: error in running COAWST
Replies: 1
Views: 506

Re: error in running COAWST

can you do me a favor and post this as a ticket here:
https://coawstmodel-trac.sourcerepo.com ... WST/report

-john
by jcwarner
Tue Jun 12, 2018 3:25 pm
Forum: ROMS Discussion
Topic: Problems with SWAN in the COAWST system
Replies: 6
Views: 1439

Re: Problems with SWAN in the COAWST system

this looks like 2 different errors. 1) Error : invalid boundary point Error : ( -99 -99) is outside computational grid so one of your boundary points (maybe you are using tpar files?) is not in the grid. SWAN thinks the point is at -99 -99, which is south of the south pole. So maybe just comment out...
by jcwarner
Wed May 23, 2018 5:31 pm
Forum: ROMS Discussion
Topic: Turn off tracer advection
Replies: 6
Views: 742

Re: Turn off tracer advection

how about
#define TS_FIXED

-j
by jcwarner
Tue May 22, 2018 11:58 am
Forum: ROMS Discussion
Topic: My constant wind-induced case seems strange
Replies: 4
Views: 891

Re: My constant wind-induced case seems strange

THETA_S == 8.0d0 ! surface stretching parameter THETA_B == 0.5d0 ! bottom stretching parameter TCLINE == 20.0d0 ! critical depth (m) use something like Gout=get_roms_grid('ocean_his.nc',tidx); this will give you a field Gout.s_rho that you can avg to get the vertical locations of your grid. how did ...
by jcwarner
Mon May 21, 2018 2:45 pm
Forum: ROMS Tools and Techniques
Topic: error as reading my atmospheric forcing files
Replies: 2
Views: 722

Re: error as reading my atmospheric forcing files

I think you have the 'coordinates' info backwards. here is some information: ROMS/Utility/get_varcoords.F ! This routine reads the spatial positions of any NetCDF variable ! ! having the attribute "coordinates", as specified by CF rules. ! ! For example, in CDL syntax: ! ! ! ! float my_var(time, lat...
by jcwarner
Thu May 10, 2018 4:46 pm
Forum: ROMS Problems
Topic: Coastal blow-ups issues: WRF typhoon to ROMS [via COAWST]
Replies: 5
Views: 1390

Re: Coastal blow-ups issues: WRF typhoon to ROMS [via COAWST

so this sounds like a roms coastline problem (not necessarily coawst). To be sure, what do the su/vstr look like? How about the surface heat fluxes?

-john
by jcwarner
Sun May 06, 2018 5:16 pm
Forum: ROMS Bugs
Topic: sustr and svstr Value has not been read
Replies: 11
Views: 1723

Re: sustr and svstr Value has not been read

Gridded: LonMin = -13.6075 LonMax = 35.9397 LatMin = 97.6335 LatMax = 146.0208 Model: LonMin = 0.0000 LonMax = 0.0000 LatMin = 0.0000 LatMax = 0.0000 your forcing file does not have any grid data. each field needs an attribute of lon and lat. see here https://www.myroms.org/forum/viewtopic.php?f=1&t...
by jcwarner
Sun May 06, 2018 3:54 pm
Forum: ROMS Bugs
Topic: sustr and svstr Value has not been read
Replies: 11
Views: 1723

Re: sustr and svstr Value has not been read

you basically have this set: /* Forcing */ # ifdef WRF_MODEL ... # else /*# define BULK_FLUXES/* # define ANA_SMFLUX # endif so you have ana_smflux on. that will look for analytical values and they are most likely set to 0. so you would want to use #undef ANA_SMFLUX and then it will look for su/vstr...
by jcwarner
Wed Apr 25, 2018 5:33 pm
Forum: ROMS Problems
Topic: Nesting Run Problem
Replies: 8
Views: 3766

Re: Nesting Run Problem

no. you need N= the same number for each nested grid.
by jcwarner
Sun Apr 22, 2018 10:59 pm
Forum: ROMS Discussion
Topic: Coupling any Atmospheric model to ROMS
Replies: 3
Views: 962

Re: Coupling any Atmospheric model to ROMS

we have coupled roms to wrf and several wave models using mct as part of the COAWST modeling system. You can use that as an example, but my experiences tell me that every model is different and there is no set specific way that will work for all models. so you would need to learn that other model co...
by jcwarner
Sun Apr 22, 2018 10:54 pm
Forum: ROMS Discussion
Topic: WRF quilting support in the WRF-ROMS setup in COAWST
Replies: 1
Views: 687

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

when the wrf + other-model job is submitted, the total number of cpus goes through an mpi_comm_split to allocate the procs for wrf. Then each wrf proc does an MCT init so all the processors for wrf will look to communicate through mct. I dont think the current configuration will allow quilting, but ...
by jcwarner
Mon Mar 26, 2018 11:49 am
Forum: ROMS Adjoint
Topic: question about the caculating of om_p, om_u, on_p, on_u
Replies: 5
Views: 1319

Re: question about the caculating of om_p, om_u, on_p, on_u

the metrics are computed in a way that is consistent in the code, and needs to be this way for mass conservation properties. Here is a post from 2003 on the similar issue:
viewtopic.php?f=19&t=150&view=next
by jcwarner
Mon Feb 26, 2018 2:58 pm
Forum: ROMS Problems
Topic: excess mixing in k-epsilon GLS mixing; bug or misconfigured?
Replies: 8
Views: 2055

Re: excess mixing in k-epsilon GLS mixing; bug or misconfigu

i have been out of touch for a awhile with some family losses but trying to get back to work. When I coded GLS into ROMS in 2003 (or whenever) I spent a considerable effort to compare my code to GOTM. I do remember some small differences to GOTM, but essentially our implementation gave the same resu...
by jcwarner
Mon Feb 05, 2018 2:18 pm
Forum: ROMS Discussion
Topic: Problem compiling using new version
Replies: 6
Views: 1163

Re: Problem compiling using new version

edit ROMS/Utility/def_quick.F and look near lines 1576-1577. IF you see this: IF (FoundError(exit_flag, NoError, __LINE__, & & cccc __FILE__)) RETURN change it to IF (FoundError(exit_flag, NoError, __LINE__, & & __FILE__)) RETURN see if that works. -john
by jcwarner
Fri Jan 26, 2018 10:26 pm
Forum: ROMS Discussion
Topic: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles
Replies: 8
Views: 1219

Re: Curvilinear grid, (Uwind,Vwind) and (sustr,svstr) angles

do you have curvgrid defined? i think the ana_wind is doing the right thing because it looks like your grid is curved. So a 'north' wind on a curved grid is in the eta direction. The set_data routine does not seem to perform any rotation for calls to either ana_wind or ana_smflux. So i think the iss...
by jcwarner
Thu Jan 11, 2018 2:51 am
Forum: ROMS Problems
Topic: problem of running coupled SWAN+ROMS, Error of Inlet_test
Replies: 1
Views: 1128

Re: problem of running coupled SWAN+ROMS, Error of Inlet_tes

typically we dont modify the makefile. make changes to the build.bash or coawst.bash file, and run that bash file.
-john
by jcwarner
Wed Jan 10, 2018 12:41 pm
Forum: ROMS Sediment
Topic: Blowing-up when using COAWST to couple ROMS and SWAN
Replies: 4
Views: 2172

Re: Blowing-up when using COAWST to couple ROMS and SWAN

the water levels go from roms to swan, so wet dry will affect swan. when you increase the min depth for swan to 1 m, then there will be locations that will turn dry and this will change the wave dynamics in the very shallow parts. suggest you do 2 runs with the different depths, and then pcolor some...
by jcwarner
Fri Dec 29, 2017 4:29 pm
Forum: ROMS Sediment
Topic: Blowing-up when using COAWST to couple ROMS and SWAN
Replies: 4
Views: 2172

Re: Blowing-up when using COAWST to couple ROMS and SWAN

Your stdout says " 17974 62 09:50:00 9.874450E-02 2.791750E+02 2.792738E+02 2.775544E+13 (095,245,20) 0.000000E+00 1.108514E-02 8.481048E+00 6.956400E+01 DIAG speed trouble 69.5640000538570 Blowing-up: Saving latest model state into RESTART file DIAG speed ijk 98.3770376595520 95 246 1 1" So you nee...
by jcwarner
Fri Dec 22, 2017 1:22 pm
Forum: ROMS Discussion
Topic: Analytical OBC Files?
Replies: 2
Views: 1248

Re: Analytical OBC Files?

Sometimes you just need to run the model to see what the impact of the boundary conditions are. if you created a domain with a boundary that is "purposefully a large distance away from my (nearshore) area of interest" that sounds like a good thing. The values you impose are not necessarily erroneous...
by jcwarner
Mon Dec 11, 2017 1:18 pm
Forum: ROMS Problems
Topic: problem of running coupled SWAN+ROMS under wind forcing
Replies: 1
Views: 874

Re: problem of running coupled SWAN+ROMS under wind forcing

you have 2 wave-current formulations defined:
WEC_MELLOR Mellore wave current interaction terms.
WEC_VF Vortex Force wave current interaction.
you need to select only one. I suggest you use WEC_VF.
-john
by jcwarner
Fri Dec 08, 2017 11:09 pm
Forum: ROMS Problems
Topic: CPU time and MPI issues with nested grids
Replies: 16
Views: 5187

Re: CPU time and MPI issues with nested grids

the child grid should be listed with the -2 in each direction as well.
try to compile in debug=on and run oceanG to see if you get more info.
-j
by jcwarner
Fri Dec 08, 2017 6:37 pm
Forum: ROMS Discussion
Topic: input forcing files
Replies: 4
Views: 934

Re: input forcing files

check the values of Tcline in your ocean.in and the media/Ahumada2016/HyCom_Tehua_Scar/GT_inic.nc file.
these all need to be the same.
by jcwarner
Fri Dec 08, 2017 1:50 am
Forum: ROMS Discussion
Topic: input forcing files
Replies: 4
Views: 934

Re: input forcing files

try this NFFILES == 3 ! number of unique forcing files FRCNAME == /media/Ahumada2016/HyCom_Tehua_Scar/GT_tidfrc.nc \ /media/Ahumada2016/HyCom_Tehua_Scar/GT_genfrc.nc \ /media/Ahumada2016/HyCom_Tehua_Scar/GT_wind2003.nc depending on the version of roms, the newer version wants the tidal forc to be li...
by jcwarner
Mon Nov 27, 2017 6:27 pm
Forum: ROMS Discussion
Topic: unrealistic pressure gradient forming along bathymetry
Replies: 10
Views: 3196

Re: unrealistic pressure gradient forming along bathymetry

this can occur for many reasons. One could be an imbalance of the initial density-velocity-bathymetry. Where did you get the density and velocity fields from? did you interpolate them to the roms grid? If you let it run for a while, do they go away?
-john
by jcwarner
Wed Nov 15, 2017 8:52 pm
Forum: Job Opportunities
Topic: USGS Woods Hole numerical modeler position
Replies: 0
Views: 1336

USGS Woods Hole numerical modeler position

The US Geological Survey, Woods Hole Coastal and Marine Science Center is seeking to fill a term position to conduct research pertaining to impacts of storms on coastal systems. The Researcher will help to develop and apply all aspects of a coupled hydrodynamic and wave modeling system, including gr...
by jcwarner
Wed Oct 11, 2017 11:37 am
Forum: ROMS Problems
Topic: floating invalid Couple ROMS and SWAN
Replies: 2
Views: 1166

Re: floating invalid Couple ROMS and SWAN

difficult to tell. maybe start out by looking at line 634 ocean_coupler.f90
do you have wettt/dry on? is this a shallow water application?

can you post a ticket here:
https://coawstmodel-trac.sourcerepo.com ... el_COAWST/

thanks
john
by jcwarner
Tue Oct 03, 2017 11:42 am
Forum: ROMS Sediment
Topic: Aout sediment.in
Replies: 1
Views: 1079

Re: Aout sediment.in

it might be the bottom roughness is using the grain size and that is creating a smaller stress.
we have some test cases like lake_signell or test_head that might be of help to get things set up.
-john
by jcwarner
Wed Sep 20, 2017 12:39 pm
Forum: ROMS Discussion
Topic: Error of Inlet_test
Replies: 1
Views: 792

Re: Error of Inlet_test

it is telling you "Number assigned processors: 002 not equal to spawned MPI nodes: 001" So for the inlet test, it is setup to use 1 proc for roms, and 1 for swan, so that means you need to submit the job with mpirun -np 2. Of course you can change the number of proces for each model, we just have th...
by jcwarner
Mon Sep 11, 2017 7:46 pm
Forum: ROMS Installation
Topic: gfortran -ffast-math experience
Replies: 2
Views: 1360

gfortran -ffast-math experience

Dear all- I want to strongly encourage removal of gfortran compiler flag: -ffast-math My recent experience shows this compiler option is not a good choice. The documentation is stating: [https://gcc.gnu.org/wiki/FloatingPointMath] " ... the -ffast-math flag which is a shortcut for several options, p...
by jcwarner
Mon Sep 11, 2017 2:21 pm
Forum: ROMS Installation
Topic: compling COAWST error
Replies: 1
Views: 890

Re: compling COAWST error

a few other people had similar issues. they were using older version of gfortran that did not allow a value that was a single integer to be allocatable. we suggested to them to upgrade their gfortran and have not herd back so that most likely worked for them. can you try the same? the posts were at ...
by jcwarner
Fri Sep 08, 2017 5:46 pm
Forum: ROMS Discussion
Topic: ROMS-WRF coupling problem
Replies: 3
Views: 1116

Re: ROMS-WRF coupling problem

looks odd. if they are same grid, make sure you dont define MCT_INTERP_OC2AT.
did you try any of the cases distributed with the code?
-john
by jcwarner
Wed Sep 06, 2017 4:29 pm
Forum: ROMS Problems
Topic: Error in sed_test1
Replies: 2
Views: 1154

Re: Error in sed_test1

made a few slight changes to ass SOURCES(ng)% in front of the Qbar, and modified a reshape call.
see if this works.
if so then you or i can submit as a ticket to be fixed.
thanks,
john
by jcwarner
Fri Aug 25, 2017 12:20 pm
Forum: ROMS Discussion
Topic: How do you access GRID variables that aren't called in tile?
Replies: 1
Views: 579

Re: How do you access GRID variables that aren't called in t

ana_vmix.h has an example for z_r and z_w. Notice that z_w is slightly different as it has a range from 0:N, as compared to z_r that goes from 1:N.

-j
by jcwarner
Wed Aug 16, 2017 6:01 pm
Forum: ROMS Discussion
Topic: How to hot-start swan?
Replies: 5
Views: 2746

Re: How to hot-start swan?

there is information in the coawst user manual, page 19. "5.3.7) SWAN restart files. SWAN has the command 'HOTFILE' that is used to write out the model state of two-dimensional spectra to a file at the end of the simulation. This file can be used as a restart for a future simulation. We have recentl...
by jcwarner
Mon Jul 24, 2017 5:05 pm
Forum: ROMS Bugs
Topic: Difference between trunk/ROMS and COAWST/ROMS with Qout
Replies: 3
Views: 2045

Re: Difference between trunk/ROMS and COAWST/ROMS with Qout

i just made a change in mod_ncparam to add:
< CASE ('idsurT(itemp)')
< idsurT(itemp)=varid
< CASE ('idsurT(isalt)')
< idsurT(isalt)=varid

Ivicia found this out. can you do an update and see if that fixed it?

-john
by jcwarner
Sun Jul 02, 2017 11:57 pm
Forum: ROMS Installation
Topic: COAWST compiling Problem
Replies: 5
Views: 1557

Re: COAWST compiling Problem

In ROMS/Utility/mp_exchange.F the only occurrence of mpi_proc_null is on lines 113-117 as: # if defined MPI Null_Value=MPI_PROC_NULL # else Null_Value=-1 # endif so if you defined MPI to compile roms only, and defined MPI to compile roms+swan, then this section of code Null_Value=MPI_PROC_NULL shoul...
by jcwarner
Sun Jul 02, 2017 1:26 pm
Forum: ROMS Installation
Topic: COAWST compiling Problem
Replies: 5
Views: 1557

Re: COAWST compiling Problem

the problem you are showing is roms code compilation. the coupling would have nothing to do with this. So if you can get roms to compile by itself, then use those same mpi settings to compile the coupled system.
by jcwarner
Sun Jul 02, 2017 11:56 am
Forum: ROMS Installation
Topic: COAWST compiling Problem
Replies: 5
Views: 1557

Re: COAWST compiling Problem

this is an issue during the build of roms. looks like an mpi could not find a correct lib or something. Can you build an application with just roms? this might be an openmpi problem.
by jcwarner
Sun Jun 25, 2017 3:20 pm
Forum: ROMS Discussion
Topic: Blended SST combined data in COAWST
Replies: 1
Views: 638

Re: Blended SST combined data in COAWST

WRF first uses SST data from its wrflow files. then it gathers sst from all roms grids and the scrip weights are used to determine the final sum of all the sources. all this should be taken care of. all you need to do is the scrip weights.
by jcwarner
Mon Jun 19, 2017 11:44 am
Forum: ROMS Sediment
Topic: A problem about coupled SWAN and ROMS model
Replies: 4
Views: 2191

Re: A problem about coupled SWAN and ROMS model

For both choices (uncoupled or coupled) you still have 2 choices: you can save Ub from swan and force roms with that or use the SSW_CALC_UB option and have the code compute Ub. For the bottom roughness, we have a minimum value of 0.05 set in MAster/mct_roms_swan.h. you can edit that file and change ...
by jcwarner
Fri Jun 16, 2017 2:57 pm
Forum: ROMS Sediment
Topic: A problem about coupled SWAN and ROMS model
Replies: 4
Views: 2191

Re: A problem about coupled SWAN and ROMS model

typically when you couple roms with swan, then you should use the bottom orbital velocity from swan. so we typically only use ssw_calc_znot with roms+swan. If you added ssw_calc_ub, then Ubot is being computed in sswbbl as: ! Set bed wave orbital velocity and excursion amplitude. Use data ! from wav...
by jcwarner
Fri Apr 28, 2017 6:46 pm
Forum: ROMS Problems
Topic: cannot read swan .in file, newest coawst code
Replies: 1
Views: 972

Re: cannot read swan .in file, newest coawst code

it says
"formatted read / unit=0 / READ not allowed for write-only file. File name = stderr "
so the file named "stderr" is not being allowed open for write.
That file is created by your operating system by PBS or Slurm or whatever. try to "rm stderr" and see if it works after that.
-j
by jcwarner
Fri Apr 28, 2017 12:23 am
Forum: ROMS Discussion
Topic: Weird Inf in the first time step, coupled with ROMS and SWAN
Replies: 3
Views: 704

Re: Weird Inf in the first time step, coupled with ROMS and

can you turn off any cpp options,
like SSWBBL, WEC, Sediment, etc.
by jcwarner
Tue Apr 25, 2017 1:37 pm
Forum: ROMS Discussion
Topic: Weird Inf in the first time step, coupled with ROMS and SWAN
Replies: 3
Views: 704

Re: Weird Inf in the first time step, coupled with ROMS and

there could be a lot of different things going on here. We need more info to provide some guidance. Can you run the grids with just swan and then just roms successfully? then try to couple them and see what happens. Are the Inf values at the wet/dry interface, along an open boundary, ... If you turn...
by jcwarner
Mon Apr 17, 2017 11:59 am
Forum: ROMS Problems
Topic: The model does not read the Forcing file for the Child grid
Replies: 6
Views: 2266

Re: The model does not read the Forcing file for the Child g

not sure where you are getting the forcing data from, but when you made the forcing files, did you interpolate the raw data to the grid resolutions? If you keep the raw data at the same resolution it came from, then the model will interpolate to the parent and child resolutions. Can you just make a ...
by jcwarner
Thu Apr 13, 2017 5:34 pm
Forum: ROMS Problems
Topic: The model does not read the Forcing file for the Child grid
Replies: 6
Views: 2266

Re: The model does not read the Forcing file for the Child g

how do you know that "[the model] interpolates data from the parent forcing file into the child grid."?
I would need to see the std out text to see what is happening.

-john
by jcwarner
Fri Mar 31, 2017 11:37 pm
Forum: ROMS Sediment
Topic: Issues with settling fluxes
Replies: 6
Views: 2382

Re: Issues with settling fluxes

sorry for taking a long time to respond. I received your email and was able to run the test case. as you said it was basically a shorter version of the distributed case. you can get different solutions depending on how many sed classes you choose, but the effect you are talking about is still there....
by jcwarner
Tue Mar 28, 2017 2:17 am
Forum: ROMS Usage
Topic: DOGBONE in COAWST
Replies: 1
Views: 1019

Re: DOGBONE in COAWST

looks like the nesting is limited to 3d (?). i thought it had worked in 2d at one point, but i have not tried a 2d nesting lately. see:

viewtopic.php?f=17&t=4475

try to
#define SOLVE3D
and see if it works with that.
-john
by jcwarner
Sun Mar 26, 2017 3:55 pm
Forum: ROMS Sediment
Topic: Issues with settling fluxes
Replies: 6
Views: 2382

Re: Issues with settling fluxes

so the sedbed_toy might have several variations over the years. the one i have ramps the forcing up and down to look at several cycles of resuspension and deposition. So to be sure that we are using the same cases, can you send the .h ocean.in and any ana files to me? jcwarner@usgs.gov
thanks
john
by jcwarner
Thu Mar 23, 2017 12:59 pm
Forum: ROMS Sediment
Topic: Issues with settling fluxes
Replies: 6
Views: 2382

Re: Issues with settling fluxes

I haven't looked at this in a while, but i can dig back into it. I suppose I am not clear what you plotted in fig 1. You said that each line is a different wsed, but then what are the different panels? One thing to check is that the bottom stress gets to be the same, especially as you change dt. You...
by jcwarner
Thu Mar 16, 2017 2:21 pm
Forum: ROMS Bugs
Topic: WET-DRY points initialization, how they go from DRY to WET?
Replies: 7
Views: 3235

Re: WET-DRY points initialization, how they go from DRY to W

can you put that rmask_wet in pre_step3d.F that i showed above?
-j
by jcwarner
Wed Mar 15, 2017 5:56 pm
Forum: ROMS Bugs
Topic: WET-DRY points initialization, how they go from DRY to WET?
Replies: 7
Views: 3235

Re: WET-DRY points initialization, how they go from DRY to W

we have also spent an effort to resolve similar issues. Are you using bulk_fluxes or imposing your own stresses? In bulk_fluxes we have the wet/dry masks multiplied so that the dry cells are not heating up. Also, this needs to be set in pre_step for the light penetration. we added: ! Add in incoming...
by jcwarner
Tue Mar 14, 2017 12:14 am
Forum: ROMS Bugs
Topic: Sediment bedload output issue
Replies: 2
Views: 1450

Re: Sediment bedload output issue

it is good to see that you can figure out the code, but i really dont see why those BC's would be needed. The bedldu is at u points, and the bedldv is at v points. So the call DO j=JstrR,JendR DO i=Istr,IendR bedldu(i,j,ised)=FX(i,j)*(pn(i-1,j)+pn(i,j))*cff would do all the u points. there is no nee...
by jcwarner
Thu Mar 02, 2017 1:46 am
Forum: ROMS Sediment
Topic: Increasing sediment volume in domain - a known issue?
Replies: 4
Views: 2351

Re: Increasing sediment volume in domain - a known issue?

if the domain is N-S periodic, are the results alongshore uniform? how long does this take to run? is it something i can grab and look at locally.
-j
by jcwarner
Wed Mar 01, 2017 4:20 am
Forum: ROMS Sediment
Topic: Increasing sediment volume in domain - a known issue?
Replies: 4
Views: 2351

Re: Increasing sediment volume in domain - a known issue?

looking at this quickly, i dont see a reason why the sediment should increase. there is no 'known issue' so it should conserve mass. where in the domain is the cross section?
by jcwarner
Fri Jan 27, 2017 7:25 pm
Forum: ROMS Sediment
Topic: Activate the SSW_BBL but the current magnitudes too small!
Replies: 1
Views: 1420

Re: Activate the SSW_BBL but the current magnitudes too smal

you need to look at the bottom stresses from the waves and currents separately. How big are the waves? what sediment characteristics did you set? IF the sd50 is too large, then the stresses will be very large and the currents small. On another note - did you post a question about 1/tau_ce or somethi...
by jcwarner
Tue Jan 24, 2017 3:04 am
Forum: ROMS Sediment
Topic: Bed load transport all zeros!
Replies: 1
Views: 1155

Re: Bed load transport all zeros!

the bedload does not move cohesive sed (NCS). It only moves the non cohesive (NNS).
-j
by jcwarner
Fri Dec 30, 2016 3:57 pm
Forum: Job Opportunities
Topic: USGS Woods Hole, MA, Numerical Modeler
Replies: 0
Views: 1902

USGS Woods Hole, MA, Numerical Modeler

The U.S. Geological Survey Woods Hole Coastal and Marine Science Center is looking for a numerical modeler to support nearshore and coastal processes studies. The candidate will work with the COAWST modeling system (WRF, ROMS, SWAN, InWave, and Sediment Transport Routines) to investigate storm dynam...
by jcwarner
Fri Nov 11, 2016 2:16 pm
Forum: ROMS Installation
Topic: Data exchange when coupling swan
Replies: 2
Views: 1506

Re: Data exchange when coupling swan

if you are using roms+swan coupling, you really should use our coawst system. wee have updated the swan version, implemented grid refinement in swan, have vortex force, etc. we have just been updating the coupling in the coawst system, and staying very close to the roms distribution. send me an emai...
by jcwarner
Mon Nov 07, 2016 1:23 pm
Forum: ROMS Installation
Topic: Netcdf downloading via Cygwin
Replies: 64
Views: 288371

Re: Netcdf downloading via Cygwin

I see you went through a lot of effort to build netcdf. so i dont want to send you astray, but I wanted ot make sure you knew that you can download pre-built netcdf in the cygwin environment, as well as the gcc/gfortran compiler, as well as the openmpi utilities so that you can just build roms and r...
by jcwarner
Fri Oct 21, 2016 12:20 pm
Forum: ROMS Installation
Topic: output files of SWAN model
Replies: 5
Views: 2547

Re: output files of SWAN model

to the Compilers/Linux-ifort.mk, you will need to add ifdef USE_SWAN FFLAGS += -assume byterecl endif i would suggest that if you intend on doing a lot of research with swan+roms, that you switch to coawst. We have updated swan to more recent versions, include vortex force wave effects on currents, ...
by jcwarner
Thu Oct 20, 2016 8:24 pm
Forum: ROMS Installation
Topic: output files of SWAN model
Replies: 5
Views: 2547

Re: output files of SWAN model

are you using Rutgers version of swan coupling or COAWST? In the past, this error was due to a byte offset needed for matlab, and it depending on the compiler. So we have added some flags , such as in Compilers/CYGWIN-ifort.mk ifdef USE_SWAN FFLAGS += /noextend_source -assume:byterecl endif what com...
by jcwarner
Thu Oct 20, 2016 6:55 pm
Forum: ROMS Discussion
Topic: Very cold water forming at boundaries
Replies: 10
Views: 3110

Re: Very cold water forming at boundaries

is the cold water forming at the boundaries of the model grid, or the land/sea boundary?

what are the Tnudg values at those locations? The climatology temp values may be = 0 because you have a finer grid than the clim data (just a guess).
by jcwarner
Thu Oct 20, 2016 12:19 pm
Forum: ROMS Sediment
Topic: Sediment from river!
Replies: 19
Views: 14123

Re: Sediment from river!

thanks for sending the whole thing. But you did not understand my question, or in fact, i didn't ask my question in a good way. I asked "where is the value of ised set" I see ised is declared as an integer, but what value is set for the variable "ised". not the value of Tsrc, but ised itself. You ne...
by jcwarner
Tue Oct 18, 2016 7:27 pm
Forum: ROMS Installation
Topic: Finished Compiling COAWST - MAC OSX Leopard 10.6.8
Replies: 5
Views: 4556

Re: Finished Compiling COAWST - MAC OSX Leopard 10.6.8

to build MCT, you need to do a ./configure in the MCT root dir. This will create a Makefile.conf, you should look in that to see if the settings are correct. Then do a make.
-j
by jcwarner
Mon Oct 17, 2016 12:17 pm
Forum: ROMS Sediment
Topic: Sediment from river!
Replies: 19
Views: 14123

Re: Sediment from river!

well, i would need to see the whole routine. but a first look
"SOURCES(ng)%Tsrc(is,k,ised)=20.0_r8"
where is the value of ised set? looks like it needs a do ised=1,NST loop around it.
by jcwarner
Sat Oct 15, 2016 8:02 pm
Forum: ROMS Sediment
Topic: Sediment from river!
Replies: 19
Views: 14123

Re: Sediment from river!

the csed is just the init value of sediment in the water column. to introduce sediment as a psource, you need to set that in ana_psource.
by jcwarner
Fri Oct 14, 2016 2:15 pm
Forum: ROMS Sediment
Topic: Cannot introduce river sediment in ROMS (svn 713)
Replies: 3
Views: 2037

Re: Cannot introduce river sediment in ROMS (svn 713)

great. keep us posted if you any other issues.
by jcwarner
Fri Oct 14, 2016 1:22 pm
Forum: ROMS Sediment
Topic: Cannot introduce river sediment in ROMS (svn 713)
Replies: 3
Views: 2037

Re: Cannot introduce river sediment in ROMS (svn 713)

There could be a couple of issues, so lets try to figure it out. - Did you set NCS (or NNS) to 1 at top of ocean.in? - how are you adding sediment to the river> netcdf file? If you look at the stdout, does it show that it is trying to read the file and getting mud_01 values? - you could just make th...
by jcwarner
Mon Oct 03, 2016 1:36 am
Forum: ROMS Problems
Topic: Segmentation Fault when executing oceanM
Replies: 11
Views: 4449

Re: Segmentation Fault when executing oceanM

seeing that you are pointing coawstm directly to the coupling.in file, i will assume this issue is related to the paths of all the files. the coupling file will have paths to ocean.in and swan.in, and they will have paths to grids etc. make sure all the files are in the correct paths.
by jcwarner
Fri Sep 23, 2016 1:32 pm
Forum: ROMS Installation
Topic: issues about inlet_test case
Replies: 6
Views: 3100

Re: issues about inlet_test case

i dont see the error in that build. try ./coawst.bash >& foo or script build.log ./coawst.bash exit and that makes a build.log file. but the issue is taht the build shows (for example) cd /home/chenghao/inlet_test/Build; /usr/local/bin/mpif90 -c -heap-arrays -fp-model precise -ip -O3 -I/opt/intelsof...
by jcwarner
Fri Sep 23, 2016 12:43 pm
Forum: ROMS Installation
Topic: issues about inlet_test case
Replies: 6
Views: 3100

Re: issues about inlet_test case

can i see the part of the build output that caused this issue? so it would say
mpif90 -c - blah blah -I/home/chenghao/soft/include ..
and see if the build line has all the correct information.

do you have correct permissions on those mod files?
by jcwarner
Fri Sep 23, 2016 12:21 pm
Forum: ROMS Installation
Topic: issues about inlet_test case
Replies: 6
Views: 3100

Re: issues about inlet_test case

the first set are just warnings about write statements. these should be changed by swan. the second issues is an error you need to correct: "waves_coupler.f90(17): error #7002: Error in opening the compiled module file. Check INCLUDE paths. [M_MCTWORLD] USE m_MCTWorld, ONLY : MCTWorld_init => init" ...
by jcwarner
Wed Sep 21, 2016 5:32 pm
Forum: ROMS Tools and Techniques
Topic: wind vector direction in NAM and NARR data
Replies: 5
Views: 3030

Re: wind vector direction in NAM and NARR data

John- thanks for spotting this. Yes, you are correct. These velocities from NAM/NARR via THREDDS need to be rotated from Lambert Conformal to Earth long lat. I found a set of code in a fortran file called nam2model.F and some corresponding code at the bottom of this page: http://ruc.noaa.gov/RUC.faq...
by jcwarner
Sun Jul 31, 2016 3:27 pm
Forum: ROMS Discussion
Topic: the latest COAWST version compile problem
Replies: 1
Views: 3196

Re: the latest COAWST version compile problem

for coawst issues, can you post them on
https://coawstmodel-trac.sourcerepo.com ... discussion

for this issue, we chose to require netcdf 4 to red the scrip weight files. we updated the scrip to write out the weights using netcdf groups, and that requires netcdf 4.
-john
by jcwarner
Wed Jul 13, 2016 12:17 pm
Forum: ROMS Installation
Topic: Will it be convenient (in the long run) to use Windows?
Replies: 7
Views: 3162

Re: Will it be convenient (in the long run) to use Windows?

running roms in cygwin these days is actually a lot easier than it used to be. When you download cygwin, you can also get gfortran compiler (nice but slower than ifort and pgi) and you can also get the netcdf libraries all precompiled with that version of gfortran. you need to make sure you select o...
by jcwarner
Sat Jun 18, 2016 2:32 am
Forum: ROMS Bugs
Topic: Infinity=one/zero
Replies: 6
Views: 4818

Re: Infinity=one/zero

is this running one of the distributed cases?
by jcwarner
Fri Jun 17, 2016 11:45 am
Forum: ROMS Bugs
Topic: Infinity=one/zero
Replies: 6
Views: 4818

Re: Infinity=one/zero

you can remove this code:

#ifdef SWAN_COUPLING
Infinity=one/zero
#endif
by jcwarner
Fri Jun 17, 2016 12:53 am
Forum: ROMS Problems
Topic: ROMS&SWAN coupling, Segmentation Fault
Replies: 6
Views: 4419

Re: ROMS&SWAN coupling, Segmentation Fault

can you post this on
https://coawstmodel-trac.sourcerepo.com ... del_COAWST

i think i will need to see your .h file.
did you compile with nesting and not have 2 swan grids?
-j
by jcwarner
Fri Jun 03, 2016 3:34 pm
Forum: ROMS Discussion
Topic: Losing tracer as time progresses
Replies: 9
Views: 3157

Re: Losing tracer as time progresses

i posted a reply, and it did not seem to show up so i will try this again. We have conducted tests in the past with sediment to make sure the total amount of tracer is conserved, and the tests worked with MPDATA and other schemes. I suggest you make your calculations on the history files as the aver...
by jcwarner
Fri Jun 03, 2016 12:49 pm
Forum: ROMS Discussion
Topic: How to get zeta var without Hycom data for COAWST forecast
Replies: 5
Views: 2538

Re: How to get zeta var without Hycom data for COAWST foreca

yes, we have developed coawst as a separate entity from the roms code. please send me an email and we can give you access:
jcwarner@usgs.gov
thanks,
john
by jcwarner
Thu May 12, 2016 8:03 pm
Forum: ROMS Discussion
Topic: How to get zeta var without Hycom data for COAWST forecast
Replies: 5
Views: 2538

Re: How to get zeta var without Hycom data for COAWST foreca

So your question is specific to ROMS to get ssh data for a boundary file. I am not sure what i said, but it was probably a mention that HYCOM data does not have tides (?). So if you want sea surface height data with tides, there are other sources. Where is your domain?

-john
by jcwarner
Fri Apr 08, 2016 12:35 pm
Forum: ROMS Problems
Topic: How to control sediment resuspension in ROMS(version783)?
Replies: 3
Views: 2209

Re: How to control sediment resuspension in ROMS(version783)

still could be a few things. I saw in your .h file

#define SEDIMENT
#if def SEDIMENT
# define SUSPLOAD
# define SED_MORPH
# undef SED_DENS
#endif

that second line looks funny. make it be
#ifdef SEDIMENT

also, can you post your stdout?
what it is the bed thickness?
by jcwarner
Thu Apr 07, 2016 12:39 pm
Forum: ROMS Problems
Topic: How to control sediment resuspension in ROMS(version783)?
Replies: 3
Views: 2209

Re: How to control sediment resuspension in ROMS(version783)

I think you should add a bunch of variables to your history file and see what is happening. activate write out of bottom stresses, ssc, bed thickness, bed mass, etc etc. then you can look at computed bottom stress as compared to tau_crit erosion and bedmass etc to get a better picture of what is hap...
by jcwarner
Wed Apr 06, 2016 12:22 pm
Forum: ROMS Problems
Topic: Post subject: Error while processing Climatology multi-file
Replies: 10
Views: 4799

Re: Post subject: Error while processing Climatology multi-f

these m files may have been updated at some time. to solve the problem, it would simply be to load the ocean_time into maltab using ot=ncread(file,'ocean_time'); use ocean_time=ocean_time/3600/24 ( I think), then ncwrite(file,'ocean_time',ocean_time) and we can move on. these issues happen a lot bec...
by jcwarner
Tue Apr 05, 2016 11:46 pm
Forum: ROMS Problems
Topic: Post subject: Error while processing Climatology multi-file
Replies: 10
Views: 4799

Re: Post subject: Error while processing Climatology multi-f

the error says CHECK_MULTIFILE - Error while processing Climatology multi-files: data does not include initialization time = 56228 00:00:00 ***** **:**:** ***** **:**:** Projects_My/scs_1/3D_BC/6_days/scs_1_clm.nc so the model is reading your netcdf file and coming up with: ***** **:**:** ***** **:*...
by jcwarner
Thu Mar 10, 2016 8:09 pm
Forum: ROMS Discussion
Topic: angle in grid file
Replies: 8
Views: 3083

Re: angle in grid file

Kate- yes I agree. I think the more important part of my post was that ROMS uses the angle, not the units of it.
-j
by jcwarner
Thu Mar 10, 2016 2:21 pm
Forum: ROMS Discussion
Topic: angle in grid file
Replies: 8
Views: 3083

Re: angle in grid file

I don't want to cause any confusion on this thread, but the Rutgers version of ROMS does internally use the variable 'angle' at several instances, so the units do matter. For example, if a user has a forcing file of Uwind_east and Vwind_north, and they impose that on a rotated ocean grid, then the c...
by jcwarner
Thu Mar 03, 2016 1:45 pm
Forum: ROMS Problems
Topic: Another Problem with SWAN
Replies: 1
Views: 1755

Re: Another Problem with SWAN

I answered you on the coawst forum. no need to also ask here.
-john
by jcwarner
Fri Feb 26, 2016 6:26 pm
Forum: ROMS Tools and Techniques
Topic: scrip mask boundary
Replies: 2
Views: 2727

Re: scrip mask boundary

so two answers here: 1) for lbressan, we are in the process to release a new version of coawst soon. then you can test the new routines, and it would be good to know if the masking problem is corrected. this can be tricky if wrf and roms have different resolutions and land/sea masking at the coast. ...
by jcwarner
Wed Feb 17, 2016 3:45 pm
Forum: Meetings/Workshops
Topic: COAWST Training August 15-19, 2016
Replies: 0
Views: 4864

COAWST Training August 15-19, 2016

We are announcing a COAWST Training workshop to be held: August 15-19, 2016, in room 507 Clark Building, Woods Hole Oceanographic Institution, Quissett Campus, Woods Hole, MA. There will be a $75 registration fee to cover room charges, snacks, and cookout. If you would like to attend the Training, p...
by jcwarner
Sun Feb 07, 2016 10:48 pm
Forum: ROMS Discussion
Topic: treatment of rivers
Replies: 1
Views: 649

Re: treatment of rivers

there is a lot of information on this forum and on the wiki here:
https://www.myroms.org/wiki/River_Runoff

rivers do impart momentum and you can have them impose tracers.
-j
by jcwarner
Thu Feb 04, 2016 1:48 pm
Forum: ROMS Installation
Topic: build fail
Replies: 1
Views: 2148

Re: build fail

the EW and NS periodic calls have been updated a while ago. Perhaps you are using an older ana_* function with the new code, and it is complaining about the periodic exchange calls. do a diff with your ana functions to the new ones.

-j
by jcwarner
Wed Feb 03, 2016 1:50 pm
Forum: ROMS Usage
Topic: Headland Test Case
Replies: 4
Views: 2329

Re: Headland Test Case

if you want to use roms+swan coupling, then please use the formulations and coupling schemes that are in coawst. roms is distributed with an older version of swan, and it has a coupling infrastructure. but over the years we have updated swan, added wrf, updated the coupling scheme, added regridding ...
by jcwarner
Tue Feb 02, 2016 3:02 pm
Forum: ROMS Usage
Topic: Headland Test Case
Replies: 4
Views: 2329

Re: Headland Test Case

Tom-
if you are interested to use the roms+swan coupling, then I suggest you use the COAWST modeling system. we are using that system for advancements in the model coupling. please send an email to jcwarner@usgs.gov and ask for that code.
thanks,
john
by jcwarner
Mon Jan 18, 2016 7:05 pm
Forum: ROMS Discussion
Topic: COAWST sediment problem
Replies: 1
Views: 741

Re: COAWST sediment problem

the problem is most likely something to do with SED_DENS. So i did a grep -nir SED_DENS ROMS/* and the most important location for the forward model is in ROMS/Nonlinear/rho_eos.F: # if defined SEDIMENT && defined SED_DENS SedDen=0.0_r8 DO ised=1,NST itrc=idsed(ised) cff1=1.0_r8/Srho(ised,ng) SedDen...
by jcwarner
Sat Dec 19, 2015 11:12 pm
Forum: ROMS Discussion
Topic: Coupling test: swan output cannot combine?
Replies: 1
Views: 664

Re: Coupling test: swan output cannot combine?

if you are running coupled roms+swan, i suggest you use coawst. if you need access to coawst, please send me an email to jcwarner@usgs.gov. if you are using coawst and posting an issue, please post it on https://coawstmodel-trac.sourcerepo.com/coawstmodel_COAWST/ or just send me an email. in short f...
by jcwarner
Tue Dec 08, 2015 2:17 pm
Forum: ROMS Problems
Topic: Checkerboard instability
Replies: 7
Views: 3420

Re: Checkerboard instability

i cant download the attachment because my computer is being overly secure. but did you activate
#define N2S2_HORAVG

-j
by jcwarner
Tue Dec 01, 2015 2:09 pm
Forum: ROMS Usage
Topic: Model blowing up using elevation at the boundary
Replies: 12
Views: 4057

Re: Model blowing up using elevation at the boundary

our coawst forecast is for the US east coast. so I suggest you use the HYCOMM data as a start, and set the correct hycomm url in roms_master_climatology_coawst_mw.
-j
by jcwarner
Thu Nov 19, 2015 3:01 pm
Forum: ROMS Usage
Topic: INLET TEST: Why are initial wave directions 81 deg?
Replies: 2
Views: 1608

Re: INLET TEST: Why are initial wave directions 81 deg?

Amy- To run the inlet_test, we set the init conditions for SWAN to be Hwave=0. This case does not have any wind forcing, so there is really no 'real' init conditions for Dwave. The value of DWAVE=81 is computed by SWAN, but there are really no waves present at start so this is just some value that r...
by jcwarner
Tue Nov 10, 2015 2:20 pm
Forum: ROMS Problems
Topic: model blow-up due to abnormal velocity
Replies: 8
Views: 4901

Re: model blow-up due to abnormal velocity

very difficult to determine what is happening from that image.
Maybe you need to use a tide ramp:
#define RAMP_TIDES
and see if that helps.
-j
by jcwarner
Mon Oct 26, 2015 12:16 am
Forum: ROMS Discussion
Topic: unreasonable wave simulation result using COAWST
Replies: 2
Views: 792

Re: unreasonable wave simulation result using COAWST

since this is a wave related issue, I am not sure this is the best place for this discussion. can you submit something on the coawst forum https://coawstmodel-trac.sourcerepo.com/coawstmodel_COAWST/discussion/2 just to add some info, i mapped the swan wave height from our forecast (see attached imag...
by jcwarner
Fri Sep 18, 2015 12:58 pm
Forum: ROMS Discussion
Topic: unreasonable results in COAWST
Replies: 2
Views: 1071

Re: unreasonable results in COAWST

I think the last post would translate to something like: "Hello , I am also being used to simulate COAWST Western Pacific typhoon. And has the results , do some analysis. Perhaps we can discuss." If you have questions specific to ROMS, then this would be a great place to post those messages. For exa...
by jcwarner
Mon Sep 14, 2015 2:34 pm
Forum: ROMS Problems
Topic: ROMS-SWAN Coupling in COAWST
Replies: 2
Views: 1991

Re: ROMS-SWAN Coupling in COAWST

If you have an issue with COAWST, can you please post the message on the trac here:
https://coawstmodel-trac.sourcerepo.com ... discussion
or send me an email at jcwarner@usgs.gov
thanks,
john
by jcwarner
Mon Sep 07, 2015 8:21 pm
Forum: ROMS Problems
Topic: Problem with SWAN
Replies: 6
Views: 2825

Re: Problem with SWAN

well, the problem you posted is exactly the reason we put those extra flags in. the file ocpcre.F has USE OCPCOMM1 40.41 and that is the error you are getting. the '40.41' is beyond the 72 char limit, and so you need to tell the compiler $(SCRATCH_DIR)/ocpcre.o: FFLAGS += -nofree so you need to chec...
by jcwarner
Mon Sep 07, 2015 2:21 pm
Forum: ROMS Problems
Topic: Problem with SWAN
Replies: 6
Views: 2825

Re: Problem with SWAN

Are you using Rutgers roms with coupling or my COAWST coupled system? If you are coupling roms+swan, I suggest you use COAWST. This error is most likely due to the Compilers/Linux-ifort.mk (I am guessing here) does not have the correct flags for the special cases to compile swan. swan files are some...
by jcwarner
Tue Aug 18, 2015 4:19 pm
Forum: ROMS Installation
Topic: Finished compiling COAWST Installation - Ubuntu 14.04
Replies: 26
Views: 13441

Re: Finished compiling COAWST Installation - Ubuntu 14.04

Please don't post COAWST questions on the roms forum. I saw your post on the coawst trac but did not get a chace to answer it yet. In reality, this is an issue with your compilation of WRF. See the COAWST user manual : Set your environment variables, for example as shown here: setenv MCT_INCDIR /PAT...
by jcwarner
Mon Aug 17, 2015 8:13 pm
Forum: ROMS Discussion
Topic: mat2roms with modern matlab
Replies: 1
Views: 988

Re: mat2roms with modern matlab

We use the attached file: mat2roms_mw.m that uses Matlab inherent netcdf calls to create a roms grid file. An example usage would be: *************************** theWRFFile='wrfinput_d01'; theROMSFile='roms_grid.nc'; % open the wrf file rho.lat=ncread(theWRFFile,'XLAT'); rho.lon=ncread(theWRFFile,'X...
by jcwarner
Sat Aug 15, 2015 1:21 pm
Forum: ROMS Problems
Topic: Three nested grid (WRF-ROMS)
Replies: 1
Views: 1431

Re: Three nested grid (WRF-ROMS)

since this is a coawst issue, and not roms, lets post this in the coawst trac forum and i will respond thru that site.
https://coawstmodel-trac.sourcerepo.com ... el_COAWST/
login and post a new ticket.
thanks,
john
by jcwarner
Fri Jul 24, 2015 5:38 pm
Forum: ROMS Discussion
Topic: question follows after splitting bry and wind force file
Replies: 12
Views: 2867

Re: question follows after splitting bry and wind force file

does "data include initialization time = ***** 00:00:00"?
what are the actual values of ocean_time in your mysimulatin_bry.nc
by jcwarner
Thu Jul 23, 2015 3:19 pm
Forum: ROMS Discussion
Topic: question follows after splitting bry and wind force file
Replies: 12
Views: 2867

Re: question follows after splitting bry and wind force file

there was just a recent update for this. is your code up to date?
by jcwarner
Mon Jul 20, 2015 1:17 pm
Forum: ROMS Problems
Topic: forrtl: severe (174): SIGSEGV, segmentation fault occurred
Replies: 24
Views: 8175

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

it looks like it took 2 time steps 1 0 00:05:00 6.479416E-13 6.585677E+02 6.585677E+02 3.884376E+11 (01,01,01) 1.068554E-09 1.950202E-06 0.000000E+00 9.216771E-06 2 0 00:10:00 9.311375E-12 6.585677E+02 6.585677E+02 3.884376E+11 (01,02,08) 6.735275E-08 1.793582E-06 1.433645E-05 3.661905E-05 forrtl: s...
by jcwarner
Fri Jul 17, 2015 12:39 pm
Forum: ROMS Problems
Topic: Four Nested Grid
Replies: 1
Views: 1421

Re: Four Nested Grid

my experience with the nesting is that the child grids do not have to step on a time scale factor equal to the spatial scale factor. That is a recommendation, but not a requirement. I have not seen that error before, but you should be able to use dt of the child as the same as the parent. You could ...
by jcwarner
Wed Jul 15, 2015 6:28 pm
Forum: ROMS Problems
Topic: Problem with netcdf file for Dogbone example
Replies: 6
Views: 3200

Re: Problem with netcdf file for Dogbone example

can you try this file?
by jcwarner
Wed Jul 15, 2015 5:45 pm
Forum: ROMS Problems
Topic: Nesting Run Problem
Replies: 8
Views: 3766

Re: Nesting Run Problem

your input files has a lot of 'tabs' in it.
No tabs allowed.
please remove all tabs and try again. only use the space bar. not sure if there is a way to show all the tab characters in emacs or someother software.

-j
by jcwarner
Fri Jul 10, 2015 3:13 pm
Forum: ROMS Problems
Topic: Tide forcing model error
Replies: 10
Views: 2249

Re: Tide forcing model error

try to recompile with
#define ANA_FSOBC
#define ANA_M2OBC
by jcwarner
Fri Jul 10, 2015 1:51 pm
Forum: ROMS Problems
Topic: Tide forcing model error
Replies: 10
Views: 2249

Re: Tide forcing model error

what are the times in your lateral BC file?
CHECK_MULTIFILE - Error while processing Lateral Boundary multi-files:
data does not include initialization time = ***** 00:00:00
by jcwarner
Fri Jul 10, 2015 1:22 pm
Forum: ROMS Problems
Topic: Tide forcing model error
Replies: 10
Views: 2249

Re: Tide forcing model error

try this:
DSTART = 734323.0d0 ! days
TIDE_START = 734690.0d0 ! days
TIME_REF = -1.0d0 ! yyyymmdd.dd
by jcwarner
Fri Jul 10, 2015 12:35 pm
Forum: ROMS Problems
Topic: Tide forcing model error
Replies: 10
Views: 2249

Re: Tide forcing model error

what is your init time?
CHECK_MULTIFILE - Error while processing Lateral Boundary multi-files:
data does not include initialization time = ***** 00:00:00

compared to DSTART?
by jcwarner
Fri Jul 10, 2015 12:24 pm
Forum: ROMS Discussion
Topic: forrtl: severe (174): SIGSEGV, segmentation fault occurred
Replies: 3
Views: 1018

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

since these are issues with coawst and not rutgers roms, lets move this discussion.
can you send to me the read_coawst_par.f90
thanks
john
jcwarner@usgs.gov
by jcwarner
Fri Jul 10, 2015 12:22 pm
Forum: ROMS Problems
Topic: Tide forcing model error
Replies: 10
Views: 2249

Re: Tide forcing model error

should
NFFILES == 1 ! number of forcing file
FRCNAME == INPUT/mymodel_tide.nc |
INPUT/mymodel_frz.nc

be

NFFILES == 2
FRCNAME == INPUT/mymodel_tide.nc \
INPUT/mymodel_frz.nc
by jcwarner
Tue Jun 30, 2015 4:23 pm
Forum: ROMS Usage
Topic: Model blowing up using elevation at the boundary
Replies: 12
Views: 4057

Re: Model blowing up using elevation at the boundary

Warner, J.C., Haas, K.H ., Defne, Z., and Arango, H.G.(2013). A wetting and drying scheme for ROMS, Computers and Geosciences, 58, 54-61.
by jcwarner
Mon Jun 22, 2015 1:44 pm
Forum: ROMS Usage
Topic: Model blowing up using elevation at the boundary
Replies: 12
Views: 4057

Re: Model blowing up using elevation at the boundary

you can also use our COAWST forecast to get open bc's for roms. We run romw+swan coupled daily, and have output available on a thredds server model_url='http://geoport.whoi.edu/thredds/dodsC/coawst_4/use/fmrc/coawst_4_use_best.ncd' One way to access this is to download nctoolbox from here: https://g...
by jcwarner
Fri Jun 19, 2015 1:45 pm
Forum: ROMS Usage
Topic: Model blowing up using elevation at the boundary
Replies: 12
Views: 4057

Re: Model blowing up using elevation at the boundary

i am not sure that using free surface from HYCOM is the best choice. it really depends on your domain. i would think adcirc/osu tides would be a better place to start. for the error, it is difficult to tell. You could try to run the modeling system with just roms, and that would limit where to look....
by jcwarner
Wed Jun 17, 2015 12:56 pm
Forum: ROMS Discussion
Topic: Instabilities with SWAN_COUPLING and WET_DRY together
Replies: 8
Views: 3702

Re: Instabilities with SWAN_COUPLING and WET_DRY together

i try my best to keep a version of roms in coawst that is consistent with the rutgers web site. however, we do make modifications in coawst that may take some time to get back to rutgers. the same is for wrf and swan - we make changes that work for the coupled system and sometimes these features mak...
by jcwarner
Tue Jun 02, 2015 1:16 am
Forum: ROMS Discussion
Topic: Publications on M2 diagnostics
Replies: 2
Views: 902

Re: Publications on M2 diagnostics

Kumar, N., Voulgaris, G., Warner, J.C., and M., Olabarrieta (2012). Implementation of a vortex force formalism in the coupled ocean-atmosphere-wave-sediment transport (COAWST) modeling system for inner-shelf and surf-zone applications. Ocean Modeling 47, pp 65-95. Warner, List, Schwab, Voulgaris, Ar...