Ocean Modeling Discussion

yang

This problem has been reported and fixed, see the following links. Update your code to the latest version and this problem will disappear.

viewtopic.php?f=19&t=4717

https://www.myroms.org/projects/src/ticket/756

yang

Yes, i have the cycle_length attribute in the time variable (river_time). I agree with you that

arango wrote: we can have the ELSE conditional that works for multi files and time cycling

.

yang

My case is that the river_time is at the middle date of each month, for example [15, 45, 75, 105, 135, 165, 195, 225, 255, 285, 315, 345], and the river_time has a cycle_length attribute(cycle_length=360 days). Thus, when we start to run the model from 0.0 day, we need the river discharges at 0.0 ...

yang

In the version of svn 884, when i run the model with the climatology forcing data, i get an error as follows: Found Error: ** Line: 2240 Source: ROMS/Modules/mod_netcdf.F, netcdf_get_fvar NETCDF_GET_FVAR_0D - error while reading variable: river_time in input file: ../00_ECS_SCS/ECS_SCS_riv_clim.nc ...

yang

def_info.F is what writes the names of the forcing files into the global attributes. There's a string variable which has length 1024 - it could be that you're overflowing that in the call to join_string for all the forcing files of one type. Try changing the size of string to 2048 or so. Thanks ...

yang

Not yet. I have not solved this problem. If you figure this out, let me know. Thanks in advance.

yang

In my application in Eash China Sea using ROMS src 766, there are weird results for tracer mixing (see the attached picture), where CLIMA_TS_MIX is turned on for salinity and temperature, but turned off for tracer. The bug resulted in this wrong result is that the tracer mixing are wrongly ...

yang

Good work. very nice!
It works well, now.

yang

Matlab-2012b

The tide data i used is from ftp://ftp.oce.orst.edu/dist/tides/Globa ... tcdf.tar.Z

yang

OK, the attached map is corresponding to the following code.
TPXO=readTPXOdata(TPXOfile_elev,TPXOfile_vel,TPXOfile_grid,lonR,latR);
zz=squeeze(TPXO.h.z(:,:,8));
contourf(abs(zz))
colorbar

yang

I think that there is something wrong. For example, when i use this code to get tide, the Cmax is often larger than 100.[imgh=][/imgh]

yang

The code that loads up data structures with the 12 or 13 forcing file names is in inp_par.F, in a routine called load_s2d. I don't see anything there that should be a problem at 13 files, but Fname is dimensioned by 100, so more than 100 forcing files would overflow. Do you have 13 each of Uwind ...

yang

Wind forcing files in *.in file are: <<*********************** ../00_ECS_hindcast_forcing_2000-2013/ECS_2000_Wind_stress@ERA_interim_12hourly.nc | ! forcing file 1, grid 1 ../00_ECS_hindcast_forcing_2000-2013/ECS_2001_Wind_stress@ERA_interim_12hourly.nc | ! forcing file 1, grid 1 ../00_ECS_hindcast ...

yang

Thanks. Now, it works.

yang

I have download the files, but the .zip file can not be unzipped correctly. Can you give me the .m files ?

yang

How to set the time attribute 'cycle_length' when using the multiple input files? I am running the climatology case(for example, run 10 model years for spin-up). When the climatology file is given by one file(the cycle_length=360),the model can be nudged correctly. But when the climatology file is ...

yang

With or Without VOLCONS
the compution time is very different.
Without VOLCONS is 6 times faster than with VOLCONS.

yang

If you want to use nudging in M2radiation, you use cpps : East_FSchampman , East_M2Radiation and East_Nudging together. There is an error that the 'zeta_east' can't be writed !! this error can be overcomed by adding the following code in globaldefs.h #ifdef OBC_M2RADIATION # if defined SOUTH ...

yang

I think there is a mistake of "f1" at the line 734 in the lmd_skpp.F: f1(i,j)=5.0_r8*MAX(0.0_r8,Bfsfc(i,j))* vonKar / & & (Ustar(i,j)*Ustar(i,j)*Ustar(i,j)*Ustar(i,j)+eps) I think the vonKar is not neccessary? f1(i,j)=5.0_r8*MAX(0.0_r8,Bfsfc(i,j))/ & & (Ustar(i,j)*Ustar(i,j)*Ustar(i,j)*Ustar(i,j ...

yang

Warner is right.
In the latest simulation, i use 64 processor and it is ok.

yang

The ROMS run is instructed to run for 86400 time steps at 360 s per time step = 360 days The SWAN run is instructed to run from 1/1/04 to 12/10/04 = 344 days This will cause SWAN to stop coupling and ROMS to sit and wait. These times need to match, for both coupling and the length of your SWAN ...

yang

4 version of initialisation file Delft University of Technology name of institute 3 command file ref. number INPUT command file name 4 print file ref. number PRINT print file name 4 test file ref. number test file name 6 screen ref. number 200 highest file ref. number hi,nganju When i changed the ...

yang

ROMS/TOMS: DONE... Thursday - January 7, 2010 - 5:48:26 AM PGFIO/stdio: Too many open files PGFIO-F-/OPEN/unit=1046/error code returned by host stdio - 24. File name = ./WPO_SWAN_output_x5/tmbot.mat-129 In source file swanparll.f90, at line number 3764 ...

yang

I want to simulate the diffusion and advection of nutrient in the middle of open sea , and i find out that if i use point source, it will include the water mass. I can not give the initial value of the dye point, because the dye is released contiously in time. Is there a way to deal with this case ...

yang

Why is there no information about wind field in swan_inlet_test.in?
SWAN is forcing by What??
SWAN is forcing by the same wind field with ROMS??

yang

Thanks
Give my best wish to you.

yang

nganju wrote:I think the problem is the time step and duration of the ROMS and SWAN runs not matching. Can you show us the number and size of the ROMS time steps?

The number and size of the ROMS time steps is:

! Time-Stepping parameters.

NTIMES == 86400
DT == 360.0d0
NDTFAST == 90

yang

something else to try, more advanced: Edit SWAN/Src/waves_coupler.F and comment out line 648 (or thereabouts, depending on code version) from !----------------------------------------------------------------------- ! Create a restart file ...

yang

When i run ROMS+SWAN,i always get error in 16th model day in src433. Then i download 436, and get the warning information of compilation as following . Has it something with my error? PGF90-W-0129-Floating point overflow. Check constants and constant expressions (ocean_coupler.f90: 436) PGF90-W-0129 ...

yang

something else to try, more advanced: Edit SWAN/Src/waves_coupler.F and comment out line 648 (or thereabouts, depending on code version) from !----------------------------------------------------------------------- ! Create a restart file ...

yang

something else to try, more advanced: Edit SWAN/Src/waves_coupler.F and comment out line 648 (or thereabouts, depending on code version) from !----------------------------------------------------------------------- ! Create a restart file ...

yang

are there PRINT-0* files? What do they look like, at the end? Time of computation -> 20040116.000000 in sec: 1296000. accuracy OK in 97.01 % of wet grid points ( 98.00 % required) Time of computation -> 20040116.020000 in sec: 1303200. accuracy OK in 83.13 % of wet grid points ( 98.00 % required ...

yang

I have tried that,and the same error appears again in 16th day.
Does that means : SWAN input currents according the grid as following
INPGRID CURRENT CURVILINEAR 0 0 432 334 EXC 9.999000e+003 &
NONSTAT 20040101.000000 30 DAY 20041210.000000

yang

Thanks. The following is my configuration for SWAN. SET DEPMIN 5.00 INRHOG 1 NAUTICAL COORDINATES SPHERICAL CGRID CURVILINEAR 432 334 CIRCLE 30 0.04 1.0 30 READGRID COORDINATES 1. './WPO_wave_input/WPO_SWAN_coord.grd' 3 1 1 FREE INPGRID BOTTOM CURVILINEAR 0 0 432 334 READINP BOTTOM 1 './WPO_wave ...

yang

When i give 1 processor to SWAN, it will be stuck at 16th model day. It have been stuck at 16th day for 4hours.

yang

I run ROMS+SWAN from time 0(that is to say, i does not run it from restart file). 64 processors for SWAN. The same error apears again: PGFIO-F-217/formatted read/unit=116/attempt to read past end of file. File name = ./WPO_wave_output/depth.mat-001 formatted, sequential access record = 150081 In ...

yang

Now i use 16 processor for SWAN,and the same error appears again at the same time(16th model day.) the following give the content of the Errfile from SWAN. Errfile-001 Severe error : All free units used the error from SWAN-ROMS: - ROMS Export: vertically integrated v-momentum component (Min= -5 ...

yang

105E-141E;17N-42N;
433x335x21

yang

Thanks
128 processors for swan.
Now i change it to 1.

yang

In addition, i have test ROMS and SWAN alone. ROMS can run ten years,and SWAN one year alone. When they are coupled , they are ok in 16 model days. An error apears in the 16th model day. OCN2WAV - (06) imported and (05) exported fields, t = 3255 18:00:00 - ROMS coupling exchanges wait clock (s ...

yang

When i run ROMS-SWAN coupling model,i get a error at 16th model day? is that a bug ? Should there be 'close(UNIT)' in swanparll.f90. 781390 3255 19:00:00 1.011380E-02 2.056170E+04 2.056171E+04 1.379430E+16 781400 3255 20:00:00 1.002353E-02 2.056333E+04 2.056334E+04 1.379487E+16 781410 3255 21:00:00 ...

yang

I have run ROMS successfully in two way couplping with SWAN in my realistic case.
I have a question : if i give 144 processers to SWAN , it works well , and it does not work if i give more than 144 processers .
are there a limit number of processeres with SWAN??

yang

Thanks
best wish to you.

yang

Thanks Now i am using the configuration (#define TS_PSOURCE,#define UV_PSOURCE) accoording your sugestion. But i don't find the Warner's presentation in the folloing links? http://www.myroms.org/index.php?page=Workshops Can you give me a link in detail,or send me a copy of Warner's PPT? my address ...

yang

Thanks.
I have downded that 3 days ago and it is helpful for me. But there are not information about combination of open boundary CCP.

yang

If i want to include river in my model. Is the following configuration right? #define UV_PSOURCE /*use to turn ON or OFF point Sources/Sinks*/ #define Q_PSOURCE /*use to turn ON or OFF mass point Sources/Sinks*/ #define TS_PSOURCE /*use to turn ON or OFF point Sources/Sinks*/ I don't know the ...

yang

Kate，please, accept my very best wishes for the holiday season. According your advice, it is ok now that the tide simulation can be run coupled with current for a whole year. Another question, are there some other boundary CPPs which can be used in the combination form(for exmple: for M3 ,use both # ...

yang

Thank you a lot. Kate In the above runing , I have changed OBCFAC from 0 to 2.0d+1.It always blows up in the second day. OBCFAC == 2.0d+1 ! nondimensional Ydz changed from 0.0d0 In addition, i have make a configurating for nudging as following. TNUDG == 2*3.0d+1 ! days Ydz changed from 0.0d0 ZNUDG ...

yang

I am runing a basin scale climatology simulation in west Pacific ocean with grid 5'x5'. It will be ok for runing 5 years with the following configuration without tide. #define EAST_FSCHAPMAN #define EAST_M2FLATHER #define EAST_M3RADIATION #define EAST_TNUDGING If included tide,(that is to say ...

yang

I have two question about your configuration? 1. Can we use both M2RADIATION and M2NUDGING options together for 2D momentum in east open boundary?? #define EAST_M2RADIATION #define EAST_M2NUDGING I have not seen such configuratin in all of the *.h in "ROMS\Include\" . 2. How did you run the ...

yang

i am learning the ROMS-swan coupled system to study how sea wave effect the sea upper mixing layer.

yang

Thanks ,m.hadfield
Now, i am using OpenSUSE instead of Cygwin to run ROMS. and i succeed to run the Upwelling test case with OpenMP.

yang

Thank you very much.Kate
I installed Intel complier and set USE_MPI= , USE_MPI90=on and FORT=ifort, then the script woks.
if i set USE_MPI=on , USE_MPI90=on and FORT=ifort, then the script will don't work.

yang

When i use G95 and MPICH2-1.0.7 to compile ROMS, i get the following error. cd /root/roms/Projects/Upwelling/Build; g95 -c -fno-second-underscore -O2 -ffast-math -ffree-form mod_strings.f90 cd /root/roms/Projects/Upwelling/Build; g95 -c -fno-second-underscore -O2 -ffast-math mod_parallel.f90 Error ...

yang

i have used the same version of g95 to complie the netcdf and ROMS,but i get two different results. When i complie in x86-32, i succeed and roms can be run. however,i get the same error ,when i complie the ROMS at x86-64(Business vista 64 bit)

yang

when i following "Tutorial: Installing ROMS in Windows XP / Vista (via Cygwin)" to compile ROMS, i get the following error: cd /home/dd/roms/Projects/Upwelling/Build; g95 -c -fno-second-underscore -O3 -ffast-math -ffree-form mod_netcdf.f90 Fatal Error: While reading module 'netcdf' found module ...

Search found 56 matches