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://github ...

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 IF ...

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. =====cff1 ...

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=Istr ...