Ocean Modeling Discussion

Hi,

I modified pyroms smoothing code.
(bath_tools.RoughnessMatrix => RoughMtrx,
bathy_smoothing.smoothing_Laplacian_rx0 => smthng_laplacian_rx0,
bathy_smoothing.smoothing_PlusMinus_rx0 => smthng_PlsMns_rx0)

They worked faster than the originals.
If you like, you can pick them up from
https ...

Hi,

I tried a test case using DEFINE nesting.

I think there might be bugs about fine2coarse.
Please let me know if these are bugs, if you know well about fine2coase.

I think that line 291-293 in main3d.F

IF ((istep.eq.Nsteps).and.(RefineScale(ng).gt.0)) THEN
CALL nesting (ng, iNLM, n2way)
END ...

Hi,

I tried a test case using nesting model + mpi parallel.
But, Blowing up at the first step.

I found that it may be a bug in Nesting.F of ROMS revison 674, if parallel application is on.
It is around calculating cff1 for mass flux at the physical boundaries in SUBROUTINE put_refine2d ...

Dear ROM people,

I tried to calculate a passive tracer with define KPP option.
The tracer was calculated to be not 0, although I. C. and all B. C. of the tracer were set to 0.d0.

I read lines form 763 to 777 of pre_step3d.F as follows;
DO itrc=1,NT(ng)
ltrc=MIN(NAT,itrc)
DO k=1,N(ng)-1
DO i ...