Ocean Modeling Discussion

pooran

Hi everybody,

I have downloaded the tidal constituents from OTPS. I got the negative value for phases in some constituents. I do not this problem when using TMD. Could you please how to fix this problem?

Thank you in advance for any suggestions

pooran

Hi everybody, I need to provide the tidal constituents for an enclosed gulf located in the Persian Gulf that its bathymetry is fine (dx= dy = 250 meters). could you please tell me which dataset is appropriate for this region? I think the TMD toolbox has a coarse resolution and it is not good for my ...

pooran

Hello everyone,
I have thousand of NetCDF files that must be merged over the time dimension,
I would be grateful if you could tell me the best way to do that.
I am looking forward to any idea.

Kind regards,
Pooran

pooran

Maybe you should start by telling us why you think you need a relaxation scheme? What problem are you trying to solve? Dear Kate, I am going to use the relaxation scheme in the open boundary, I have two options to define salinity and temperature in the open boundary 1- using hycom or Noaa data 2- u...

pooran

Dear all,
I am going to use relaxation scheme, but I don't know the meaning of it exactly, I was wondering if you could give me some necessary information, or if you could introduce me a good reference about that.
I would really appreciate your help in advance.

kind regards,
Pooran

pooran

kate wrote:Well, it depends. I would use pyroms, others would use Matlab. pyroms comes with a HYCOM interpolation example using scrip (which has to be compiled).

Many thanks for your help

pooran

Well, the steps you need to take might well depend on what you are planning to use the interpolated fields for. For ROMS, we typically use something like the scrip or ESMPy remapping tools for the horizontal remapping, but sometimes the land masks don't match and we need to do some lateral extrapol...

pooran

I don't understand what you are saying. What exactly is the problem? What tools are you using? Dear Kate, Actually, I don't want to use Roms model, I am going to find a good way to interpolation hycom data ( with resolution 1/18 degree) on my grid ( with resolution 2 min). I used cdo to do it, but ...

pooran

Hi everybody,
I interpolated hycom data on a new grid, but I have a problem and the new interpolation file does not coincide on a grid in the boundaries,
could you please explain to me how to solve this problem?
looking forward to any suggestions.

Best,
Pooran

pooran

Hi
I want to save simulation results in different history files in a run.
for example for :for 3_year simulation period,3 history files created that each one contains results of one year.
Many thanks for your help
Pooran

pooran

Hi
I want to save simulation results in different history files in a run.
for example for :for 3_year simulation period,3 history files created that each one contains results of one year.
Many thanks for your help
Pooran

pooran

many thanks for your kind attention

pooran

many thanks for your attention.
it is an excellent article that i find.
would you Please if you know any more case send me their information
i cant find any more using Google Search
cheers
pooran

pooran

Dear
I want to examine ROMS applications in Lake by reviewing of published studies. Can anyone introduce any published references on Lake simulation topic by ROMS?
Cheers
pooran.

pooran

dear
would you please tell me what is the benefit of FV(vortex force ) compare to Mellor08 formulation?
i know that Mellor (2008) that is more appropriate for sigma coordinates. what cause FV to be better that Mellor?
cheers
pooran

pooran

i am using roms in domain with one open boundary(in south) in tidal region for first time. before, i used rome in closed body without tide. i don't know how can add tidal force to roms :oops: by some posts it seems that tide is activated by SSH_TIDES and UV_TIDES options and requires a tide harmonic...

pooran

Dear There is ROMS current result over Caspian Sea in fig.1. As you can see current speed along red curve, which is consist with where there is great steep in bathymetry gradient (dash line in fig.2), is greater than other region (such as red curve's interior) . I did not know its result. Also I did...

pooran

Dear
is surface elevation equal to zeta?

pooran

Dear friends i faced with open boundaries for first time. My domain has 3 open boundaries (east,west and north) and one closed (south). I think that it needs T,S,U,v and zeta data for each boundary. i become grateful if you help me: 1- would you please introduced the best way (especially Matlab rout...

pooran

Can you provide incoming LONGWAVE instead of net

do you mean downwelling longwave radiation flux?

pooran

dear i turn on follow options: #define BULK_FLUXES # ifdef BULK_FLUXES # define EMINUSP # define EMINUSP_SSH # define ALBEDO # define CLOUDS # define COOL_SKIN !# define ANA_SRFLUX !# define LONGWAVE #define SHORTWAVE #endif and there are cloud fraction, surface air pressure (mb),air temperature at ...

pooran

many thanks for your nice reply 1- i am using net solar radiation from ECMWF_Interim database. I don't know if they applied the albedo correction to all incoming shortwave or not. Do you know about it? Do you mean that if they applied albedo correction, then we don't require AlBEDO option? 2- CLOUDS...

pooran

So if it acts as a nudging to that SST:
1-does it required to define TCLIMATOLOGY and TCLM_NUDGING or QCORECTION is enough?
2- will this nudging do in each step or base Nudging/relaxation time scale (TNUDG)? Or TNUDG will applied whenever TCLIMATOLOGY and TCLM_NUDGING are defined?

pooran

I am trying to use bulk_flux. My question is about ALBEDO and COOL_SKIN options. 1-please tell me that Since there is not ice cover in my domain, ALBEDO should not be defined? 2-what about COOl_SKIN and CLOUDS? When we should define them? Sorry but base on description in cppdefs.h i cant understand ...

pooran

Dear Kobashi i have studied mentioned article and will chose -32 as dqdsst (my studding area is located between 36N and 47N which has seasonal variation of dQdSSt). i become appreciate if you help me to understand more about it by answering follow question: 1-i found that using SST and dQdSST will h...

pooran

Dear kobashi many thanks for your reply. I have checked those databases before but they don not have enough data in Caspian sea region :oops: because of this i was computing data from different formula and functions. So, IS there any recommended constant value for dqdsst (especially in in Caspian se...

pooran

1 – does dqdsst's relation with sst and wins speed is directly and Reverse, respectively. 2- Pair (mb),Tair (c), qsea (kg/kg), air density (kg/m3), sst (c)( which 2 first one created by d_ecmwf2roms.m and third one by first post's equation and fix value (1.29) as a air density) have benn used for co...

pooran

would you please guide me that It is correct to use 2 meter dewpoint temperature (c) and surface pressure (mb) or 2 meter temperature (Tair_era.nc (c)) and mean seal level pressure (Pair_era.nc (mb)) as Ta and Pa for computing saturate specific humidity (in above formula)? Because Their calculated s...

pooran

dear i want to create dqdsst for first time by my own data. Base on different posts get_dqdsst.m is a good idea for this aim which needs sea surface temperature (Celsius), sea surface atmospheric temperature (Celsius), atmospheric density (kilogram meter-3), wind speed (meter s-1)and sea level speci...

pooran

Dear kate
thanks for your reply.
If it is possible please tell me what are the upper and lower bound for TNU2 and VISC2.
Cheers
pooran

pooran

My question is about horizontal diffusion and Viscosity coefficient. As you know TNU2 and TNU4 are same in all test (equal to zero) and VISC2 value has been changed in different tests.my grid cells are about 3300 m and changed to 31 m for child grid. How can recognize the appropriate values for thos...

pooran

It has been mentioned that it currently only supports one parent with one child grid. No grandchildren as yet.

is it true for COAWST updated code, too?

pooran

Sorry but I become grateful to know that for figure in second post is it true to set Ngrids, nestlayers and gridinlayer flags (in ocean.in file) like bellow? Ngrids= 4 nestlayers=4 gridsinlayer=1 2 3 4 Sorry but I can’t find any test to find my answer therefore I need your experience for my set up(i...

pooran

sorry I cant completely catch nestlayers and gridsinlayer. Would you please help me with second post's figure as an example: should it be set as: Ngrids= 4 nestlayers=4 gridsinlayer=1 2 3 4 2-if there is two child grids and one parent (in mentioned figure) then is it true to set: Ngrids= 3 nestlayer...

pooran

WOW, model is running with new wec_stokes.F

Sorry but I am worried to continue using COAWST recent version for next approach. For example nesting. I am worried that their sources need update, too. And it gives long time

What is your recommendation?

pooran

Dear warner sorry for sending post in roms forum again. I prefer to continue this subject here but for another one I will use COAWST forum. I used #undef WEC_VF #undef WDISS_WAVEMOD #undef ROLLER_RENIERS #undef UV_KIRBY #define SSW_BBL #define SSSW_CALC_ZNOT #undef UV_QDRAG model could work and past...

pooran

Dear I have tried another states which in all of them GRA chosen for all lateral WEC boundary condition (just defined options have been written and others are undef) a) #define WEC_VF and #define ROLLER_RENIERS b) #define WEC_VF and #define ROLLER_SVENDSEN c)#define WEC_VF , #define WAVE_MIXING d)#d...

pooran

I have tried
#define WEC_VF
#define WDISS_THORGUZA
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY

and

#define WEC_VF
#define WDISS_CHURTHOR
#undef WDISS_THORGUZA
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY

but same blow up appeared.

pooran

thay are in attached fig

pooran

exactly
it works with out any swan output

pooran

roms itself and swan itselp have ran without problem

pooran

Dear Kate and Warner coupled model has been ran for different scenarios: a) decreasing dt of ROMS from 90s to 1s, b) adding WEC_VF, WDISS_WAVEMOD,ROLLER_RENIERS,UV_KIRBY to header file step by step,c) changing of Wec boundary conditions types (Clo to Gra), d) using ana_initai instead of created init...

pooran

1 304 00:01:30 NaN NaN NaN NaN (145,323,26) 2.828234E-06 2.260028E-06 3.989284E-03 4.544180E-03 Blowing-up: Saving latest model state into RESTART file as i told there are just u and v for initial time step (equal to zero) in rst file (becau...

pooran

thanks for rapid reply. My domain in Caspian sea (closed basin)which i want to do nesting (4 times) in coupled model to investigate rip current near the coast. 1- this must be a different grids than we distribute yes, i used Lm == 194 ! Number of I-direction INTERIOR RHO-points Mm == 344 ! Number of...

pooran

Base on your suggestion to use new version of COAWST, I have updated my coupled test and then each model (ROMS and SWAN) ran ,separately without any problem . But their coupling faced with blow up very very soon, in the first step :cry: and therefore all data in rst and his files has their initial v...

pooran

1-Do you mean: for river we should turn of flags like this?

Code: Select all

LuvSrc == T ! horizontal momentum transport
LwSrc == F ! volume vertical influx
LtracerSrc == T T  ! temperature, salinity, inert

2-I can't completely caught what it “volume vertical influx”.
When should turn it on?

pooran

Dear I am updating my test to COAWST new version. 1- river file need #define UV_PSOURCE and #define TS_PSOURCE in old version. It seems that new version does not need them. Does it need another options? 2- we add river file's name in FRCNAME flag for old version. but it seems its names should be add...

pooran

Did you plunk a new SWAN version into the SWAN directory?

no
i did not know that just some test cases get up to the new version. I can run inlet-test

.
What is SWAN's version in new COAWST? is swan's current version able just to nest in one-way (as old versions) ?
Cheers
pooran

pooran

My system move to Ubuntu .After a while I can compile new version of COAWST by some changes: 1- changes in linux.gfortran.mk from LIBS := $(shell $(NC_CONFIG) –flibs) to LIBS := -L$(NETCDF_LIBDIR) -lnetcdf -lnetcdff 2-and in coawst.bash use: USE_MPI=on, USE_MPIF90=on, USE_NETCDF4=on and export COMPI...

pooran

In wikiroms you said:

Notice the duality of a ROMS nested grid: data donor in one contact region and data receiver in its conjugate contact region. The exchange of information is always two-way.

Is that confirmed even if using #define ONE-WAY option?
cheers
pooran

pooran

would you please tell me that in Multiple Refinement Sub-Class state , there is just one-way nesting? It meas: does coarser grid is just as a donor and finer grid as a receiver?
all the best
pooran

pooran

Sorry :oops: i would like to change my question 1-does it important that child's masking and bathymetry to be consistent with parent's masking and bathymetry, respectively, just at the boundaries points:? it means: is it possible to change masking and bathymetry of new refined grid in overlapping re...

pooran

Dear in my test, nesting will occur between near-shore and offshore. and i have Multiple Refinement Sub-Class grids. would you please tell me that 1-is it possible to define another minimum depth (hmin) for first Refinement Grid from parent grid? If yes, how? 2-is it possible to define another minim...

pooran

Dear with regard to matlab routines which create child grid for nesting (such as create_nested_grid.m) it is possible to define new masking and bathymetry. Would you please tell me when creating child grid it is possible to: 1-change number of layers in child grid compared to parent grid? 2-change m...

pooran

Dear kate many thanks for your attention you're right. dqdsst is discussed in forum and i have read almost all posts about it 8) . It can correct model's SST. But my question was about how i can provide it :?: I follow forum which has suggested:1- get_dqdsst.m routine and ECMWF database for providin...

pooran

dear friend
i would be so appreciate to receive any reply to two above questions
is DIGNISTIC equal to DIGNISTIC_UV and DIGNISTIC_TS?or DIGNOSTIC just define if no diagnostic file to be created or written (therefore is different from DIGNISTIC_UV and DIGNISTIC_TS)?
And about dqdsst?
Cheers
pooran

pooran

Dear i have seen that some tests have : #define DIGNISTIC and some other have both #define DIGNISTIC_UV and also #define DIGNISTIC_UV. Would you please tell me that #define DIGNISTIC is equal to #define DIGNISTIC_UV and #define DIGNISTIC_UV. Also ROMS has ana_dqdsst.h but COAWST does not have it. An...

pooran

Dear i expect that by increasing number of nodes, model run time decrease linearly. For example: by choosing NTilei=1 and NTilej=1 and then Ntilei=2 and Ntilej=2 , run time decrease to 1/4 times. But it decreases just 1/2 times :?: . Would you please tell me that if there are some things wrong in mp...

pooran

Before applying CDO command on force file, when i run model it gives following error: NETCDF_INQ_VAR - requested variable: lat sms_time not found in input file: /home/forcing/CS_sms_era.nc call from: get_varcoords.F after doing CDO command on force file, roms was able to read force file (which creat...

pooran

dear friends sorry to repeat my question but because i am a newbie with what ROMS need, i am worried that changed force file with mentioned header (in above post) have been damaged (after doing CDO command) and changed to incomplete format. Would you please guide me that its new header is complete a...

pooran

dear after doing cdo selvar command to change lon and lat's names in Dimensions 's part to x and y, force file's header change to: ncdisp('shflux.nc') Source: /home/forcing/shflux.nc Format: 64bit Global Attributes: CDI = 'Climate Data Interface version 1.5.9 (http://code.zmaw.de/projects/cdi)' Conv...

pooran

sorry
i hope to receive your answer. i cant find my answer not in manual nor in cppdefs.h.
cheers
pooran

pooran

Thanks. Are you mean that whenever CHARNOk and ZOH_HSIG options be defined in header file then CHARNOK_ALPHA and ZOS_HSIG_ALPHA will be used, respectively? Sorry but i cant fined any thing about and CRGBAN_CW and SZ_ALPHA?are they related to CHARNOk and ZOH_HSIG options , respectively? Sorry for con...

pooran

Dear warner thanks. With regard that article it seems that those lines are important when we have coupling with roms and swan in VF formulation (is it true?). But it is my fault that i keep those lines when i run just roms model. Would you please tell me that it cause error in result or because i ju...

pooran

Dear i have seen COAWST example. Just 2 of them have following lines (rip current and inlet) !Constants used in surface turbulent kinetic energy flux computation. CHARNOK_ALPHA == 1400.0d0 ! Charnok surface roughness ZOS_HSIG_ALPHA == 0.5d0 ! roughness from wave amplitude SZ_ALPHA == 0.02d0 ! roughn...

pooran

Dear friends i need to use your experience to correct my results. I used wind stress, shflux, swflux and swrad as input force files which provided from ECMWF-Interim. After running model for about 20 days, model 's currents speed and measurement did not have same trends :cry: that is while meansurme...

pooran

Dear i have used set_depth.m and then checked z values in bottom layer with amount of h which are in grid file. There were different. For example for one grid point which its maximum depth is grid file is 227.1217 m , computed z from surface to bottom are, -0.0927,-.2868,....,-10.1409,-44.2666,....,...

pooran

Dear As you know by using set_depth.m, number of z levels will be same as number of sigma layers. In the other hand each sigma layer has thickness in vertical direction. So please tell me that computed z values for each layer are corresponded to maximum or mean depth in each sigma level. Cheers pooran

pooran

dear friends
many many thanks
i will try them
cheers
pooran

pooran

Dear actually, my first post was about that . I want to know how i can compute each layer depth by ROMS output :?:. would you please tell me that the correct way to compute depth (by ROMS output) is to multiplying each layer's S_coor (or CS_curve :?: ) with h values which can be gotten from h matrix...

pooran

Sorry, but i have used rnt_setdepth.m (which located in RNT/matlib/rnt subdirectory) and depths.m (which is in ROMS/MATLAB/utility) for computing depth at rho points. Results were different :!: Actually, i my aim is to extract u and v values in specific depth.if you have another recommend i become a...

pooran

hi
How can compute each sigma level's depth at rho point? By multiplying rho point's depth by value of each level's S_coor Or CS_curve ?is there another way?
Thanks

Search found 72 matches