Ocean Modeling Discussion

mariafattorini

Dear all, Please find below the announcement for the JONSMOD 2018 Conference that will be held from 17-19 October 2018 in Florence , Italy. Aims and scope JONSMOD means Joint Numerical Sea Modelling Group and started in 1981. The general theme of JONSMOD conferences is ‘Latest modeling techniques f ...

mariafattorini

Hello to everyone again, I have not found yet where the problem is. :roll: I have looked for it at the lines of the files indicated by debugging message: 1. line 2674 of distribute.f90: DO j=LB2,UB2 DO i=LB1,UB1 np=np+1 WRITE(*,*) 'i =', i ! MARIA WRITE(*,*) 'j =', j WRITE(*,*) 'A(i,j) =', A(i,j ...

mariafattorini

Hello Everyone, I have tried to run the FTE double-gyre test, as distributed plus USE-DEBUG on, and I have got a segmentation fault: ... 119 4 23:00:00 -3.169893E-11 3.953715E-07 3.953398E-07 2.000000E+15 120 5 00:00:00 -3.120319E-11 4.231900E-07 4.231588E-07 2.000000E+15 TL_WRT_HIS - wrote tangent ...

mariafattorini

Hello all, I have a problem with a model on a horizontal jet. - I have run a model with UV_PSOURCE and TS_PSOURCE activated, and I have defined in ana_psource Tsrc=T0 (T0=20°Celsius). The resulting hydrodynamic fields are good. - Then I have run another model with a horizontal cold jet with again UV_ ...

mariafattorini

Hello jivica, of course I agree what you noted. In other words,if the horizontal velocity field is zero (i.e. the ambient water is at rest) ROMS can not simulate any vertical source point. Is it correct? In priciple any vertical velocity (let say acceleration) violate the hydrostatic assumption but ...

mariafattorini

Hello Kate, I did roms model with a dimension grid about 10m and discharge features given by a "near field" model (done with a other software), as suggested by jivica in the other post. I have an other question: in ana_psource, how is it distinguished a sink from a source? Thank you very much, Maria

mariafattorini

Hello all, Studying the roms, I have found the other following doubts: 1. What is the difference between R0 (density for state equation) and RHO0 (mean density)? 2. Usually what are common values of mixing coefficients (both vertical, Akt and Akv, and horizontal, VISC2 and TNU2)? 3. There are many ...

mariafattorini

Thank you for your replies. I am studing a submerged downward cold sea water jet. The jet comes from a port of diameter=1.9m that is located at 12m under the free surface. I used the Q_PSOURCE and TS_PSOURCE options. I set a horizontal grid spacing of about 1.7m in x and y directions, then, in the ...

mariafattorini

Dear all and thank you zhaoqian, still concerning the Q_PSOURCE option, I just noted the post at the page "https://www.myroms.org/forum/viewtopic.php?f=31&t=867&hilit=point+source", where jcwarner pointed out that, althought Q_PSOURCE is actually a vertical source term (see also http://www.myroms ...

mariafattorini

Dear all, I'm wondering if I can use roms in order to simulate the heat exchange between the ambient sea water and a costant discharge of cold water that is released at a point -10 m below the free surface and is downward directed (i.e. the horizontal velocity components are zero but there is a ...

mariafattorini

Hello, I need an advice: I have to define a discharge characterized by zero horizontal velocity (v=u=0m/s) and vertical velocity of 2m/s toward the bottom (w=-2m/s). How can I set this case? In the ana_psource.h it is possible set only u- or v- velocity components; is there nothing for w- vertical ...

mariafattorini

Hello Angelolemos, model inputs (as forcing, grid, boundary conditions, initial conditions, ect) can be defined analytically or through external files (netCDF format). I think that the easiest way is to define analytically them. To do this you have to define the appropriate CPP options in 'myfile.h ...

mariafattorini

hello kate,

- with dt=1 sec, the model blows up after 6 hours and 9 minutes;
- with dt=0.5 sec, the model blows up after 6 hours and 45 minutes;
- with dt=0.25 sec, the model blows up after 6 hours and 7 minutes;
- with dt=0.1 sec, the model blows up after 6 hours and 27 minutes.

mariafattorini

Hello, to build the example previously posted, I have defined the following obc and cpp to define inflow on West boundary: SOUTHERN_WALL NORTHERN_WALL EAST_FSCHAPMAN - EAST_M2FLATHER - EAST_M3RADIATION - EAST_M3NUDGING WEST_FSCLAMPED - WEST_M2CLAMPED - WEST_M3CLAMPED ANA_FSOBC - ANA_M2OBC - ANA ...

mariafattorini

Thank you very much Kate; I am beginning to understand how to set the boundary conditions. I should define a rectangular channel with the following characteristics: - the North and South boundaries are the lateral walls of the channel; - the West and East boundaries are open; - the inflow is equal ...

mariafattorini

Thanks Kate, your help is very useful for me. But I don't know which is the meaning of 'nudging time scales' and I am not sure where I have to set them: are they the parameters in the ocean.in called M3NUDG (timescale for attive condition) and OBCFAC (factor between timescales for attive and passive ...

mariafattorini

Hello, I proved that in absent of user definition of boundary conditions, the Roms model runs. I wonder which default boundary conditions are used in the model run. Moreover how can I set a classical "open boubdary"? I mean, a bc that let the flux to go through the boundary. Thanks very much for any ...

mariafattorini

Hello Kate and thank you, I have tried different setting: - just inflow on W boundary ( with u_west in ana_m3obc.h); - either inflow on W and equivalent outflow on E (with u_west and u_east in ana_m3obc.h). but in either cases the result is not good; the water oscillates between O and E boudaries ...

mariafattorini

Hello all, I am trying to create a rectangular channel with walls on the boundaries N and S. I'd like set an inflow on the W side that come out of the channel through the E side: to do this I have used ana_m3obc.h where I have set the following values: #if defined CHANNEL # ifdef WEST_M3OBC IF ...

mariafattorini

Hello Shchepet and thank you very much for your clear explanations. I have another question: I tried to simulate a tank with flat bathymetry and a hump in the middle of the base of the tank without forcing (not waves, nor tide, no wind, no inflow, ..) and with zero initial conditions, but the result ...

mariafattorini

Hello all, I am setting my first simple application and I have need to know some things: 1. is Akv the classic eddy viscosity or is it VISC? 2. is there a condition to choice the size of grid cells and the size of time step to maintain numerical stability, eg. a condition on the Courant number? 3 ...

mariafattorini

Hello all, I would like to see how it changes the Sed_test1 test-case if I add a wind forcing. To do this, I have to add the following cpp options in the header file: bulk_fluxes, ana_winds, ana_pair, ana_tair, ana_humidity, ana_rain, ana_cloud, longwave and I have to eliminate ana_smflux and ana ...

mariafattorini

Hello, I have corrected the values in the boundary ascii file (a my distraction!): 9.700 42.4 1 9.700 42.9 0 9.700 43.4 1 10.45 43.4 0 11.20 43.4 1 11.20 42.9 0 11.20 42.4 1 10.45 42.4 0 and now I get the following errors (in Matlab, using Seagrid): theMaskTool = [] theMaskToolFlag = [] ??? warning ...

mariafattorini

Hello all, I am trying to create the ascii boundary file for Seagrid. The boundary file contains as following: 9.700 4.24 1 9.700 4.29 0 9.700 4.34 1 1.045 4.34 0 1.120 4.34 1 1.120 4.29 0 1.120 4.24 1 1.045 4.24 0 where the first column is the longitude, the second one is the latitude and the last ...

mariafattorini

hello Duckweed,
if you load the grid.nc in Matlab, which variables apper in the workspace?

mariafattorini

Hi, Nacholibre, and thank you for your reply. I tried to do how you suggested me but perhaps I have not understood what I have to do. I loaded coastline and bathymetry data in Seagrid and I set my grid; then I have set "MinimumDepth" (in "Seagrid Setup" window) to 2 and I saved Seagrid.mat file. But ...

mariafattorini

Dear Wilkin, thank you very much for your help. Actually I used the old version of the routine. At the present I am trying to plot the vector field related to the "Sed_test1" test case. I have setted the following script: >> ncload 'ocean_his.nc' >> vplot=v(11,1,:,:); >> uplot=u(11,1,:,:); >> grd ...

mariafattorini

Hello all, I too have a problem with Seagrid use. I am preparing the ROMS input grid file using the source data arising from "http://topex.ucsd.edu/WWW_html/srtm30_plus.html" web site and the Matlab utility named Seagrid. Seagrid has produced the "s" matlab structured variable where "s.gridded ...

mariafattorini

Dear all, I am experiencing a similar error trying to plot the upwelling model results. I am using the grd=roms_get_grid('ocean_his.nc') but the command stops because the variables long_rho, lat_rho,.........are absent in the input file. Actually, by using ncload('ocean_his.nc') I have checked that ...

mariafattorini

Hello, I have a similar error message but I don't know where the error is. When the application runs I get: ... T Hout(idBott) Write out bottom property 18: dep_net Output/Input Files: Output Restart File: ocean_rst.nc Output History File: ocean_his.nc Output Averages File: ocean_avg.nc READ_PHYPAR ...

mariafattorini

yes, removing Lwave in ana_wwave.h the compiler runs correctly.
I want not being tedious but I would like to know when Lwave is used;

many thanks again, Maria

mariafattorini

Ok, it's clear now, I understand. thank you very much. I have added SSW_CALC_UB in the header file. Now only this error remains: s -O3 -ffast-math -ffree-form -ffree-line-length-none analytical.f90 analytical.f90:618.28: Lwave(i,j)=20.0_r8 1 Error: Unexpected STATEMENT FUNCTION statement at (1) make ...

mariafattorini

Hello jcwarner, you are right: in my mod_forces.f90 (in the Build directory) there are "Ub_swan" "Dwave" and "Pwave_bot" and there are not Hwave and Lwave. In my test I defined the following CPP options: #define UV_ADV #define UV_COR #define UV_LDRAG #define UV_VIS2 #define MIX_S_UV #define TS ...

mariafattorini

Hello Kate and Jcwarner, thank you much for your help. 1. in my "mod_scalars.f90" there is deg2rad; 2. the Lake signell test case is correctly compiled; 3. in my "mod_forces.f90" (in Build directory), any variable about wave (Hwave and friends) are not definited; 4. in the "globaldefs.h" (in the ...

mariafattorini

Thanks Jcwarner, I have tried "dos2unix -U ana_wwave.h", without good result. In the ana_wwave.h I have added the following lines about my application "PROVA1": ... #if defined BL_TEST wdir=210.0_r8*deg2rad DO j=JstrR,JendR DO i=IstrR,IendR Hwave(i,j)=0.5_r8 Dwave(i,j)=wdir Pwave_bot(i,j)=8.0_r8 END ...

mariafattorini

Hello all and many thanks to Kate for reply. now, when compiling I have the following error: s -O3 -ffast-math -ffree-form -ffree-line-length-none analytical.f90 analytical.f90:463.27: Hwave(i,j)=2.0_r8 1 Error: Unexpected STATEMENT FUNCTION statement at (1) analytical.f90:464.36: Dwave(i,j)=90.0_r8 ...

mariafattorini

Hello All, I would like to know if it is possible to build a 3D grid with Seagrid. I have downloaded and opened Seagrid but I seem to see only tools to build a 2D grid (for x- e y- directions but not for z- direction (depth)). Is the third dimension set by N variable in ocean_application.in? Thank ...

mariafattorini

Hello, I am trying to build a my simple application, but when I type "./build.bash" to compile it, the following error appears: ROMS/Bin/cpp_clean /home/maria/roms/Projects/prova1/Build/mp_exchange.f90 cd /home/maria/roms/Projects/prova1/Build; /usr/bin/gfortran-4 -c -frepack-array s -O3 -ffast-math ...

mariafattorini

Hello Kate and thank you very much for your reply.
I found the lines
*/TS_MDATA:
*/
in my header file and I removed them.
Now, I know that if I write something in header file, this something will be writen in the make_macro.mk

Thank you again,
maria

mariafattorini

Hello Kate and Arango, in my make_macro.mk (in the directory roms/trunk/Compilers)there are these lines: */TS_MDATA: */ USE_ADJOINT := USE_TANGENT := USE_REPRESENTER := USE_SEAICE := USE_WRF := USE_SWAN := USE_REFDIF := USE_ARPACK := USE_MCT := USE_ESMF := This content is as created from the ...

mariafattorini

Hello all, I have a similar messagge error: maria@PC-Satellite ~/roms/Projects/prova1 $ ./build.bash rm -f -r core *.ipo /home/maria/roms/Projects/prova1/Build Compilers/make_macros .mk Compilers/make_macros.mk:2: *** missing separator. Stop. I am working on the Windows Vista PC by CYGWIN with ...

mariafattorini

Hello all, I have compiled OTPSnc by gfortran-4, modifying makefile as my above post described. So, I have found some errors in the OTPSnc code that I reported below since they could be useful to other users: - File extract_HC Errors at row 231 ORIGINAL FRAGMENT 231) write(11,fmt)lat(k),lon0(k), 232 ...

mariafattorini

Hello all, thanks to Kate for your reply. Trying to create tidal forcing input file using OSU Tidal Prediction, I have found some problems. I am working on my Window Vista PC by CYGWIN. Following the instructions in the READ ME file: - Open the makefile e change it in the following red lines: ARCH ...

mariafattorini

Hello all,

If I have to set a constant wind (constant in the space and in time), can I use the ANA_WIND?

And, if ANA_WIND is defined (as in some test cases), where I set velocity and direction values?

I hope in your help,
Maria

mariafattorini

Hello, sorry again I discovered that I already had the package of gmake, or at least I think. Gmake is inlcuded in the ‘make’ package which is found in the CYGWIN's setup, is it right? So I don't understand why the CYGWIN doesn't find gmake command when it is called from the command line. I sent a ma ...

mariafattorini

Well, I perhaps understood the source of my problem:
I simply haven't the gmake command!!!
So I have downloaded the gmake package from the web
page "http://mirtchovski.com/p9/gmake/", and I am trying to install it.

thanks and bye,
Maria

mariafattorini

Hello, Thank you very much to Carlos and F.Daryabor. I have tried to use OSU Tidal Prediction, as suggested by Carlos and by Wikipage, but I found some problems. I am working on my Window Vista PC by CYGWIN. I have saved the OTPSnc package and the Mediterranean Model within the ‘home’ directory. Th ...

mariafattorini

Hello all, I am studying to build my first application. About my application, I have tide data at two points of the domain. The wikipage ‘Tidal Forcing Files in ROMS’ say that I should interpolate the tidal constituent database to the desired ROMS grid. I don’t understand how I have to do this inter ...

mariafattorini

Hello, many thanks for the above replies. I have deleted the previous installation and I installed the Plotting Package again. I installed the binary's NCL 5.0.1 and I downloaded the Plotting Package in Mydir directory. I followed the indications of the Tutorial (web page: https://www.myroms.org ...

mariafattorini

Dear all, I have a problem with the installation of the Plotting Package with Cygwin (my PC is a Window Vista). First I have installed a precompiled binary NCL (version 5.1.0) following the procedure on the web page http://www.ncl.ucar.edu/Download/cygwin.shtml.During this operation I don’t find the ...

mariafattorini

Dear Robertson, thank you very much for your quick answer. I followed the updated tutorial but again the nf_test.exe file is not created. However I went forward in the installation process. I opened the file ‘build.bash’ and changed some lines according to the instructions in the tutorial. In the ‘bu ...

mariafattorini

Hi everyone, I'm trying to install ROMS on my Windows Vista PC using Cygwin, but I got some errors. May you help me? I followed the installation procedure provided on Wikiroms (https://www.myroms.org/wiki/index.php/ROMS_Cygwin). After installing Cygwin the root folder has the following folders: bin ...

Search found 52 matches