Ocean Modeling Discussion

Inside ROMS/Nonlinear/Biology/fennel_mod.h, the 15 biological tracers (here I'm assuming that the user has activated RIVER_DON,PO4,CARBON,OXYGEN) are attributed the following identification index: 1-NO3,2-NH4,3-chlorophyll,4-phytoplankton,5-zooplankton,6-LdetritusN,7-SdetritusN,8-RdetritusN,9-PO4,10 ...

Line #1 of ROMS/Nonlinear/Biology/red_tide.h appears to have a typo:
#undef NEWMAN -> #undef NEUMANN

Information about the vertical discretization is to be provided in the text file roms_*.in. I think the idea is that if one uses an analytical initial condition, then they would be able to change `N' (number of vertical levels) or the stretching parameters inside roms_*.in, without having to modify ...

A quick update on this---"artificial ponding" is the proper term for what I was referring to. E.g., https://doi.org/10.1016/j.coastaleng.2022.104246 Since the initial post, I came up with a satisfyingly-fast function to handle the ponds in postprocessing. Happy to share if someone is interested (psl ...

When using realistic coastal topography and inundation (WET_DRY), it is not uncommon for wetdry_mask_rho to feature a few wet `puddles' that are disconnected from the main body of water, at least temporarily. I'm looking for a fast `puddle drying' algorithm, i.e. one that is provided a logical mask ...

The Department of Meteorology and Atmospheric Science at The Pennsylvania State University is seeking a postdoctoral scholar in the area of numerical modeling of coastal waters . The postdoctoral scholar will work with Raymond Najjar, Professor of Oceanography ( https://www.ems.psu.edu/directory ...

Since you mention it, there was a short paper circulating around back in ~2011 about running ROMS on a NVIDIA GPU. I don't remember what it said, but I attached it to this post.

The line: has a typo; it should be a upper case letter o (oxygen) but the line has a zero.

A few typos in the trunk version. In file User/External/roms.in: ! Input lateral boundary conditions and climatology file names. The ! USER has the option to split input data time records into several ! NetCDF files (see prologue instructions above). If so, use a single ! line per entry with a ...

I'm running ROMS with the trunk version, BIOLOGY activated, and I have "multi-file" BRY files (one BRY file per calendar year). The code is compiled in serial mode with ifort, USE_DEBUG and USE_NETCDF4. At the end of the year, when ROMS is ready to switch to the next BRY file, it enters this block ...

should be setenv

I accidentally ran ROMS with the default values of Linux-gfortran.mk, meaning ffast-math and O3. Since I had this unexpected data point in hand, I thought I should conduct a few more tests and then report on this thread, in support of the original post (by John Warner) recommending the removal of ...

Good afternoon, I have a little question related to ROMS/Utility/get_nudgcoef.F, first introduced in the post: https://www.myroms.org/projects/src/ticket/627 As described in the post, get_nudgcoef.F is designed to handle values in either day-1 or s-1. Its logic follows the "units" metadata ...

Good afternoon, I have two suggestions regarding Nonlinear/Biology/fennel.h (Rutgers trunk version): (1) On line 516, the code splits into two cases according to surface PAR: IF (PARsur(i).gt.0.0_r8) THEN followed on line 714 by: ELSE I'll refer to those two cases as "day-time" (PAR>0) and "night ...

Good morning, I apologize in advance if I'm wrong about what follows. I could be misinterpreting the code. Early on in step3d_t.F (lines 334-350), oHz is computed as: IF (Lmpdata) THEN DO k=1,N(ng) DO j=Jstrm2,Jendp2 DO i=Istrm2,Iendp2 oHz(i,j,k)=1.0_r8/Hz(i,j,k) END DO END DO END DO ELSE DO k=1,N ...

Just to complete the earlier post---here is the relevant part of ROMS/Utility/read_phypar.F (as of rev.1008): CASE ('Hadvection') IF (itracer.lt.(NAT+NPT)) THEN itracer=itracer+1 ELSE itracer=1 ! next nested grid END IF itrc=itracer Npts=load_tadv(Nval, Cval, line, nline, itrc, igrid, & & itracer, 1 ...

Good afternoon, I think there is a typo in the GLOSSARY section of roms.in. It says: ! or in nested applications ! ! Hadvection == U3 \ ! temperature, Grid 1 ! U3 \ ! temperature, Grid 2 ! U3 \ ! temperature, Grid 3 ! HSIMT \ ! salinity, Grid 1 ! HSIMT \ ! salinity, Grid 2 ! HSIMT ! salinity, Grid 3 ...

A Post-Doctoral Scholar position is available at The Pennsylvania State University to study the transport and fate of plastics in the Chesapeake Bay using numerical models and to collaborate with experimentalists. The Post-Doctoral Scholar will work with an interdisciplinary team at two institutions ...

Hi,
There is a typo in ad_set_vbc.F on line 1165: should be (needs a "l" between "f" and "x").
Pierre

Good afternoon, I think there's a tiny typo in the recent versions of mod_scalars.F (revision >=919): #ifdef T_PASSIVE IF (.not.allocated(inert)) THEN allocate ( inert(NPT) ) END IF #endif Since "inert" is a pointer, I think it should be "IF(.not.associated(inert)" rather than "IF(.not.allocated ...

Maybe you can figure out what the parameters should be - probably a very small but finite value not zero WTYPE=9 (dark coastal water) is absent from the original papers (e.g. Paulson & Simpson JPO 1977) but the following publication (attached) from: https://publications.europa.eu/en/publication ...

I would like to report a division by zero that occurs when WTYPE=9 and SOLAR_SOURCE are selected. WTYPE=9 is defined in mod_scalars.F as: real(r8), dimension(9) :: lmd_mu1 = & & (/ 0.35_r8, 0.6_r8, 1.0_r8, 1.5_r8, 1.4_r8, & & 0.42_r8, 0.37_r8, 0.33_r8, 0.0_r8 /) Note that lmd_mu1 is 0.0_r8 for WTYPE ...

Just a suggestion---sensitivity to the choice of the advection scheme may possibly indicate that your grid stretching is more than it should be. Could be worth trying: (1) decreasing the number of vertical levels, (2) use conservative stretching parameters, (3) increase the horizontal resolution ...

Yes, that would work for my purpose. I see the option in my own version of the code (v.~3.3) and I will give it a try. Thanks for the reply.
Pierre

Hi, I have high-resolution outputs from ROMS and I would like to advect+diffuse an arbitrary 3-D scalar field using an `offline' implementation of ROMS advection-diffusion schemes. In other words, I'm looking for a Fortran/C/whatever program that takes as input the netcdf files, interpolates them in ...

I updated my code and it works great!
Thank you very much for providing this fix.

Best regards,

Pierre

I would be very grateful if anybody that does not encounter this issue could describe what is happening in his roms code during the reading of boundary input data such as zeta_north(:,:). It can be simply a matter of inserting a couple of write(*,*) in your code to see what roms does. Just (please ...

I get the runtime error with: gfortran (gcc version 4.3.2, from Linux Debian Lenny). I also get the error on a different machine that uses: f95: Sun Fortran 95 8.3 SunOS_i386 Patch 127002-04 2008/04/16 and it says ****** FORTRAN RUN-TIME SYSTEM ****** An array expression does not conform : dim 1 is ...

Thanks for this reply! Unfortunately, I checked with ncdump -h, and by opening the file in Matlab, and zeta_north correctly has two dimensions (zeta_time and xi_rho), with one record in time. Numerous users of ROMS use boundary files and this is the first time that this is reported. I agree that ...

Perhaps a minimal testcase would make things clearer: program testcase implicit none integer, parameter :: UBi = 403, LBi = -2, UBj = 1 ! Dummy values. integer :: ng, model, ncid, Trec, Iend, Istr character(400) :: ncfile, Vname real, dimension((UBi-LBi+1)*UBj) :: A ! Dummy values for testcase. ng ...

Thanks for your help! ROMS is trying to get zeta_north(:,:) from my ocean_bry.nc file. My zeta_north(:,:) array has only one record in time, and the other dimension is space (size Lm+2). What happens then is that net_cdf_get_fvar is called over line 400 of get_ngfld: CALL netcdf_get_fvar (ng, mode ...

Hello, I have a small question regarding subroutine netcdf_get_fvar. This routine is called within get_ngfld.F around line 400 (SVN rev.491M, very recent) to obtain a one-dimensional array called 'A(:)', yet the 'start' and 'total' arguments each have two elements (as if A(:) was two-dimensional ...