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 12149 for NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW – NEMO

Ignore:
Timestamp:
2019-12-10T15:03:24+01:00 (5 years ago)
Author:
ayoung
Message:

Updated trunk to 12072

Location:
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/context_nemo.xml

    r9930 r12149  
    66<context id="nemo"> 
    77<!-- $id$ --> 
     8    <variable_definition> 
     9    <!-- Year of time origin for NetCDF files; defaults to 1800 --> 
     10       <variable id="ref_year" type="int"   > 1800 </variable> 
     11       <variable id="rau0"     type="float" > 1026.0 </variable> 
     12       <variable id="cpocean"  type="float" > 3991.86795711963 </variable> 
     13       <variable id="convSpsu" type="float" > 0.99530670233846  </variable> 
     14       <variable id="rhoic"    type="float" > 917.0 </variable> 
     15       <variable id="rhosn"    type="float" > 330.0 </variable> 
     16       <variable id="missval"  type="float" > 1.e20 </variable> 
     17    </variable_definition> 
    818<!-- Fields definition --> 
    919    <field_definition src="./field_def_nemo-oce.xml"/>   <!--  NEMO ocean dynamics                     --> 
     
    1828     
    1929    <axis_definition> 
    20       <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" /> 
    21       <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" /> 
    22       <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" /> 
    23       <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" /> 
    24       <axis id="nfloat" long_name="Float number"      unit="-"                 /> 
    25       <axis id="icbcla"  long_name="Iceberg class"      unit="1"               /> 
    26       <axis id="ncatice" long_name="Ice category"       unit="1"               /> 
    27       <axis id="iax_20C" long_name="20 degC isotherm"   unit="degC"            /> 
    28       <axis id="iax_28C" long_name="28 degC isotherm"   unit="degC"            /> 
     30      <axis id="deptht"  long_name="Vertical T levels" unit="m" positive="down" /> 
     31      <axis id="depthu"  long_name="Vertical U levels" unit="m" positive="down" /> 
     32      <axis id="depthv"  long_name="Vertical V levels" unit="m" positive="down" /> 
     33      <axis id="depthw"  long_name="Vertical W levels" unit="m" positive="down" /> 
     34      <axis id="nfloat"  long_name="Float number"      unit="-"                 /> 
     35      <axis id="icbcla"  long_name="Iceberg class"     unit="1"                 /> 
     36      <axis id="ncatice" long_name="Ice category"      unit="1"                 /> 
     37      <axis id="iax_20C" long_name="20 degC isotherm"  unit="degC"              /> 
     38      <axis id="iax_28C" long_name="28 degC isotherm"  unit="degC"              /> 
    2939    </axis_definition> 
    3040  
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_cen-ahm1000_cfg

    r10075 r12149  
    210210   !                       ! coefficients 
    211211   rn_avm0     =   1.e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
    212    rn_avt0     =   0.0e      !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     212   rn_avt0     =   0.0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
    213213   nn_avb      =    0        !  profile for background avt & avm (=1) or not (=0) 
    214214   nn_havtb    =    0        !  horizontal shape for avtb (=1) or not (=0) 
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_ubs_cfg

    r10075 r12149  
    210210   !                       ! coefficients 
    211211   rn_avm0     =   1.e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
    212    rn_avt0     =   0.0e      !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     212   rn_avt0     =   0.0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
    213213   nn_avb      =    0        !  profile for background avt & avm (=1) or not (=0) 
    214214   nn_havtb    =    0        !  horizontal shape for avtb (=1) or not (=0) 
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_sco_FCT4_flux_ubs_cfg

    r10075 r12149  
    210210   !                       ! coefficients 
    211211   rn_avm0     =   1.e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
    212    rn_avt0     =   0.0e      !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     212   rn_avt0     =   0.0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
    213213   nn_avb      =    0        !  profile for background avt & avm (=1) or not (=0) 
    214214   nn_havtb    =    0        !  horizontal shape for avtb (=1) or not (=0) 
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_zps_FCT2_flux_ubs_cfg

    r10075 r12149  
    211211   !                       ! coefficients 
    212212   rn_avm0     =   1.e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
    213    rn_avt0     =   0.0e      !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     213   rn_avt0     =   0.0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
    214214   nn_avb      =    0        !  profile for background avt & avm (=1) or not (=0) 
    215215   nn_havtb    =    0        !  horizontal shape for avtb (=1) or not (=0) 
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_flux_ubs_cfg

    r10075 r12149  
    3434   cn_exp      =   "OVF_zps_FCT4_flux_ubs"  !  experience name 
    3535   nn_it000    =       1   !  first time step 
    36    !nn_itend    =    6120  ! here 17h of simulation  (=6120 time-step)  
    37    nn_itend    =    5760   ! here 16h of simulation  (=5760 time-step) abort after 5802 for zps: pb of physics conditions 
     36   nn_itend    =    6120   ! here 17h of simulation  (=6120 time-step)  
    3837   nn_istate   =       0   !  output the initial state (1) or not (0) 
    39    nn_stock    =    1080   !  frequency of creation of a restart file (modulo referenced to 1) 
    40    nn_write    =    1080   !  frequency of write in the output file   (modulo referenced to nn_it000) 
     38   nn_stock    =    6120   !  frequency of creation of a restart file (modulo referenced to 1) 
     39   nn_write    =    6120   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    4140/ 
    4241!----------------------------------------------------------------------- 
     
    255254&namzdf        !   vertical physics manager                             (default: NO selection) 
    256255!----------------------------------------------------------------------- 
     256   !                       ! adaptive-implicit vertical advection 
     257   ln_zad_Aimp = .true.       !  Courant number dependent scheme (Shchepetkin 2015) 
     258   ! 
    257259   !                       ! type of vertical closure (required) 
    258260   ln_zdfcst   = .true.      !  constant mixing 
     
    279281   !                       ! coefficients 
    280282   rn_avm0     =   1.e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
    281    rn_avt0     =   0.0e      !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     283   rn_avt0     =   0.0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
    282284   nn_avb      =    0        !  profile for background avt & avm (=1) or not (=0) 
    283285   nn_havtb    =    0        !  horizontal shape for avtb (=1) or not (=0) 
     
    291293!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    292294!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    293 !!   namflo       float parameters                                      ("key_float") 
    294 !!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
    295 !!   namdct       transports through some sections                      ("key_diadct") 
     295!!   namflo       float parameters                                      (default: OFF) 
     296!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
     297!!   nam_diadct   transports through some sections                      (default: OFF) 
    296298!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
    297299!!   nam_dia25h   25h Mean Output                                       (default: OFF) 
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_vect_een_cfg

    r10075 r12149  
    210210   !                       ! coefficients 
    211211   rn_avm0     =   1.e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
    212    rn_avt0     =   0.0e      !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     212   rn_avt0     =   0.0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
    213213   nn_avb      =    0        !  profile for background avt & avm (=1) or not (=0) 
    214214   nn_havtb    =    0        !  horizontal shape for avtb (=1) or not (=0) 
  • NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/MY_SRC/usrdef_nam.F90

    r10074 r12149  
    3939CONTAINS 
    4040 
    41    SUBROUTINE usr_def_nam( ldtxt, ldnam, cd_cfg, kk_cfg, kpi, kpj, kpk, kperio ) 
     41   SUBROUTINE usr_def_nam( cd_cfg, kk_cfg, kpi, kpj, kpk, kperio ) 
    4242      !!---------------------------------------------------------------------- 
    4343      !!                     ***  ROUTINE dom_nam  *** 
     
    5151      !! ** input   : - namusr_def namelist found in namelist_cfg 
    5252      !!---------------------------------------------------------------------- 
    53       CHARACTER(len=*), DIMENSION(:), INTENT(out) ::   ldtxt, ldnam    ! stored print information 
    5453      CHARACTER(len=*)              , INTENT(out) ::   cd_cfg          ! configuration name 
    5554      INTEGER                       , INTENT(out) ::   kk_cfg          ! configuration resolution 
     
    5756      INTEGER                       , INTENT(out) ::   kperio          ! lateral global domain b.c.  
    5857      ! 
    59       INTEGER ::   ios, ii   ! Local integer 
     58      INTEGER ::   ios   ! Local integer 
    6059      !! 
    6160      NAMELIST/namusr_def/ ln_zco, ln_zps, ln_sco, rn_dx, rn_dz 
    6261      !!---------------------------------------------------------------------- 
    6362      ! 
    64       ii = 1 
    65       ! 
    6663      REWIND( numnam_cfg )          ! Namelist namusr_def (exist in namelist_cfg only) 
    6764      READ  ( numnam_cfg, namusr_def, IOSTAT = ios, ERR = 902 ) 
    68 902   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namusr_def in configuration namelist', .TRUE. ) 
     65902   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namusr_def in configuration namelist' ) 
    6966      ! 
    70       WRITE( ldnam(:), namusr_def ) 
     67      IF(lwm)   WRITE( numond, namusr_def ) 
    7168      ! 
    7269      cd_cfg = 'OVERFLOW'           ! name & resolution (not used) 
     
    7976      ! 
    8077      !                             ! control print 
    81       WRITE(ldtxt(ii),*) '   '                                                                          ;   ii = ii + 1 
    82       WRITE(ldtxt(ii),*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg'   ;   ii = ii + 1 
    83       WRITE(ldtxt(ii),*) '~~~~~~~~~~~ '                                                                 ;   ii = ii + 1 
    84       WRITE(ldtxt(ii),*) '   Namelist namusr_def : OVERFLOW test case'                                  ;   ii = ii + 1 
    85       WRITE(ldtxt(ii),*) '      type of vertical coordinate : '                                         ;   ii = ii + 1 
    86       WRITE(ldtxt(ii),*) '         z-coordinate flag                     ln_zco = ', ln_zco             ;   ii = ii + 1 
    87       WRITE(ldtxt(ii),*) '         z-partial-step coordinate flag        ln_zps = ', ln_zps             ;   ii = ii + 1 
    88       WRITE(ldtxt(ii),*) '         s-coordinate flag                     ln_sco = ', ln_sco             ;   ii = ii + 1 
    89       WRITE(ldtxt(ii),*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters'   ;   ii = ii + 1 
    90       WRITE(ldtxt(ii),*) '      vertical   resolution                    rn_dz  = ', rn_dz, ' meters'   ;   ii = ii + 1 
    91       WRITE(ldtxt(ii),*) '      OVERFLOW domain = 200 km x 3 grid-points x 2000 m'                      ;   ii = ii + 1 
    92       WRITE(ldtxt(ii),*) '         resulting global domain size :        jpiglo = ', kpi                ;   ii = ii + 1 
    93       WRITE(ldtxt(ii),*) '                                               jpjglo = ', kpj                ;   ii = ii + 1 
    94       WRITE(ldtxt(ii),*) '                                               jpkglo = ', kpk                ;   ii = ii + 1 
     78      WRITE(numout,*) '   ' 
     79      WRITE(numout,*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg' 
     80      WRITE(numout,*) '~~~~~~~~~~~ ' 
     81      WRITE(numout,*) '   Namelist namusr_def : OVERFLOW test case' 
     82      WRITE(numout,*) '      type of vertical coordinate : ' 
     83      WRITE(numout,*) '         z-coordinate flag                     ln_zco = ', ln_zco 
     84      WRITE(numout,*) '         z-partial-step coordinate flag        ln_zps = ', ln_zps 
     85      WRITE(numout,*) '         s-coordinate flag                     ln_sco = ', ln_sco 
     86      WRITE(numout,*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters' 
     87      WRITE(numout,*) '      vertical   resolution                    rn_dz  = ', rn_dz, ' meters' 
     88      WRITE(numout,*) '      OVERFLOW domain = 200 km x 3 grid-points x 2000 m' 
     89      WRITE(numout,*) '         resulting global domain size :        jpiglo = ', kpi 
     90      WRITE(numout,*) '                                               jpjglo = ', kpj 
     91      WRITE(numout,*) '                                               jpkglo = ', kpk 
    9592      ! 
    9693      !                             ! Set the lateral boundary condition of the global domain 
    9794      kperio = 0                    ! OVERFLOW configuration : close basin 
    9895      ! 
    99       WRITE(ldtxt(ii),*) '   '                                                                          ;   ii = ii + 1 
    100       WRITE(ldtxt(ii),*) '   Lateral boundary condition of the global domain'                           ;   ii = ii + 1 
    101       WRITE(ldtxt(ii),*) '      OVERFLOW : closed basin                  jperio = ', kperio             ;   ii = ii + 1 
     96      WRITE(numout,*) '   ' 
     97      WRITE(numout,*) '   Lateral boundary condition of the global domain' 
     98      WRITE(numout,*) '      OVERFLOW : closed basin                  jperio = ', kperio 
    10299      ! 
    103100   END SUBROUTINE usr_def_nam 
Note: See TracChangeset for help on using the changeset viewer.