Ocean Modeling Discussion

shiyf

Hello, everyone, There is a problem that I am very confused about, I have been running for a month, ready to output a his file every day, the Settings of.in are as follows, but my output file is like this(roms_his_****.nc), and only one file is output, I did not find similar problems in the forum ...

shiyf

Have you solved this problem? I have the same problem

shiyf

dear all,

I want to change some program in ROMS for my study,

and want to find the program which contains the calculation of CW,

could you help me?

shiyf

Dear all, my roms model is running from 2016.01.01-2016.12.31, the frc data get from ECMWF in 3h step, the bry data get from HYCOM data in 1day step, but when the model running 201 days it is blow up, i check the temperature,and find that the pattern of temperature becomes more and more wrong , and ...

shiyf

As the model running, The temperature calculation will be abnormal around the land, And the area is getting bigger and bigger(the flowing figures are the 0002 and 0012 days ); But these are all water areas, i have checked the ini and frc data which are all ok( i think It has nothing to do with the ...

shiyf

I have down load the 3-hours step data of slhf, ssr,str,strd (W m-2 s)from ECMWF,

if the shflux(W m-2)= -(ssr+str+sshf+slhf) / (3*3600)?

could you help me, thank you very much!

shiyf

kate wrote:This is the next step. Learn to view the ROMS output.

Thank you very much kate!
The problem has been solved which is caused by steep depth in some special regions(see figure below),

I was even more smooth the deep, and the model runs smoothly!

shiyf

My model is blowing and I don't know why. This problem seems to have been discussed many times. Reference before discussion,I changed my time step, cheaked the grid file, DxMIN,DYmin,DZmin,and boundary conditions. However, the model is still blowing, I don't know how to solve this problem, could you ...

shiyf

You should recompile with USE_DEBUG turned on to see if you get more useful messages. It will be oceanG instead of oceanM. Have you tried starting with one of the pre-configured cases like UPWELLING? Thank you very much kate, I have found the mistake,This is a syntax error in my ocean.in files ...

shiyf

I build successfully,and has been submitted task several times,but today when I submited the task, there are some errors,i am very confuse,could you help me? The errors are: forrtl: severe (59): list-directed I/O syntax error, unit -5, file Internal List-Directed Read Image PC Routine Line Source ...

shiyf

hi,everyone, I am preparing the forcing data from NCEP dataset which are shflux,swflux,swrad, But i am not sure it is correct, could you help me,the data are: swrad get from Downward Solar Radiation Flux at surface minus Upward Solar Radiation Flux shflux get form -1.0* (net shortwave flux plus net ...

shiyf

Dear ROMS users- my forcing data is NCEP long term mean data which were averaged from 1981-2010; Uwind: surface u wind Vwind: surface v wind Tair : air temperature at sigma level 0.995, (degC) Qair: Relative Humidity at sigma level 0.995,(%) Pair :Surface Pressure,(millibars) Swrad: net solor ...

shiyf

hi, everyone, I have struggle of this problem so long,it is that the ocean temperature getting from model is inaccurate, my data is monthly mean.I have change forcing data source,or bulk option, and try kpp or MY2 mixing option, but the result from model is not much difference,As shown in the below ...

shiyf

hi, guys my result from ROMS run sames to have some problem, see:https://www.myroms.org/forum/viewtopic.php?f=17&t=3290,It is more and more deviated from the actual value from the first year to 15 years. I have check my forcing file,It seems to be no mistake,But What else can I do? Can you give me ...

shiyf

hi, I try to modify the kpp mixing scheme through adding a coefficient to AKT and AKv, I modified some fortran program to read the bv coefficient as a forcing file,as: nc = netcdf('E:\bvdata\bv_hrho1.nc', 'noclobber'); if isempty(nc), return, end %% Global attributes: %% (none) %% Dimensions: nc('xi ...

shiyf

hi, I have learn ROMs some time. I calculate 15 years,and get some results, when I check the 15th year result which are shown in the below figure, I am confused, I want to know Whether my result is reasonable,if not,how to improve my calculation? I use force,ini,bry and grid files, my cpp option ...

shiyf

Hi, I have run ROMS smoothly,but when i check the rerult,i find that the temperature which get from ROMS is much higher in some regions,and much lower in others, as is shown in the figure along with the latitude line.i have calculate 15 years,it is the 15th resut! could you give me some advice ...

shiyf

Hi, I have run ROMS smoothly,but when i check the rerult,i find that the temperature which get from ROMS is much higher in some regions,and much lower in others, as is shown in the figure along with the latitude line.i have calculate 15 years,it is the 15th resut! could you give me some advice ...

shiyf

kate wrote:Did you add a check for WAVE_MIXING in checkdefs.F? You can see how it's done there.

Thank you very much! it is ok!

shiyf

Did it attempt to recompile mod_forces first? If the makefile dependencies are correct, it should have. See if your build.bash is trying a "make clean" first for a fresh go at it. You have only declared the variables in mod_forces. You also need to allocate and initialize them. Copy the example of ...

shiyf

Hi, I have run roms smoothly at the first time,Now,i try add a certain coefficient to kpp scheme, the coefficient was stored as a forcing file which have variable,s_rho,time,lon_rho and lat_rho. I have modified five .f program and varinfo.dat file, as follows: mod_forces.F : # ifdef WAVE_MIXING real ...

shiyf

kate wrote:Yes, that's right, try it out.

Thanks a lot,it is ok!

shiyf

You have told it to run some NTIMES steps from the beginning. It knows it already did that and doesn't need to run longer. I suggest you add to NTIMES for more timesteps. thanks a lot! I do not know if I understand correctly .when i run the first time making nrrec =0 for one year, NTIMES steps * dt ...

shiyf

hi, I have run roms successfully first time, then i make nrrec = -1,and use ocean_rst file instead of inifile, it stop,and do not output data,when i check the log message, i can not find any error! I am very confused,Can you give me some advice? log file: Model Input Parameters: ROMS/TOMS version 3 ...

shiyf

hi， when i run roms with my own application, it blows up.i get the error message from the log file, NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00 (File: shiyf_ini.nc, Rec=0001, Index=1) - free-surface (Min = 0.00000000E+00 Max = 0.00000000E+00) - vertically integrated u-momentum c ...

shiyf

hi, I want to check if my results are right,My original data WOA2009 is z coordinates, NOW I want to know how to convert S coordinates in result data to z coordinates ? do you have good method to check the correctness of results? i am a new user,i am very looking forward to your any suggestions ...

shiyf

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME C => (i,j,k) Cu Cv Cw Max Speed 0 0 00:00:01 1.567195E-06 2.153599E+04 2.153599E+04 3.828818E+17 (166,001,05) 0.000000E+00 3.603793E-03 0.000000E+00 9.407347E-01 DEF_HIS - creating history file: ocean_his.nc WRT_HIS - wrote history ...

shiyf

hi, i am a new mode guy,when i run my application first,there are some errors: error while writing variable: zeta; the log file: Model Input Parameters: ROMS/TOMS version 3.6 Thursday - November 21, 2013 - 7:42:55 PM ----------------------------------------------------------------------------- SHIYF ...

shiyf

You should not need to download ROMS all over again. When you build ROMS, it creates intermediate files in some $SCRATCH_DIR location. Do you know where that is? You want to find that location and investigate what those files look like when there are problems like your bottom drag problem. Then ...

shiyf

Checking the read_phypar.F code shows that it will only try to read the stations file if STATIONS is defined. You can check the read_phypar.f90 intermediate file to check for the lines inside the #ifdef: #ifdef STATIONS fname=sposnam IF (.not.find_file(ng, fname, 'SPOSNAM')) GO TO 30 WRITE (out,230 ...

shiyf

That looks scarily like one of my cppdefs.h files. I see you have both floats and stations on. Do you have the floats and stations input files for it to read? I would debug this in a debugger to see where it's failing because it's not all that obvious from what you sent. hi, I delete the cpp ...

shiyf

That looks scarily like one of my cppdefs.h files. I see you have both floats and stations on. Do you have the floats and stations input files for it to read? I would debug this in a debugger to see where it's failing because it's not all that obvious from what you sent. In fact,I only want to ...

shiyf

It says you can't have Tnudg(2) set to zero, yet you have it set to zero. Have you tried changing it in your ocean.in file? Thank you very much！the error has been solved,But when I run again,it is shown that Output/Input Files: Output Restart File: ocean_rst.nc Output History File: ocean_his.nc O ...

shiyf

I would agree with Kate that linux is better os to install and use Pyroms. You will encounter enough problems in installing pyroms on linux taht i would not recommend winxp. For using pyroms you will have to install support packages such as netcdf 4.2, netcdf-cxx, netcdf-fortran etc. Detailed steps ...

shiyf

Hello everyone, I have got an error as follow: ( I don't know the reason, I am confused. Anyone can tell me how to fix it?) Model Input Parameters: ROMS/TOMS version 3.6 Friday - November 15, 2013 - 3:50:05 PM ----------------------------------------------------------------------------- SHIYF ...

shiyf

Hi ,
I have download pyroms and the Enthought Python Distribution on my winXP system,but i do not know how to install pyroms on EPD,
could you help me， Thank you very much!
best,
shiyf

shiyf

I have download pyroms and the Enthought Python Distribution on my winXP system,but i do not know how to install pyroms on EPD,
could you help me， Thank you very much!

shiyf

Now, I try to create my application.when i create initial filed using d_intial.m, there are some errors, I suspect that some of the parameters set wrong, the GRDname is my grid file name(seagrid.nc)? the OAname is my temp and salt file name(tempsalt.nc)? and I do not know what is according to set ...

shiyf

In my experience, the model behaves better when provided with consistent initial (and boundary) conditions from a global ocean model. We have python scripts to extract initial conditions from SODA. I'm sure the Matlab community has something similar. I want to know where can load python scripts and ...

shiyf

Now, I try to create my application.when i create initial filed using d_intial.m, there are some errors, I suspect that some of the parameters set wrong, the GRDname is my grid file name(seagrid.nc)? the OAname is my temp and salt file name(tempsalt.nc)? and I do not know what is according to set ...

shiyf

fereshteh wrote:hi
if you want to create Rectangular grid ,i suggest you use roms-agrif version. it is so easy to use program.
http://www.romsagrif.org/
be Successful
fereshte

thank you very much! Now,i try to use seagrid to create grid!

shiyf

hello,everyone, I am a new user, i will be very grateful for your any Suggestions. I want create grid from 5S-70N and 30E-255E. In the beginning, i used easy_grid.m script,and setting: " % Geographical and Grid parameters -------- lat = -5.0; % Latitude (degrees) of the bottom-left corner of the ...

shiyf

Hello everybody, I'm trying to set up the "Fjord tidal case" . So first I made the grid.nc and ini.nc files. Both files were created with easygrid.m and ROMS ran without problems. Afterwards, I build some files,eg,fjord.h,ana_fsobc.h,and so on according to the instructions on the web. ROMS ran ...

shiyf

shiyf wrote:Hello,everyone,

Now,I want to simulate the ocean temperature.First,i need to prepare the grid,initial field and forcing field.But,i do not know where can load the ROMS_tools of MATLAB, could you tell me ? thank you very mach!

shiyf

The problem has been resolved!!!

shiyf

Hello,everyone,

Now,I want to simulate the ocean temperature.First,i need to prepare the grid,initial field and forcing field.But,i do not know where can load the ROMS_tools of MATLAB, could you tell me ? thank you very mach!

shiyf

shiyf

You would need to have both Matlab itself as well as the ROMS matlab toolbox. Thank you very much for helping me!Problem has been solved. Now,there are some questions about Latte_C case.Instruction shows that"We have placed the CPP options file latte_c.h, standand input ocean_lattec.in, and a ...

shiyf

if i have Matlab tools in my cumputer,I do not need download, is it?

shiyf

I have run UPWELLING case successfully, and generate four nc files, but I can not open it,and i do not know which is the file we need that contains the result data, eg U,V,TEM,SAIT ?

could you tell me?
Thank you very much!

shiyf

shiyf

Thank you very much for your help, just contact mode, no experience for errors, please forgive me. the error have solved,but it have anther error: cd /lustre/home/shiyongfangWKJ/roms/projects/upwelling/Build; /lustre/inspur/mpi/openmpi-intel/bin/mpif90 -c -heap-arrays -fp-model precise -g -check ...

shiyf

Hello, everyone, When I first compile build.bash for the ROMS according to the instructions-Installing and Running ROMS for First Time Users , I come across a problem as below: ROMS/Bin/cpp_clean /lustre/home/shiyongfangWKJ/Projects/upwelling/Build/mp_exchange.f90 cd /lustre/home/shiyongfangWKJ ...

shiyf

Thank you for your help, I have added mpif90 path,but it have anther error: cd /lustre/home/shiyongfangWKJ/Projects/upwelling/Build; /lustre/inspur/mpi/openmpi-intel/bin/mpif90 -c -heap-arrays -fp-model precise -g -check uninit -ftrapuv -traceback -free analytical.f90 analytical.f90(396): error ...

shiyf

Hello, everyone, When I first compile build.bash for the ROMS according to the instructions-Installing and Running ROMS for First Time Users , I come across a problem as below: cd /lustre/home/shiyongfangWKJ/Projects/upwelling/Build; /lustre/inspur/mpi/openmpi-intel/bin/mpif90 -c -heap-arrays -fp ...

Search found 52 matches