Ocean Modeling Discussion

Hello everyone
is some time I try to solve a problem.
I have an output of a meteorological model,
and I have interpolated from a grid of 500m grid spacing.
dimensions

meteo :
GLAT 115x50
GLON 115x50

Roms: the name and size lat and lon

lat_rho 576x480

lon_rho 576x480

I started in this way, but I ...

Hello everyone
I would be interested to know if there is a way to use the output of WRF running on my PC to initialize the ROMS ...
Generates WRF grib files but I did not understand how to manage them for roms ...
There are tools??
can someone help?


thanks a lot


Best regards

Hello everyone

I'm using lake_signel as a test of starting.
I created a file of tidal forcing.
on my PC works perfectly but does not work on the cluster.
this is the problem and .h files
...
...
...
31 106.50 212.50 195.50 208.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1 ...

oookkk!!!

I solved by NCO


I extracted lat_rho and lon_rho from file grid.nc.
then I copied forcing.nc.
1)ncks -C-v lat_rho, lon_rho grid.nc coo.nc
2)ncks coo.nc forcing.nc (this is not the best command)

of course I modified the file .cdl in this way

float Uwind(wind_time, eta_rho, xi_rho ...

hi
I have the same problem.
this is my netcdf initial "template":

netcdf roms_wind_frc {
dimensions:
wind_time = 49 ;
sms_time = 49 ;
eta_rho = 156 ;
xi_rho = 114 ;
variables:
float Uwind(wind_time, eta_rho, xi_rho) ;
Uwind:_FillValue = 1.e+35f ;
Uwind:units = "m s-1" ;
Uwind:long_name ...

first problem

I solved by setting consistently the number of cpus for each model , number of node and the number of tile partitions to use, remembering that "the product of NtileI and NtileJ
! must be EQUAL to the number of parallel nodes specified during
! execution with the "mprun" or "mpirun ...

ok...thanks for your advice
I changed the start time and now it works,
but I have a persistent error, this error occurs even with the test case ...
depends on what you know??


...
86397 13462 23:59:57 7.221366E-04 2.485874E+02 2.485881E+02 1.711143E+08
86398 13462 23:59:58 7.221378E-04 2.485874E ...

Hello everyone
some days I'm trying to work with 2 models coupled.
test cases (inlet_test , test_head) are executed with a good outcome.
There are problems when using a realistic case.
The following attached file and error.
I hope very much that can help me.



NTIMES == 1440
DT == 30.0d0
NDTFAST ...

1433 0 11:56:30 1.357294E-01 9.617819E+01 9.631391E+01 1.952436E+10
1434 0 11:57:00 1.352699E-01 9.617825E+01 9.631352E+01 1.952438E+10
1435 0 11:57:30 1.348071E-01 9.617832E+01 9.631313E+01 1.952439E+10
1436 0 11:58:00 1.343405E-01 9.617839E+01 9.631273E+01 1.952441E+10
1437 0 11:58:30 1 ...

what is the compiler that you used
and the flags for MCT??
mpif90 is not a good idea
at your suggestion I prefer to use the same compiler for all.
now I have used
for
mpich: ifort
MCT: mpif90
make: gfortran

clima@perico:~/a/ROMSanto> mpirun -np 16 ./oceanM ROMS/External/coupling_test_head.in
WARNING: Can't read mpd.hosts for list of hosts, start only on current
Coupled Input File name = ROMS/External/coupling_test_head.in

Model Coupling Parallel Threads:

Ocean Model MPI nodes: 000 - 007
Waves ...

okkk resolved
thanks

I've compiled so MCT

. / configure FC = mpif90

Currently I have a new (and known) problem (I think it is caused by gfortran)



cd Build; /opt/intel/ict/3.0/mpi/3.0/bin64/mpif90 -c -frepack-arrays -O3 -ffast-math -I/home/clima/a/ROMSanto/MCT/mct -I/home/clima/a/ROMSanto/MCT ...

hi jw

thanks for the fast
you think is a good idea to use gfortran???

I compiled mct but I received this error

Machine=x86_64
OS=Linux
Fortran Compiler is Intel ifort

Output Variables: {CC=cc} {ALLCFLAGS=-DFORTRAN_UNDERSCORE_ -DSYSLINUX -DCPRINTEL -O} {FPP=cpp} {FPPFLAGS=-P -C -N -traditional ...

Hello everyone

I tried to run intel_test on a cluster, with MPI and gfortran.
I compiled the libraries MCT (. / configure, make, with this Makefile.conf)
makefile.conf:


# COMPILER, LIBRARY, AND MACHINE MAKE VARIABLES

# FORTRAN COMPILER VARIABLES #

# FORTRAN COMPILER COMMAND
FC = mpif90 ...

i can do it with one core (in lake_signell test case , and in test_head received error write on),
but i can't with two core...is it possibile run it???



There are other ways that do not know??


thanks

Hello everybody

I can't use ROMS very well , but till now I have run in a correct way several Test-case bot idealized and realistic case .
Now I to work whit ROMS and SWAN coupled but for personal need I must set the work on my laptop first (intel dual core processor / Ubuntu 9.04 ).
At present I ...

I used this configuration

%--------------------------------------------------------------------------
dims = { 'eta_rho'; 'xi_rho'};
nc{ 'angle'} = ncdouble(dims);
nc{ 'angle'}(:) = angle; % angle of conterclock-wise grid rotation (degrees)
nc{ 'angle'}.description = 'angle of conterclock-wise ...

dat=load('boxgeog.dat');
xbathy=dat(:,1);
ybathy=dat(:,2);
zbathy=dat(:,3);
>> save bat.mat xbathy ybathy zbathy
where
xbathy = lat
ybathy = lon
zbathy = z

z = neg. in boxgeog.dat

I use Global Mapper.
impot my data
I export ascii file (first column lat second column lon third column z)
z ...

ok...
i solved (for previous error , not for clima.nc file)


rnt for president

I also use rtn but I have several problems.
for example when I create analytical grid:








% EDIT: set a name to recognize the grid
nameit='test'; % 8 km grid

% EDIT: set path where you want to store the grid.
grdfile=[nameit,'-grid.nc'];

dx=500; % m in X direction
dy=500; % m in Y ...

hi
I have multibeam data, xyz format (any format), I want to create the grid file.
what tools should I use?? (MATLAB please )

gracias

hi
What is the most simple and complete tools for MATLAB to create these files???
tutorials???

gracias

I can not extract a bathymetry.
I decided to follow this path: https: / / www.myroms.org / wiki / index.php / easygrid

following this path I extract the coastline: https: / / www.myroms.org / wiki / index.php / seagrid

the coastline is ok!


I have to extract the bathymetry
or used 3 methods
1 ...