GLS y MY25 in ROMS

[oartal]

Hi, I'm a student of physics in the University of Concepción Chile.
I'm working now in the test example of Upwelling. I want to see how the process of upwelling varies with different vertical mixing schemes. My problem is that my code not have GLS_mixing and MY25_mixing, where can I download or as I get it?

thanks in advance

Osvaldo Artal
University of Concepción, Chile

[ggerbi]

The mixing schemes are in gls_corstep.F and my25corstep.F, which are in the Nonlinear directory of your source code. You can get the essential features of Mellor-Yamada by choosing kkl in GLS.

[jacopo]

your code?!? did you get ROMS via SVN or what?

Let's assume you did (If you didn't, go here https://www.myroms.org/wiki/index.php/Subversion)

In ROMS/Include/upwelling.h you have

Code: Select all

#if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
#else
# define ANA_VMIX
#endif

this means that you have to define either GLS_MIXING or MY25_MIXING otherwise the model will use ANA_VMIX.

In case you define GLS_MIXING, you have to select actually which scheme are you using, whether k-e, k-w, k-kl or gen specifying correct values in your ROMS/External/ocean_upwelling.in.

At the end of the file you find all the info you need

Code: Select all

!------------------------------------------------------------------------------
! Generic length-scale turbulence closure parameters.
!------------------------------------------------------------------------------
!
!  GLS_P       Stability exponent (non-dimensional).
!
!  GLS_M       Turbulent kinetic energy exponent (non-dimensional).
!
!  GLS_N       Turbulent length scale exponent (non-dimensional).
!
!  GLS_Kmin    Minimum value of specific turbulent kinetic energy
!
!  GLS_Pmin    Minimum Value of dissipation.
!
! Closure independent constraint parameters (non-dimensional):
!
!  GLS_CMU0    Stability coefficient.
!
!  GLS_C1      Shear production coefficient.
!
!  GLS_C2      Dissipation coefficient.
!
!  GLS_C3M     Buoyancy production coefficient (minus).
!
!  GLS_C3P     Buoyancy production coefficient (plus).
!
!  GLS_SIGK    Constant Schmidt number (non-dimensional) for turbulent
!              kinetic energy diffusivity.
!
!  GLS_SIGP    Constant Schmidt number (non-dimensional) for turbulent
!              generic statistical field, "psi".
!
! Suggested values for various parameterizations:
!
!              k-kl         k-epsilon    k-omega      gen
!
!      GLS_P = 0.d0          3.0d0       -1.0d0        2.0d0
!      GLS_M = 1.d0          1.5d0        0.5d0        1.0d0
!      GLS_N = 1.d0         -1.0d0       -1.0d0       -0.67d0
!   GLS_Kmin = 5.0d-6        7.6d-6       7.6d-6       1.0d-8
!   GLS_Pmin = 5.0d-6        1.0d-12      1.0d-12      1.0d-8
!
!   GLS_CMU0 = 0.5544d0      0.5477d0     0.5477d0     0.5544d0
!     GLS_C1 = 0.9d0         1.44d0       0.555d0      1.00d0
!     GLS_C2 = 0.52d0        1.92d0       0.833d0      1.22d0
!    GLS_C3M = 2.5d0        -0.4d0       -0.6d0        0.1d0
!    GLS_C3P = 1.0d0         1.0d0        1.0d0        1.0d0
!   GLS_SIGK = 1.96d0        1.0d0        2.0d0        0.8d0
!   GLS_SIGP = 1.96d0        1.30d0       2.0d0        1.07d0

set the values you want where you find

Code: Select all

! Generic length-scale turbulence closure parameters.

       GLS_P == 3.0d0                           ! K-epsilon
       GLS_M == 1.5d0
       GLS_N == -1.0d0
    GLS_Kmin == 7.6d-6
    GLS_Pmin == 1.0d-12

    GLS_CMU0 == 0.5477d0
      GLS_C1 == 1.44d0
      GLS_C2 == 1.92d0
     GLS_C3M == -0.4d0
     GLS_C3P == 1.0d0
    GLS_SIGK == 1.0d0
    GLS_SIGP == 1.30d0

hope this helps

[oartal]

hi, what happens is that I have the version of ROMS AGRIF. How do I get the schematics or how can I adapt gls_corstep.F and my25corstep.F for this code?

Thanks

[kate]

ROMS_AGRIF is a different animal. Since each of these ROMS versions uses different timestepping algorithms, you can't just mix and match without a lot of careful rewriting, checking and testing. You'd best ask the AGRIF gang what vertical mixing schemes they support.

[oartal]

I donwload ROMS via SVN and ran the model for the test example of upwelling.
Modify the file upwelling.h and check ANA_VMIX,GLS and MY25. My new question is
how modify upwelling.h for LMD_MIXING and BVF_MIXING.
For LMD define LMD_RIMIX, LMD_CONVEC,LMD_SKPP and LMD_BKPP but I get error with ocean_frc.nc
For BVF I have not done anything. How can I run these two schemes. Thank you very much

[oartal]

I achieve running the schemes BFV, but not LMD, the following is the massage of error.

READ_PHYPAR - could not find input file: ocean_frc.nc

Elapsed CPU time (seconds):

Thread # 0 CPU: 0.012
Total: 0.012

Nonlinear model elapsed time profile:

Total: 0.000 0.0000

All percentages are with respect to total time = 0.012

ROMS/TOMS - Output NetCDF summary for Grid 01:

ROMS/TOMS - I/O error ............... exit_flag: 4


ERROR: I/O related problem.


thanks in advance

[kate]

The LMD scheme is looking for shortwave radiation. See these clauses in globaldefs.h:

Code: Select all

 ( defined LMD_SKPP     && !defined ANA_SRFLUX)   || \
:
#  define FRC_FILE

and

Code: Select all

#if defined LMD_SKPP     || defined SOLAR_SOURCE   || \
    defined BULK_FLUXES  || defined BIOLOGY
# define SHORTWAVE
#endif

You might try ANA_SRFLUX.

I don't know about the BVF scheme, but it has no such clause in globaldefs.

[oartal]

thank kate, your suggestion worked correctly.