Changeset 12551


Ignore:
Timestamp:
2020-03-13T15:51:01+01:00 (8 months ago)
Author:
smasson
Message:

dev_r12472_ASINTER-05: update to trunk@12550, see #2156

Location:
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback
Files:
1 deleted
12 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg

    r12495 r12551  
    2828&namdom        !   time and space domain 
    2929!----------------------------------------------------------------------- 
    30    rn_Dt      = 5400.     !  time step for the dynamics and tracer 
     30   rn_Dt       = 5400.     !  time step for the dynamics and tracer 
    3131/ 
    3232!----------------------------------------------------------------------- 
     
    381381   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
    382382   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
     383 
     384   cn_dir      = './'      !  root directory for the iwm data location                                                                            
     385   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     386   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     387   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     388   sn_mpb      = 'mixing_power_bot'      , -12               , 'field'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
     389   sn_mpp      = 'mixing_power_pyc'      , -12               , 'field'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
     390   sn_mpc      = 'mixing_power_cri'      , -12               , 'field'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
     391   sn_dsb      = 'decay_scale_bot'       , -12               , 'field'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
     392   sn_dsc      = 'decay_scale_cri'       , -12               , 'field'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
    383393/ 
    384394!!====================================================================== 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/cfgs/SHARED/namelist_ref

    r12499 r12551  
    7272   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    7373   ! 
    74    rn_Dt      = 5400.     !  time step for the dynamics and tracer 
     74   rn_Dt       = 5400.     !  time step for the dynamics and tracer 
    7575   rn_atfp     =    0.1    !  asselin time filter parameter 
    7676   ! 
     
    217217   nn_ice      = 0         !  =0 no ice boundary condition 
    218218      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
    219       !                    !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    220       !                    !          except in AGRIF zoom where it has to be specified 
     219      !                    !  =2 or 3 for SI3 and CICE, respectively 
    221220   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
    222221      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
     
    272271      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
    273272      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
    274       rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity  
     273      rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity 
    275274      !                     !  used to calculate the wind stress 
    276275      !                     ! (0. => absolute or 1. => relative winds) 
     
    292291   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    293292   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    294    sn_hpgi     = 'NONE'                       ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
    295    sn_hpgj     = 'NONE'                       ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
     293   sn_hpgi     = 'NOT USED'                   ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
     294   sn_hpgj     = 'NOT USED'                   ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
    296295   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    297296   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     
    999998      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    1000999         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed 
    1001          nn_e      = 30          ! =F : the number of sub-step in rn_Dt seconds 
     1000         nn_e         = 30          ! =F : the number of sub-step in rn_Dt seconds 
    10021001      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F) 
    10031002/ 
     
    11931192   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
    11941193   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    1195 / 
    1196  
     1194 
     1195   cn_dir      = './'      !  root directory for the iwm data location                                                                            
     1196   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     1197   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     1198   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     1199   sn_mpb      = 'NOT USED'              , -12               , 'mixing_power_bot' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1200   sn_mpp      = 'NOT USED'              , -12               , 'mixing_power_pyc' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1201   sn_mpc      = 'NOT USED'              , -12               , 'mixing_power_cri' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1202   sn_dsb      = 'NOT USED'              , -12               , 'decay_scale_bot'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1203   sn_dsc      = 'NOT USED'              , -12               , 'decay_scale_cri'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1204/ 
    11971205!!====================================================================== 
    11981206!!                  ***  Diagnostics namelists  ***                   !! 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/ABL/sbcabl.F90

    r12495 r12551  
    7575      !!--------------------------------------------------------------------- 
    7676 
    77       REWIND( numnam_ref )              ! Namelist namsbc_abl in reference namelist : ABL parameters 
     77      ! Namelist namsbc_abl in reference namelist : ABL parameters 
    7878      READ  ( numnam_ref, namsbc_abl, IOSTAT = ios, ERR = 901 ) 
    7979901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_abl in reference namelist' ) 
    80       ! 
    81       REWIND( numnam_cfg )              ! Namelist namsbc_abl in configuration namelist : ABL parameters 
     80      ! Namelist namsbc_abl in configuration namelist : ABL parameters 
    8281      READ  ( numnam_cfg, namsbc_abl, IOSTAT = ios, ERR = 902 ) 
    8382902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_abl in configuration namelist' ) 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/ICE/icectl.F90

    r12495 r12551  
    725725       
    726726      CALL prt_ctl_info(' ') 
    727       CALL prt_ctl_info(' - Heat / FW fluxes : ') 
    728       CALL prt_ctl_info('   ~~~~~~~~~~~~~~~~~~ ') 
    729       CALL prt_ctl(tab2d_1=sst_m  , clinfo1= ' sst   : ', tab2d_2=sss_m     , clinfo2= ' sss       : ') 
    730       CALL prt_ctl(tab2d_1=qsr    , clinfo1= ' qsr   : ', tab2d_2=qns       , clinfo2= ' qns       : ') 
    731       CALL prt_ctl(tab2d_1=emp    , clinfo1= ' emp   : ', tab2d_2=sfx       , clinfo2= ' sfx       : ') 
    732        
    733       CALL prt_ctl_info(' ') 
    734727      CALL prt_ctl_info(' - Stresses : ') 
    735728      CALL prt_ctl_info('   ~~~~~~~~~~ ') 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/BDY/bdydta.F90

    r12396 r12551  
    9292      INTEGER ::  ii, ij, ik, igrd, ipl               ! local integers 
    9393      INTEGER,   DIMENSION(jpbgrd)     ::   ilen1  
    94       INTEGER,   DIMENSION(:), POINTER ::   nblen, nblenrim  ! short cuts 
    9594      TYPE(OBC_DATA)         , POINTER ::   dta_alias        ! short cut 
    9695      TYPE(FLD), DIMENSION(:), POINTER ::   bf_alias 
     
    108107         DO jbdy = 1, nb_bdy 
    109108            ! 
    110             nblen    => idx_bdy(jbdy)%nblen 
    111             nblenrim => idx_bdy(jbdy)%nblenrim 
    112             ! 
    113109            IF( nn_dyn2d_dta(jbdy) == 0 ) THEN  
    114                ilen1(:) = nblen(:) 
    115110               IF( dta_bdy(jbdy)%lneed_ssh ) THEN  
    116111                  igrd = 1 
    117                   DO ib = 1, ilen1(igrd) 
     112                  DO ib = 1, idx_bdy(jbdy)%nblenrim(igrd)   ! ssh is allocated and used only on the rim 
    118113                     ii = idx_bdy(jbdy)%nbi(ib,igrd) 
    119114                     ij = idx_bdy(jbdy)%nbj(ib,igrd) 
     
    123118               IF( dta_bdy(jbdy)%lneed_dyn2d) THEN  
    124119                  igrd = 2 
    125                   DO ib = 1, ilen1(igrd) 
     120                  DO ib = 1, SIZE(dta_bdy(jbdy)%u2d)   ! u2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init 
    126121                     ii = idx_bdy(jbdy)%nbi(ib,igrd) 
    127122                     ij = idx_bdy(jbdy)%nbj(ib,igrd) 
     
    129124                  END DO 
    130125                  igrd = 3 
    131                   DO ib = 1, ilen1(igrd) 
     126                  DO ib = 1, SIZE(dta_bdy(jbdy)%v2d)   ! v2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init 
    132127                     ii = idx_bdy(jbdy)%nbi(ib,igrd) 
    133128                     ij = idx_bdy(jbdy)%nbj(ib,igrd) 
     
    138133            ! 
    139134            IF( nn_dyn3d_dta(jbdy) == 0 ) THEN  
    140                ilen1(:) = nblen(:) 
    141135               IF( dta_bdy(jbdy)%lneed_dyn3d ) THEN  
    142136                  igrd = 2  
    143                   DO ib = 1, ilen1(igrd) 
     137                  DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 
    144138                     DO ik = 1, jpkm1 
    145139                        ii = idx_bdy(jbdy)%nbi(ib,igrd) 
     
    149143                  END DO 
    150144                  igrd = 3  
    151                   DO ib = 1, ilen1(igrd) 
     145                  DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 
    152146                     DO ik = 1, jpkm1 
    153147                        ii = idx_bdy(jbdy)%nbi(ib,igrd) 
     
    160154 
    161155            IF( nn_tra_dta(jbdy) == 0 ) THEN  
    162                ilen1(:) = nblen(:) 
    163156               IF( dta_bdy(jbdy)%lneed_tra ) THEN 
    164157                  igrd = 1  
    165                   DO ib = 1, ilen1(igrd) 
     158                  DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 
    166159                     DO ik = 1, jpkm1 
    167160                        ii = idx_bdy(jbdy)%nbi(ib,igrd) 
     
    176169#if defined key_si3 
    177170            IF( nn_ice_dta(jbdy) == 0 ) THEN    ! set ice to initial values 
    178                ilen1(:) = nblen(:) 
    179171               IF( dta_bdy(jbdy)%lneed_ice ) THEN 
    180172                  igrd = 1    
    181173                  DO jl = 1, jpl 
    182                      DO ib = 1, ilen1(igrd) 
     174                     DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 
    183175                        ii = idx_bdy(jbdy)%nbi(ib,igrd) 
    184176                        ij = idx_bdy(jbdy)%nbj(ib,igrd) 
     
    338330            DO jbdy = 1, nb_bdy      ! Tidal component added in ts loop 
    339331               IF ( nn_dyn2d_dta(jbdy) .GE. 2 ) THEN 
    340                   nblen => idx_bdy(jbdy)%nblen 
    341                   nblenrim => idx_bdy(jbdy)%nblenrim 
    342                   IF( cn_dyn2d(jbdy) == 'frs' ) THEN   ;   ilen1(:)=nblen(:) 
    343                   ELSE                                 ;   ilen1(:)=nblenrim(:) 
     332                  IF( cn_dyn2d(jbdy) == 'frs' ) THEN   ;   ilen1(:)=idx_bdy(jbdy)%nblen(:) 
     333                  ELSE                                 ;   ilen1(:)=idx_bdy(jbdy)%nblenrim(:) 
    344334                  ENDIF 
    345335                  IF ( dta_bdy(jbdy)%lneed_ssh   ) dta_bdy_s(jbdy)%ssh(1:ilen1(1)) = dta_bdy(jbdy)%ssh(1:ilen1(1)) 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/LBC/lib_mpp.F90

    r12377 r12551  
    402402# if defined key_mpi2 
    403403      IF( ln_timing ) CALL tic_tac( .TRUE., ld_global = .TRUE.) 
    404       CALL  mpi_allreduce( y_in(:), todelay(idvar)%y1d(:), isz, MPI_DOUBLE_COMPLEX, mpi_sumdd, ilocalcomm, ndelayid(idvar), ierr ) 
     404      CALL  mpi_allreduce( y_in(:), todelay(idvar)%y1d(:), isz, MPI_DOUBLE_COMPLEX, mpi_sumdd, ilocalcomm, ierr ) 
     405      ndelayid(idvar) = 1 
    405406      IF( ln_timing ) CALL tic_tac(.FALSE., ld_global = .TRUE.) 
    406407# else 
     
    469470# if defined key_mpi2 
    470471      IF( ln_timing ) CALL tic_tac( .TRUE., ld_global = .TRUE.) 
    471       CALL  mpi_allreduce( p_in(:), todelay(idvar)%z1d(:), isz, MPI_DOUBLE_PRECISION, mpi_max, ilocalcomm, ndelayid(idvar), ierr ) 
     472      CALL  mpi_allreduce( p_in(:), todelay(idvar)%z1d(:), isz, MPI_DOUBLE_PRECISION, mpi_max, ilocalcomm, ierr ) 
     473      ndelayid(idvar) = 1 
    472474      IF( ln_timing ) CALL tic_tac(.FALSE., ld_global = .TRUE.) 
    473475# else 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/SBC/cpl_oasis3.F90

    r12377 r12551  
    365365                     WRITE(numout,*) '     - Minimum value is ', MINVAL(pdata(nldi:nlei,nldj:nlej,jc)) 
    366366                     WRITE(numout,*) '     - Maximum value is ', MAXVAL(pdata(nldi:nlei,nldj:nlej,jc)) 
    367                      WRITE(numout,*) '     -     Sum value is ', SUM(pdata(nldi:nlei,nldj:nlej,jc)) 
     367                     WRITE(numout,*) '     -     Sum value is ',    SUM(pdata(nldi:nlei,nldj:nlej,jc)) 
    368368                     WRITE(numout,*) '****************' 
    369369                  ENDIF 
     
    444444                     WRITE(numout,*) 'oasis_get:   kstep', kstep 
    445445                     WRITE(numout,*) 'oasis_get:   info ', kinfo 
    446                      WRITE(numout,*) '     - Minimum value is ', MINVAL(pdata(:,:,jc)) 
    447                      WRITE(numout,*) '     - Maximum value is ', MAXVAL(pdata(:,:,jc)) 
    448                      WRITE(numout,*) '     -     Sum value is ', SUM(pdata(:,:,jc)) 
     446                     WRITE(numout,*) '     - Minimum value is ', MINVAL(pdata(nldi:nlei,nldj:nlej,jc)) 
     447                     WRITE(numout,*) '     - Maximum value is ', MAXVAL(pdata(nldi:nlei,nldj:nlej,jc)) 
     448                     WRITE(numout,*) '     -     Sum value is ',    SUM(pdata(nldi:nlei,nldj:nlej,jc)) 
    449449                     WRITE(numout,*) '****************' 
    450450                  ENDIF 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/ZDF/zdfiwm.F90

    r12495 r12551  
    2323   USE phycst         ! physical constants 
    2424   ! 
     25   USE fldread        ! field read 
    2526   USE prtctl         ! Print control 
    2627   USE in_out_manager ! I/O manager 
     
    383384      !!              de Lavergne et al. in prep., 2017 
    384385      !!---------------------------------------------------------------------- 
    385       INTEGER  ::   inum         ! local integer 
     386      INTEGER  ::   ifpr               ! dummy loop indices 
     387      INTEGER  ::   inum               ! local integer 
    386388      INTEGER  ::   ios 
    387389      REAL(wp) ::   zbot, zpyc, zcri   ! local scalars 
    388       !! 
    389       NAMELIST/namzdf_iwm/ nn_zpyc, ln_mevar, ln_tsdiff 
     390      ! 
     391      CHARACTER(len=256)            ::   cn_dir                 ! Root directory for location of ssr files 
     392      INTEGER, PARAMETER            ::   jpiwm  = 5             ! maximum number of files to read 
     393      INTEGER, PARAMETER            ::   jp_mpb = 1 
     394      INTEGER, PARAMETER            ::   jp_mpp = 2 
     395      INTEGER, PARAMETER            ::   jp_mpc = 3 
     396      INTEGER, PARAMETER            ::   jp_dsb = 4 
     397      INTEGER, PARAMETER            ::   jp_dsc = 5 
     398      ! 
     399      TYPE(FLD_N), DIMENSION(jpiwm) ::   slf_iwm                ! array of namelist informations 
     400      TYPE(FLD_N)                   ::   sn_mpb, sn_mpp, sn_mpc ! informations about Mixing Power field to be read 
     401      TYPE(FLD_N)                   ::   sn_dsb, sn_dsc         ! informations about Decay Scale field to be read 
     402      TYPE(FLD  ), DIMENSION(jpiwm) ::   sf_iwm                 ! structure of input fields (file informations, fields read) 
     403      ! 
     404      NAMELIST/namzdf_iwm/ nn_zpyc, ln_mevar, ln_tsdiff, & 
     405         &                 cn_dir, sn_mpb, sn_mpp, sn_mpc, sn_dsb, sn_dsc 
    390406      !!---------------------------------------------------------------------- 
    391407      ! 
     
    422438      IF( zdf_iwm_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'zdf_iwm_init : unable to allocate iwm arrays' ) 
    423439      ! 
     440      ! store namelist information in an array 
     441      slf_iwm(jp_mpb) = sn_mpb ; slf_iwm(jp_mpp) = sn_mpp ; slf_iwm(jp_mpc) = sn_mpc 
     442      slf_iwm(jp_dsb) = sn_dsb ; slf_iwm(jp_dsc) = sn_dsc 
     443      ! 
     444      DO ifpr= 1, jpiwm 
     445         ALLOCATE( sf_iwm(ifpr)%fnow(jpi,jpj,1)   ) 
     446         IF( slf_iwm(ifpr)%ln_tint )ALLOCATE( sf_iwm(ifpr)%fdta(jpi,jpj,1,2) ) 
     447      END DO 
     448 
     449      ! fill sf_iwm with sf_iwm and control print 
     450      CALL fld_fill( sf_iwm, slf_iwm , cn_dir, 'zdfiwm_init', 'iwm input file', 'namiwm' ) 
     451 
     452      !                             ! hard-coded default definition (to be defined in namelist ?) 
     453      sf_iwm(jp_mpb)%fnow(:,:,1) = 1.e-6 
     454      sf_iwm(jp_mpp)%fnow(:,:,1) = 1.e-6 
     455      sf_iwm(jp_mpc)%fnow(:,:,1) = 1.e-10 
     456      sf_iwm(jp_dsb)%fnow(:,:,1) = 100. 
     457      sf_iwm(jp_dsc)%fnow(:,:,1) = 100. 
     458 
    424459      !                             ! read necessary fields 
    425       CALL iom_open('mixing_power_bot',inum)       ! energy flux for high-mode wave breaking [W/m2] 
    426       CALL iom_get  (inum, jpdom_data, 'field', ebot_iwm, 1 )  
    427       CALL iom_close(inum) 
    428       ! 
    429       CALL iom_open('mixing_power_pyc',inum)       ! energy flux for pynocline-intensified wave breaking [W/m2] 
    430       CALL iom_get  (inum, jpdom_data, 'field', epyc_iwm, 1 ) 
    431       CALL iom_close(inum) 
    432       ! 
    433       CALL iom_open('mixing_power_cri',inum)       ! energy flux for critical slope wave breaking [W/m2] 
    434       CALL iom_get  (inum, jpdom_data, 'field', ecri_iwm, 1 ) 
    435       CALL iom_close(inum) 
    436       ! 
    437       CALL iom_open('decay_scale_bot',inum)        ! spatially variable decay scale for high-mode wave breaking [m] 
    438       CALL iom_get  (inum, jpdom_data, 'field', hbot_iwm, 1 ) 
    439       CALL iom_close(inum) 
    440       ! 
    441       CALL iom_open('decay_scale_cri',inum)        ! spatially variable decay scale for critical slope wave breaking [m] 
    442       CALL iom_get  (inum, jpdom_data, 'field', hcri_iwm, 1 ) 
    443       CALL iom_close(inum) 
    444  
    445       ebot_iwm(:,:) = ebot_iwm(:,:) * ssmask(:,:) 
    446       epyc_iwm(:,:) = epyc_iwm(:,:) * ssmask(:,:) 
    447       ecri_iwm(:,:) = ecri_iwm(:,:) * ssmask(:,:) 
     460      CALL fld_read( nit000, 1, sf_iwm ) 
     461 
     462      ebot_iwm(:,:) = sf_iwm(1)%fnow(:,:,1) * ssmask(:,:) ! energy flux for high-mode wave breaking [W/m2] 
     463      epyc_iwm(:,:) = sf_iwm(2)%fnow(:,:,1) * ssmask(:,:) ! energy flux for pynocline-intensified wave breaking [W/m2] 
     464      ecri_iwm(:,:) = sf_iwm(3)%fnow(:,:,1) * ssmask(:,:) ! energy flux for critical slope wave breaking [W/m2] 
     465      hbot_iwm(:,:) = sf_iwm(4)%fnow(:,:,1)               ! spatially variable decay scale for high-mode wave breaking [m] 
     466      hcri_iwm(:,:) = sf_iwm(5)%fnow(:,:,1)               ! spatially variable decay scale for critical slope wave breaking [m] 
    448467 
    449468      zbot = glob_sum( 'zdfiwm', e1e2t(:,:) * ebot_iwm(:,:) ) 
    450469      zpyc = glob_sum( 'zdfiwm', e1e2t(:,:) * epyc_iwm(:,:) ) 
    451470      zcri = glob_sum( 'zdfiwm', e1e2t(:,:) * ecri_iwm(:,:) ) 
     471 
    452472      IF(lwp) THEN 
    453473         WRITE(numout,*) '      High-mode wave-breaking energy:             ', zbot * 1.e-12_wp, 'TW' 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_cfg_orca025_like

    r12495 r12551  
    5050!----------------------------------------------------------------------- 
    5151  ln_usr      = .true.    !  user defined formulation                  (T => check usrdef_sbc) 
     52  nn_ice      = 2         !  =0 no ice boundary condition 
     53     !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     54     !                    !  =2 or 3 for SI3 and CICE, respectively 
    5255  ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    5356/ 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_cfg_orca12_like

    r12495 r12551  
    5050!----------------------------------------------------------------------- 
    5151  ln_usr      = .true.    !  user defined formulation                  (T => check usrdef_sbc) 
     52  nn_ice      = 2         !  =0 no ice boundary condition 
     53     !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     54     !                    !  =2 or 3 for SI3 and CICE, respectively 
    5255  ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    5356/ 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_cfg_orca1_like

    r12495 r12551  
    5050!----------------------------------------------------------------------- 
    5151  ln_usr      = .true.    !  user defined formulation                  (T => check usrdef_sbc) 
     52  nn_ice      = 2         !  =0 no ice boundary condition 
     53     !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     54     !                    !  =2 or 3 for SI3 and CICE, respectively 
    5255  ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    5356/ 
  • NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/BENCH/EXPREF/namelist_pisces_cfg

    r10343 r12551  
    8181/ 
    8282!----------------------------------------------------------------------- 
    83 &nampissbc     !   parameters for inputs deposition 
     83&nampisbc      !   parameters for inputs deposition 
    8484!----------------------------------------------------------------------- 
    85    ln_dust     =  .false.   ! boolean for dust input from the atmosphere 
    86    ln_solub    =  .false.   ! boolean for variable solubility of atm. Iron 
    87    ln_river    =  .false.   ! boolean for river input of nutrients 
    88    ln_ndepo    =  .false.   ! boolean for atmospheric deposition of N 
    8985   ln_ironsed  =  .false.   ! boolean for Fe input from sediments 
    9086   ln_ironice  =  .false.   ! boolean for Fe input from sea ice 
Note: See TracChangeset for help on using the changeset viewer.