Ocean Modeling Discussion

dhtjtks

kate wrote:I'm sure it was, but things change in the ROMS world. 0.0 is a valid temperature so not a very good _FillValue.

Then, it's o.k. I have to follow the rule.
Thanks Kate!

Lee

dhtjtks

You can use the script matlab/grid/c_grid.m to create an empty ROMS grid NetCDF file: > status=c_grid(Lp, Mp, 'my_grid.nc', true) where Lp=L+1=Lm+2 and Mp=M+1+Mm+2 . That is, Lp and Mp are the total number of points in the I - and J -directions ( 0:L , 0:M ). Once that the file is created, you can ...

dhtjtks

kate wrote:For messing with NetCDF files I use NCL and/or Python. I was hoping someone else would help you with the Matlab if that's the road you want to take.

Thanks Kate. I think it is the time to learn Phthon. Thanks anyway.

Lee

dhtjtks

Maybe it's about your expectations. In serial mode, the size of the arrays is smaller than in parallel mode. It might even be unsafe to be setting Istr-2 in serial mode. In parallel mode you are guaranteed to have 2 or 3 halo points around all the edges on each tile. As for the 9.9e36, that's in ...

dhtjtks

Hi all, I tried to use 'ana_mask.h' to block some area of my domain, which have 100(i)*10(j) grids in interior rho point. When I looked at the results it didn't work like I wanted. Extra areas were blocked when I used parallel mode. However, if I used a single processor, it worked fine but also gave ...

dhtjtks

You can create a file with ncgen to get the correct variables and dimensions. How did you specify where you wanted your land mask? Did you have the full array of ones and zeros in the CDL file? Many, many scripting languages have a NetCDF interface which can be used to set the values or you can use ...

dhtjtks

You can create a file with ncgen to get the correct variables and dimensions. How did you specify where you wanted your land mask? Did you have the full array of ones and zeros in the CDL file? Many, many scripting languages have a NetCDF interface which can be used to set the values or you can use ...

dhtjtks

Hi all, I have a problem to simulate a simple tidal basin case. I made a simple rectangular domain (100 km in x-direction by 10 km in y-direction, and the grid spacing is 1 km by 1km). I put the tide at the eastern boundary. I tried both using 'ana_grid.h' and 'ocean_grd.nc', and both worked fine ...

dhtjtks

I would recompile with the USE_DEBUG flag turned on, then execute the oceanG file. It should be able to tell you line numbers then, perhaps even give array-out-of-bounds info if appropriate. Hi Kate, When I turned on USE_DEBUG flag, it worked fine both in serial and parallel mode. I don't know why ...

dhtjtks

Hi all, When I tried to run ROMS (example case, 'upwelling'), I had this error. I tried to use both netcdf-3.6.3 and netcdf-4.1.1, however, both cases have same errors. What does it mean? Thanks, dhtjtks #################################################################### # Alabama Supercomputer ...

dhtjtks

arango wrote:Use the build script I have addressed this issue countless times in this forum. See build script information in WikiROMS.

Thanks!

dhtjtks

Hi all, When I try to generate executable file using 'make' I got an error message like below. Actually, NETCDF and HDF is not located in '/usr/local', but I am using 'modules' to manage softwares. So, the library is locate in /opt/asn/etc/modulefiles/netcdf. Is this the problem of this? If so, how ...

dhtjtks

Hi everyone,

Do 'u_hadv' and 'u_vadv' terms refer 'uux + vuy' and 'wuz'? If so, can I separate 'uux + vuy' term?

Thanks,

Lee

dhtjtks

Hi, Angelolemos, I don't know which version of ROMS you are using, however, you don't need to have ocean_grd.nc and ocean_ini.nc in the upwelling case. Because 'grid' and 'initial condition' is defined 'analytically'. Let's see the 'upwelling.h'. #define ANA_GRID #define ANA_INITIAL #define ANA ...

dhtjtks

Hi Hiroshi: Thanks for your interest. I just followed Wilkin's solution and it works now. Thanks anyway. Lee Hi Lee, Can I make some questions? We need more information in detail; 1) What is version of your MATLAB? 2) Why don't you use (download) the newest netCDF Library: netcdf-4.0.1 so far? 3) I ...

dhtjtks

Hi Wilkin:

Thanks, it works now. I really appreciate.

Thanks again.

Lee

dhtjtks

Hi everyone, I am using MATLAB for reading netcdf file in my pc. Now I am trying to read netcdf file in the UNIX computer. So I set path in the matlab (in UNIX computer). However, it doesn't work. I downloaded 'netcdf-3.6.2' and 'mexcdf' from website and set path to these folders. It worked in PC ...

dhtjtks

kate wrote:Maybe if you show us what you are doing in ana_grid, we can tell you how it should be. Follow the examples given and see how they work in parallel (which they should).

Thanks Kate. I attach a pdf file. So you can check my problem. I will wait for your comment.

Jung

dhtjtks

Hello ROMS users, I have a problem with parallel running. If I generate evenly spaced grids, it works good. But if I generate un-evenly spaced grids, it doesn't work for parallel running with MPI. What I figured out my problem is this. If I have 10 by 10 grids, and the grids are evenly spaced in y ...

dhtjtks

Thanks Meiw! I figured it out yesterday. It works! Thanks again. dhtjtks Hi all, First question: I am using MATLAB to read output file, 'ocean_his.nc'. If I use small grid numbers such as 100 * 100 *10, it works o.k.. But if I used many grids (not that many I think) such as 500 * 500 * 10, I cannot ...

dhtjtks

Hi all, First question: I am using MATLAB to read output file, 'ocean_his.nc'. If I use small grid numbers such as 100 * 100 *10, it works o.k.. But if I used many grids (not that many I think) such as 500 * 500 * 10, I cannot load output file in MATLAB. If I used 'ncload', it shows like below. Does ...

dhtjtks

Hi everyone! I just want to put wind stress. It is same as wind-basin test case. So I turned off every other forcing terms except wind stress term. And I gave 1dyne/cm2 (0.1 Pa) in 'ana_smflux.h'. I needed to do unit conversion [dyne/cm2] to [m2/s2] unit. So I put 0.0001 m2/s2 which is equal to 1 ...

dhtjtks

Now I am using MATLAB for reading '*.nc' file but it is too slow. So I want to use fortran or C. I read several scripts about reading netcdf file with fortran and C but still don't know how to read. Does anyone know how to do it? Any suggestion will be helpful for me. Before abandoing Matlab, you ...

dhtjtks

Hi everyone!

Now I am using MATLAB for reading '*.nc' file but it is too slow. So I want to use fortran or C. I read several scripts about reading netcdf file with fortran and C but still don't know how to read. Does anyone know how to do it? Any suggestion will be helpful for me.

Thanks.

JW

dhtjtks

Hi all, I have a question about open boundary conditions. There are many options to use in ROMS. However, I still don't know which combination of open boundary conditions is good for some general cases. For example, let's think about an estuary which has a river inflow from west and other boundary ...

dhtjtks

Try: NORTH,SOUTH,EAST _TRADIATION (instead of gradient) and show us the results. By the way, thanks for so clearly providing plots of your interim results. This makes it much easier to help. John. Hi John, I did what you said but the result was same as before. Do you have any other suggestion ...

dhtjtks

Thanks John! I will try it and show you result.
Thanks again.

Jungwoo

dhtjtks

Hi all, I am trying to run simple case with 100*100*10 (x,y,z) grids. It has river inflow from west side and open boundary at north, east and south. I just want to let the flow go out through open boundary because I don't have any information at those boundaries. So I used 'radiation' condition at ...

dhtjtks

Thanks Paul. However, I cannot find anything wrong. I am using CDL file included in ROMS. But, it works when I retype whole sentences in unix computer.
Thanks again.

jungwoo

dhtjtks

Hi, all.

I have a problem with ncgen. When I try to make a *.nc file with 'ncgen -b *.cdl', I have an error message, cyntax error. I am using netcdf 3.6.2. I don't know how to fix it. Please help me.

Thanks,

Jungwoo

dhtjtks

First of all, I am sorry for late response. Now I just copy 'estuary_test' case and closed east boundary. So it looks like just closed_estuary. And what I want to do is that just 1m amplitude with M2 tide on west open boundary. In this case I set up 1m amplitude at 'ana_fsobc'. Usually, if I give ...

dhtjtks

No I didn't tried without SOLVE_3D. I think if I set amplitude equals 2m at western boundary, then it must have 2m amplitude. But it didn't have 2m amplitude. Result showed around 1m amplitude at western boundary and 2m at eastern wall. How can I fix it?

dhtjtks

I just want to simulate simple bay type problem. I closed north, east and south. I just open west boundary. And I gave 2m amplitude free surface elevation. But it looks wierd at west boundary. What is the problem and solution? This is the part of my cpp file. #define EASTERN_WALL #define SOUTHERN ...

dhtjtks

I am sorry! I just have a difficulty to find solution. That's why I wrote the message. I do not want to make other forum ^^.

dhtjtks

Hi all!!!
Let's begin to run ROMS together!!! Let's try to do simple example file first, make our own and share everthing together! Please email me if you want to join me.
email: jungwoo33@gmail.com

dhtjtks

Hi all!
I just began to learn ROMS but I can't find a manual. Is there a manual? Actually, I can run example cases included in the model. But I have a problem to make my own grid and run. How to make it? Where can I find some information?

Ocean Modeling Discussion

Search found 36 matches

Re: Question for 'ana_mask.h'

Re: Problem with 'ocean_grd.nc'

Re: Problem with 'ocean_grd.nc'

Re: Question for 'ana_mask.h'

Question for 'ana_mask.h'

Re: Problem with 'ocean_grd.nc'

Re: Problem with 'ocean_grd.nc'

Problem with 'ocean_grd.nc'

Re: forrtl: severe (174): SIGSEGV, segmentation fault occurr

forrtl: severe (174): SIGSEGV, segmentation fault occurred

Re: make & netcdf problem

make & netcdf problem

The meaning of terms in Diagnostic file

Re: Input files of Downwelling/Upwelling test case

Re: reading netcdf file in matlab

Re: reading netcdf file in matlab

reading netcdf file in matlab

Re: Grid generating failure with parallel running

Grid generating failure with parallel running

Re: 'ncload' problem in MATLAB

'ncload' problem in MATLAB

How to setup wind forcing?

Re: netcdf file reading with FORTRAN

netcdf file reading with FORTRAN

open boundary conditions

Re: Radiation Open Boundary Condition

Re: Radiation Open Boundary Condition

Radiation Open Boundary Condition

Thanks!

ncgen error!!!

Just simple tidal boundary problem!

Hmmm...

Let's begin to run ROMS together!!!

How to make grid?