New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 15260 for NEMO/branches/NERC/dev_r4.0.6_GEOMETRIC/cfgs – NEMO

Ignore:
Timestamp:
2021-09-14T18:31:24+02:00 (3 years ago)
Author:
acc
Message:

Untested, initial port of GEOMETRIC changes. Includes a bug fix to eken diagnostic in diawri.F90 which, if confirmed, is also relevant to current 4.0-HEAD. #2722

Location:
NEMO/branches/NERC/dev_r4.0.6_GEOMETRIC/cfgs/SHARED
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/NERC/dev_r4.0.6_GEOMETRIC/cfgs/SHARED/field_def_nemo-oce.xml

    r13648 r15260  
    138138        <field id="cfl_cv"       long_name="v-courant number"   unit="#" /> 
    139139        <field id="cfl_cw"       long_name="w-courant number"   unit="#" /> 
     140 
     141         <!-- GEOMETRIC fields (requires nn_aei_ijk_t = 32)  --> 
     142         <field id="eke"               long_name="total EKE (EKE+EPE)"                             unit="m3/s2" /> 
     143         <field id="trd_eke_adv_ubt"   long_name="ubt advective trend of EKE (LHS)"                unit="m3/s3" /> 
     144         <field id="trd_eke_adv_wav"   long_name="wav advective trend of EKE (LHS)"                unit="m3/s3" /> 
     145         <field id="trd_eke_lap"       long_name="diffusive trend of EKE (RHS)"                    unit="m3/s3" /> 
     146         <field id="trd_eke_peS"       long_name="PE to EKE source trend (RHS)"                    unit="m3/s3" /> 
     147         <field id="trd_eke_keS"       long_name="KE to EKE source trend (RHS)"                    unit="m3/s3" /> 
     148         <field id="trd_eke_dis"       long_name="dissipation trend of EKE (RHS)"                  unit="m3/s3" /> 
     149 
    140150 
    141151      </field_group> <!-- grid_T --> 
     
    545555        <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                unit="s-2" /> 
    546556 
     557         <!-- GEOMETRIC fields (requires nn_aei_ijk_t = 32)  --> 
     558         <field id="aeiv_geom"      long_name="3D w-EIV coefficient from GEOMETRIC param."                       unit="m2/s" /> 
     559         <field id="rossby_rad"     long_name="internal Rossby defromation radius"                               unit="m"    grid_ref="grid_W_2D"/> 
     560         <field id="bn2"            long_name="squared Brunt-Vaisala frequency"                                  unit="s-1"  /> 
     561         <field id="c1_vert"        long_name="1st baroclinic mode phase speed"                                  unit="m/s"  grid_ref="grid_W_2D"/> 
     562         <field id="c_ros"          long_name="long Rossby phase speed"                                          unit="m/s"  grid_ref="grid_W_2D"/> 
     563 
    547564      </field_group> 
    548565       
  • NEMO/branches/NERC/dev_r4.0.6_GEOMETRIC/cfgs/SHARED/namelist_ref

    r14372 r15260  
    836836      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    837837      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     838      !                             !   = 32 F(i,j,t)  = GEOMETRIC parameterization        (=> fill namldf_eke) 
    838839      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le  
    839840      rn_Ue        = 0.02           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
     
    841842      ! 
    842843      ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     844   ln_eke_equ    =.true.   ! switch on update of GEOMETRIC eddy energy equation            (=> fill namldf_eke) 
     845                           ! forced to be true if nn_aei_ijk_t = 32 
     846/ 
     847!---------------------------------------------------------------------------------- 
     848&namldf_eke !   GEOMETRIC param. (total EKE equation)                           (nn_aei_ijk_t = 32) 
     849!---------------------------------------------------------------------------------- 
     850   rn_ekedis      =  180.      ! dissipation time scale of EKE [days] 
     851      nn_eke_dis  =    0       ! dissipation option 
     852      !                             !   =  0  constant in space 
     853      !                             !   =-20  read in geom_diss_2D.nc file 
     854   rn_geom        =  0.05      ! geometric parameterization master coefficient (>0 & <1) 
     855   rn_eke_init    =  1.e-1     ! initial total EKE value 
     856   rn_eke_min     =  1.e+0     ! background value of total EKE 
     857   rn_ross_min    =  4.e+3     ! tapering of aeiv based on min Rossby radius [m] 
     858   !                           !   set to zero to not taper it 
     859   rn_eke_lap     =   500.     ! Laplacian diffusion coefficient of EKE 
     860   !                           !   this is in all options below, so set it to zero and nothing is done 
     861   rn_aeiv_min    =  1.e+1     ! minimum bound of eiv coefficient 
     862   rn_aeiv_max    =  1.5e+4    ! maximum bound of eiv coefficient 
     863   rn_SFmin       =  1.0       ! minimum bound of Structure Function 
     864   rn_SFmax       =  1.0       ! maximum bound of Structure Function 
     865   nn_eke_opt     =  1         ! options for terms to include in EKE budget 
     866   !                                !   =  0  PE->EKE conversion, dissipation only 
     867   !                                !   =  1  as 0 but with advection 
     868   !                                !   =  2  as 1 but with additional KE->EKE conversion 
     869   !                                !   for testing purposes: 
     870   !                                !   = 88  only advection by depth-averaged flow 
     871   !                                !   = 99  only Laplacian diffusion 
     872   ln_adv_wav     =  .false.  ! include advection at long Rossby speed 
     873   ln_beta_plane  =  .false.  ! beta plane option for computing long Rossby speed (default: sphere option) 
    843874/ 
    844875!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.