Ocean Modeling Discussion

javadi

Hi, I have taken some actions but am still having problems. My data is for NEMO model in the whole Mediterranean Sea with grid [1016 380 141]but my grid is the east med sea with Lm=803 and Mm=629 and 20 layers. if directly use some code like d_inital. I faced some errors mentioned that the second ...

javadi

Hi all,

In the past, I used a global dataset to make initial conditions.

Now I try to use Mediterranean reanalysis data for ROMS as the initial condition, how can I do this? is there any code for other model output for ROMS?
For example, d_initial in MATLAB for this type of work?

Best,

javadi

Hi, I have run ROMS many times without any problems. No changing .in and.h, configuration .... Only I try to run for the increased domain. So, only I changed .nc files I have faced a memory corruption error. When I checked my log file, I saw that ROMS had faced NLM Horizontal, harmonic mixing ...

javadi

Here is my domain!

javadi

Thanks for your response!

The problem is the open boundary!

Up to now, I had one open boundary in the south but now I have two open boundaries both North and west.

In this case, what should I do as both are close together in the attached Fig.

javadi

Hi all,

I ran the model but I have faced an error. This error can be a result of low Ram or CPU or the problem is in input file.

Best,

javadi

Dear John, Now it is the first time I have tried to extract data with a domain that is both a negative and positive longitude(-5 to 5). I remember that your Matlab code worked only for positive or negative separately in the past. Your updated code now works for both negative and positive longitude ...

javadi

Hi John

Thanks for the update!

It works perfectly now!

Best,

javadi

Hi, I use the roms_wikin-master code to download ERA5 data from NCAR. It worked nicely. I have used it many times. Now I have faced this error; Maybe because of some changes in the NCAR accessibility. Please guide me to solve it Reading variable 1 of 9: mean sea-level pressure e5.oper.an.sfc.128_151 ...

javadi

Hi Thanks, kate for your response! I am a new user of Python so it is a little difficult for me. In the past, I used Matlab for plotting ROMS. I have two different files: My grid file is in this address in my system :/home/jabmubuntu/ROMS_OUTPUT/roms_grd.nc My average file is in this address in my ...

javadi

Hi I have installed Pyroms using kate's branch https://github.com/ESMG/pyroms. I try plotting figures on the model's output but cannot understand how to use it. There are some codes in pyroms_toolbox for plotting. For example zview. This is for depth I think. For plotting, where should I put Roms's ...

javadi

Hi wikin Thanks for your response. I was waiting for running to see outputs. Now there are many new products in average file. I am a little confused. The first thing is that the outputs is based on layers. So, I should sum of them for all layers. and then for transects I should sum along longitude(j ...

javadi

Excuse me.

This is true in my comments.

#define DIAGNOSTICS_TS
#define DIAGNOSTICS_UV
#define AVERAGES

not

# undef DIAGNOSTICS_TS and
# undef DIAGNOSTICS_UV
# undef AVERAGES

javadi

Thanks for your reply. I have done some action for it. 1) .h # undef DIAGNOSTICS_TS and # undef DIAGNOSTICS_UV # undef AVERAGES 2)NAVG = NDIA; should be the same 3).in, I activated this parameter Dout(iTrate) == T F ! temp_rate, ... time rate of change Dout(iThadv) == T F ! temp_hadv, ... horizontal ...

javadi

Hi all, In my domain, I need to calculate some parameters enters to the my basin from the strait. v is northward. v.dx.dz v.T.dx.dz v.S.dx.dz v.rho.dx.dz As I did not understand how can I calculate it after running, I decide to calculate them in running time. My first question is Hvom. It gives me ...

javadi

Hi all, I have model outputs and I have to calculate transport that enters my basin from the strait. I think that I can achieve it when defining some options before running. As and now I have not activated them. How can calculate it? Are there some codes to calculate transport or how can do this ...

javadi

Hi Thanks for your response. So how can I access them? I should run the model, and after that, I should use another code to calculate it separately. Or can we activate some options for example ocean.in the file to give use salt transport or temperature transport from one transect in history file or ...

javadi

Hi all, I try to calculate how much water enters from the strait to our basin. My box is the surface while bounding in the depth of 200 m. How can calculate how much heat enters temperature or salinity flux enters our basin when defining a transect? There are some codes for it to extract after ...

javadi

Vorticity from the history file and pvorticity from the average file!

javadi

Hi all I try to analyze model outputs. The model gives me vorticity and pvorticity. The vorticity values range from (1.5-2) *0.0001. This is relative vorticity and I think that it is rational. But, Pvorticity is probably potential vorticity (1/ms unit) but here the model gives me the value in some ...

javadi

Hi all, I have one question about jerlov water type WTYPE. This is in the .in file and we should define it 1 (clear ocean)to 7. Here is my question. In my domain, it should be better to have two water types in the western and eastern coast. The western is not clear and the eastern is clear. In this ...

javadi

Ok. It has been solved!

I try to explain it as others maybe use it!

There is a limitation in terms of the number of characters!

should use this
/scratch/f85015jb/MedSea/Ad/2012/201202.nc|

instead of:
/scratch/f85015jb/MedSea/Ad/2012/frc_MY_APPLICATION_ERA5_bulkflux_201202.nc|

Thanks!

javadi

I have faced an unusual error. when defining forcing files with multiple files, I have faced this error 0 2012-01-01 00:00:00.00 0.000000E+00 3.029579E+03 3.029579E+03 9.673655E+13 (000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 DEF_HIS - creating history file, Grid 01: roms_his.nc ...

javadi

Hi Based on this method, I have made it for 7 years(84 months). Now, I define it in .in file like this NFFILES == 1 ! number of unique forcing files FRCNAME == /scratch/f85015jb/MedSea/Adriatic/2012/frc_MY_APPLICATION_ERA5_bulkflux_201201.nc| /scratch/f85015jb/MedSea/Ad/2012/frc_MY_APPLICATION_ERA5 ...

javadi

Hi Wikin Thanks! This problem has been solved! But, Another point is that this error when applying on ROMS: REGRID - input gridded data does not contain model grid: Gridded: LonMin = 12.0000 LonMax = 20.5000 LatMin = 39.2500 LatMax = 45.7500 Model: LonMin = 12.0000 LonMax = 20.8000 LatMin = 38.9000 ...

javadi

Hi I try to make a forcing file using era5. I use roms_get_era5_NCARds633_bulkflux in this code which is for wikin. I can use it. The problem is bbox. latitude is work! Ok! My domain along longitude is 12-20.8 E. But, the forcing gives me -348-339.5. I put -348 -339.5, the problem has not been ...

javadi

Hi all I try to extract data from era5 using https://github.com/trondkr/ERA5-ROMS When running ECMWF_tools.py, I have seen roughly 10 variables(instead of 13) that are extracted and converted. But, I have faced this error File "/usr/lib/python3.8/multiprocessing/pool.py", line 771, in get raise self ...

javadi

After checking the outputs, the model tries to show too high a temperature, particularly in summer. I think that it would be a flux problem, a bad configuration. Are there any suggestions for it? here is my .h file. #define UV_ADV #define UV_COR #define UV_QDRAG #define UV_VIS2 #undef MIX_GEO_UV ...

javadi

Hi I have run the model with all files. But, the model shows the temperature a little higher than roughly 2 degrees. In this condition what we should do it? I have added some options to reform flux which is in CPP but the problem hasn't been solved. Another issue is that model has problems in the ...

javadi

Dear users I have run the model for all inputs. After one year, I check the outputs. The model shows higher temperature and salinity. 4 degrees higher temperature and 3 PSU more salinity near the surface particularly in hot seasons. When plotting the profile, I find the model has a problem ...

javadi

Hi I have not understood one point! We have three files output including history, average, and diagnostics! I am going to save the monthly mean file. For this reason, I put the average file, for example, NAVG=24*3600*30/100=25920(100 is the time step) to average one month for example from January ...

javadi

Hi Thanks! I guess that I have solved the problems of installation with your suggestions! However, I cannot use it! For example, file = 'E:\roms_his.nc' >> grd=roms_get_grid(file,file); I will face: Error using ncread Too many input arguments. Error in roms_zslice (line 56) t = ncread(file,time ...

javadi

Dear wikin I use some packages for plotting and can see the history file. But, honestly, I have problems with your package ! For a few days, I am working on it but I can't use it, and all of the time the error is g=roms_get_grid(file,file); You need to add stretching.m to your matlab path Error ...

javadi

I try to define a profile thus I weigh in different layers. Here is it

I considered s_rho=1 to s_rho=40 is not right?

Do you mention that I should s_rho=0 to s_rho=40? It will be 41 layer?

javadi

Hi

Thanks!

I have a dilemma about one point.

In my case, I have 40 layers.

I should define V_shape (s_rho)from 1 to 40 or 0 to 39?

Another dilemma would be the least number would be for the bottom!

Best,

Javad

javadi

Hi Thanks for your reply! I have checked u velocity. It is here. I guess that the flow is towards the left. The transport is positive so this can not be the problem. With regard to i,j, I consider i=15 and j=295 in my file(for example). is it right? the result is based on this approach! now should i ...

javadi

Dear all I run the model but I have faced a problem. I have 4 rivers and define monthly discharge with salinity=3 and Temperature=12. Here is the location of river points in my grid. As my river comes from left to right, I define river transport as positive as mentioned here https://www.myroms.org ...

javadi

Thank you so much for your response

I change only NAVG and I estiamte

NAVG=30.2*24*3600/400=6523

I put this value but the result would be:
15 Jan,16 Feb,17 Aug...

Excuse me!

Should I change other variables NDEFAVG ,NTSDIA,NDIA...?

javadi

Dear users, I am going to record the average of variables at the end of each month. My DT=400 and I guess that I should change this part of .in. By default, it records 8 hours. LDEFOUT == T NHIS == 72 NDEFHIS == 0 NQCK == 0 NDEFQCK == 0 NTSAVG == 1 NAVG == 72 NDEFAVG == 0 NTSDIA == 1 NDIA == 72 ...

javadi

Dear users

I have activated rho on .in and I can see it in my history file. But, the rho varies from 24.5 to 28.5 where it should be 4.5 to 8.5 in my domain based on my estimation.

Why the model does not calculate it correctly?

How to correct it?

javadi

Dear users I have a problem downloading data from ecmwf. I have downloaded data from this website and then use ecmwf2roms to make forcing files. The problem is that the first one timestep finishes and in the middle of the file, another time step begins. In other words, the downloaded data is not ...

javadi

Dear wikin Thank you so much for your suggestions. Excuse me! I am completely confused about this error! I applied your points but the error is the same. I add some extra months to the interval but the error is the same. I know why this occurrence. You mentioned that my forcing files can lead to the ...

javadi

Thanks for your response. It was my problem. How should I change my .h file for trunk code to be compatible with trunk codes. Now my .h is: #define UV_ADV #define UV_COR #define UV_LDRAG #define UV_VIS2 #undef MIX_GEO_UV #define MIX_S_UV #define SPLINES_VDIFF #define SPLINES_VVISC #define DJ_GRADPS ...

javadi

Dear kate I ask you one question. Based on your suggestion, I use this cpp (bulk fluxes)for my model. Does this work correctly in the trunk branch? or should I use your branch to use this option? #ifdef SOLVE3D # define CORE_FORCING # define BULK_FLUXES # define CCSM_FLUXES # if defined BULK_FLUXES ...

javadi

Dear users I used d_ecmwf2roms and bulk fluxes. I have faced a problem to use this script. Based on codes and this forum, I should use both 00:00:00 12:00:00 and all-time step (0 3 6 9 12) in this link https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/ to download data. Here is my ...

javadi

Here are my lrf_time and lhf_time variables from the input. It is clearly shown that why this error occurs.lrf_time which is related to lwrd do not cover the whole month. lrf_time = "2012-12-31 22:30", "2013-01-01 10:30", "2013-01-01 22:30", "2013-01-02 10:30", "2013-01-02 22:30", "2013-01-03 10:30 ...

javadi

Dear all,

I have checked. All inputs initiate from the same time but the model faces this error. I could not understand why the tmin is not the same for variables.

I use d_ecmwf2roms to mak forcing. How should I the same tmin and tmax?

javadi

Thank you for your response. The model runs for 30 days and then the model has faced the error. For example, I have started from January the error is 31 January. I have started on 1 February, the error occurs on 28 Feb. In other words at the end of each month, this error occurs. I have tested for ...

javadi

Dear users I ran Roms when considering all inputs for 6 years. I have faced an error when applying forcing file: SET_2DFLD - current model time exceeds ending value for variable: lwrad_down TDAYS = 30.3773 Data Tmin = -0.1250 Data Tmax = 364.8750 Data Tstr = 27.8750 Data Tend = 28.3750 TINTRP1 = 30 ...

javadi

No, this is not matter. I have corrected it. I think that i have different variables with different tmin. It creates problmes.For this reason٫ at the end of each month there is less input forcing٫leading issue. This problem is as a resualt of new .h. I ran the model for 6 years without any problem w ...

javadi

Dear Kate Thank you so much for your points. I have changed the time based on your sugestions but I have faced the same error. DSTART = 4749 ! days TIDE_START = 0.0d0 ! days TIME_REF = 2000101.0d0 GET_2DFLD - solar shortwave radiation flux, 2013-01-31 13:30:00.00 (Grid=01, Rec=245, Index=1, File ...

javadi

Excuse me! I am a little confused about how you estimate the time to start. Now, I have defined DSTART = 662214.4375 ! days TIDE_START = 0.0d0 ! days TIME_REF = 2000101.0d0 But I have faced this error SET_2DFLD - current model time exceeds ending value for variable: lwrad_down TDAYS = 662215.4398 ...

javadi

This is ok. I have checked it.I am surprised as in the previous flags the model ran without any problems. But, with this .h we have this problem cloud_time:long_name = "cloud fraction time" ; cloud_time:units = "days since 2000-01-01 00:00:00" ; cloud_time:calendar = "gregorian" ; float cloud(cloud ...

javadi

I use ecmwf2roms Matlab manuscript to make forcing. I use DSTART = 0.0d0 ! days TIDE_START = 0.0d0 ! days TIME_REF = 20130101.5d0 ! yyyymmdd.dd I do not know but we have different times for each variable. GET_2DFLD - solar shortwave radiation flux, 2013-01-31 13:30:00.00 (Grid=01, Rec=245, Index=2 ...

javadi

No. I think that this is ok. I have downloaded from ECMWF and merge months. Like this, is it ok? maybe, the time is a problem? NFFILES == 13 ! number of unique forcing files FRCNAME == /media/java/G/2013/gom_cloud_era1.nc | /media/java/G/2014/gom_cloud_era2.nc | /media/java/G/2015/gom_cloud_era_3.nc ...

javadi

I have defined #defined MIX_GEO_TS and compile the model successfully. The model works correctly for one month but at the end of each month, we have faced an error. For example, when the model starts in January, the model has faced problems on 31 January. While I started in February, the model has ...

javadi

Dear Kate

Thanks!

When omitting #define MIX_S_TS, I have faced this error :

18 | T3DMIX: must define one of MIX_S_TS, MIX_GEO_TS, MIX_ISO_TS
| 1
Error: Unclassifiable statement at (1)
t3dmix.f90:16:20:

javadi

Dear Kate Thanks for your response. I try to model the temperature, salinity, density, and water circulation(eddies) of my domain, not biology work . Based on your response, I have reviewed your application at Github(Apps-master). I have changed my upwelling.h file. here is the new one. Is it ok ...

javadi

I used river runoff when activating #undef ANA_PSOURCE. Also, I think that the model has a problem in flux too as the water is a little cooler than observation. I think that the model has a problem with vertical forcing and fluxes.

javadi

Dear Users, I ran the model for 6 years with this .h. The result does not show proper vertical temperature and salinity.Here I compare with observation. I try to model the temperature, salinity, and circulation of my domain. How to solve this problem? Please help me with your interesting comments ...

javadi

Hi all, I have a question. I ran the ROMS model successfully for one year with all input files. I made one year forcing and I decided to repeat for four years. In other words, after finishing one year: how the model uses this forcing file again for the second, third…? I use rst file as an initial c ...

javadi

Dear Kate

Thanks for your response!

My domain is the Caspian Sea. What is your suggestion to add ice to my model? What is the easiest way to add ice to my model?

I try to understand ice coverage and ice fraction and the effects of ice on the circulation of my area.

javadi

Dear users, I have prepared all grides, my initial, forcing, river runoff … and run the model when using the trunk for one year. However, in some areas of my domain, we have ice with a thickness of 0.5 m between December to February. I have reviewed this form and I understood that I should use a ...

javadi

Hi

Thanks(!)

I have selected tempearure=5 and salinity=1 PSU for the model. I think that your points can be true as the maximum freshwater is recorded in May which leads to blow-up. However, I have divided this freshwater into three cells too. The problem is not solved.

javadi

Thanks for your responses. I have repeated many times. Unfortunately, the model has faced blow up after 6 months. I have selected i=97 and j=101 for my model. The blow-up is (096,101,01). For this reason, I think that the runoff can lead to blow up. This point has selected as for example i=96, the ...

javadi

Dear users I have made my river run-off NC file based on this link https://www.myroms.org/wiki/River_Runoff. I think that my river position is correct as the model show the inflow when running the riverplum1 test case. However, after 6 months the model has faced blowup because of Maxdensity. What is ...

javadi

Thanks(!) I do not understand what you mean ''which field is blowing up''. Initially, I guessed that this point is related to the masking process as I have edited grids(masking water or land when making a grid.). For this reason, I have checked again this point. This grid is land and I have not ...

javadi

Dear users I have applied all files including grid, initial, forcing, rivers... to my model. After that, the model began to run. The model ran for 26 days without any problem but after that, I have faced a blow-up. Based on my previous experiences and this forum, I have done some actions: 1) I ...

javadi

Dear Kate Thanks! No. I think that I should explain more about my problem. I am a little confused although I have reviewed manual. I have changed my rupwelling.in as I have one open boundary in my southern part. After that I have faced this error there is not roms_bry.nc. Then, I have made roms_bry ...

javadi

My region is a part of the sea. The northern part is coast similar to West and south. In the southern part, we have an open boundary with a depth between 100-120 m. Here is my upwelling.in and upwelling.h. I have not changed anything. Which open boundary is suitable for me? How should I define the ...

javadi

Dear users, I try to run a upwelling_test for my domain. I have applied my grid, initial, climatology, and forcing the file to the model. The model runs without any problem. In my test case, we have one open boundary in the southern part. I think that I should define the sponge layer for the ...

javadi

Dear Robertson

Yes. Your opinion is completely right.

Thank you so much for your help

javadi

Dear Robertson When it comes to compiling without root, I have faced some permission errors. However, when I run ‘touch junk’, a file junk is made in my output directory without errors. The error when compiling without root: svn: E155007: '/home/javadbabagoli/Documents/roms' is not a working copy sv ...

javadi

Dear Robertson Thank you o much for your response. I usually use the trunk to run the model and I ran my test case without any problem. However, I try to add ice in my model; thus, I add kate branch to my system. I download from git when using clone. After that, I have faced this error. I use the ...

javadi

Dear users I have faced an error when running upwelling_test case which is related to Netcdf. Here are the details of this error. Please help me with your interesting comments. Wind-Driven Upwelling/Downwelling over a Periodic Channel Operating system : Linux CPU/hardware : x86_64 Compiler system ...

javadi

Dear kate

I do not change anything in the code only varinfo.dat. I add some variables at the end of the file when your greb .... does not work in my system. Here is my varinfo.dat. I try to run upwelling test case. It may be the source of the problem?

javadi

Dear kate I add some variables in vainfo.dat thanks to your help. I solve these errors but I have faced a new one. It is not related to my system as I have checked with the trunk. 0 0001-01-01 00:00:00.00 0.000000E+00 6.585677E+02 6.585677E+02 3.884376E+11 (00,00,00) 0.000000E+00 0.000000E+00 0 ...

javadi

Dear Kate

Thank you so much for your help

I understand how to add it in my varinfo.dat. I solve the previous error when looking your comments in this forum. Now I have faced the new error about id3dPVI
READ_PHYPAR - Variable index not yet loaded, id3dPVI

I do not find anything about this.

javadi

Dear kate Yes. Thanks for your help. I compile the upwelling test case using your branch both single and parallel. However, I have faced a new error in running time. Here is the detail of that. In the trunk, there is not any problem Model Input Parameters: ROMS/TOMS version 3.9 Sunday - January 24 ...

javadi

Dear Kate No. My problem is not related to the version of my netcdf and gfortran. Based on your suggestion, I have checked with the trunk and the model compile without any problem both in single and parallel. Thus, it shows that my packages install correctly. I guess that my problem is installing ...

javadi

Dear kate
thank you for your response. where can download compatible NetCDF and version? I use gfortran to compile.

javadi

Dear users

I have problems in compiling as follow:

svn: E155007: '/home/jav/Documents/roms' is not a working copy

/import/home/kshedstrom/bin/nf-config: Command not found

please help me with your interesting comments

javadi

Dear users I have run my test case without any problem when applying ROMS(trunk). Now I run the kate codes(https://github.com/kshedstrom/roms) but I have the following error. It is unusual for me as my Linux environment is the same for ROMS and roms_kate_svn. I have reviewed this forum and WIKIROMS ...

javadi

Dear kate I have run without any problem ROMS with the trunk. Based on your suggestion, I have installed your branch in my system BUT I have faced this error when compiling the model. svn: E155007: '/home/jav/kate/roms' is not a working copy svn: E155007: '/home/jav/kate/roms' is not a working copy ...

javadi

You mean I should copy and paste it from the site(!). is it correct as we install roms from command on the terminal? I have downloaded but there are many files one of them is Sea ice. How to use it?

javadi

Dear kate
Thank you for your response
How can I add your branch in my Roms through svn ? Then How can I add the files in my test case?i I am a little confused.
Please explain more about it

javadi

Dear users I try to run ROMS in my domain. I run ROMS with rivers, fluxes, and also the initial conditions in the upwelling test case. Now I try to add ice into the model but I do not understand how to do it although I review this forum. In this case, I have some questions. Firstly, ROMS have this ...

javadi

Dear kate and wikin Thank you so much for your suggestions. ROMS is sensitive to order in arguments and it was my problem. Now my model accepts my river file. I use the new version of ROMS(Rutgers). Now my model runs with my river files when I try to run the river_plume1 test case. However, I have ...

javadi

Dear Kate

I have changed my river_runoff.nc based on your suggestion. Here is the new one. However, I have faced this error which is related to river variable. Although I have input river, the model cannot find this variable. Please note that my ROMs is a new version.

javadi

Dear Kate Thank you so much for your response. Based on your comments, I conclude that I should make vshape variable. It should be the number of my vertical level. The sum of all values should be one. Here the river_transport variable is discharge(?) and v_shape divides river_transport into ...

javadi

Dear users

I have three rivers in my domain. I have made my river runoff file but I have faced this error. I am not sure about my input file as I do not define v-shape; is it essential? Here are my runoff.nc and this error.

please help me with your interesting comments

javadi

Dear wikin

Thank you so much for your response. Is there any Matlab script to make runoff?

javadi

Dear users, Happy new year. ................................ I have three rivers in my domain. For this reason, I have successfully tested the riverplume1 example in Roms when using Ana_psource.h as my river plumes. After that, I try to test in a real case. I have made river files including Qbar ...

javadi

Dear users

I think that this problem originates from a bad understanding of OA and Zout. I could not correctly define zout; hence, I have face this error. I need some information about Zout.

javadi

Dear users I try to make initial conditions using c_initial and d_inital which are based on MATLAB. Although I understand some parts of the codes, I have faced some problems. Firstly, I cannot interpolate correctly after using this script. Secondly, although my temperature and salinity are not zero ...

javadi

It can lead to instability as we do not consider ice in the model(?)The water converts to ice when the temperature is approximately zero. The observation supports this; however, in the model, I do not consider ice. If this leads to instability, how can I add ice to the model?

javadi

Dear wikin

Yes. It is a real area. The values for temperature and salinity are correct based on my previous experiences.

javadi

Dear wikin

Thank you so much for your help.I have reformed my .h based on new conditions and I have turned ANA_VMIX off but the model is unstable too.

javadi

Dear Kate Thank you so much for your cooperation. I have applied your suggestion and I have changed my upwelling.h. However, I have faced this error: CHECKDEFS - only one vertical closure scheme is allowed. Found Error: 05 Line: 1255 Source: ROMS/Utility/inp_par.F Found Error: 05 Line: 112 Source ...

javadi

Dear kate
Thank you so much for your response. I have another question in this case. As there is not GLS or LMD flags in my upwelling.h, I should copy this flag in cppdefs and add this flag in my upwelling.h like this:
#define LMD_MIXING
Or should I another action.
Please help me in this case

javadi

Dear kate and hezi I try to run upwelling example and here is my upwelling.h file. I guess that the problem comes from my initial condition. Without initial condition, the model ran without any problem using ANA file. Surprisingly, when I apply my grid file, the minimum z grid and minimum volume are ...

javadi

Dear users, I have faced blowing up when running a model based on my domain. I have made my all grid, initial, and forcing files. In the upwelling example, I added my input files and I started running the model. The model has faced blowing up after 5-time steps without any changing in upwelling.in ...

javadi

Dear users, I have made my grid and initial condition files. The model runs after applying both files. However, after applying the forcing files, I have faced the following error. To make forcing files, I have used decmwf2roms Matlab manuscript and my forcing files(13 files) are for one month from 1 ...

javadi

Dear users I have faced some problems when using parallel mode. I review all comments in this forum but I am not able to solve it. My system cannot find the directory; thus, I am not able to compile. I use build_roms.sh and I know that I should add some paths in my_build_paths.sh. I find the path ...

javadi

Dear Kate, To be honest, I am confused a little. Can you explain more. Should I use make clear in build_roms.sh? If yes,please address in this section of file. Excuse me. #-------------------------------------------------------------------------- # Compile ...

javadi

Dear KATE, Thank you so much for your cooperation. Yes, I used my own grid file and Roms run without any problem. Then, I made my initial file and then I changed ANA_INI to undef in my upwelling.in as I tried to use my own ini_file which is named roms_ini.nc. After that, I compiled again and then ...

javadi

Dear users, I try to run ROMS for my test case. I have faced an error when applying my initial file. Firstly, I have changed define Ana_grid in my ocean.in to undef. Then I put my grid file. In this step, fortunately roms accept my file and the model run without any problem. However, when I use my ...

javadi

Thank you so much Kate.
Yes there is one line.
To be honest, I do not understand exactly how to turn on flag. Can you correct my upwelling.h? I try to use smflux. Now i should use my own grid named roms_grd.nc and forcing files. How active them in upwelling.h?
please help me

javadi

HI,
when it comes to your question, Here is my file. You mean that the error comes from ....h file?

javadi

Hi users, I have faced a FORTRAN runtime error in running time. Firstly, I made the grid for my region. The model run without any problem when applying my real grid. Then, I created my forcing file using d_ecmwf2roms script. I have defined 13 files in upwelling.in file and active smflux in upwelling ...

javadi

Hi users I made my grid file for my region. Then I made forcing using d_ecmwf2roms. Now I have witnessed gom_cloud_era, gom_lwrad_era, gom_shflux_era ,…..files. My question is that what is the next step to make forcing file for Roms? Can Roms accept gom files as these files are based on long and l ...

javadi

Dear users Thank you so much for your help. I have faced an error when using d_ecmwf2roms.m . I downloaded data from ECMWF and then I tried to use d_ecmwf2roms.m. I ran this code when checking line by line to detect error of this code in my system. The first time I have faced an error in line 396 ...

javadi

Dear wikin Thank you so much for your help. I learned many things and I did some steps with new strategies. I download data from ECMWF and then I tried to use d_ecmwf2roms.m. I ran this code when checking line by line to detect error of this code in my system. The first time I have faced an error in ...

javadi

Hello users, Thank you so much for your help. It means that I should download all of files http://apps.ecmwf.int/datasets/data/interim_full_daily/ and then put all of them in one folder. Then I should define direction to this folder. However, there are some questions in this case: 1- The d ...

javadi

Dear users, I do not understand a point when making forcing. I try to use d_ecmwf2roms.m but I have a problem. Which of the following scenarios can be true? 1. When I run d_ecmwf2roms.m , these codes automatically download data from ECMWF and then convert to ROMS. Before that, I log in ECMWF. 2 ...

javadi

Hello
I am new user of ROMS and Matlab. The version of my matlab is 2014a. What is solution? you mean i install new version of matlab?

javadi

Hi, Thank you so much for your cooperation. Based on the previous comments, I do not understand how to do this Remove directories including netcdf functions added in matlab paths and call nc functions from the toolbox that is available in ROMS repository , please explain more. By the way,After ...

javadi

I am a user of ROMS. I have faced an error to make forcing when using d_ecmwf2roms.m. I tested many methods to fix the following error; however, I am not able to fix this error. Please help me to solve this error. Error using fclose Invalid file identifier. Use fopen to generate a valid file ...

Search found 118 matches