Search found 125 matches

by jivica
Sat Dec 07, 2019 1:21 am
Forum: ROMS Discussion
Topic: Parallel IO not scaling over Serial IO for ROMS
Replies: 13
Views: 173

Re: Parallel IO not scaling over Serial IO for ROMS

First of all you should read your log file (at the end) where you have profiling of your run. There you can see what portion of time is spent in IO and what in real computing kernel etc. In some cases when I need often IO I get bad ratio but this is to be expected. Remember Sasa once told he had ver...
by jivica
Fri Dec 06, 2019 5:00 am
Forum: ROMS Discussion
Topic: Parallel IO not scaling over Serial IO for ROMS
Replies: 13
Views: 173

Re: Parallel IO not scaling over Serial IO for ROMS

I'm using CRAY XC40 and must admit, quite happy with performance (dragonfly net topology), 24 core nodes. Made some benchmarking and it really depends on your application and grid size. Note that depending on users load on scratch you can have totally different results (and you can tweak parallel wr...
by jivica
Tue Dec 03, 2019 7:44 am
Forum: ROMS Wish List
Topic: Rotated 2-way nested grids ?
Replies: 0
Views: 55

Rotated 2-way nested grids ?

Hernan and the rest, are there any plans to make nesting scripts (contact.m) capable of ingesting rotated grids? Have feeling that all the bits of the engine are already available to make more general contacts between the grids. I am currently doing 1 way downscaling of arbitrary nested grids which ...
by jivica
Fri Jun 21, 2019 4:51 am
Forum: ROMS Adjoint
Topic: depth of assimilated data
Replies: 2
Views: 424

Re: depth of assimilated data

ROMS is not black box model, and it is very well documented. In that sense you should dig into the source and, for example, have a look at the Utility/extract_obs.F There you will find relevant part of the code: p2=REAL(i2-i1,r8)*(Xobs(iobs)-REAL(i1,r8)) q2=REAL(j2-j1,r8)*(Yobs(iobs)-REAL(j1,r8)) p1...
by jivica
Sat Jun 08, 2019 10:34 am
Forum: ROMS Discussion
Topic: bulk flux advice from clivar community
Replies: 3
Views: 872

Re: bulk flux advice from clivar community

Just to add on top of that, ERA-5 is amazing source of atmo forcing. From 1979 -> recent at 0.25 deg with HOURLY (!) fields. It was couled with ocean and used advanced data assimilation. Even for strong tropical cyclones it gets decent results, think it is better than JRA55 and definitely better the...
by jivica
Fri Mar 22, 2019 4:56 am
Forum: ROMS Discussion
Topic: Which tool is the stablest tools to create ROMS grid?
Replies: 20
Views: 2774

Re: Which tool is the stablest tools to create ROMS grid?

Just adding comment when Sasa is not going to ;) I find extremely useful Alex's tool for generating curvlinear orthogonality error free grids and suggest you to have a look at the older posts: https://www.myroms.org/forum/viewtopic.php?f=23&t=3878&p=14908 and https://www.myroms.org/forum/viewtopic.p...
by jivica
Sun Mar 10, 2019 4:07 am
Forum: ROMS Bugs
Topic: Pair and MPI version has problem with regridding
Replies: 10
Views: 908

Re: Pair and MPI version has problem with regridding

I've sent you the link privately on the email..

Ivica
by jivica
Sat Mar 09, 2019 12:09 pm
Forum: ROMS Bugs
Topic: Pair and MPI version has problem with regridding
Replies: 10
Views: 908

Re: Pair and MPI version has problem with regridding

Not sure if I am right, but think the problem is in new version of regrid and new variable "myxout" After compiling in debug mode (with mpi) I managed to trap it for 6x6 tile configuration (4x4 and 8x8 works?!): At line 155 of file regrid.f90 At line 155 of file regrid.f90 At line 155 of file regrid...
by jivica
Sat Mar 09, 2019 6:20 am
Forum: ROMS Bugs
Topic: Pair and MPI version has problem with regridding
Replies: 10
Views: 908

Re: Pair and MPI version has problem with regridding

It is getting even more interesting; it is working for certain tile (2x2, 3x2) configuration and crashing for 6x4, and then working for 32, 36, 48 but having wrong pressure as I wrote. I recompiled roms with debug option and gfortran + openmpi and bomb turns to be in inp_par.f90 which in line 77 is ...
by jivica
Sat Mar 09, 2019 5:09 am
Forum: ROMS Bugs
Topic: Pair and MPI version has problem with regridding
Replies: 10
Views: 908

Re: Pair and MPI version has problem with regridding

John, tried mpirun -np 1 and it is working OK, gives identical result as serial run. It is southern hemisphere system, NO BULK_FLUX, only storm surge case with surface wind stress and pressure. Reading of original data is OK as well (in all ceases, MPI or serial), reasonable values within range. Wil...
by jivica
Fri Mar 08, 2019 7:09 am
Forum: ROMS Bugs
Topic: Pair and MPI version has problem with regridding
Replies: 10
Views: 908

Re: Pair and MPI version has problem with regridding

Hernan, I know you are quite busy and that this bug is on your todo list, this post was more to have others aware of the problem in the latest version of ROMS code (if they use MPI and Pair as I am). Using in addition (to ATM_PRESS) option for bry pressure correction --> PRESS_COMPENSATE doesn't cha...
by jivica
Fri Mar 08, 2019 4:49 am
Forum: ROMS Bugs
Topic: Pair and MPI version has problem with regridding
Replies: 10
Views: 908

Pair and MPI version has problem with regridding

Having problem with application where I am using Pair on atmo model native grid and *ONLY* MPI-parallel version of the latest code. In serial version is regridding and working OK, but parallel version still has problem.. Trying to nail that, and by digging I think I've found bug in ./Modules/mod_for...
by jivica
Tue Feb 26, 2019 7:52 am
Forum: ROMS Bugs
Topic: Error compiling ROMS with new flag SINGLE_PRECISION (svn928)
Replies: 7
Views: 1349

Re: Error compiling ROMS with new flag SINGLE_PRECISION (svn

If you are using SINGLE_PRECISION and wet-dry option then you can't compile the latest version of ROMS. There is DSIGN problem (at least with gfortran), so after changing DSIGN -> SIGN for Build/u2dbc_im.f90 Build/v2dbc_im.f90 Build/wetdry.f90 Build/step2d.f90 and ./build_roms.bash -noclean you can ...
by jivica
Tue Jan 08, 2019 4:44 am
Forum: ROMS Discussion
Topic: Doubts regarding strong surface currents
Replies: 12
Views: 2580

Re: Doubts regarding strong surface currents

I suggest you first to read a few papers related to the wind stress formulations in extreme events like TC, many of them are John's papers using COAWST and many Hurricane cases etc. Not to just blindly run ROMS model and then expect to get result that is perfect. For example: Bryant & Akbar in JMSE ...
by jivica
Tue Nov 27, 2018 7:08 am
Forum: ROMS Discussion
Topic: Set a land mask on the sea in SeaGrid
Replies: 22
Views: 3969

Re: Set a land mask on the sea in SeaGrid

Dear Sasha, really like you answers, always getting to the bottom of the problem. Regarding mesh, I can see that many wants to nest inside global models (Hycom or NEMO) and do not want/need to download the whole grid domain just to get stripes along boundary. In other words, getting stripe along con...
by jivica
Wed Aug 08, 2018 4:42 am
Forum: ROMS Discussion
Topic: River bottom salinity too high
Replies: 9
Views: 1299

Re: River bottom salinity too high

Actually the best option would be to dig a uniform depth/width channel inland (a few pixels should be enough) and have mixing processes happening "over the land", so once entering the ocean it is as should be.

Cheers,
Ivica
by jivica
Mon Aug 06, 2018 4:50 am
Forum: ROMS Discussion
Topic: River bottom salinity too high
Replies: 9
Views: 1299

Re: River bottom salinity too high

Adding my 2c to the topic; I am running "many rivers scenario" without specifying temperature (only salinity) of the river sources and all is working fine (this option is using ambient ocean temperature for river sources as I don't know exactly temperature and it is shallow and really similar to the...
by jivica
Thu Jun 07, 2018 6:33 am
Forum: ROMS Discussion
Topic: Nonbreaking wave-induced mixing in ROMS?
Replies: 7
Views: 1532

Re: Nonbreaking wave-induced mixing in ROMS?

Thanks for all responses! Farhid in which version of ROMS you've implemented it? Looking in Rutger's version and lmd_vmix.F there is only diffusivity due to the internal waves and not because of surface wind induced waves? lmd_vmix.F ! Compute interior diffusivity due to wave breaking (Gargett and !...
by jivica
Thu May 31, 2018 5:27 am
Forum: ROMS Discussion
Topic: Nonbreaking wave-induced mixing in ROMS?
Replies: 7
Views: 1532

Nonbreaking wave-induced mixing in ROMS?

Dear ROMS users, before duplicating already done job I am wondering if someone already did it? I am looking into the additional mixing (via k-e turbulence) with NONBREAKING(!) waves as additional term in the shear production term (along with shear and bouy2) as a function of depth, Hs and frequency ...
by jivica
Thu May 24, 2018 4:56 am
Forum: ROMS Problems
Topic: ROMS is blowing up after just one time step
Replies: 15
Views: 2073

Re: ROMS is blowing up after just one time step

Have you checked your initial condition's file values? This line looks really suspicious: 1 2018-01-01 12:15:00.00 NaN NaN NaN NaN Do you have missing values in the ini.nc? What about bry.nc? atmo.nc? you have to check all those files to be sure you are not giving NaN to ROMS model... Cheers I.
by jivica
Mon May 21, 2018 7:24 am
Forum: ROMS Bugs
Topic: calendar in dateclock.F
Replies: 0
Views: 821

calendar in dateclock.F

I think that gregorian_proleptic is not a standard name according to https://www.unidata.ucar.edu/software/thredds/current/netcdf-java/CDM/CalendarDateTime.html Should be proleptic_gregorian Seems semantic but ncview is picky and some other software as well. Easy fix in Utility: dateclock.F:809: cal...
by jivica
Mon May 21, 2018 7:02 am
Forum: ROMS Bugs
Topic: River forcing LtracerSrc and setting temp to 0 ?!
Replies: 3
Views: 984

Re: River forcing LtracerSrc and setting temp to 0 ?!

I think I've fixed it. It is working now as I would expect (MP_DATA) and the problem was in step3d_t.F and pre_step3d.F (note that "_original" is current version in roms trunk) Shortly, indexing was (I think) wrong in the step3d_t.F and in pre_step3d.F part with source direction == 1 was missing etc...
by jivica
Sun May 20, 2018 4:17 am
Forum: ROMS Bugs
Topic: River forcing LtracerSrc and setting temp to 0 ?!
Replies: 3
Views: 984

Re: River forcing LtracerSrc and setting temp to 0 ?!

Hi Kate
I am using MPDATA in old 3.4 and the latest code
And with old all is as should be... looking into the code
to figure what/when went wrong
Cheers Ivica
by jivica
Sun May 20, 2018 3:06 am
Forum: ROMS Bugs
Topic: River forcing LtracerSrc and setting temp to 0 ?!
Replies: 3
Views: 984

River forcing LtracerSrc and setting temp to 0 ?!

Think there is a "feature" introduced in the recent ROMS related to river source temperature set to 0C even you chose not the alter ocean ambient temperature. What I am doing is simulating only fresh water inflow using river forcing and setting source salinity but do not want to change ocean tempera...
by jivica
Tue May 15, 2018 1:12 am
Forum: ROMS Adjoint
Topic: Normalization coefficients
Replies: 4
Views: 1426

Re: Normalization coefficients

Regarding cores; I've made some benchmark for CRAY XC40 machine (dragonfly topology, ~70G/s infiniband), which agrees well with what long time ago Kate wrote, it is close to up to 20x20 tiles. In other words, you get reasonable speedup/CPU_cost if your tiles are up to 20x20. For your 600x400 that wo...
by jivica
Mon May 14, 2018 2:42 am
Forum: ROMS Adjoint
Topic: Normalization coefficients
Replies: 4
Views: 1426

Re: Normalization coefficients

Hi Ivan, I think that nothing is wrong, it just takes time. Your grid is not that small (and depending on your de-correlation scales) and given small number of cores it takes quite a lot. My recipe is to make exact solution only (!) for zeta (2D so should be faster) and then compare that with random...
by jivica
Fri Mar 09, 2018 5:00 am
Forum: ROMS Problems
Topic: Contact point problem with salinity in Two Way Nesting run
Replies: 6
Views: 1740

Re: Contact point problem with salinity in Two Way Nesting r

Actually, you might be able to run your child grid with time step the same as in your parent grid (large time step, efficient, as for example I have the case in coastal nested domain where is quite shallow region). If your child grid is in shallow region and parent has some deep parts then interplay...
by jivica
Mon Jan 15, 2018 5:18 am
Forum: ROMS Discussion
Topic: problem in tide simulation
Replies: 3
Views: 1020

Re: problem in tide simulation

This is really not enough info to give you any decent help. Checking ROMS output at the boundary and comparing it with tidal boundary forcing is just to make sure you haven't messed with creating forcing. What does your grid looks like, how big domain? do you resolve shelf and consequently resonance...
by jivica
Tue Nov 28, 2017 5:41 am
Forum: ROMS Discussion
Topic: unrealistic pressure gradient forming along bathymetry
Replies: 10
Views: 3229

Re: unrealistic pressure gradient forming along bathymetry

I think there was already a quite a lot discussion about HPGE and smoothing techniques in the forum. In short, I would make pure HPGE simulation using only vertical profile of T and S (using ANA_INITIAL), and no forcing, closed brys (or clamp to the same profile of T,S using again ANA_). You can com...
by jivica
Fri Nov 24, 2017 5:26 am
Forum: ROMS Problems
Topic: Boundary condition type definition at river point sources
Replies: 2
Views: 1076

Re: Boundary condition type definition at river point source

Seems to me that tidal forcing doesn't match :?: :!: at your northern region either, this is where your BC are. You can have amplification in the the bays, like resonance, and if you are having say 10% wrong BC than that can result with large errors elsewhere, even more in bays like that, sensitive ...
by jivica
Tue Jul 25, 2017 2:10 pm
Forum: ROMS Bugs
Topic: Difference between trunk/ROMS and COAWST/ROMS with Qout
Replies: 3
Views: 2062

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

Hey both Johns,

Actually it is ivica ;) and not ivicia

and if you remember we met at Canberra last year at IMOS/ACOMO when I showed you some ROMS data assimilation for WA...

Cheers,
Ivica@UWA
by jivica
Fri Jul 07, 2017 7:46 am
Forum: ROMS Usage
Topic: COAWST problem when using 2way WRF-ROMS
Replies: 6
Views: 2066

Re: COAWST problem when using 2way WRF-ROMS

Just to update you, I've managed to trap and fix it .. not nice and full fix but working on it. John's script to plot weights helped in a way to see that wav2atm and ocn2atm differ at the narrow strip along the coast i.e. landmask. I just plotted diffs btw them and then you can see the region clearl...
by jivica
Thu Jul 06, 2017 11:26 am
Forum: ROMS Usage
Topic: COAWST problem when using 2way WRF-ROMS
Replies: 6
Views: 2066

Re: COAWST problem when using 2way WRF-ROMS

Well it is strange, ROMS is crashing at the first step, or better to say between init and 1 step while as it didn't saved 1st his record. WRF is happily running till next coupling exchange even after ROMS died. I've put WRF output every time step and fields looks OK. So, I am confident that problem ...
by jivica
Wed Jul 05, 2017 12:05 pm
Forum: ROMS Usage
Topic: COAWST problem when using 2way WRF-ROMS
Replies: 6
Views: 2066

Re: COAWST problem when using 2way WRF-ROMS

I don't think it is the cause, as I wrote in my initial email WRF is bigger and encompassing ROMS (and SWAN) grids; My feeling is that mapping btw WRF and ROMS is the problem (scrip/weights), as they are not defined on the same grid (WRF is on 10km grid and ROMS 2km, so I have to use #define MCT_INT...
by jivica
Mon Jul 03, 2017 4:30 am
Forum: ROMS Usage
Topic: COAWST problem when using 2way WRF-ROMS
Replies: 6
Views: 2066

COAWST problem when using 2way WRF-ROMS

Dear all, I'm having problem with setting/using COAWST (rev 1169) for 2-way WRF and ROMS. The problem is with WRF sending unrealistic forcing fields to ROMS which is crashing at the 1st ROMS time step (can't save his as not ending 1st step). Possibly someone had similar problem and could shed some l...
by jivica
Fri Sep 23, 2016 5:21 am
Forum: ROMS Tools and Techniques
Topic: How to decrease nc file bulk using mexnc
Replies: 2
Views: 2117

Re: How to decrease nc file bulk using mexnc

Just a hint you can try saving variables as short using add_offset and scale_factor.
This is useful if you have already created netCDF files but want to make them smaller.
With nco you have already ncpdq doing that for you...
I.
by jivica
Thu Aug 18, 2016 5:16 am
Forum: ROMS Installation
Topic: upwellingcase
Replies: 8
Views: 4358

Re: upwellingcase

It is exactly telling you what is missing! In your .in file, at the top, you have line where you define location of varinfo.dat So, you can put right link or even copy varinfo.dat into the folder where you run your model, but it has to be correct. example: my ocean.in file has line with: ! Applicati...
by jivica
Thu Aug 11, 2016 4:56 am
Forum: ROMS Discussion
Topic: atmospheric forcing time problem
Replies: 3
Views: 2325

Re: atmospheric forcing time problem

If you are into bulk flux then you probably need heat flux etc. Here is mine forcing file ncdump -h to give you idea about structure (and you could see that in the $ROMS_SOURCE/Data/ROMS/CDL as well) ncdump -h forcing.nc netcdf forcing { dimensions: time = UNLIMITED ; // (56 currently) dim1 = 220 ; ...
by jivica
Fri Dec 04, 2015 4:53 am
Forum: ROMS Problems
Topic: problem with mpirun
Replies: 4
Views: 2883

Re: problem with mpirun

not sure if you made typo but when running with MPI you do not use "<" to feed in .in file
this is wrong: mpirun -np 8 ./oceanM<tidebasinforcing.in
and should be: mpirun -np 8 ./oceanM tidebasinforcing.in > log &

Cheers
I.
by jivica
Sat Jun 27, 2015 2:15 pm
Forum: ROMS Discussion
Topic: weird problem with forcing file
Replies: 10
Views: 1976

Re: weird problem with forcing file

This is because of your meteo data are flipped up-down, and lower left coordinate is not correct. You have to make sure that lower left coordinate is (1,1), because of search/interpolate algorithm within ROMS. Easy to check this is to make your lat coords (or lon, depending where you model and switc...
by jivica
Mon Jun 22, 2015 4:38 am
Forum: ROMS Source
Topic: ROMS code modernization
Replies: 1
Views: 2092

Re: ROMS code modernization

Should be easy if you are using git.
What speedup you got when compared to "standard" version?
I am interested into making benchmark for my application on my cluster (using intel compilers + openmpi 1.8.6)

Cheers,
Ivica
by jivica
Tue Dec 02, 2014 7:29 am
Forum: ROMS Problems
Topic: Cool SSTs over Congo's submarine canyon; is this realistic?
Replies: 4
Views: 1465

Re: Cool SSTs over Congo's submarine canyon; is this realist

do you interpolate atmo forcing (i.e. T2) before run or you leave that to ROMS linear interpolation? If you use ROMS then it might use atmo values over the land :!: and interpolate them as they were above the channel, hence used as forcing. To test that you have to interpolate atmo forcing before :!...
by jivica
Sun Nov 30, 2014 8:13 am
Forum: ROMS Discussion
Topic: one problem in run time
Replies: 7
Views: 1317

Re: one problem in run time

Well you can get better interconnection btw nodes to get better efficiency; For example I setup cluster last week (4 x Dell r620 E5 with 40 cores) using infiniband Mellanox FDR 56G/s (!!) and I have really nice low latency and high speed communication. So, in short you hit the wall, you have to chan...
by jivica
Thu Oct 23, 2014 6:44 am
Forum: ROMS Problems
Topic: When using the adjoint model, Rad+Nud boundary can't be used
Replies: 2
Views: 1434

Re: When using the adjoint model, Rad+Nud boundary can't be

This is not a problem in forward nonlinear model setup because you can radiate outside, however if using adjoint back in time then you are in trouble. That is basically what it is, so use clamped boundary and then you can get boundary increments from data assimilation fixing that if needed... Have a...
by jivica
Thu Oct 23, 2014 6:40 am
Forum: ROMS Problems
Topic: wetting and drying problem
Replies: 22
Views: 10301

Re: wetting and drying problem

Just wondering would it be possible to use one if statement in a way if cell is dry then use tracer old state (from the previous time step) without applying anything i.e. heatflux etc?
And keep doing that until it get wet again.
probably it is like that already in new version ;)
Ivica
by jivica
Thu Oct 23, 2014 6:33 am
Forum: ROMS Discussion
Topic: Digital filter initialization and negative time steps
Replies: 1
Views: 604

Re: Digital filter initialization and negative time steps

I had problem with data assimilation and shocks when applying all increments to the initial state at t0.
As far as I know Andy is working on it and in the near? future it should be available.
This is "normal" approach to use in the atmospheric DA.
Ivica
by jivica
Thu Oct 23, 2014 6:30 am
Forum: ROMS Adjoint
Topic: how to write different type of data into one obs file
Replies: 1
Views: 1644

Re: how to write different type of data into one obs file

Create 1 file using obs_merge.m or something similar. This will give you idea what to do.
$ROMS/matlab/4dvar/obs_merge.m

Ivica
by jivica
Sat Jul 12, 2014 9:20 pm
Forum: ROMS Problems
Topic: Near-bottom numerical instability over the slope
Replies: 16
Views: 5506

Re: Near-bottom numerical instability over the slope

viewtopic.php?f=17&t=1080#p12021

see at the bottom -> my last response, hope it helps.
cheers
I.
by jivica
Tue Jun 17, 2014 7:17 am
Forum: ROMS Discussion
Topic: tidal analysis of large output files
Replies: 22
Views: 6581

Re: tidal analysis of large output files

If you decide to go the way of python or matlab then you can use pp (parallel python) or parallel toolbox in matlab and distribute the node_loop along the cluster (cores); doing that you are minimizing memory footprint and speedup calcs. Note that MPI is reducing memory footprint and not openmp (cor...
by jivica
Thu Apr 10, 2014 11:59 am
Forum: ROMS Usage
Topic: Construction of forcing file (stationary forcing)
Replies: 5
Views: 2549

Re: Construction of forcing file (stationary forcing)

glad it helped
the best option is to use ana_ if you can, it is then, during compilation, included directly into the code, so no reading at all.
Cheers,
Ivica
by jivica
Thu Apr 10, 2014 7:33 am
Forum: ROMS Sediment
Topic: Mediterranean Sea Model setup advises
Replies: 5
Views: 4166

Re: Mediterranean Sea Model setup advises

OK, looking into .h more deeply I can notice that you use central 2nd order for advection, why not standard upwind 3rd order is there some specific reason for C2? you can add diffusion and/or viscosity within ocean.in... If you want more/less temp/salt conservative then use MPDATA, I do not use smag...
by jivica
Thu Apr 10, 2014 7:25 am
Forum: ROMS Discussion
Topic: A benchmark case for tide
Replies: 1
Views: 803

Re: A benchmark case for tide

I believe it is too simple to have it ;) You can setup it easily without any grid or forcing, just using ana_* For example use ana_grid.h to make rectangular basin, in that case you can increase resolution within ocean.in file where you specify number of cells, no need for recompile (specify some an...
by jivica
Thu Apr 10, 2014 7:14 am
Forum: ROMS Usage
Topic: Construction of forcing file (stationary forcing)
Replies: 5
Views: 2549

Re: Construction of forcing file (stationary forcing)

If you have static field you can create forcing file with actually 2 time steps sms_time 0 and 1000000 because ROMS is linearly interpolating to your model time. Having those limits encompassing your model run period you are sure to have it correct within the model interval. You can use ROMS CDL exa...
by jivica
Thu Apr 10, 2014 6:59 am
Forum: ROMS Sediment
Topic: Mediterranean Sea Model setup advises
Replies: 5
Views: 4166

Re: Mediterranean Sea Model setup advises

From your .h file I can see that you do not have realistic atmo forcing, you do not have river forcing, etc ... It is hard to expect model to reproduce physical features if not using physical assumptions. To get correct dynamics in the Adriatic I use 39 river sources, realistic as possible atmo forc...
by jivica
Thu Apr 10, 2014 6:53 am
Forum: ROMS Discussion
Topic: potential temperature
Replies: 1
Views: 796

Re: potential temperature

You can use http://www.teos-10.org/software.htm which is "official"
There is matlab suite of routines and is easy to follow.
Cheers,
Ivica
by jivica
Wed Feb 26, 2014 8:51 am
Forum: ROMS Discussion
Topic: Input Sources/Sinks forcing file
Replies: 4
Views: 1479

Re: Input Sources/Sinks forcing file

What I did for one test case of LNG (temp) impact to the environment was: run the model with source but neutral with respect to temp (as suggested switch F for temp). Before that switch option, many years ago :), I had to run model without point source and saved station file with temp at the positio...
by jivica
Wed Jan 29, 2014 7:23 am
Forum: ROMS Usage
Topic: river points
Replies: 6
Views: 3161

Re: river points

This topic is recycled over and over, many posts about river forcing, position, etc have been answered and all you need is to look *search* and scratch a bit under the forum, after all this is the main goal of the forum. Rich Signell wrote matlab function to plot you river sources for a given grid a...
by jivica
Thu Jan 09, 2014 8:36 am
Forum: ROMS Discussion
Topic: courant
Replies: 3
Views: 4751

Re: courant

From my point of view it is waste of CPU time to have CFL 0.1, but should look at the dynamics and make a few tests to be sure how far with dt you can go. I believe Kate once wrote that for Rutgers version CFL should be less than 0.85. Quick and dirty way you can compute it by yourself, for barotrop...
by jivica
Wed Jan 08, 2014 8:07 am
Forum: ROMS Tools and Techniques
Topic: wind force
Replies: 2
Views: 1440

Re: wind force

This was asked before :!: so have a look at the
viewtopic.php?f=17&t=2238

basically you do not put number of forcing files == 0 but leave it empty

Ivica
by jivica
Mon Dec 16, 2013 9:08 am
Forum: ROMS Discussion
Topic: some questions about the WET_DRY
Replies: 5
Views: 1270

Re: some questions about the WET_DRY

Use the source Luke ;)

ROMS/Nonlinear/set_depth.F

there is all you need, you can invent even new function ;)

Cheers,
Ivica
by jivica
Tue Dec 10, 2013 8:33 am
Forum: ROMS Problems
Topic: Stormsurge simulation problem
Replies: 7
Views: 3237

Re: Stormsurge simulation problem

Got some spare time and quickly fixed it; so the steps were: 1) for your problem you can use 2D, no need to mess with 3D and is by far faster stromsurge.h: /* Momemtun equation */ #define UV_VIS2 #define UV_ADV #define UV_COR #undef UV_QDRAG #define UV_LDRAG #define CURVGRID #define AVERAGES #define...
by jivica
Mon Dec 09, 2013 10:25 am
Forum: ROMS Problems
Topic: Stormsurge simulation problem
Replies: 7
Views: 3237

Re: Stormsurge simulation problem

I have looked at your files and there are many things you should reconsider; you want to get storm surge but configure application to be nor full atmo driven nor only storm surge which is a bad idea. You want bulk flux but put constant swrad and lwrad through the whole simulation. Make up your mind ...
by jivica
Sat Dec 07, 2013 9:23 am
Forum: ROMS Problems
Topic: Stormsurge simulation problem
Replies: 7
Views: 3237

Re: Stormsurge simulation problem

It is hard to debug your application from the info you gave.
If you can put your grid file and forcing so I can grab it (when got time)
and see what is going on.
Ivica
by jivica
Thu Nov 21, 2013 7:56 am
Forum: ROMS Problems
Topic: Sea cliffs
Replies: 8
Views: 4591

Re: Sea cliffs

If you care only about storm surge possibly you are happy with 2D barotropic case which has only N=! !! What is killing you (blowup) is not rx0 but rx1, this is important to notice (i.e. density fields, number of vertical nodes, theta_b, ... all of them play a role). For wet_dry one have to be caref...
by jivica
Fri Nov 08, 2013 7:59 am
Forum: ROMS Benchmarks
Topic: Experience with new Intel Xenon Phi
Replies: 1
Views: 2229

Experience with new Intel Xenon Phi

Dear community,
does any of you have experience with new Intel Xeon Phi arch and ROMS performance?

Thanks in advance for your time,
cheers
Ivica
by jivica
Fri Nov 08, 2013 7:56 am
Forum: ROMS Tools and Techniques
Topic: TPXO to ROMS EXPRESS (tide extraction software)
Replies: 40
Views: 21131

Re: TPXO to ROMS EXPRESS (tide extraction software)

I think tidal ellipses, to be more precise inclination, is not rotated to the ROMS grid within forcing file. It is rotated during applying BC see angler(i,j) in the last line of code , in Nonlinear/set_tides.F: # if defined UV_TIDES ! !----------------------------------------------------------------...
by jivica
Fri Nov 08, 2013 7:33 am
Forum: ROMS Adjoint
Topic: Observations for normalization?
Replies: 3
Views: 2266

Re: Observations for normalization?

You are right, you do not need obs file just to compute norms, they depend only on the grid. Whenever you change something with the grid you have to recompute norms. They are expensive, so be sure the grid is final before kick. To be more precise you need to define in addition error decorrelation sc...
by jivica
Fri Mar 01, 2013 7:48 am
Forum: ROMS Discussion
Topic: stiffness number for rx0 and rx1
Replies: 11
Views: 9926

Re: stiffness number for rx0 and rx1

You want to have as real as possible bathy (tides sensitive) so smooth as less as possible. In one of my real applications, there are many big and small passages for tides and currents to in and out which usually with big slope on depth. The problem is i can not smooth these deep-narrow patches as ...
by jivica
Thu Feb 28, 2013 2:41 pm
Forum: ROMS Discussion
Topic: stiffness number for rx0 and rx1
Replies: 11
Views: 9926

Re: stiffness number for rx0 and rx1

I am afraid I do not follow you; If you for example copy your original roms grid to say new_grid.nc (in order to have original to compare with) and than smooth variable h and put it back inside that new file. Than in ocean.in you have to change GRDNAME to point to the new_grid.nc; run roms and check...
by jivica
Wed Feb 27, 2013 9:23 am
Forum: ROMS Discussion
Topic: stiffness number for rx0 and rx1
Replies: 11
Views: 9926

Re: stiffness number for rx0 and rx1

Well it is really important topic, so here is my 2c. All is your grid, problem to solve related question. For example if you want to do internal tides or tides only than you will look at the some spectral lines and currents that are generated because of HPGE (horizontal pressure gradient errors) are...
by jivica
Fri Feb 22, 2013 7:49 am
Forum: ROMS Usage
Topic: SWAN input - offline coupling
Replies: 2
Views: 1271

Re: SWAN input - offline coupling

just hit search for SWAN in the forum and got :!: +371 posts :!: related to it
by jivica
Fri Feb 22, 2013 7:44 am
Forum: ROMS Usage
Topic: SWAN input - offline coupling
Replies: 2
Views: 1271

Re: SWAN input - offline coupling

Well that is one of the " many time commented/discussed topic " so I really do not get it why not to look what other people wrote; strongly would suggest to search at the forum posts first. With really small effort, if you look at the roms source folder then you will find src/User/Functionals/ana_ww...
by jivica
Tue Nov 20, 2012 8:52 am
Forum: ROMS Information
Topic: CUBIC_MASKED INTERPOLATION
Replies: 0
Views: 1367

CUBIC_MASKED INTERPOLATION

Dear developers,
is there working version of Utility/interpolation.F that is using options CUBIC_MASKED along with LIMTR capable of taking care LANDMASK?

Thanks,
Ivica
by jivica
Wed Jan 18, 2012 10:20 am
Forum: ROMS Documentation
Topic: imposing tidal boundary conditions in parent and child 1 way
Replies: 5
Views: 5185

Re: imposing tidal boundary conditions in parent and child 1

Well I have to say that I do not use the latest version :oops: of ROMS (using version 3.4 from Aug 2011)so I am doing all in a classical 1-way nesting. For that I found the easiest and by far fastest and memory efficient is to use "station file" approach. I create child grids (as much as I want) and...
by jivica
Tue Jan 17, 2012 8:29 pm
Forum: ROMS Documentation
Topic: imposing tidal boundary conditions in parent and child 1 way
Replies: 5
Views: 5185

imposing tidal boundary conditions in parent and child 1 way

Dear all, in the past months there was never-ending story about tides, tidal dynamics BC inside ROMS model, hence think that some of you could benefit from paper we have recently published. Analysis of tidal dynamics as a function of BC that you use in a constellation of parent-child (or even only p...
by jivica
Tue Jun 21, 2011 8:12 pm
Forum: ROMS Discussion
Topic: Dwave
Replies: 3
Views: 1311

Re: Dwave

OK figured it out, It is in the degrees and during the loading phase it is converted into radians (varinfo.dat, the last line) 'Dwave' ! Input 'wind-induced wave direction' 'degrees' ! [radians] 'Dwave, scalar, series' 'wave_time' 'idWdir' 'r2dvar' 0.017453292519943295 ! pi/180 This resolves analyti...
by jivica
Mon Jun 20, 2011 10:20 pm
Forum: ROMS Discussion
Topic: Dwave
Replies: 3
Views: 1311

Re: Dwave

Thanks John on swift replay! Nevertheless, I am not using SWAN but other wave model (WWM), the most important is to have the same convention DEGREES not RADIANS (Nautical SWAN convention) calculating direction from the North clockwise and that direction, as you said, is pointing direction "from wher...
by jivica
Mon Jun 20, 2011 7:26 am
Forum: ROMS Discussion
Topic: Dwave
Replies: 3
Views: 1311

Dwave

Dear users, I have question regarding NEARSHORE_MELLOR05 . In the varinfo.dat I have found that Dwave (Direction for waves) should be in degrees, so I have used forcing file where I defined Dwave in degrees from North CW. Should I rotate the values into grid coordinate, I mean to rotate for -angler?...
by jivica
Sat May 14, 2011 12:56 am
Forum: ROMS Problems
Topic: MORE SEAGRID PROBLEMS
Replies: 6
Views: 3550

Re: MORE SEAGRID PROBLEMS

in the old matlab "finite" was OK but in new version it is replaced with "isfinite"
so go into the code and find/replace "finite" with "isfinite"

Cheers,
Ivica
by jivica
Fri May 13, 2011 9:39 pm
Forum: ROMS Problems
Topic: Seagrid Error
Replies: 26
Views: 8478

Re: Seagrid Error

You can do landmasking by a) using landmask utility in roms matlab toolbox or b) if you have coastal line in polyline (closed lines stored in matlab struct ie poly(1).lon poly(1).lat and so on for other islands like poly(n).lon poly(n).lat ...) you can do it easy in matlab by checking roms grid poin...
by jivica
Fri May 13, 2011 12:48 am
Forum: ROMS Problems
Topic: Seagrid Error
Replies: 26
Views: 8478

Re: Seagrid Error

It is simply not true that MASKED points are not used in computations! They are just zeroed out at the end of time step (!) we have been through that zillion times. ROMS is different animal than say SWAN which excludes those points from calculations. Navidad if you are still struggling with your gr...
by jivica
Thu May 12, 2011 10:14 pm
Forum: ROMS Problems
Topic: SOURCE POINT FOR RIVER RUNOFF (topic locked)
Replies: 115
Views: 24073

Re: SOURCE POINT FOR RIVER RUNOFF

ncdump is a application that is "made" when you installed netcdf library and is located inside netcdf bin folder, so in your case should be: /opt/netcdf/3.5.1-intel/bin/ncdump It is veeeery useful command and you will use it a lot. So if you can show us how does this grid file looks like: /opt/netcd...
by jivica
Wed May 11, 2011 8:43 pm
Forum: ROMS Problems
Topic: SOURCE POINT FOR RIVER RUNOFF (topic locked)
Replies: 115
Views: 24073

Re: SOURCE POINT FOR RIVER RUNOFF

give us:
ncdump -h your_grid_file.nc

hope you used seagrid2roms.m because seagrid is giving you .mat file
cheers,
Ivica
by jivica
Tue May 10, 2011 12:52 am
Forum: ROMS Problems
Topic: SOURCE POINT FOR RIVER RUNOFF (topic locked)
Replies: 115
Views: 24073

Re: SOURCE POINT FOR RIVER RUNOFF

If you do not want to use any forcing file than simply leave it clear:
NFFILES == ! number of forcing files

FRCNAME ==
And that should do, do not put any number like 0.
This means that all needed forcing should be done through ANA_

Cheers,
Ivica
by jivica
Fri Mar 11, 2011 12:45 am
Forum: ROMS Adjoint
Topic: Backwards trajectory
Replies: 7
Views: 3905

Re: Backwards trajectory

Think UCLA group has offline float tracking code for ROMS that is capable of going back in time. So, basically you feed it with ROMS precomputed history files and use similar input file for initial float positions. There is a option for going frwd or back in time.. Capet and Sasa did that long time ...
by jivica
Fri Jan 28, 2011 12:17 am
Forum: ROMS Bugs
Topic: possible bug in t3dbc_im.F, u3dbc_im.F
Replies: 0
Views: 935

possible bug in t3dbc_im.F, u3dbc_im.F

Think that there is a possible bug :?: in routines for applying boundary conditions for tracers and velocity when used NUDGING option. Issue is regarding time step in the last term t(Istr-1,j,k, nstp ,itrc) should be t(Istr-1,j,k, nout ,itrc) to be consistent. # ifdef WEST_TNUDGING t(Istr-1,j,k,nout...
by jivica
Fri Jan 14, 2011 12:53 am
Forum: ROMS Usage
Topic: Problems with basic model
Replies: 9
Views: 3879

Re: Problems with basic model

You should check your setup, think there is a problem 1. we need more info to help you out - do you use "tide.nc" file with tidal info for zeta and tidal ellipses or pure ANA? (if so you have specified ANA_FSOBC as well, hope you know what you are doing) 2. you used ndtfast==1 ?! if not in 2D does n...
by jivica
Wed Oct 13, 2010 1:20 am
Forum: ROMS Discussion
Topic: zeros in bathymetry
Replies: 8
Views: 3810

Re: zeros in bathymetry

As Hernan suggested run the case with some zero forcing and vertical (stable) density profile to see HPGE (Horizontal Pressure Gradient Error) effect. Just 10cents :idea: to the topic; you can afford higher rx0 in shallow region than in deeper (I do smoothing with different rx0 target in different r...
by jivica
Thu May 06, 2010 1:11 am
Forum: ROMS Discussion
Topic: model horizontal velocities error
Replies: 4
Views: 2269

Re: model horizontal velocities error

I should say that this is really good to do at the begining :wink: , just to make sure you have all setup right (!) and can explain "ghost currents" if you have them :roll: What I do is run model similar like you but with analitical (for my region specific) vertical T & S; just to be sure that there...
by jivica
Fri Feb 12, 2010 12:50 am
Forum: ROMS Discussion
Topic: offline-nested BCs vanish!
Replies: 2
Views: 1095

Re: offline-nested BCs vanish!

if you are adding tides to the open boundaries then you have to do something like: /* if you are adding tides on turn on SSH_TIDES and UV_TIDES */ # define SSH_TIDES # define UV_TIDES # ifdef SSH_TIDES # define ADD_FSOBC # define WEST_FSCHAPMAN # else # define WEST_FSCHAPMAN # endif # if defined UV_...
by jivica
Wed Feb 10, 2010 12:42 am
Forum: ROMS Bugs
Topic: Problem with scaling of non-gridded data
Replies: 5
Views: 1757

Re: Problem with scaling of non-gridded data

Are you sure that you are using the last version of ROMS? This was corrected recently, we had the problem with atmo forcing scaled to short int (to save diskspace)! Think there are still some issues if trying to read boundary files using scales (short int), but those files are small anyway and we ar...
by jivica
Thu Jan 21, 2010 12:53 am
Forum: ROMS Discussion
Topic: how set a vertical discharge?
Replies: 8
Views: 1931

Re: how set a vertical discharge?

Well be sure that you know what are the limits of the model that you use... In your example you have vertical velocities order of m/s which is comparable with horizontal ones.... hydrostatic/non-hydrostatic models? From my experience I use "near-field model" for close to the source dynamics and afte...
by jivica
Tue Jan 19, 2010 6:53 am
Forum: ROMS Discussion
Topic: how set a vertical discharge?
Replies: 8
Views: 1931

Re: how set a vertical discharge?

Be careful, the last time I used it was not working properly. Application was: Q_SOURCE with sediments and after when all sediment was on the bottom there was no mass conservation. On the other hand while using UV_ there it was exactly. Do some tests to see if it is working OK, the code was running ...
by jivica
Fri Aug 28, 2009 12:33 pm
Forum: ROMS Discussion
Topic: FSCHAPMAN/M2FLATHER/M3RADIATION not fully radiative
Replies: 5
Views: 1669

Re: FSCHAPMAN/M2FLATHER/M3RADIATION not fully radiative

Well I would not be happy with that phrase "... a few cells inside things in velocity field is OK." We have seen this in our application for the Adriatic Sea when one can get strong along boundary current that is due only to boundary conditions; like yours.. Imagine that you have some tracer like T ...
by jivica
Fri Mar 20, 2009 12:32 pm
Forum: ROMS Bugs
Topic: possible bug when restarting from history file (missing val)
Replies: 1
Views: 1430

possible bug when restarting from history file (missing val)

Well I have trouble with restarting from history file, think the problem could be when reading initial fields like zeta, temp, ubar ,.... then mask is not applied like it should be... restarting his file that I am using have for example zeta (ncdump output:) like this: zeta = _, _, _, _, _, _, _, _,...
by jivica
Thu Mar 12, 2009 11:24 am
Forum: ROMS Discussion
Topic: Coupling ROMS and WRF
Replies: 3
Views: 2302

Re: Coupling ROMS and WRF

Well for fully coupling I can not help you :( , but if you need info for preparing WRF forcing for ROMS (variables, and so on) I can.. :wink: One thing that you can do is run ROMS forced with WRF (1 way) and than save in time SST from ROMS and then feed it into next loop for WRF... in a few iteratio...
by jivica
Fri Feb 06, 2009 10:16 am
Forum: ROMS Usage
Topic: struggling for sharp halocline
Replies: 9
Views: 6756

Re: struggling for sharp halocline

Jacopo, sure I do agree that hybrid vertical could be thing that improve many thinks, BUT you should know how rx* looks like (for the ALL POINTS of your domain!) not just max value. So you can pcolor it and see where it is 'insane' and believe me, at the regions HPGE is devel! you'll get your real s...
by jivica
Fri Feb 06, 2009 7:54 am
Forum: ROMS Usage
Topic: struggling for sharp halocline
Replies: 9
Views: 6756

Re: struggling for sharp halocline

Jacopo, have you checked your sigma layer angle with the respect to "z", and hydrostatic stability among rx0 and rx1 at the region? We have made some experiments with regions that have big slope of sigmas and result was plaque by HPG a lot, giving unrealistic velocity which after some time went to a...
by jivica
Wed Oct 29, 2008 7:56 am
Forum: ROMS Documentation
Topic: ROMS Presentaions
Replies: 1
Views: 2478

ROMS Presentaions

Dear users, I have to talk to the "non physical oceanographers" in our ministry of env. protection (decision makers and possible funding source) and make intro to the ROMS. Idea is ti show them quality of the model, capabilities in as many aspects as possible, from bio, sediment, chem to the pure ph...
by jivica
Tue Jul 08, 2008 7:16 am
Forum: ROMS Discussion
Topic: Problem installing plotting package
Replies: 3
Views: 1110

Re: Problem installing plotting package

Safe bet is if you get binary for NCL (compiled with g77 ... gcc) and compile roms_plotting package as well with g77 see if that is working *get netcdf compiled with g77 as well*.. compilation of ncl with ifort take some time and can be tricky.. Just one quick question for plotting; say temp field a...
by jivica
Fri Jun 20, 2008 7:17 am
Forum: ROMS Discussion
Topic: How to set variables in Time Stamps section
Replies: 5
Views: 1251

Re: How to set variables in Time Stamps section

TIDE_START is used in the code for tidal phase lag (check it by tour self) in simple manner. so If you have tidal phase lags calculated (using say t_tide or some other tools) for specific day say 20080610 and you MODEL TIME is for current time (say 20080620) thanyou have to specify the two times in ...
by jivica
Thu Apr 24, 2008 6:31 am
Forum: ROMS Information
Topic: server recommendations?
Replies: 2
Views: 1580

From my point of view of building small linux cluster for ROMS + WRF I have to say that I am not happy with performance on quad cores connected with 1Gb switch... I have HP Blade 460c (2x quad cores 2.8G) with 4G of RAM and HP procurve switch.. well it turns out that running on 2 blades (2x8CPUs) is...
by jivica
Thu Mar 27, 2008 10:20 am
Forum: ROMS Problems
Topic: Trouble in shallow water
Replies: 5
Views: 1889

What have you used for optimization of compiler?
We had problems when using too aggressive -03 option (default) and had to put -O2
just idea...
Ivica
by jivica
Wed Mar 26, 2008 7:27 am
Forum: ROMS Bugs
Topic: SWAN_COUPLING: Trunk-r165 vs Trunk-r168
Replies: 2
Views: 2207

I am just posting my comment on the WRF and ROMS; Discussing with J. Wilkin I got to quite happy usage of WRF and ROMS without messy cloud fraction at all. I am using net shortwave at surf., only incoming longwave and the rest for bulk (T2, mslp, U10,V10, Rh) and let ROMS compute longwave outgoing ....
by jivica
Fri Mar 14, 2008 12:08 pm
Forum: ROMS Tools and Techniques
Topic: NCL
Replies: 9
Views: 5124

Or just simply specify startdate in the first line of .in file (second variable) for ccnt?
by jivica
Fri Mar 14, 2008 11:35 am
Forum: ROMS Tools and Techniques
Topic: NCL
Replies: 9
Views: 5124

model time

I have problem using plotting package when I use TIMEREF=-2,
so I have offset of datenum(1968,23,5) which is state in @units and not in the add_offset ....
Should I change opencdf.F where are attributes checked to get right offset?
Thanks
Ivica
by jivica
Fri Mar 07, 2008 8:49 pm
Forum: ROMS Tools and Techniques
Topic: NCL
Replies: 9
Views: 5124

Thanks Kate; I find Hernan's ploting stuff very usefull and have installed it without a problem. More over, I would like to change vector's color and so on, play little bit with contour @res and other options, add something and put it to run via cron (without DISPLAY which matlab req.).. so I think ...
by jivica
Wed Mar 05, 2008 7:21 am
Forum: ROMS Tools and Techniques
Topic: NCL
Replies: 9
Views: 5124

Thanks for the reply, I have installed week ago :wink: , plotting stuff from roms svn, I am using it but seems to me it's easier to use scripts like .ncl where one can change more options like res@ and so on on the fly... At think that on the wiki you should put ncargwww and ncarg2gif commented that...
by jivica
Tue Mar 04, 2008 4:17 pm
Forum: ROMS Tools and Techniques
Topic: NCL
Replies: 9
Views: 5124

NCL

Dear users,
I am thinking to move on NCL to get ROMS figures instead of Matlab
so I am collecting scripts .ncl for vecplots and other fields..
If you are willing to share you'll save me a lot of time
Thanks in advance!
Cheers, Ivica
by jivica
Wed Feb 27, 2008 8:36 am
Forum: ROMS Discussion
Topic: open boundary condition problems
Replies: 6
Views: 3282

Well I can confirm that model blows up is one have quite large values (or low) for forcing.. What I do is I find nans and replace them with nanmean(field) rather than putting just 0 or some missing_value
... After that all is working oky..
Cheers,
Ivica
by jivica
Fri Jan 18, 2008 8:29 am
Forum: ROMS Benchmarks
Topic: ROMS performance characteristics and tuning
Replies: 2
Views: 3990

dual quad performance

Hej guys have any of you tried to load ROMS on that structure? I have build 6 HP blades 460c with dual quads so at the end 6x8=48 cores, 2G ram / blade... I have tried ROMS with ifort_10.1.x and pgi_7.0.7 as well mich1-27 and mpich2 latest one, have GB net with HP procurve switch 1800 G24... at the ...
by jivica
Tue Oct 16, 2007 7:00 am
Forum: ROMS Information
Topic: WRF and ROMS
Replies: 3
Views: 2402

Well...

Hej Kate (If you do not recall I'm the one with funny shirt at Gordon RC :D ), I am not able to run the same grids for WRF and ROMS, I use NCEP GFS 0.5deg for init of WRF at 20km and than I use nest to my domain at 5km and this is CPU disaster (on my hand made cluster) ... I can not afford it to run...
by jivica
Mon Oct 15, 2007 6:46 am
Forum: ROMS Information
Topic: WRF and ROMS
Replies: 3
Views: 2402

WRF and ROMS

Dear community, I am running WRF as forcing counterpart for ROMS and I am wondering what are impressions of any of you in the same business... I have question about forcing variables for ROMS, what do you use (some .m scripts?) for computing (or better saying preparing) forcing fields like relative ...
by jivica
Tue Jul 10, 2007 1:09 pm
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 6341

Thnks Rich, I have found it and test it as well. It is working the way I wanted :D this bug is fixed. (we were talking about that at the GRC; in old way with Q_PSOURCE I got negative sand_01 and/or mud_01 (it used UV_PSOURCE) and because I put it on the top of the column and in the middle of domain....
by jivica
Mon Jul 09, 2007 8:00 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 6341

In which version is that? official roms with svn or SED_CODE?
and BTW what is policy of that now, is there merging of SED stuff done occasionally as before, and we should look for the last SEC CODE or is it done in no time & we can find it in the last official roms version?
Thanks for the info!
Ivica
by jivica
Thu Jul 05, 2007 2:32 pm
Forum: ROMS Information
Topic: Vertical flux for pointsource
Replies: 1
Views: 1328

Vertical flux for pointsource

Dear all, I am wondering if some of you have used option in model for sediments by specifying vertical flux like: #undef UV_PSOURCE #define Q_PSOURCE (Q_PSOURCE should allow a vertical flux tracer instead of a horizontal flux) I have made test grid of 60x30x20 (grid cell size is about 50m in both di...
by jivica
Wed Jul 04, 2007 3:29 pm
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 6341

It is exactly what I needed ... John please if you can provide little bit of info about Q_PSOURCE option .... Like what direction notation is used (1 -1 downward or upward), and if it is specified at the surface like sand storm (just for testing how it works) at some rho points, than it actually act...
by jivica
Wed Jul 04, 2007 7:56 am
Forum: ROMS Bugs
Topic: Minus values of sediment next to sediment point sources
Replies: 18
Views: 6341

Well, I had the same situation when I played with sand storms applying sediment flux at the surface in the middle of the basin :D .... Point is that I (and think you) want to have VERTICA L flux downward like rainfall and not horizontal like is implemented in the ROMS code... If using horizontal flu...
by jivica
Thu Apr 12, 2007 7:09 am
Forum: ROMS Installation
Topic: a problem in modeling ROMS with SWAN by using MCT
Replies: 18
Views: 7779

MCT compiling

Hej, just wondering if there is any gain in compiling MCT using mpicc?! (What do you use for netcdf lib compilation?) We are using SWAN+ROMS but using for fortran mpif90 and for preprocessor cpp. One more thing, if using ifort and -O3 we got blowups quite often, and by reducing so aggressive optimiz...
by jivica
Tue Mar 06, 2007 8:03 am
Forum: ROMS Installation
Topic: start file
Replies: 2
Views: 1835

Well think that you are using symbolic toolbox in matlab where you have defined sub function...
what do you get by typing in matlab prompt: which sub
you can disable it by rmpath /your_matlab_path_to_symbolic_toolbox

Hope it helps
Ivica
by jivica
Thu Nov 30, 2006 4:27 pm
Forum: ROMS Usage
Topic: Open Boudary Conditions: Inflow
Replies: 11
Views: 7779

Tides with REDUCED OBCS

Well I just wonder
what to #def and #undef for specifying only TIDE_SSH
and not TIDE_UV?? but to be able to use FLATER OBC?
Thnks!
ivica
by jivica
Thu Nov 16, 2006 9:13 am
Forum: ROMS Bugs
Topic: MPI code bug in interpolating station locations
Replies: 10
Views: 4938

station locations for parallel run...

Well, For my application (Adriatic Sea in spherical) it is not working, your patch I ment. Still hindices does not find corresponding I,J in tiles. The simplest solution for me is to run serial code for 1 iteration with stations defined in station.in like 1 1 lon,lat then get from out_sta.nc Ipos an...
by jivica
Wed Oct 11, 2006 7:09 am
Forum: ROMS Bugs
Topic: MPI code bug in interpolating station locations
Replies: 10
Views: 4938

Yip, there is a problem when using more than 1 tile (*MPI*) at the connecting nodes with Lon,Lat , so you have to find locations in I and J before (*NOT LON,LAT) and put them in the station file. You can test this with serial run and then grab i,j model coordinates and feed them into station file an...
by jivica
Mon Oct 10, 2005 7:10 am
Forum: ROMS Installation
Topic: ifort 9.0 with Red Hat Linux EM64T
Replies: 11
Views: 7572

libs without root passwd

You can use LM_LIBRARY_PATH variable to specify directory which should be incluced in search path.
like (bash):
export LM_LIBRARY_PATH=$LM_LIBRARY_PATH:/usr/your_path_to_spec_lib

should work..
Cheers, Ivica
by jivica
Thu Sep 16, 2004 7:40 am
Forum: ROMS FAQ
Topic: Tidal forcing at the open boundary
Replies: 23
Views: 42532

tidal open boundary conditions

Well. I saw the example where you used ADCIRC tidal fields but you have prescribed tidal fields at all grid points instead only at the open boundaries. If I only want to have tidal forcing open boundary, how do I specify this in the NetCDF file? Do I need to set it at all other interior points to a ...