Search found 71 matches

by ymamoutos
Tue Feb 27, 2024 4:02 pm
Forum: ROMS Discussion
Topic: Issue with Output File Generation in ROMS
Replies: 3
Views: 684

Re: Issue with Output File Generation in ROMS

Greetings John

Many thanks for the correction on
my post. I misinterpreted the
description on roms*.in file.

Kind Regards

Giannis
by ymamoutos
Tue Feb 27, 2024 8:38 am
Forum: ROMS Discussion
Topic: Issue with Output File Generation in ROMS
Replies: 3
Views: 684

Re: Issue with Output File Generation in ROMS

Greetings

you need to set NDEFHIS and NDEFAVG not equal to zero.
roms*.in file has an extensive description about these switches

Kind Regards

Giannis
by ymamoutos
Mon Oct 25, 2021 12:10 pm
Forum: ROMS Problems
Topic: Analytical surface shortwave radiation flux - RBL4DVAR algorithm
Replies: 0
Views: 2250

Analytical surface shortwave radiation flux - RBL4DVAR algorithm

Greetings Recently I change the atmospheric dateset that I was using with another from a local operational weather model because of its high spatial resolution. Model's operational output doesn't include shortwave and longwave radiation data and includes cloud fraction. So I am trying to compile the...
by ymamoutos
Tue Aug 10, 2021 4:18 pm
Forum: ROMS Problems
Topic: ROMS BLOWUP MaxDensity
Replies: 3
Views: 3106

Re: ROMS BLOWUP MaxDensity

Greetings

Is it possible to post your log file?
What are the options you are using
at the *.h? What is your dt and dtfast?

Giannis
by ymamoutos
Wed Mar 03, 2021 8:45 pm
Forum: ROMS Discussion
Topic: ROMS+SWAN - WAM cycle 4.5 boundary conditions
Replies: 1
Views: 3008

ROMS+SWAN - WAM cycle 4.5 boundary conditions

Greetings I am trying to setup a coupled ROMS-SWAN coastal application using COAWST. ROMS setup is OK and running smoothly but SWAN - coupled and not - gives the following error when I run it using debugger : forrtl: error (182): floating invalid - possible uninitialized real/complex variable. Image...
by ymamoutos
Mon Feb 01, 2021 3:52 pm
Forum: ROMS Problems
Topic: Error while executing ROMS for a self created grid region
Replies: 37
Views: 118318

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

Greetings from the log file you attached it seems that you have inconsistent calendar 18485.000 dstart 755666.00 time_ref dstart is the initialization time and your value (18485) refers to 01-January-2019 00:00:00 if your are using modified Julian days calendar and time_ref = -2.0 (see roms.in secti...
by ymamoutos
Fri Jan 29, 2021 3:37 pm
Forum: ROMS Tools and Techniques
Topic: Extraction of ERA 5 data for custom grid
Replies: 34
Views: 68913

Re: Extraction of ERA 5 data for custom grid

Greetings

nc_constant.m is part of ROMS matlab toolbox.
You can download it with svn using your ROMS credentials
from here:

https://www.myroms.org/svn/src/matlab/

Giannis
by ymamoutos
Sat Dec 12, 2020 7:43 am
Forum: ROMS Discussion
Topic: evaluating tidal constituents, amplitude and phases from water level data
Replies: 5
Views: 4476

Re: evaluating tidal constituents, amplitude and phases from water level data

Greetings

check this link https://github.com/wesleybowman/UTide
is the implementation of Matlab's UTIDE package.
According to author is still under development so
use it with some caution. Another option is TAPPY
(https://sourceforge.net/projects/tappy/)

Giannis
by ymamoutos
Wed Sep 02, 2020 4:58 pm
Forum: ROMS Discussion
Topic: swrad and lwrad not reading properly in ROMS
Replies: 6
Views: 6739

Re: swrad and lwrad not reading properly in ROMS

Hi Matheus check here https://confluence.ecmwf.int/display/CKB/ERA5%3A+data+documentation#ERA5:datadocumentation in Temporal frequency section. Also if you download, let's say, recent reanalysis data for example July 2020 some fields are not complete (i.e. total precipitation). There is a ~ 2 month ...
by ymamoutos
Wed Sep 02, 2020 10:47 am
Forum: ROMS Discussion
Topic: swrad and lwrad not reading properly in ROMS
Replies: 6
Views: 6739

Re: swrad and lwrad not reading properly in ROMS

Greetings, Are you using 6hr timestep ERA5 data? if yes your conversion from J/m^2 to W/m^2 (1/(6*3600)) is valid. If you are using a different timestep then you must change it accordingly i.e 1hr 1/3600, 3hr 1/(3*3600) etc. Also are you downloading from ECMWF site directly or from CDS (https://cds....
by ymamoutos
Wed Jul 15, 2020 9:09 am
Forum: ROMS Adjoint
Topic: Adjoint model blow up
Replies: 13
Views: 16154

Re: Adjoint model blow up

Greetings once more Try to run using the same options for boundary conditions, i.e. set NLM obc to be the same with adjoint. Another suggestion is to unset the relaxation options (QCORRECTION, SRELAXATION) and run the nonlinear model with BULK_FLUXES enable. Also which of the available 4DVAR schemes...
by ymamoutos
Tue Jul 14, 2020 12:19 pm
Forum: ROMS Adjoint
Topic: Adjoint model blow up
Replies: 13
Views: 16154

Re: Adjoint model blow up

Greetings

You are using nested boundary conditions in NLM model
but it seems you didn't enable the NESTING cpp option.
Also why you are using mixed radiation conditions in
NLM model? In general the choice of the boundary conditions
in the NLM model I think is completely inconsistent.

Giannis
by ymamoutos
Tue May 12, 2020 8:21 pm
Forum: ROMS Adjoint
Topic: Adjoint model blow up
Replies: 13
Views: 16154

Re: Adjoint model blow up

Greetings once more Raquel, Hernan and Joonho I am more than grateful for your suggestions. Finally I was able to run without any problem the 4DPSAS model :D . Raquel you had a point about my decorrelation scales. I found some literature about using variograms in order to estimate decorrelation scal...
by ymamoutos
Fri May 08, 2020 12:06 pm
Forum: ROMS Adjoint
Topic: Adjoint model blow up
Replies: 13
Views: 16154

Re: Adjoint model blow up

Greetings once more

After several test runs with small to very small dt (30sec to 1 sec and Maximum barotropic Courant Number = 3.56216612E-03)
the problem persists.

Any help or suggestion is more than welcome.

Kind Regards

Giannis
by ymamoutos
Thu May 07, 2020 5:48 pm
Forum: ROMS Adjoint
Topic: Adjoint model blow up
Replies: 13
Views: 16154

Re: Adjoint model blow up

Greetings after several weeks i was able to solve the blowup at the first step of the adjoint model. To be honest I had to redesign completely my approach in order to run - at some point - the PSAS model. At the moment I am assimilating only SST. The model again blows up at the first inner loop. I a...
by ymamoutos
Fri Apr 17, 2020 8:17 pm
Forum: ROMS Adjoint
Topic: Adjoint model blow up
Replies: 13
Views: 16154

Adjoint model blow up

Greetings, the last few weeks I am trying to set up a stable 4DVAR setup but I think that there is something fundamentally wrong on my approach. My application has relative high resolution - approximately 1.7km - with 30 vertical levels. The free runs - with or without tides - are completely stable ...
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: 3913

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: 5488

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: 3913

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: 12781

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: 7394

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: 5679

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: 5679

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: 5679

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: 12841

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: 12841

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: 12841

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: 12841

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: 6384

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: 15742

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: 15742

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: 15742

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: 15742

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: 15742

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: 15742

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: 15742

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: 15742

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: 9711

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: 6839

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: 4482

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: 4482

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: 8141

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: 8141

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 ...
by ymamoutos
Tue Dec 22, 2015 7:05 pm
Forum: ROMS Installation
Topic: mpirun problem?
Replies: 7
Views: 10872

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: 2380

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: 8003

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: 8003

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: 8003

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: 8003

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: 9845

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: 2577

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: 8074

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: 8074

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: 3098

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: 5726

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: 5726

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: 15991

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: 20560

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: 20560

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: 4170

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: 3782

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: 3782

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: 6052

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: 3619

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: 3619

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: 2603953

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: 2603953

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: 7182

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: 7182

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: 7182

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: 7182

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.