Ocean Modeling Discussion

susonic

Hello

For the Horizontal and vertical mixing combination, please try U3 and C4, or both A4.

You set the temperature and salinity closed condition at Northern boundary. Try RadNud and set the proper value with Tnudg, M3nudg and OBCFAC.

susonic

Seems like your initial condition does not match with the boundary condition.
Could you plot your initial condition(u,v, temp, salt) here?

susonic

Please check out your lateral boundary conditions. I saw many 0 value at eastern and western boundaries. Also, when it blows up, could you plot horizontal u, v, temp, salt and upload here? It might not the major cause but I doubt about your unit at the humidity. Are you sure that it is relative ...

susonic

hello

1. The bathymetry was not smoothed. The value of rx0, rx1 are too high.
2. Please check out the unit of the Sea level pressure.
3. Please check out the humidity(specific, or relative) unit and match the unit in varinfo.dat

-JH

susonic

Hello, We are currently using the ROMS-SWAN coupled model. The COAWST version we're using is 3.8, and there are no issues when running ROMS or SWAN individually. However, when attempting the coupling, especially with #define WEC_VF, we are encountering the following message. For your information, we ...

susonic

I am using COAWST model (V3.8) with ROMS and SWAN couple. I tried to compile the RIP_CURRENT case and faced with this error. I had no problem compile with ROMS only. I only get this error when I compile with SWAN only or ROMS & SWAN. I am using ifort(intel18). Netcdf(4.6.1) Total log file from ...

susonic

#define TS_C4HADVECTION #define TS_U3VADVECTION needs to be written as #define TS_U3HADVECTION #define TS_C4VADVECTION Seems like you provide 6 hourly wind stress while the heat flux is coming from mean value. I would not recommend diurnal_srflux if this is the case. From what I experienced, it ...

susonic

Well if that is the case, First of all, I would plot every snapshot of temperature, salinity, current and sea surface height variable. And also I would run the model without surface forcing and lateral boundary condition(closed). If the model is stable without forcing(surface/lateral), then I would ...

susonic

Dear Sankar,

I just saw that your E-P file shows large value which is possibly filled with blank or NaN value.
Please check out fresh water flux file.

-Joonho

susonic

1. I see your nothern zeta,vbar and ubar value are 0. 2. Your rx0 value of the nesting model is too high. 3. It is not possible reason for blowing up, but I was wondering setting the boundary condition of the northern area as gradient for tracers? Also western is set as a gradient instead of close ...

susonic

I would suggest use Cha for zeta (no Che) Fla for 2d momentum (no Sch) Cla for all 3d value including Nonlinear model use also Same mixing/viscosity coefficient for NLM ADM. I saw zero bottom drag coefficient in your ocean.in. Is there a reason for 0 ? I don't know the reason but I suspect that the ...

susonic

Hello everyone,

In ROMS 4dvar, we can assimilate the observed ssha to the model. When assimilate the ssha, I was wondering how assimilated ssha influence the other variable such as temperature and salinity? I was trying to find the code regarding this issue, but I couldn't.

-Joonho

susonic

I think your set up for boundary condition is not common set up.

Why don't you give below boundary condition and see if the result change.

For free surface, Chapman
2d momentum(Ubar, Vbar), Flather
3d momentum(Uvel, Vvel), Clamp
MTKE, Gradient
Tvar (TS), Clamp

-JH

susonic

Thank you for your answer Dr.Arango. I've got more question. You mentioned that If the 1D variable arguement to mp_bcastf is singe precison, it will usee routine mpbcastf_1d in distribute.F. But my understanding looking at the code is if the single precision is defined, ROMS will use mp_bacastf_1dp ...

susonic

Thank you for your report. I also faced with same situation as you mentioned and was about to report.

susonic

First of all, check out the barotropic time step as Kate mentioned above. Secondly, I saw several weird option with your boundary condition. I saw gradient condition for temp and salt at the northern boundary. For West, I saw gradient condition for temp and salt(Not closed). I was wondering whether ...

susonic

Regarding this issue,
please give definition to your tide file with all available point when you build a tide file(say u,v,rho,psi dimension).

-JH

susonic

Would you show us the log file?

-JH

susonic

Kates's FAQ would be helpful to you.

For me, I would save every snapshot as possible as I can(here give NHIS==1), and then I plot every field

along all levels. You will find something.

-JH

susonic

For your information, you'd better run the PGE test with homogeneous(no stratification).

susonic

Hi,
The name in the code 'fac' means factor.

Therefore fac determines that the largest value of the viscosity(or diffusivity) in your model boundary.

For example, if you set 4.0_r8 and your viscosity is 10 m2/s(in ocean.in), you will have maximum 40 m2/s along your model boundary.

-JH

susonic

Oh! I saw wrong log file(I should have read third one). Sorry for the wrong post. Yes, like Kate mentioned, I also recommend you to try shorter time step. I also have a question about your boundary condition. As far as I know, the conventional boundary condition for the sea surface and barotropic ...

susonic

Hi,
You are not allowed to have zero value in your depth.
Check out more information on the below website.

https://myroms.org/forum/viewtopic.php? ... w=previous
viewtopic.php?f=23&t=2546

susonic

What is your minimum depth?
What if you run the model without WET_DRY?

susonic

My suggestions are
1.Try different horizontal advection scheme (say U3 instead of C4)
2.Try quadratic drag instead of linear drag(UV_QDRAG instead of UV_LDRAG)
3.Try different vertical turbulent scheme(say LMD or MY2.5 instead of GLS)
-JH

susonic

Minor stuff: I guess you are using 4 times daily surface forcing and I don't think you need to use diurnal srflux in your cpp. Wind data is critical for your mixed layer depth calculation as Duran mentioned above. BTW, is your SST reasonable compare with OBS? My suggestion is 1. Use Akima 4th order ...

susonic

Would you show us the log file? Right now I doubted the vertical turbulence model, K-epsilon which is pretty sensitive. My recommandation is that use different mixing (say, K-kl, k-w). Log file can give us more information about the causes of blowing up.
-JH

susonic

Check out the cyclone strength(Zonal and Meridonal wind) from WRF output and see if it's reasonable. You are giving very large bottom drag coefficient now. I was wondering whether this large value influence the thermocline depth. Try different vertical turbulent closure(say MY2.5, GLS) What is the ...

susonic

How did you give the surface forcing time interval(here, wind will be the critical factor)?

What are the horizontal and vertical resolutions?

Which scheme did you use for the vertical turbulent closure? (say, MY2.5, LMD, GLS, etc)

susonic

1. What if you run the model without tide? Does it still blow up?

2. What if you reduce the time step? Say, half from current time step.

-JH

susonic

I checked out the log file you attached several times but it is difficult to diagnose the causes of NaN value.

By the way, your Haney(rfactor) and Beckman & Haidvogel number is quiet high compare with known value(say, rfactor < 0.2) which means you did not smooth the bottom topography much.

susonic

Also, I noticed that, when the model only includes the SLA data or the SST data in the observation file, the model rejected most of them, but if the model includes both the SLA data and the SST data, all the observation data can be accepted by the model. For above situation, try to use super ...

susonic

I've never done SANITY_CHECK before. But I compiled with SANITY_CHECK and faced with same error message as you noted above(with gfortran). Maybe code developer can helped this out to solve the problem. I saw the log file you attached and found TLM increments are zero even though ROMS continuing the ...

susonic

Hi,
Does your domain include equator?
If it's ok for you, I want to see your model domain.

-JH

susonic

Ok, Dr.Arango. I will retry that model again. Thank you for your reply.

-JH

susonic

Dear ROMS developers, I know that there are lots of to do list ahead of you, but I'd like to add one more wish list. I'm currently using a two way nesting in ROMS and this is very useful and powerful tool to see the regional phenomena with one time calculation. While running the model, I'm also ...

susonic

Dear Developers,

I have a simple question.
Can ROMS simulate bio-physical (like NPZD model) using two way coupling?

-JH

susonic

Dear developers and users, In nesting application, if I want to use triply nested refiment grid and one additional nested grid (just like below figure) what should be the NestLayers and GridsInLayer be like? I read users should set NestLayers = 3 GridsInLayer = 1 2 1 in the ROMS wiki. https://www ...

susonic

Kunal & Zhang,

Thank you for your reply.

I know that the Fennel type bio model has a Chlorophyll variable. But I just want to start the bio-physical model as simple as possible.
Thanks for the reference.

Best,

-JH

susonic

Hi Luke,

How can I get the matlab script?
Sould I contact to Willis or
Would you send the script to my e-mail or upload to this post?
ocean0510@gmail.com

-JH

susonic

I also agree with Luke's opinion.
I have been trying to find a way to remove the steric signal from SSH.
But since I'm a novice user of ROMS, I do not know how to do it.
Are there any recommendation or tips to do it?
Any answer would be appreciated.

JH

susonic

Dear ROMS Users,

I am trying to calculates chlolophyll concentration from the output of NPZD model.

I was wondering how can I calculate the chlolophyll concentration.

Would you show me the way how to calculate the chlolophyll?

Formulae or matlab mfile would be very helpful.

Thank you.

-JH

susonic

Did you run the model without tide?

Did you run the model without holland wind input and put constant or 0 wind?

-JH

susonic

Hi Dr.Arango, Thank you very much for your reply. Domain averaged volume can get through the log file from ROMS, right? I can't log in the computer for some reason, right now. I will plot domain volume with a time series as soon as I reach there. By the way, I also ran the model with same domain and ...

susonic

Yes. It looks reasonable kludge. Thank you for your post, Dr.Lanerolle. The original ROMS code in bulk_flux.F cff=(0.7859_r8+0.03477_r8*TairC(i))/ & & (1.0_r8+0.00412_r8*TairC(i)) e_sat=10.0_r8**cff ! saturation vapor pressure (hPa or mbar) vap_p=e_sat*RH ! water vapor pressure (hPa or mbar) cff2 ...

susonic

ROMS can get either specific or relative humidity as a forcing input. And ROMS has three ways to calculate longwave radiation. When we calculate net longwave radiation internally (with LONGWAVE defined) ROMS requests humidity(RH, 412 line in bulk_flux.F) and this time ROMS does not know wheather ...

susonic

Maybe my question was too naive and also I did not understand the theory well. So I change my question. In Pacific basin, I ran ROMS for decadal order. The forcing is derived from AR4 and AR5(IPCC). After I finished the calculation, I saw the zeta in ROMS was becoming larger as time went by. First ...

susonic

As far as I know, ROMS is a Bousnessq approximated model which the desnsity field is fixed.

But as Dr.Arango mentioned, the zeta in ROMS seems varying with respect to density variation.

How do I understand about the zeta in ROMS?

susonic

Can I see your domain and the area problem occuring? Is it far from the boundary?

Have you ever changed the RDRG2 value? Mine is (0.0025 and the many test cases in ROMS are 0.003)

-JH

susonic

Why don't you switch off the M2NUDGING?
How did you set the TNUD, M3NUD and OBCFACTOR in ocean.in file?

This one is little different issue but does it work LONGWAVE_OUT with undefined BULK_FLUX?
It may not harm the lateral boundary value but I'm just curious.

-JH

susonic

Dear code developer and users, Recently I was told that there is a branch of ROMS which can add steric effect. That model is developed by Song and it is not a open source code. These days the climate variability is becoming more important and I think that it is necessary to consider the steric ...

susonic

I sincerely appreciate your file.
It will be very helpful.
Thank you so much.

Best,
-JH

susonic

Hi,

Is there a toolbox to generate attached figure like domain?

susonic

1. What about the depth from the river source point to the river mouth?
Does it vary?
2. How about change the river starting point?(ie. masking above the curving point and starting from it)

susonic

Well, actually I experienced that kind of problem many moons ago. 1. At that time I was giving only salinity on and the problem was fixed when I give both temperature and salinity on. But in your case, you already tried that before right? 2. Is your model including tide? If so, what happen if you ...

susonic

What happen if you give sources at whole vertical layer?

susonic

I believe that you have used the climatological monthly windstress data.
Why don't you use daily mean or 6 hourly wind data?
I'm sure that you will get a better result.
-JH

susonic

Hi, Kate. Thanks for the reply.

My ROMS case is calculating the heat flux via the bulk formulae.

Attached figures are the yearly mean, surface and bottom level avergaged temperature variation for 41 years.

The model domain is North pacific basin (E98~W76, S20~N65).

susonic

Dear users, In the climatology run, either using Levitus T,S data or longterm averaged value from OGCM (SODA, ECCO) for the ocean initial state and COADS or the ECMWF mean data for the surface forcing, I found that the temperature sharply rises on first year in most cases. And also many literature ...

susonic

Hi Kate,

Thank you for your reply and the useful paper link.

I want to see the result of the differences between the regional (Pacific ocean) which used different model boundary condition and the Global ocean model that has no boundary condition.

-JH

susonic

Dear Users,

Are there any grid generator program which is useful for the global domain?

-JH

susonic

Did you fix the problem?

I've just noticed that there is no boundary condition for M3 in your application.

I would recommend you to check out the current right before the model blow up.

Cheers,

-JH

susonic

Dear users, I had run a model with using 4Dvar data assimilation. Temperature and salinity profile are used to assimilate. I ran 2 cases, January and August. Assimilated temperature on January seems pretty reasonable compare with obs data but the problem occurs on August. The temperature in August ...

susonic

Thank you for your detailed comment and adding an entry of WikiROMS.
That is kind of you to do so.
And I'm sure that users can get this fundamental information via the WikiROMS.

Regards,

-JH

susonic

Thnank you for your reply, Dr. Arango. Good clue.

So, we have to set non-zero value as a minimum depth in order not to make the vertical level thickness zero.

-JH

susonic

Dear ROMS users,

I was wondering why we set minimum depth in ROMS ?

susonic

Hi,

Did you check out the variables in varinfo.dat file whether they matches
with the variables in your climatology file?

susonic

Thank you for your reply, Kate. Yes, you are exatly right about what you said. The problem occurred what I put the number of forcingfile did not match what really exist. My mistake. Now I've learned more about how to find a solution from error. I appreciate your instruction Kate. Thank you very much ...

susonic

Thank you for your reply, Kate. 3007 !----------------------------------------------------------------------- 3008 ! Load I/O information into structure. 3009 !----------------------------------------------------------------------- In inp_par.f90. There is Load I/O section. 3044 ! Initialize and ...

susonic

Hi Kate, Thank you for your reply. With debugger, I got below message. Model Input Parameters: ROMS/TOMS version 3.5 Thursday - July 21, 2011 - 10:13:41 AM ----------------------------------------------------------------------------- 0: Subscript out of range for array s%files (inp_par.f90: 3063 ...

susonic

Dear ROMS users, I encountered this problem only with specific case. I ran same domain with three different dimension in new pc, Quad core, 16GB RAM, SUSE 11.2, 3.6.3 netcdf, PGI fortran, mpich2-1.2.1, svn 564. case 1 dimension 88*96*20 case 2 dimension 177*194*20 case 3 dimension 292*332*20 The ...

susonic

There are two more. One from UCLA and the other from AGRIF.

As far as I know, AGRIF can do two way nesting.

Check this out.

http://roms.mpl.ird.fr/

susonic

Then the monthly ERA-40 data is probably averaged value from 6 hourly data.

Thank you very much for your replies.

-JH

susonic

Hi back, Thank you for your note. I tried to search in google and got an answer like below. http://stommel.tamu.edu/~baum/ecmwf.html 2. The "Surface solar radiation" in "W m**-2 s" needs to be converted into a "shortwave radiation flux" in "W m**-2". shortwave radiation flux (W/m^2) = surface solar ...

susonic

Any suggestion, please?

susonic

Hi,

I'm trying to use ERA40 forcing. The solar radiation unit is w/m^2/s.

I need to convert this to w/m^2 which is unit for ROMS.

What should I do?

Sorry for asking basic question.

-JH

susonic

Are there any specific reasons why you want to use tide information from
large domain result?

-JH

susonic

Before you think about the wrong use of horizontal/vertical diffusion cpp option,
did you check out the value of surface forcing whether it's reasonable?

Sometimes we find some weird value in the forcing file due to mistake
from the building script.

susonic

I can't see the images.

Would you attach them?

susonic

I mixed up with your answer. I edited my answer.
So, both of your model(large and nested domain) are initialized with climatology T and S? And the nested domain blows up? Can I see figures of your horizontal and vertical
speed at the blowing time?

susonic

Are you doing one-way nesting ?

If so, is your initial file also from larger domain?

susonic

Thank you for your answer, Mr.Warner. I downloaded sea level pressure and surface pressure from NCEP reanalysis version 1. I plotted same time step of both dataset(see attached) The general features of the ocean looks about the same although there is a small wave like motion in the surface pressure ...

susonic

Thank you for your reply, Mr.Hadfield.

So, if I can get the sea level pressure, then I'd better use it.

But if I can not get the sea level pressure, I can also use the surface pressure (or surface air

pressure) as an alternative since it is expected showing negligible effect.

susonic

Dear Users,

There are two different surface pressure variable from atmospheric model such as NCEP.

The one is surface pressure and the other is sea level pressure.

I was wondering which one should I use when I use 'Pair' in the bulk formulae.

Regards,

-Peter

susonic

Thanks for your message and link, David.

I'll follow your suggestion.

Regads,

-JH

susonic

Thanks for your notes, Ivica and Rodrigo.

So, I've got two options.

1. Modify the step_floats.F and interp_floats.F. and run the ADJOINT model with using floats.
2. Use history file to trace back.

susonic

Dear ROMS users,

I was wondering whether ROMS can track back the orginal source place using adjoint model with floats?

-JH

susonic

Hi folks, I've encountered a problem when I use d_ssh_obs.m First of all, when I executed d_sst_obs.m, there was no any problem. But when I executed d_ssh_obs.m, it gave me the below error message. h has the following dimensions (input order): 1) eta_rho = 55 2) xi_rho = 56 h loaded into an array of ...

susonic

Hi, Dr. Arango. Thank you for your reply.

At first, I used original(wrong) ARPACK but later I used a right one which you are describing above (I wrote about that).

You are right about the mpif.h. It does not exist in the original ARPACK library.

Thank you for your note.

-JH

susonic

I'm back to here to report the solution which I was stuck in before. I was trying to run the example, WC13 from recent assimilation Workshop. Thanks to the detailed explanation, I was able to run the I4DVAR. To run the I4DVAR, one need to compile the ARPACK tools before they run the model. I ...

susonic

Hi all, I have encountered the compiling error when I tried to use the R4DVAR. I used build.bash and got below result /opt/mpich2/bin/mpif90 -O3 -tp core2-64 /home/ocean/model/trunk32/project/WC13/R4DVAR/Build/esmf_roms.o /home/ocean/model/trunk32/project/WC13/R4DVAR/Build/master.o /home/ocean/model ...

susonic

Hi Mr.Ban, Did you fix the problem ? I had faced with similar problem which you mentioned above. Recently, I tried with updated ROMS svn 514 and it shows better performance than previous ROMS version. Dr.Arango amended some bugs with parrallel and I believe that the amending brought better result ...

susonic

I had faced with similar situation before. I gave only salinity at that time. It was solved with both temperature and salinity on with my river forcing. See this disscussion, https://www.myroms.org/forum/viewtopic.php?f=31&t=867&hilit=susonic Now I don't have a good solution but I doubt whether you ...

susonic

Hi,

Did you give salinity only to the river forcing file or both temperature and salinity?

-Peter

susonic

Dear users, I was able to see the Dr. Arango's presentation at the latest ROMS workshop. In that ppt file, there is a slide which compares about the bathymetry smoothing methods. And I was wondering whether users can get the linear programing method for the bathymetry smoothing through downloading ...

susonic

Dear ROMS Users,

I have a simple question.
Does ROMS calculate sea level variation with respect to the temperature variation?
If it does, would you point out where the code is?

Thanks in advance,

-Peter

susonic

Dear Users,
What should I do if I want to write the streamfunction?

Thanks in advance

-Joonho

susonic

Dear Users, I have some problem with plotting. The ncl was successfully complied but I got an error message when I run the idt gmeta. The error message is, 'reading stdout of translator: Resource temporarily unavailable. Translator aborted' Have you ever seen this message? and if so, what would be ...

susonic

Have you ever tried to run the model without tidal forcing under same condition?

susonic

Hi, Kate. Thanks for the reply. My rst file size is only 504mb when it stops. I just came to know that it's concerned with vertical mixing scheme. When I ran the model with LMD and PERFECT_RESTART activated, it couldn't save more than 3 step. However, when I define different vertical mixing scheme ...

susonic

Hi all, I found possible bug with perfect restart option in ROMS. without perfect restart option, we can write restart nc file continusouly(say up to 5th or 10th). But with defining perfect restart option, ROMS cannot write more than 3 step in rst nc file. The error message is 'WRT_RST - error while ...

susonic

CHECKDEFS - no bottom stress formulation is selected. See above error message. You need to give the formulation which caculates the bottom stress. In ROMS, we have logarithmic, Linear and Quadratic bottom friction formulation. Those are UV_LOGDRAG, UV_LDRAG and UV_QDAG respectively in cppdefs.h ...

susonic

Thanks Shchepetkin and Kate, Storing 3D data in restart file can be an alternative solution but it is only a snapshot(not a average). And we don't have split option for the restart file as the history or average file does. So it will store all variables in one file which will give us one huge file ...

susonic

Hi Mathieu, Thanks for the reply. Sorry. I might give you a wrong idea. What I want to do is I want to store 2d variable(zeta, ubar,vbar) frequently(say hourly or daily) without 3d variable(temperature, salinity, u, v) on history file. And I want to store 3d vaiable with longer time step(say 10day ...

susonic

Dear ROMS users, I was told that ROMS is using a 4DVAR data assimilation scheme which takes 20 times longer than normal calculation. But I think that 3DVAR is also powerful and less cost effective. My question is Are there any plans that making ROMS also be capable of doing 3DVAR data assimilation ...

susonic

Hi Mr.Liu,

Thank you for your reply.

What do I need to do if I want to consider the sea ice effects?

Do you mean that I should turn on the ICESHELF cpp option if I want to consider sea ice effects?

Thanks in advance

-Peter

susonic

Hi all, So far in ROMS we can store the variables on history file and average file with turning on/off the switches in ocean.in. And both output files are operated by only one switch. I think it would be better if we have each switches on each file(history and average respectively) so that we can ...

susonic

Dear ROMS users, I've got a same problem which Mr.Lanerolle had mentioned above. My domain is including the Okhotsk basin which is cold area and the air temperature goes down below -10 degree celcius in winter(especially on February). And then I got a unphysical sea surface temperature value which ...

susonic

Check it out below address and download the first ppt.

http://www.myroms.org/index.php?page=ROMS_2005_agenda

-Peter

susonic

Hi, yang
please check out below address and download nobuhitomori's pdf file.
There is a boundary code information in that file.

viewtopic.php?f=17&t=930&p=4200#p4200

Hope that it helps to you.

-Peter

susonic

I believe that you'd better use #define UV_PSOURCE #define TS_PSOURCE or #define Q_PSOURCE #define TS_PSOURCE Detailed explanation is given in Warner's presentation in 2005 ROMS workshop. You can get a ppt file at below address. http://www.myroms.org/index.php?page=Workshops Hope it helps to you ...

susonic

Ah... I should have looked at the code carefully.
Thank you for the advices, prof. Wilkin.

susonic

Thank you prof. Wilkin. You're very kind. Here is another question. The explanation which you linked to me is related with lmd mixing and Jerlov water type. Then, what if I want to use the other vertical mixing. Say, GLS vertical mixing scheme. My model result is telling me that different water type ...

susonic

Hi, Would you give me more information of the Jerlov water type. The water type 1 is representing the clearest water and 5 corresponding to the shortest penetration of solar radiation. This is all I know. Am I understanding right? Then, what is the criteria determining the water type? Any ...

susonic

First of all, I doubt the cpp option. In cpp they are warning try not to use 'TS_SADVECTION' unless it's a shallow and high vertical resolution application. If it's not a solution, then what about changing the ROMS version? Say, go back to svn 390.

susonic

Hi everyone,

It sounds basic stuff, but I was wondering which variable should I
use if I want to plot surface elevation gradient force and
baroclinic pressure gradient force term in the diagnostic file?

Thanks in advance

-Peter

susonic

Thanks for your quik reply, Dr. Warner. In your previous article which is right above my question, you mentioned this but i recommend to use the values in the ocean*.in files. and you just wrote that From the paper, the gls_cmuo value should be 0.5544 when used with the KC stability function for any ...

susonic

Hi Dr. Warner, I've got some questions about various GLS vertical mixing schemes. According to Warner et al(2005) on table no.2, 1. if I use k-epsilon with the Kantha and Clayson stability function, is nessasary to change the GLS_CMU0 value from 0.5477(ocean.in) to 0.5544(from the paper)? 2. What ...

susonic

It took a while for me to fix the high/low temperature with point source erorr. Now, I think I don't get the weird value anymore. The solution is turn on both the temperature and salinity flag. I had been doing all the experiments with only salinity value defined because I didn't know the ...

susonic

Dear ROMS users, I'm driving a regional model with a river source included. I want to use a passive tracer in order to visualize more about the distribution of the river plume. At first, I was able to add passive tracer in my initail and also boundary condition file. The model ran well and the ...

susonic

Hi ggerbi,

Thank you for your suggestion.
I'm running the model with time reduced(one third of the original time step). And now
it's woring without any problem.
It has still a little bit unstable speed near the boundary, though.

Again, thank you so much for your advice.

Best,

-Peter

susonic

Hi Dr.Warner, I already tested my application with using LMD and MY 2.5 vertical scheme. It ran well without any problem. And I have just tested k-omega and gen vertical mixing and it also ran well without any problem. For the k-epsilon, my final history file gives me that there is an increasing ...

susonic

Dear all,

I wanted to try to use the passive tracer in river forcing as mentioned above.
Would anyone help me out and give specific informations to me about
'add passive tracer in river forcing' ?
Thanks in advance,

-Peter

susonic

Hi, I'm trying to use GLS vertical mixing with k-epsilon. First of all, I ran under Kantha clayson with k-epsilon. It ran until 9th model day and blew up. Then I changed stability function to CANUTO_A. The model ran until 20th model day. I also changed the viscosity value. But it didn't help me out ...

susonic

Hi Dr.Wilkins, I wanted to try to use the passive tracer. And I was able to add dye_01 in my initail condition and also boundary condition file. The model was run well. Now I want to run with passive tracer in my river forcing file. I have been searching out the answer for this question in ROMS ...

susonic

How often is your BC value given in your domain?
Have you ever changed the BC code from RADIATION to CLAMPED for M3 and T ?
From my experience, CLAMPED BC code might be more promising.

-Peter

susonic

Thank you for the answer, Schepet. That's valuable information for me.
And thank you, Gcreager. That's another important information. I'm going to try that one.
I really hope that this new opteron would get better.

Best,

-Peter

susonic

Thank you for the reply. I was able to find an information about PASCAL.
Yes like you said, it is not a fortran compiler.
Is there a way to improve the performance with this AMD based cpu?

Any suggestion would be appreciated.

Best,

-Peter

susonic

Hi all, Recently I had bought a new system which is based on AMD cpu. And I was wondering how fast it is. I had run under same condition such as same OS(SUSE 10.3), PGI compiler, same NETCDF version, same MPI, same SVN ROMS. But the new version of system did not show good performance. Therefore I ...

susonic

Hi, all
In set_data.F line 729, unnecessary parenthesis is included.
Check that one out please.

-Peter

susonic

It seems that the error is about the same error which I experienced a couple of months ago with
eariler svn in ROMS 3.x.
Did you update with current svn ?
Try to check out current ROMS svn.
Maybe the number with current svn is 347.

susonic

Thank you Kate and Warner. I'm running the model again under temperature "TRUE" and staion with river point. I'm also testing with same discharge but more river point source. We'll see what happens. I'll report the result soon. I took a few snapshot and attached. I just hope that this is what you ...

susonic

Thank you for your reply, crode and kate.
I'm going to try it again based on your information and will report the result here.

Best,

-Peter

susonic

where are these points on the grid and which way is 'xi'? My domain is rotated so the xi direction would be rotated anticlockwise 40 degree from its original direction. And also if you see the above figure, named river_point, you can see the three blue points. Those are river flux points. Did you ...

susonic

What is the text on that line, and what is wrong? Well, it's not a big deal. I underlined the character in figure(see below). Even though I'm using a mpdata, I'm still getting weird value. This time, it has both hot water and cold water at the same time. Actually, I had run my model several time ...

susonic

Hi warner,

Thank you for your reply. I'm using default uv_advection scheme and u3h c4v for TS.
I attached again with png format.(sorrry and thanks)
I'm going to use the mpdata and will report the result here.

FYI, it has wrong spell in mpdata_diff.F on 19th line.

Regards,

susonic

Hi, all I'm trying to represent the river plume model which has a large discharge with uv_psource & ana_psource defined on cpp flag. The discharge varies 10000 upto 60000 ms-1. The size of delta x is approximately 11 km. At first, when I put the tidal forcing together with river, every thing was ...

susonic

Dera user,

I'm trying to install ROMS in NEC-SX system.
Is there a person who uses this system?
If so, would you show me a complie sctipt(sxf90 compiler)?

Please help me out.
Thanks in advance,

-Peter,

susonic

Hi all, While I was trying to add tidal forcing, I came to see that the ROMS is using both surface hieght and depth-averaged velocity boundary conditions in tidal forcing. I think that using both the height and velocity boundary conditions might lead into trouble due to the inconsistence between the ...

susonic

Hi, Mr.Arango. First of all, happy new year(sorry, it's kind of late greeting) Thank you so much for your clear explanation. It's a really specific information and also thank you for the link :D . It helped me out to understand better about those two method. Actually, I've been researching with ...

susonic

As far as I know in ROMS, if you want to do a prognostic running you need to run the model with undef xx_CLIMATOLOGY in contrast, if you want do a diagnostic running you need to run the model with define xx_CLIMATOLOGY in your header file. Am I right? Is there anyone who can tell me so that I can ...

susonic

Hi all It might sound basic thing but I can't find a definition of prognostic run and diagnostic run. Is there anyone who can define those two methods and explain what should I do if I want to run by those two separate methods with ROMS? For example, turn on specific cpp flag. Thank you so much! Best

susonic

The first one is visualization tools and the second one is for both visualization and preprocessing
tools http://eros.eas.gatech.edu/RNT/ here is the direct address.

Best,

susonic

viewtopic.php?t=916
viewtopic.php?t=1059
These two forums are helpful for you.

Regards,

susonic

Hi, there.
The tools which you are trying to use is from a different version of ROMS.
It's ROMS_Agrif from brest(France).
You'd better ask those questions
on http://roms.mpl.ird.fr/ which you downloaded the tools.

Good luck!

susonic

Hi all. First of all, thanks to Arango's advice, finally I made a forcing file for representing certain typhoon with using Netcdf file. Well, I was poor at programming. But I'm sure that I learnt a lot throughout this opportunity. Anyways, I put wind stress and pressure as a forcing for typhoon. It ...

susonic

Thank you so much for your reply, Mr.Wilkin and Mr.Johnluick. And yes I will try to make a netcdf file to force the tide. Before that I wanted to use the analytical method(I drove just one boundary-thanks to Johnluick) and it works! Now I will try to make a netcdf file. And I have a question. How ...

susonic

Hi all, I want to use tidal effect using analytical expression in my domain. Since my domain is too small to use large scale data such as TPXO6. So, I think it's better to use the ana_fsobc.h to make tadal effect. Hence I put this m2=0.757_r8 s2=0.337_r8 k1=0.245_r8 o1=0.177_r8 ttime=time(ng)+162.0 ...

susonic

Hi arango, I'm still analyzing Netcdf file to make a forcing file for typhoon. :cry: But I think I can make the forcing file soon. :roll: Right now I'm trying to simulate the effect of typhoon in the ocean such as 'cold wake'. Would you let me know what I need to define in my cppflag ? The problem ...

susonic

Thank you Kate. That's very nice of you to upload here so that anyone who
is wondering about the function can downdload and use it.
Well I need to analyze the script which you uploaded and I will apply in my case.
Thanks again. Have a great day!

Regards,
-Peter

susonic

Hi Kate,
I've got one more question to make sure I understood right or not.
So, does 'ng' mean each grid?

Best,
-Peter

susonic

Hi Arango.
Thank you so much for your early reply. And I'm sorry for the lack of understanding
and also made you surprise.

I'll take your tips and promise that I'll try to study more about ROMS.

Regards,
-Peter

susonic

Hi all, I'm trying to force both wind and pressure by using an analytical expression. It works but I got a weird value in my output file. So, is there anyone who can check out what the problem is? I put like this # elif defined NWPACIFIC OPEN(103,file='pres001.p',status='old') DO j=JstrR,JendR read ...

susonic

Hi Kate. Thank you so much for your information. It really helps me to understand more about passive tracer. Would you mind if I request your file which has passive_tracer or matlab xxx.m file which is used for making passive boundary or initial file? Or your ana_passive.h file would be also great ...

susonic

Hi all I'm trying to use inert passive tracer by ANA_PASSIVE. But I don't know for sure how to use it. Is there anyone who can give me an example file of ANA_PASSIVE so that I can reference it in my case. Right now I'm going to use those tracers as a diffusion of oil spill. ! Logical switches (TRUE ...

susonic

Thank you Kate. I really appreciate your reply.
Have a good weekend.

susonic

Hi all
I've been wondering that what does 'BOUNDARY(ng)%' mean?
Would you let me know that each of 'BOUNDARY' '(ng)' '%' meaning?

BOUNDARY(ng)%ubar_west(j)=my_flux/my_area
BOUNDARY(ng)%vbar_west(j)=0.0_r8

I appreciate any of your lessons.

Best,
-Peter

susonic

Hi! linzhenhua Thank you for the reply. I just solved the error a few minutes ago. :lol: It was the exact right error which you mentioned. I was so naive about Linux system therefore I didn't know which netcdf I should download. (I just downloaded netcdf-3.6.2_linux64.tar.gz) And I got a wrong one ...

susonic

Hi~ I've just installed Linux fedora8 and set up Intel Fortran 10.1.015_x86_64 together with netcdf 3.6.2_x96_64 I also set up roms-3.0 However whenever I try to compile I got below message. What seems to be the problem? What else I should install more about Linux system. I don't know much about ...

susonic

Thank you for your quik answer. I really appreciate that cecilia.
And I'll try~ Have a good day.

susonic

Hi everyone !
I've been trying to find out how to download ROMS 3.0 but I couldn't
have done yet. So, could you explain how/where to download ROMS 3.0 please?
And I don't know much about ROMS 3.0 so I'm wondering wheather I should pay for downloading ROMS 3.0 or not.

Thanks in advance

-Peter

susonic

Hi all I'd like to ask three questions. First, we can input river plume with salinity of zero in the model. I'm wondering if it is possible to dicribe the nutrient density of the river plume to see the influence of the river plume. Second, is there a way to calculate the diffusion of spilled oil ...

susonic

What a detail answer! Thanks for a Gooooooood explanation.

I'll apply the shorter time step and let you know how it goes.

Thanks again! And have a GOOD day~

susonic

Yes, It stopped while calculating a time step.Thanks Kate~ And I'll try to reduce the time step~ By the way, I am just wondering why should I change the time step? Is the time step change concerned with blow up ? I think the time step is just the time step. I don't understand.... Would you mind ...

susonic

Hi~ I'm using ROMS2.1 and I've got some problem when I put tidal effect. My case is North Pacific case(North-west Pacific region) and my cppdefs option is #define UV_ADV #define UV_COR #define UV_VIS2 #undef UV_VIS4 #define TS_DIF2 #undef TS_DIF4 #define MIX_GEO_TS #undef MIX_GEO_UV #define MIX_S_UV ...

susonic

I really appreciate your answer! Kate~

susonic

Hi All~I'm using ROMS 2.1 and I'm goning to ask 3 questions. the 1st question is I'm trying to make new grid file, but the bathymetry is different with each rtarget value. Is there an appropriate rtarget value when I'm making the grid? the 2nd one is when I'm making a grid for a very zoomed area ...

susonic

Thanks for your detailed reply. I really appreciate it !

susonic

ETOPO5 is 5minite grid... Then what does 5 minute grid mean?
(It might be a basic question but I don't know...

)
And is ETOPO2 a confirmed data? Does it reliable?
And the last but most important question, I'm using ROMS2.1, so can you tell me where the menual is?

Many thanks... Kate!

susonic

Hi Everyone! I'm a begginer with ROMS. And I wanna know about the topography which is used when we make grid. My friend said to me that use etopo5. He just said "You have to use etopo5 becasuse etopo5 has better depth data than etopo2" Are there any special difference between etopo2 and etopo5? I ...

Search found 170 matches