Search found 55 matches

by ymamoutos
Fri Aug 23, 2019 4:25 pm
Forum: ROMS Discussion
Topic: Barotropic model Blowing-up with high courant and dt = 1.
Replies: 6
Views: 212

Re: Barotropic model Blowing-up with high courant and dt = 1

oops I missed that

I am apologizing. Please ignore my post

Giannis
by ymamoutos
Fri Aug 23, 2019 4:23 pm
Forum: ROMS Discussion
Topic: Abnormal Cooling at Shallow Water Area
Replies: 7
Views: 396

Re: Abnormal Cooling at Shallow Water Area

Greetings Are you downloading ERA5 data from ECMWF site or from Copernicus CDS (https://cds.climate.copernicus.eu)? If your are downloading from 2nd source there is no need to use d_ecmwf2roms.m script because the fields that need to be accumulated are already accumulated according to ERA5 CDS manua...
by ymamoutos
Fri Aug 23, 2019 4:02 pm
Forum: ROMS Discussion
Topic: Barotropic model Blowing-up with high courant and dt = 1.
Replies: 6
Views: 212

Re: Barotropic model Blowing-up with high courant and dt = 1

Greetings,

how many barotropic steps you set (ndtfast)?

if ndtdast=1 try more i.e 10 for beginning and
check the value of maximum barotropic Courant number.
In any case - if i recall correctly - for stability
reasons it cannot exceed 0.8

Giannis
by ymamoutos
Fri Jun 22, 2018 9:59 am
Forum: ROMS Problems
Topic: ROMS is blowing up after just one time step
Replies: 15
Views: 2001

Re: ROMS is blowing up after just one time step

Greetings Do you have any specific reason to use the 3rd-order upstream split advection scheme for momentum and tracers? I am asking because I recall that whenever it tried to use these options in horizontal resolutions similar or higher than yours - 2km to 1km - I had stability issues. Undefine the...
by ymamoutos
Thu Sep 28, 2017 8:29 am
Forum: ROMS Ecosystem
Topic: boundary for biogeochemical-ecological modeling
Replies: 5
Views: 2528

Re: boundary for biogeochemical-ecological modeling

Greetings At the last record (328) - when your model is blowing up - you have a maximum speed at point (i,j,k)(006,023,20) of 48.6 m/s... It's clear that the problem is on your physical model as already mentioned. Check your restart file to see where this point is and what is going on. Try to reduce...
by ymamoutos
Tue Jun 06, 2017 10:15 am
Forum: ROMS Discussion
Topic: sst bad result near coast
Replies: 6
Views: 1662

Re: sst bad result near coast

Greetings Kate and sorry for my late reply. Yes it is possible to run it both ways. Currently I am using the default option. I'll probably use the non periodic option to check if there are any significant differences to my results. I don't know if it's right but my domain is relatively small and I w...
by ymamoutos
Sun Jun 04, 2017 7:27 pm
Forum: ROMS Discussion
Topic: sst bad result near coast
Replies: 6
Views: 1662

Re: sst bad result near coast

Thanks a million Kate. From the above script I understand that you
use the default option - when you don't set explicitly -p (periodicity)
flag - which is east-west periodicity with zero (0) overlapping points.

Giannis
by ymamoutos
Sun Jun 04, 2017 2:09 pm
Forum: ROMS Discussion
Topic: sst bad result near coast
Replies: 6
Views: 1662

Re: sst bad result near coast

Hi Kate Sorry for writing to this old post but I don't think there is a point to start a new thread. When you using sosie to extrapolate MERRA2 dataset are you using east-west periodicity (ewper flag) for the input grid? I am searching to MERRA2 site and manual to find information about MERRA2 grid ...
by ymamoutos
Mon May 15, 2017 9:03 pm
Forum: ROMS Discussion
Topic: Tidal Forcing File
Replies: 16
Views: 3522

Re: Tidal Forcing File

Greetings Miguel thanks for your reply. I am using two different regional models (Global Mercator with 1/12 degrees resolution and MFS 1/16 degrees resolution) to create the input data for mine and to check for any significant differences on my results. According to myocean Copernicus service docume...
by ymamoutos
Sun May 14, 2017 11:53 pm
Forum: ROMS Discussion
Topic: Tidal Forcing File
Replies: 16
Views: 3522

Re: Tidal Forcing File

John The sum of M2 and S2 harmonics is a bit bigger than 0.15m... :shock: I suppose that I can't use this dataset to validate my results. Once more any suggestions??? Kind Regards Giannis P.S. Sorry editing the post multiple times. I had to double check the harmonics because I was looking on the wro...
by ymamoutos
Sun May 14, 2017 11:01 pm
Forum: ROMS Discussion
Topic: Tidal Forcing File
Replies: 16
Views: 3522

Re: Tidal Forcing File

Greetings John first of all many thanks for your reply. I am using the latest version of ROMS code and time_ref = -2 still uses the same calendar with the older versions. Yes my domain is inside the Mediterranean sea. I forgot to mention to my previous post that I had already used the regional solut...
by ymamoutos
Sat May 13, 2017 4:17 pm
Forum: ROMS Discussion
Topic: Tidal Forcing File
Replies: 16
Views: 3522

Re: Tidal Forcing File

Greetings, I am also having - let's say - a similar problem. When I am applying tides on my model I get almost 2 times higher amplitude than the actual (I am comparing with a tide gauge data, see attached figure). During my tests I change several things and rerun but the results in all of the cases ...
by ymamoutos
Wed Apr 26, 2017 10:42 am
Forum: ROMS Ecosystem
Topic: Compiling error: dateclock?
Replies: 5
Views: 1761

Re: Compiling error: dateclock?

Greetings

FYI I am also experiencing the same
problem after updating to version 841.
I am using gfortran 4.9.2 on Debian.

Giannis
by ymamoutos
Fri Dec 02, 2016 12:52 pm
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Re: Surface and Bottom Boundary Layer combination

Hi Kate and sorry for the late response

I run it at the debugger using -g -fbounds-check flag
(I am using gfortran 4.8.4) just for a day (480 timesteps)
and it didn't have any problem. The 2D and 3D fields as far
I can see looks ok. Any other suggestion?

Giannis
by ymamoutos
Wed Nov 30, 2016 10:40 am
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Re: Surface and Bottom Boundary Layer combination

Yes I was defining SOLAR_SOURCE with the other vertical
mixing schemes (MY25 and GLS) and the water
type variable was read from the grid file - the
figure that I previously post.

Giannis
by ymamoutos
Wed Nov 30, 2016 12:13 am
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Re: Surface and Bottom Boundary Layer combination

Kate you are absolutely right. I change the values of
land - just set 2 and 3 all over the grid
- and the model run without blowing up. But I am still
curious. Why that happened with LMD and not with the
other schemes. Once more, many thanks Kate.

Giannis
by ymamoutos
Tue Nov 29, 2016 11:37 pm
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Re: Surface and Bottom Boundary Layer combination

Hi Kate

I am setting it through the wtype_grid
variable in my grid file. I had already used
the same grid with MY25 and GLS with
SOLAR_SOURCE and WTYPE_GRID options and
the model run without any problem.
I am attaching a figure with my wtype_grid
variable to see it.
by ymamoutos
Tue Nov 29, 2016 11:11 pm
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Re: Surface and Bottom Boundary Layer combination

Even more interesting. Turning off WTYPE_GRID
makes the model to run smoothly. Any ideas why
that is happening?

Giannis
by ymamoutos
Tue Nov 29, 2016 11:03 pm
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Re: Surface and Bottom Boundary Layer combination

I run it on the debugger (I am using gfortran) and gave me that: At line 123 of file lmd_swfrac.f90 Fortran runtime error: Index '0' of dimension 1 of array 'lmd_mu1' below lower bound of 1 That's interesting because, I am using WTYPE_GRID option. I 'll switch it of and see what is going to happen G...
by ymamoutos
Tue Nov 29, 2016 10:27 pm
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Re: Surface and Bottom Boundary Layer combination

Hi Kate and greetings from Greece,

no I am not using WET_DRY option.
It is also weird because in another
setup - only the central and north
part of the grid I am using with
almost 1km resolution - BKPP and
SKPP works without any problem.
Any suggestions?

Again thanks for the answer.

Giannis
by ymamoutos
Tue Nov 29, 2016 4:26 pm
Forum: ROMS Problems
Topic: Surface and Bottom Boundary Layer combination
Replies: 14
Views: 6152

Surface and Bottom Boundary Layer combination

Greetings, When I am trying to use LMD_SKPP and LMD_BKPP at the same time I am experiencing a weird blowup. The model blows up at the first step and according to log file, the point that the model blows up it's on the mask. I also use PERFECT_RESTART option and I notice that my mask at rho points is...
by ymamoutos
Fri Jun 10, 2016 12:09 pm
Forum: ROMS Installation
Topic: Blowing-up: Saving latest model state into RESTART file
Replies: 11
Views: 4250

Re: Blowing-up: Saving latest model state into RESTART file

Greetings, from your logfile, it seems that you have a fairly large vertical courant number (Cw) and maximum speed. From step 86 and after at point 078,153,32/31 the maximum speed is almost 8.0+(!!!) m/s. Is this one of your point sources (river)? Is river runoff mass transport vertical profile (riv...
by ymamoutos
Wed Mar 23, 2016 11:46 am
Forum: ROMS Discussion
Topic: SST relaxation
Replies: 7
Views: 3600

Re: SST relaxation

Greetings,

check this thread

viewtopic.php?f=14&t=3066

it's very helpful, especially Hernan's post. I think answers a part
of your question.

Giannis
by ymamoutos
Fri Feb 19, 2016 11:04 pm
Forum: ROMS Discussion
Topic: Optimal tilling once again
Replies: 4
Views: 1419

Re: Optimal tilling once again

Kate thank you very much for your immediate reply. Indeed i forgot to attach the interesting part :oops: . I am attaching now and is huge :shock: Thanks again. Nonlinear model elapsed time profile, Grid: 01 Allocation and array initialization .............. 19.252 ( 0.0001 %) Ocean state initializat...
by ymamoutos
Fri Feb 19, 2016 10:42 pm
Forum: ROMS Discussion
Topic: Optimal tilling once again
Replies: 4
Views: 1419

Optimal tilling once again

Greetings, I understand that asking once more about this matter is kind of annoying and i am apologizing in advance, but to be completely honest I am completely confused /lost after reading several posts concerning ROMS tilling. I recently gained access to a HPC system ( https://hpc.grnet.gr/en/ ) a...
by ymamoutos
Wed Jan 06, 2016 6:24 pm
Forum: ROMS Installation
Topic: Problem of abnormal cooling
Replies: 9
Views: 3565

Re: Problem of abnormal cooling

Change your OBCFAC to a bigger number (ex. 120.0)
As Kate said check your surface heat fluxes and you
should turn off WET_DRY.I suggest to run with
BULK_FLUXES activated. Use TS_U3HADVECTION and
TS_C4VADVECTION for tracers.

Giannis
by ymamoutos
Mon Jan 04, 2016 9:30 pm
Forum: ROMS Installation
Topic: Problem of abnormal cooling
Replies: 9
Views: 3565

Re: Problem of abnormal cooling

Greetings, I suggest first to change the lateral boundary conditions for zeta and ubar,vbar. Use chapman implicit if you are planning to use flather for ubar,vbar or Chapman explicit for Shchepetkin. As Kate already told you at that post https://www.myroms.org/forum/viewtopic.php?f=31&t=4072 if your...
by ymamoutos
Tue Dec 22, 2015 7:05 pm
Forum: ROMS Installation
Topic: mpirun problem?
Replies: 7
Views: 4997

Re: mpirun problem?

Greetings, recently i gain access to a HPC system and I had the same problem. System has various mpi libraries (openmpi and intelmpi). With openmpi 1.8.5 for all available compilers (gnu,intel) i had no luck but when i start using openmpi 1.8.7 had no problem at all. I suggest to change your openmpi...
by ymamoutos
Mon Mar 23, 2015 7:04 pm
Forum: ROMS Discussion
Topic: Vertical mass fluxes issues
Replies: 0
Views: 988

Vertical mass fluxes issues

Greetings to all, here at the University of the Aegean we have been using ROMS to simulate processes taking place in deep secluded sub-basins during and following dense-water formation episodes. We try to close the mass budgets during and following deep-water formation episodes. In order to understa...
by ymamoutos
Mon Oct 27, 2014 1:16 am
Forum: ROMS Problems
Topic: Unexpected(?) blow up
Replies: 7
Views: 3048

Re: Unexpected(?) blow up

Hernan, thank you very much for making clear to me all this. The script I've used for the forcing files is d_ecmwf2roms.m. Concerning Tmax, the total number of days for this run this is 5114, 1096 is the last day of the aeg_wind_era_2002.nc file. ... Input Forcing File 07: /home/yiannis/ROMS/aegean/...
by ymamoutos
Sat Oct 25, 2014 6:42 pm
Forum: ROMS Problems
Topic: Unexpected(?) blow up
Replies: 7
Views: 3048

Re: Unexpected(?) blow up

Greetings, Lanerolle thanks for your advise, I'll check H+zeta and I'll make a new grid with different hmin. Kate i don't have any problem viewing my wind file with ncview (i'll attach a figure). The minimum value for all records is -20.12764 m/s the maximum 18.314148 m/s and for the last record min...
by ymamoutos
Sat Oct 25, 2014 2:03 am
Forum: ROMS Problems
Topic: Unexpected(?) blow up
Replies: 7
Views: 3048

Re: Unexpected(?) blow up

:shock: :shock: :shock: Kate thank you so much once more. I have check everything except from that... The weird is that when i open the aeg_wind_era_2002.nc file on that record (2920) doesn't have any extreme value.How is that possible? To be honest i feel bit silly but I've never seen anything like...
by ymamoutos
Fri Oct 24, 2014 4:34 pm
Forum: ROMS Problems
Topic: Unexpected(?) blow up
Replies: 7
Views: 3048

Unexpected(?) blow up

Greetings, I am running a 14 year simulation over the Aegean sea but for some reason, which I am struggling to understand, the model blows up after 4 years. To be honest I am surprised because the current configuration of the model appeared to be ok. I have run the same region for 6 years with diffe...
by ymamoutos
Mon Dec 09, 2013 9:45 pm
Forum: ROMS Problems
Topic: ecmwf2roms
Replies: 10
Views: 4974

Re: ecmwf2roms

Dear fereshteh,

if you check the lines from 465 to 468 of d_ecmwf2roms.m you will see
that evaporation is used on case "swflux".

Have a pleasant evening
Giannis
by ymamoutos
Sun Dec 08, 2013 4:58 pm
Forum: ROMS Problems
Topic: bulk file
Replies: 1
Views: 1052

Re: bulk file

Greetings,

varinfo.dat, which is located at /<your folder>/trunk/ROMS/External
have a short description of every field that ROMS can handle (input
and output). Also check roms_metadata.m (~/trunk/matlab/netcdf/).

Have a pleasant evening.
Giannis
by ymamoutos
Thu Oct 17, 2013 3:15 pm
Forum: ROMS Tools and Techniques
Topic: Problem
Replies: 5
Views: 4513

Re: Problem

Hi Thanos, I have never used Roms_tools for grid and forcing files. I am using the last 2.5 years the matlab routines that are provided through the svn repository. If you download the latest version of Rutgers ROMS matlab tools I am sure that you are not going to have any problems with the netcdf ro...
by ymamoutos
Wed Oct 16, 2013 3:29 pm
Forum: ROMS Tools and Techniques
Topic: Problem
Replies: 5
Views: 4513

Re: Problem

Greetings, you can use geodesic_dist.m script which is located in ROMS/matlab/utility folder. It's the same with the old spheric_dist.m plus it gives you the term galpha which is the Geodesic azimuth from point 1 to point 2 clockwise from North in decimal degrees. To be honest I am using seagrid. Gi...
by ymamoutos
Wed Sep 04, 2013 7:26 pm
Forum: ROMS Discussion
Topic: Question about ECMWF data
Replies: 0
Views: 1006

Question about ECMWF data

Greetings, I have been running simulations for the north Aegean region over a month now, and I have used two different atmospheric datasets. The configuration for all simulations was the same and the only difference was the atmospheric forcing files. I used the BULK_FLUXES formula, clamped obc for m...
by ymamoutos
Fri Jul 12, 2013 5:03 pm
Forum: ROMS Discussion
Topic: Blowing up at the last timestep
Replies: 7
Views: 1864

Re: Blowing up at the last timestep

Evridiki,

if your POM outputs are in netcdf form and in mercator
projection you can use the d_mercator2roms.m script
to create your initial conditions file. Also I have
the same question with Kate about the vertical interpolation.

Giannis
by ymamoutos
Fri Jul 12, 2013 3:21 pm
Forum: ROMS Discussion
Topic: Blowing up at the last timestep
Replies: 7
Views: 1864

Re: Blowing up at the last timestep

Greetings, When I try to open your log file my editor tells me that there was a problem opening the file... At boundary conditions section I see weird characters for salt and I don't see any choice for TKE. I think a good choice for your log file is to set on your *.in file NINFO == 1 to see at each...
by ymamoutos
Tue Jun 25, 2013 12:38 pm
Forum: ROMS Installation
Topic: Error in compile
Replies: 15
Views: 9755

Re: Error in compile

fereshteh, perhaps you have download and build only the C libraries of netcdf. Check here http://www.unidata.ucar.edu/downloads/netcdf/index.jsp and you will see that if you want to use the latest version of the netcdf library you must download two zip files. The first contains the C language librar...
by ymamoutos
Thu May 23, 2013 3:15 pm
Forum: ROMS Discussion
Topic: Surface heat retention, year-round fake "upwelling"
Replies: 21
Views: 8738

Re: Surface heat retention, year-round fake "upwelling"

Greetings from Greece and many thanks again Hernan.

Giannis
by ymamoutos
Thu May 23, 2013 2:00 pm
Forum: ROMS Discussion
Topic: Surface heat retention, year-round fake "upwelling"
Replies: 21
Views: 8738

Re: Surface heat retention, year-round fake "upwelling"

Hi, I am currently using BULK_FLUXES for my application and i have a kind of silly/stupid question because i am bit confused. Which atmospheric fields are necessary, if bulk_fluxes isn't defined, except from uv surface stresses, solar shortwave radiation, surface freshwater flux, surface net heat fl...
by ymamoutos
Wed Apr 24, 2013 10:12 pm
Forum: ROMS Discussion
Topic: why warm ssts in summer months ?
Replies: 4
Views: 1518

Re: why warm ssts in summer months ?

Greetings mashinde,

I would like to ask you if you are using any point sources (forcing file) for the inflow into the Mediterranean from the Dardanelles straits, Gibraltar and Nile river. If not why? These sources are critical for the circulation and T-S in the Mediterranean.

Giannis
by ymamoutos
Wed Apr 24, 2013 9:54 pm
Forum: ROMS Discussion
Topic: Is it artifact or not?
Replies: 3
Views: 1214

Re: Is it artifact or not?

Greetings, Kate thank you very much for your answer. I would like to ask a couple more thing about what you have suggested. If i am understanding right you suggest to use RadNud boundary conditions for my 3D fields (T,S,U,V). Is there any specific/accurate(??) way to calculate/set some how optimum n...
by ymamoutos
Thu Apr 18, 2013 5:20 pm
Forum: ROMS Discussion
Topic: Is it artifact or not?
Replies: 3
Views: 1214

Is it artifact or not?

Greetings, I am running a simulation of the North Aegean region for the period 2002-2007. The general picture, concerning circulation and tracers, seems ok except an small area close to the southern boundary. A cyclone is created in that area and also the temperature field is cooler from that was ex...
by ymamoutos
Fri Mar 29, 2013 3:21 pm
Forum: ROMS Ice
Topic: CHECKDEFS - cannot activate option: TS_MPDATA
Replies: 3
Views: 2798

Re: CHECKDEFS - cannot activate option: TS_MPDATA

Hi Kee, try to compile your code by enabling the MPI option on the build.bash/build.sh script, specifically uncomment USE_MPI and USE_MPIF90 commands. Now if your machine have, for example, the mpich2 library after your compile is finished type this echo MPD_SECRETWORD=roms > ~/.mpd.conf chmod 600 ~...
by ymamoutos
Wed Jan 30, 2013 7:47 pm
Forum: ROMS Problems
Topic: Unrealistic - extremely high zeta values
Replies: 3
Views: 1182

Re: Unrealistic - extremely high zeta values

Kate thank you for your reply and for your time, i would like to ask a couple more things. From your answer and from Marchesiello's paper i understand that i have to set Radiation and Nudging boundary conditions for the tracers (salinity and temperature), for the 3D u,v momentum components and final...
by ymamoutos
Sun Jan 27, 2013 10:04 pm
Forum: ROMS Problems
Topic: Unrealistic - extremely high zeta values
Replies: 3
Views: 1182

Unrealistic - extremely high zeta values

Greetings, i am trying to simulate the circulation of the North Aegean region. My problem is that after a year run zeta value ends to be approximately 5.7( :shock: :shock: !!!!) meters. I understand that probably one of the reasons (or the main ?) that causes the problem is that i am using the river...
by ymamoutos
Thu Mar 08, 2012 1:28 pm
Forum: ROMS Problems
Topic: SOURCE POINT FOR RIVER RUNOFF (topic locked)
Replies: 115
Views: 23768

Re: SOURCE POINT FOR RIVER RUNOFF

Hi to all,

My Problem is solved. It appears
that occurred due to small
number of points (5 points) for a relatively big
volume flow.

Yiannis.
by ymamoutos
Tue Mar 06, 2012 8:02 pm
Forum: ROMS Problems
Topic: SOURCE POINT FOR RIVER RUNOFF (topic locked)
Replies: 115
Views: 23768

Re: SOURCE POINT FOR RIVER RUNOFF

Hi to everybody, I am having a problem with point source option. I am giving with an external file the flow of sea of Marmara (5 points) through to Dardanelles straits to North Aegean sea but the only thing i am getting is blow ups...I am sure that my CFL condition is ok and the values of rx1 and rx...
by ymamoutos
Tue Oct 04, 2011 3:11 pm
Forum: ROMS Discussion
Topic: what's wrong with Roms3.2?
Replies: 8
Views: 2424

Re: what's wrong with Roms3.2?

Hi Kate and greetings from Greece, My problem is solved.I modified the ana_*.h files that i wanted for my test case following your suggestion (i set the variables that i needed equal to zero) and my model configuration is running smoothly until now. Thanks again for your help and have a pleasant eve...
by ymamoutos
Sat Oct 01, 2011 7:08 pm
Forum: ROMS Discussion
Topic: what's wrong with Roms3.2?
Replies: 8
Views: 2424

Re: what's wrong with Roms3.2?

Kate hi, on line 74 get_data.f90 calls subroutine get_2dfld for surface wind stress components. You need to provide a bottom salinity flux one way or another - most of us analytically set it to zero. So if am not providing an external file i have to modify ana_btflux.h and the other files if am unde...
by ymamoutos
Fri Sep 30, 2011 10:42 pm
Forum: ROMS Discussion
Topic: what's wrong with Roms3.2?
Replies: 8
Views: 2424

Re: what's wrong with Roms3.2?

Thanks for your advise Kate. My compiler is gfortran 4.4.3 and i used -fbounds-check flag. The result for my current cpp definitions file (myapp.h) is: At line 74 of file get_data.f90 Fortran runtime error: Array reference out of bounds for array 'frc', upper bound of dimension 1 exceeded (1 > 0) If...
by ymamoutos
Thu Sep 29, 2011 6:21 pm
Forum: ROMS Discussion
Topic: what's wrong with Roms3.2?
Replies: 8
Views: 2424

Re: what's wrong with Roms3.2?

Hi to everyone,

I've got "similar" problem.My ROMS version is 3.6. Kate you suggest to recompile with bounds checking. How can i do that?

Thanks for any help and for your time.