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 1546 – NEMO

Changeset 1546


Ignore:
Timestamp:
2009-07-27T15:39:29+02:00 (15 years ago)
Author:
ctlod
Message:

small optimisation, see ticket: #495

Location:
trunk/NEMO
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/C1D_SRC/diawri_c1d.F90

    r1533 r1546  
    251251         CALL histdef( nid_T, "votkeavt", "Vertical Eddy Diffusivity"          , "m2/s"   ,   &  ! avt 
    252252            &          jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout ) 
    253          CALL histdef( nid_T, "votkeevd", "Enhanced Vertical Diffusivity",       "m2/s"   ,   &  ! avt_evd 
    254             &          jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout ) 
    255253         CALL histdef( nid_T, "votkeavm", "Vertical Eddy Viscosity",             "m2/s"   ,   &  ! avmu 
    256254            &          jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout ) 
    257255         ! 
    258          CALL histdef( nid_T, "votkeevm", "Enhanced Vertical Viscosity",         "m2/s"   ,   &  ! avmu_evd 
    259             &          jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout ) 
    260256         IF( lk_zdfddm ) THEN 
    261257            CALL histdef( nid_T,"voddmavs","Salt Vertical Eddy Diffusivity"    , "m2/s"   ,   &  ! avs 
    … …  
    338334      CALL histwrite( nid_T, "voeosbn2", it, rn2           , ndim_T , ndex_T )    ! Brunt-Vaisala Frequency 
    339335      CALL histwrite( nid_T, "votkeavt", it, avt           , ndim_T , ndex_T )    ! T vert. eddy diff. coef. 
    340       CALL histwrite( nid_T, "votkeevd", it, avt_evd       , ndim_T , ndex_T )    ! T enhan. vert. eddy diff. coef. 
    341336      CALL histwrite( nid_T, "votkeavm", it, avmu          , ndim_T , ndex_T )    ! T vert. eddy visc. coef. 
    342       CALL histwrite( nid_T, "votkeevm", it, avmu_evd      , ndim_T , ndex_T )    ! T enhan. vert. eddy visc. coef. 
    343337      IF( lk_zdfddm ) THEN 
    344338         CALL histwrite( nid_T, "voddmavs", it, fsavs(:,:,:), ndim_T, ndex_T )    ! S vert. eddy diff. coef. 

  • trunk/NEMO/OPA_SRC/DIA/diawri.F90

    r1528 r1546  
    400400         CALL histdef( nid_W, "votkeavt", "Vertical Eddy Diffusivity"          , "m2/s"   ,   &  ! avt 
    401401            &          jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout ) 
    402 # if defined key_zdftmx  
    403          CALL histdef( nid_W, "votidavt", "Tidal Vertical Eddy Diffusivity"    , "m2/s"   ,   &  ! avt_tide 
    404             &          jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout ) 
    405  
    406          CALL histdef( nid_W, "voitfavt", "Tidal ITF Vertical Eddy Diffusivity", "m2/s"   ,   &  ! avt_tide_itf 
    407             &          jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout ) 
    408 # endif 
    409          CALL histdef( nid_W, "votkeevd", "Enhanced Vertical Diffusivity",       "m2/s"   ,   &  ! avt_evd 
    410             &          jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout ) 
    411          ! 
    412402         CALL histdef( nid_W, "votkeavm", "Vertical Eddy Viscosity"             , "m2/s"  ,   &  ! avmu 
    413             &          jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout ) 
    414  
    415          CALL histdef( nid_W, "votkeevm", "Enhanced Vertical Viscosity",         "m2/s"   ,   &  ! avmu_evd 
    416403            &          jpi, jpj, nh_W, ipk, 1, ipk, nz_W, 32, clop, zsto, zout ) 
    417404 
    … …  
    525512#   endif 
    526513      CALL histwrite( nid_W, "votkeavt", it, avt            , ndim_T, ndex_T )    ! T vert. eddy diff. coef. 
    527 # if defined key_zdftmx  
    528       CALL histwrite( nid_W, "votidavt", it, av_tide        , ndim_T, ndex_T )    ! vert. mix. related to internal tides 
    529 # endif 
    530       CALL histwrite( nid_W, "votkeevd", it, avt_evd        , ndim_T, ndex_T )    ! T enhan. vert. eddy diff. coef. 
    531514      CALL histwrite( nid_W, "votkeavm", it, avmu           , ndim_T, ndex_T )    ! T vert. eddy visc. coef. 
    532       CALL histwrite( nid_W, "votkeevm", it, avmu_evd       , ndim_T, ndex_T )    ! T enhan. vert. eddy visc. coef. 
    533515      IF( lk_zdfddm ) THEN 
    534516         CALL histwrite( nid_W, "voddmavs", it, fsavs(:,:,:), ndim_T, ndex_T )    ! S vert. eddy diff. coef. 

  • trunk/NEMO/OPA_SRC/ZDF/zdf_oce.F90

    r1537 r1546  
    3737   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   avmu, avmv   !: vertical viscosity coeff. at uw- & vw-points   [m2/s] 
    3838   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   avm , avt    !: vertical viscosity & diffusivity coeff. at  w-point   [m2/s] 
    39    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   avt_evd      !: enhanced vertical diffusivity coeff. at  w-point   [m2/s] 
    40    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   avmu_evd     !: enhanced vertical viscosity   coeff. at uw-point   [m2/s] 
    41 #if defined key_zdftmx 
    42    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   av_tide      !: Tidal mixing 
    43 #endif 
    4439  
    4540   !!---------------------------------------------------------------------- 

  • trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90

    r1537 r1546  
    4545      !!      sivity coefficients when a static instability is encountered. 
    4646      !! 
    47       !! ** Method  :   avt, and the 4 neighbouring avmu, avmv coefficients 
     47      !! ** Method  :   avt, avm, and the 4 neighbouring avmu, avmv coefficients 
    4848      !!      are set to avevd (namelist parameter) if the water column is  
    4949      !!      statically unstable (i.e. if rn2 < -1.e-12 ) 
    5050      !! 
    51       !! ** Action  :   Update avt, avmu, avmv in statically instable cases 
    52       !!                and avt_evd which is avt due to convection 
     51      !! ** Action  :   avt, avm, avmu, avmv updted in static instability cases 
     52      !! 
    5353      !! References :   Lazar, A., these de l'universite Paris VI, France, 1997 
    5454      !!---------------------------------------------------------------------- 
     55      USE oce,   zav_evd  =>   ua   ! use ua as workspace 
     56      !! 
    5557      INTEGER, INTENT( in ) ::   kt         ! ocean time-step indexocean time step 
    5658      !! 
    … …  
    6567      ENDIF 
    6668 
    67       ! Initialisation of avt_evd (vertical diffusion due to convection) to avt and avmu_evd to avmu 
    68       avt_evd  (:,:,:) = avt(:,:,:)  
    69       avmu_evd (:,:,:) = avmu(:,:,:)  
     69      zav_evd(:,:,:) = avt(:,:,:)           ! set avt prior to evd application 
    7070 
    7171      SELECT CASE ( nn_evdm ) 
    7272  
    7373      CASE ( 1 )           ! enhance vertical eddy viscosity and diffusivity (if rn2<-1.e-12) 
    74          !                                                ! =============== 
    75          DO jk = 1, jpkm1                                 ! Horizontal slab 
    76             !                                             ! =============== 
     74         DO jk = 1, jpkm1  
    7775#if defined key_vectopt_loop 
    7876            DO jj = 1, 1                     ! big loop forced 
    … …  
    8381#endif 
    8482#if defined key_zdfkpp 
    85 !! no implicit mixing in the boundary layer with KPP 
     83                  ! no evd mixing in the boundary layer with KPP 
    8684                  IF( ( MIN( rn2(ji,jj,jk),  rn2b(ji,jj,jk) ) <= -1.e-12 ) .AND. ( fsdepw(ji,jj,jk) > hkpp(ji,jj) ) ) THEN 
    8785#else 
    … …  
    9795               END DO 
    9896            END DO 
    99             !                                             ! =============== 
    100          END DO                                           !   End of slab 
    101          !                                                ! =============== 
    102  
    103          ! Lateral boundary conditions on ( avt, avmu, avmv )   (unchanged sign) 
    104          ! -------------------------------=================== 
    105          CALL lbc_lnk( avt , 'W', 1. ) 
    106          CALL lbc_lnk( avm , 'W', 1. ) 
    107          CALL lbc_lnk( avmu, 'U', 1. ) 
    108          CALL lbc_lnk( avmv, 'V', 1. ) 
    109  
     97         END DO  
     98         CALL lbc_lnk( avt , 'W', 1. )   ;   CALL lbc_lnk( avm , 'W', 1. )   ! Lateral boundary conditions   (unchanged sign) 
     99         CALL lbc_lnk( avmu, 'U', 1. )   ;   CALL lbc_lnk( avmv, 'V', 1. ) 
     100         ! 
    110101      CASE DEFAULT         ! enhance vertical eddy diffusivity only (if rn2<-1.e-12)  
    111          !                                                ! =============== 
    112          DO jk = 1, jpkm1                                 ! Horizontal slab 
    113             !                                             ! =============== 
    114 !!!         WHERE( rn2(:,:,jk) <= -1.e-12 ) avt(:,:,jk) = tmask(:,:,jk) * rn_avevd   ! agissant sur T SEUL!  
     102         DO jk = 1, jpkm1 
     103!!!         WHERE( rn2(:,:,jk) <= -1.e-12 ) avt(:,:,jk) = tmask(:,:,jk) * avevd   ! agissant sur T SEUL!  
    115104#if defined key_vectopt_loop 
    116105            DO jj = 1, 1                     ! big loop forced 
    … …  
    121110#endif 
    122111#if defined key_zdfkpp 
    123 !! no implicit mixing in the boundary layer with KPP 
     112                  ! no evd mixing in the boundary layer with KPP 
    124113                  IF( ( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 ) .AND. ( fsdepw(ji,jj,jk) > hkpp(ji,jj) ) ) &           
    125114#else 
    … …  
    129118               END DO 
    130119            END DO 
    131             !                                             ! =============== 
    132          END DO                                           !   End of slab 
    133          !                                                ! =============== 
     120         END DO 
     121         ! 
    134122      END SELECT  
    135123 
    136       ! update of avt_evd and avmu_evd 
    137       avt_evd (:,:,:) = avt (:,:,:) - avt_evd (:,:,:)  
    138       avmu_evd(:,:,:) = avmu(:,:,:) - avmu_evd(:,:,:)  
    139       CALL iom_put( "avt_evd", avt_evd  )    ! T enhan. vert. eddy diff. coef. 
    140       CALL iom_put( "avm_evd", avmu_evd )    ! T enhan. vert. eddy visc. coef. 
    141  
     124      zav_evd(:,:,:) = avt(:,:,:) - zav_evd(:,:,:)      ! change in avt due to evd 
     125      CALL iom_put( "avt_evd", zav_evd )                ! output this change 
     126      ! 
    142127   END SUBROUTINE zdf_evd 
    143128 

  • trunk/NEMO/OPA_SRC/ZDF/zdftmx.F90

    r1537 r1546  
    8888      !!              Koch-Larrouy et al. 2007, GRL. 
    8989      !!---------------------------------------------------------------------- 
     90      USE oce,   zav_tide  =>   ua   ! use ua as workspace 
     91      !! 
    9092      INTEGER, INTENT(in) ::   kt   ! ocean time-step  
    9193      !! 
    … …  
    98100      IF( kt == nit000  )   CALL zdf_tmx_init      ! Initialization (first time-step only) 
    99101 
    100       !                                        ! ----------------------- ! 
    101       !                                        !  Standard tidal mixing  !  (compute av_tide) 
    102       !                                        ! ----------------------- ! 
     102      !                          ! ----------------------- ! 
     103      !                          !  Standard tidal mixing  !  (compute zav_tide) 
     104      !                          ! ----------------------- ! 
    103105      !                             !* First estimation (with n2 bound by rn_n2min) bounded by 60 cm2/s 
    104       av_tide(:,:,:) = MIN(  60.e-4, az_tmx(:,:,:) / MAX( rn_n2min, rn2(:,:,:) )  ) 
     106      zav_tide(:,:,:) = MIN(  60.e-4, az_tmx(:,:,:) / MAX( rn_n2min, rn2(:,:,:) )  ) 
    105107 
    106108      zkz(:,:) = 0.e0               !* Associated potential energy consummed over the whole water column 
    107109      DO jk = 2, jpkm1 
    108          zkz(:,:) = zkz(:,:) + fse3w(:,:,jk) * MAX( 0.e0, rn2(:,:,jk) ) * rau0 * av_tide(:,:,jk)* tmask(:,:,jk) 
     110         zkz(:,:) = zkz(:,:) + fse3w(:,:,jk) * MAX( 0.e0, rn2(:,:,jk) ) * rau0 * zav_tide(:,:,jk)* tmask(:,:,jk) 
    109111      END DO 
    110112 
    … …  
    115117      END DO 
    116118 
    117       DO jk = 2, jpkm1              !* Mutiply by zkz to recover en_tmx, BUT bound by 30/6 ==> av_tide bound by 300 cm2/s 
    118          av_tide(:,:,jk) = av_tide(:,:,jk) * MIN( zkz(:,:), 30./6. )   !kz max = 300 cm2/s 
    119       END DO 
    120  
    121       IF( kt == nit000 ) THEN       !* check at first time-step: diagnose the energy consumed by av_tide 
     119      DO jk = 2, jpkm1              !* Mutiply by zkz to recover en_tmx, BUT bound by 30/6 ==> zav_tide bound by 300 cm2/s 
     120         zav_tide(:,:,jk) = zav_tide(:,:,jk) * MIN( zkz(:,:), 30./6. )   !kz max = 300 cm2/s 
     121      END DO 
     122 
     123      IF( kt == nit000 ) THEN       !* check at first time-step: diagnose the energy consumed by zav_tide 
    122124         ztpc = 0.e0 
    123125         DO jk= 1, jpk 
    … …  
    125127               DO ji= 1, jpi 
    126128                  ztpc = ztpc + fse3w(ji,jj,jk) * e1t(ji,jj) * e2t(ji,jj)   & 
    127                      &         * MAX( 0.e0, rn2(ji,jj,jk) ) * av_tide(ji,jj,jk) * tmask(ji,jj,jk) * tmask_i(ji,jj) 
     129                     &         * MAX( 0.e0, rn2(ji,jj,jk) ) * zav_tide(ji,jj,jk) * tmask(ji,jj,jk) * tmask_i(ji,jj) 
    128130               END DO 
    129131            END DO 
    … …  
    134136      ENDIF 
    135137        
    136       !                                        ! ----------------------- ! 
    137       IF( ln_tmx_itf )   CALL tmx_itf( kt )    !    ITF  tidal mixing    !  (update av_tide) 
    138       !                                        ! ----------------------- ! 
    139  
    140       !                                        ! ----------------------- ! 
    141       !                                        !   Update  mixing coefs  !                           
    142       !                                        ! ----------------------- ! 
     138      !                          ! ----------------------- ! 
     139      !                          !    ITF  tidal mixing    !  (update zav_tide) 
     140      !                          ! ----------------------- ! 
     141      IF( ln_tmx_itf )   CALL tmx_itf( kt, zav_tide ) 
     142 
     143      !                          ! ----------------------- ! 
     144      !                          !   Update  mixing coefs  !                           
     145      !                          ! ----------------------- ! 
    143146      DO jk = 2, jpkm1              !* update momentum & tracer diffusivity with tidal mixing 
    144          avt(:,:,jk) = avt(:,:,jk) + av_tide(:,:,jk) 
    145          avm(:,:,jk) = avm(:,:,jk) + av_tide(:,:,jk) 
     147         avt(:,:,jk) = avt(:,:,jk) + zav_tide(:,:,jk) 
     148         avm(:,:,jk) = avm(:,:,jk) + zav_tide(:,:,jk) 
    146149         DO jj = 2, jpjm1 
    147150            DO ji = fs_2, fs_jpim1  ! vector opt. 
    148                avmu(ji,jj,jk) = avmu(ji,jj,jk) + 0.5 * ( av_tide(ji,jj,jk) + av_tide(ji+1,jj  ,jk) ) * umask(ji,jj,jk) 
    149                avmv(ji,jj,jk) = avmv(ji,jj,jk) + 0.5 * ( av_tide(ji,jj,jk) + av_tide(ji  ,jj+1,jk) ) * vmask(ji,jj,jk) 
     151               avmu(ji,jj,jk) = avmu(ji,jj,jk) + 0.5 * ( zav_tide(ji,jj,jk) + zav_tide(ji+1,jj  ,jk) ) * umask(ji,jj,jk) 
     152               avmv(ji,jj,jk) = avmv(ji,jj,jk) + 0.5 * ( zav_tide(ji,jj,jk) + zav_tide(ji  ,jj+1,jk) ) * vmask(ji,jj,jk) 
    150153            END DO 
    151154         END DO 
    … …  
    153156      CALL lbc_lnk( avmu, 'U', 1. )   ;   CALL lbc_lnk( avmv, 'V', 1. )      ! lateral boundary condition 
    154157 
    155       IF(ln_ctl)   CALL prt_ctl(tab3d_1=av_tide ,clinfo1=' tmx - av_tide: ',tab3d_2=avt,clinfo2=' avt: ',ovlap=1,kdim=jpk) 
     158      !                             !* output tidal mixing coefficient 
     159      CALL iom_put( "av_tide", zav_tide ) 
     160 
     161      IF(ln_ctl)   CALL prt_ctl(tab3d_1=zav_tide , clinfo1=' tmx - av_tide: ', tab3d_2=avt, clinfo2=' avt: ', ovlap=1, kdim=jpk) 
    156162      ! 
    157163   END SUBROUTINE zdf_tmx 
    158164 
    159165 
    160    SUBROUTINE tmx_itf( kt ) 
     166   SUBROUTINE tmx_itf( kt, pav ) 
    161167      !!---------------------------------------------------------------------- 
    162168      !!                  ***  ROUTINE tmx_itf  *** 
    163169      !!                    
    164170      !! ** Purpose :   modify the vertical eddy diffusivity coefficients  
    165       !!              (av_tide) in the Indonesian Through Flow area (ITF). 
     171      !!              (pav) in the Indonesian Through Flow area (ITF). 
    166172      !! 
    167173      !! ** Method  : - Following Koch-Larrouy et al. (2007), in the ITF defined 
    … …  
    178184      !! References :  Koch-Larrouy et al. 2007, GRL  
    179185      !!---------------------------------------------------------------------- 
    180       INTEGER, INTENT(in) ::   kt   ! ocean time-step 
     186      INTEGER , INTENT(in   )                         ::   kt   ! ocean time-step 
     187      REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk) ::   pav  ! Tidal mixing coef. 
    181188      !!  
    182189      INTEGER  ::   ji, jj, jk    ! dummy loop indices 
    … …  
    267274      ENDIF 
    268275 
    269       !                             ! Update av_tide with the ITF mixing coefficient 
     276      !                             ! Update pav with the ITF mixing coefficient 
    270277      DO jk = 2, jpkm1 
    271          av_tide(:,:,jk) = av_tide (:,:,jk) * ( 1.e0 - mask_itf(:,:) )   & 
    272             &            + zavt_itf(:,:,jk) *          mask_itf(:,:)  
     278         pav(:,:,jk) = pav    (:,:,jk) * ( 1.e0 - mask_itf(:,:) )   & 
     279            &        + zavt_itf(:,:,jk) *          mask_itf(:,:)  
    273280      END DO 
    274281      ! 
    … …  
    311318      !!              Koch-Larrouy et al. 2007, GRL. 
    312319      !!---------------------------------------------------------------------- 
     320      USE oce,   zav_tide  =>   ua   ! use ua as workspace 
     321      !! 
    313322      INTEGER ::   ji, jj, jk    ! dummy loop indices 
    314323      INTEGER ::   inum          ! temporary logical unit 
    … …  
    379388         ! Control print 
    380389         ! Total power consumption due to vertical mixing 
    381          ! zpc = rau0 * 1/rn_me * rn2 * av_tide 
    382          av_tide(:,:,:) = 0.e0 
     390         ! zpc = rau0 * 1/rn_me * rn2 * zav_tide 
     391         zav_tide(:,:,:) = 0.e0 
    383392         DO jk = 2, jpkm1 
    384             av_tide(:,:,jk) = az_tmx(:,:,jk) / MAX( rn_n2min, rn2(:,:,jk) ) 
     393            zav_tide(:,:,jk) = az_tmx(:,:,jk) / MAX( rn_n2min, rn2(:,:,jk) ) 
    385394         END DO 
    386395 
    387396         ztpc = 0.e0 
    388          zpc(:,:,:) = MAX(rn_n2min,rn2(:,:,:)) * av_tide(:,:,:) 
     397         zpc(:,:,:) = MAX(rn_n2min,rn2(:,:,:)) * zav_tide(:,:,:) 
    389398         DO jk= 2, jpkm1 
    390399            DO jj = 1, jpj 
    … …  
    401410 
    402411         ! control print 2 
    403          av_tide(:,:,:) = MIN( av_tide(:,:,:), 60.e-4 )    
     412         zav_tide(:,:,:) = MIN( zav_tide(:,:,:), 60.e-4 )    
    404413         zkz(:,:) = 0.e0 
    405414         DO jk = 2, jpkm1 
    406415         DO jj = 1, jpj 
    407416            DO ji = 1, jpi 
    408                zkz(ji,jj) = zkz(ji,jj) + fse3w(ji,jj,jk) * MAX( 0.e0, rn2(ji,jj,jk) ) * rau0 * av_tide(ji,jj,jk)* tmask(ji,jj,jk) 
     417               zkz(ji,jj) = zkz(ji,jj) + fse3w(ji,jj,jk) * MAX( 0.e0, rn2(ji,jj,jk) ) * rau0 * zav_tide(ji,jj,jk)* tmask(ji,jj,jk) 
    409418            END DO 
    410419         END DO 
    … …  
    429438 
    430439         DO jk = 2, jpkm1 
    431             av_tide(:,:,jk) = av_tide(:,:,jk) * MIN( zkz(:,:), 30./6. )   !kz max = 300 cm2/s 
     440            zav_tide(:,:,jk) = zav_tide(:,:,jk) * MIN( zkz(:,:), 30./6. )   !kz max = 300 cm2/s 
    432441         END DO 
    433442         ztpc = 0.e0 
    434          zpc(:,:,:) = Max(0.e0,rn2(:,:,:)) * av_tide(:,:,:) 
     443         zpc(:,:,:) = Max(0.e0,rn2(:,:,:)) * zav_tide(:,:,:) 
    435444         DO jk= 1, jpk 
    436445            DO jj = 1, jpj 
    … …  
    444453 
    445454         DO jk = 1, jpk 
    446             ze_z =                  SUM( e1t(:,:) * e2t(:,:) * av_tide(:,:,jk)     * tmask_i(:,:) )   & 
     455            ze_z =                  SUM( e1t(:,:) * e2t(:,:) * zav_tide(:,:,jk)     * tmask_i(:,:) )   & 
    447456               &     / MAX( 1.e-20, SUM( e1t(:,:) * e2t(:,:) * tmask (:,:,jk) * tmask_i(:,:) ) ) 
    448457            ztpc = 1.E50 
    449458            DO jj = 1, jpj 
    450459               DO ji = 1, jpi 
    451                   IF( av_tide(ji,jj,jk) /= 0.e0 )   ztpc =Min( ztpc, av_tide(ji,jj,jk) ) 
     460                  IF( zav_tide(ji,jj,jk) /= 0.e0 )   ztpc =Min( ztpc, zav_tide(ji,jj,jk) ) 
    452461               END DO 
    453462            END DO 
    454463            WRITE(numout,*) '            N2 min - jk= ', jk,'   ', ze_z * 1.e4,' cm2/s min= ',ztpc*1.e4,   & 
    455                &       'max= ', MAXVAL(av_tide(:,:,jk) )*1.e4, ' cm2/s' 
     464               &       'max= ', MAXVAL(zav_tide(:,:,jk) )*1.e4, ' cm2/s' 
    456465         END DO 
    457466 

  • trunk/NEMO/OPA_SRC/step.F90

    r1537 r1546  
    229229      ! 
    230230      IF( lk_ldfslp ) THEN                            ! slope of lateral mixing 
    231                          CALL eos( tb, sb, rhd, rhop )          ! before in situ density 
     231                         CALL eos( tb, sb, rhd )                ! before in situ density 
    232232         IF( ln_zps )    CALL zps_hde( kstp, tb, sb, rhd,  &    ! Partial steps: before horizontal gradient 
    233233            &                              gtu, gsu, gru,  &    ! of t, s, rd at the last ocean level 
Note: See TracChangeset for help on using the changeset viewer.