Ocean Modeling Discussion

Below are the CPP options that I used for the 2004 Global Roms, but 1st let me suggest a perhaps easier alternative for your southern boundary. Why not extending your southern bdry to Antarctica and using EW_PERIODIC in the east-west direction for the ACC? that has been sucessfully done in the past ...

I found an old hard drive where I **might** have the original cpp.h file. I will post/send it later tonight. Your cpp.h has some non-standard choices. Are you using tidal forcing?

For a Pacific domain, 3-4 years should be enough to get the main features of the upper 400-500 m of the ocean. Could you upload some figures? but not here (i.e., ROMS Results), open a new topic in "ROMS Problems" instead. Guillermo.

This run was done more than 6 years ago. I don't keep the CPP file but you will be able to find some info in Auad (2008), GRL. I'll be happy to take a look at your results. Best, Guillermo.

BAN, you will need to do some coding but technically, I don't see why it wouldn't be possible. I envision 2 ways, a) one is by modifiying the code where Hernan dealt with the variable grid size and viscosity, and perhaps importing the bottom topography if that variable is not currently included in ...

running roms on a grid with e-w periodic boundary conditions and compilation done with mpich1 works fine (32-bit). when I compile it using mpich2 (and 64 bit, on another cluster), it generates negative salinities in two spots. I ran it on yet another cluster (64 bit) without problem. Anyone had this ...

Thank you Kate,

I haven't read anywhere that the source/sink needs to be against the wall/coastline. I need to place a source/sink in the open i.e., with neighboring grid cells being all wet ones. thx once more.

Guillermo

when defining UV_Psource, what parameter defines if it's a source or a sink? the sign of the vertically mass transport? positive for rivers, negative for sinks?

thx a lot.

Guillermo.

I don't have "Intel Mpi". Do I need it to compile ROMS, or the libraries contained in ROMS/LIBS/MCT should do the job?

Any other tips related to mpi compilation with the Intel compiler are more than welcome

thx

guiyo

anyone has encountered this compilation error before? I have up to revision 110 which is the latest as of this posting....thx. Guiyo. fortcom: Error: ad_t3dbc_im.f90, line 197: This is not a field name that is defined in the encompassing structure. [T_NORTH_CE] IF (BOUNDARY(ng)%t_north_Ce(i,k,itrc ...

after typing make (upw case), I get: fortcom: Error: mod_ncparam.f90, line 692: This name does not have a type, and must have an explicit type. [HEADER] Hfile=HEADER --------------^ fortcom: Error: mod_ncparam.f90, line 692: A CHARACTER data type is required in this context. [HEADER] Hfile=HEADER ...

thx again Kate. That got ride of that error message but now it complains about the following: $ make mv ./Compilers/MakeDepend ./Compilers/MakeDepend.orig mv: cannot stat `./Compilers/MakeDepend': No such file or directory make: *** [Compilers/MakeDepend] Error 1 When I used previous versions the ...

thx Kate, GNU version is okay I think: 3.80

I paste below line 100 and the next few ones (first one is line 100):

<<<<<<< .mine
<<<<<<<mine>>>>>>> .r41
>>>>>>> .r40

I just retrieved update 108.

when I type makefile (after configuring it and after configuring Linux-ifort.mk) I get:

$ make
makefile *** missing separator. Stop.

thx for any help

guiyo

I used a wall at 80N - depending what you are looking at, it might be ok. The midlatitude and tropical features (including climatic indexes such as PDO and NINOs) were very well reproduced in the Pacific for our cuasi global ROMS. I think that Hernan posted these results in one of the sections under ...

I got the following compilation error (the compilation is succesful if I undefine ICESHELF) for the following set of choices. #define ICE_MODEL #define ICE_THERMO #define ICE_MK #define ICE_MOMENTUM #define ICE_MOM_BULK #define ICE_EVP #define ICE_ADVECT #define ICE_SMOLAR #define ICE_UPWIND #define ...

if anyone could email me the ice.in file, it would be great. It seems to be missing from my ROMS 3.0 package. Thx a lot. Guillermo.
ps: email to guillo@ucsd.edu

Afer some reasonable troubles, I managed to compile and run ROMS on my Dell laptop. This was with Linux Red Hat Enterprise Advanced Server 4, the intel compiler 9.1, only one supporting, as of May 2006, the Intel Centrino Duo (dual core) processor. Most important. The intel fortran compiler for ...

When compiling ROMS 3.0 with sea ice (PGI + Linux) I get the following message: pgf90 -u -Bstatic -fastsse master.o ocean_control.o propagator.o -o roms libRPM.a libTLM.a libNLM.a libICE.a libUTIL.a libMODS.a -L/usr/local/netcdf-3.6.1/lib -lnetcdf libICE.a(ice_elastic.o)(.text+0x4df4): In function ...

got it.

If one uses the PG linux86 compiler (on the opteron server) you get similar results but not quite what it should be.

So, things worked when I compiled with the PG linux86-64 compiler (opteron is 32-64bit).

Guillermo.

I ran the same code with identical initial conditions, grid and forcing etc (compiled both with 32bit PG compiler). The model's kinetic energy is not the same. In the dual Xeon the model's kinetic energy/volume is 0.000001 units smaller than in a quad Opteron, i.e., only difference is the type of ...

Below you may find results from a global ROMS simulation with climatological forcing started from rest. After 12 years spin-up, all main ocean currents are noticeable. Even at depth the Equatorial counter-current is seen flowing eastward. These results/grid/forcing etc are from an early spin-up test ...

Errata: where I said "area" above, it should read "diagonal" i.e., sqrt(om**2+on**2). So, it's proportional to sqrt(dx**2+dy**2) as the code is now. So, to make each grid cell contribute the same amount of vorticity to the formation and maintainance of western boundary layers, i.e., as a uniform ...

When using a variable-resolution grid, ROMS offers the advantage of making the viscosity (and diffusivity) coefficients proportional to the grid size (user specifies coefficient for the coarsest grid point, then by %define VISC_GRID in cppdefs.h, ROMS computes the visc coeff for the rest of the ...

did you try Conkright (2002)?? Not sure if it is finer than Levitus but it's becoming more popular these days (it is complemented with satellite data).

cheers,

Guiyo.

Anyone does know how to efficiently incorporate climatological sea-ice to ROMS?

I noticed that the Weddell Sea application uses sea ice but it does not vary with time throughout the year.

Thank you.

Guillermo Auad.