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 1681 for trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90 – NEMO

Ignore:
Timestamp:
2009-10-27T18:44:08+01:00 (14 years ago)
Author:
smasson
Message:

slight bugfix in iodef.xml, see ticket:566

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90

    r1546 r1681  
    5353      !! References :   Lazar, A., these de l'universite Paris VI, France, 1997 
    5454      !!---------------------------------------------------------------------- 
    55       USE oce,   zav_evd  =>   ua   ! use ua as workspace 
     55      USE oce,   zavt_evd  =>   ua          ! use ua as workspace 
     56      USE oce,   zavm_evd  =>   va          ! use va as workspace 
    5657      !! 
    5758      INTEGER, INTENT( in ) ::   kt         ! ocean time-step indexocean time step 
     
    6768      ENDIF 
    6869 
    69       zav_evd(:,:,:) = avt(:,:,:)           ! set avt prior to evd application 
     70      zavt_evd(:,:,:) = avt(:,:,:)           ! set avt prior to evd application 
    7071 
    7172      SELECT CASE ( nn_evdm ) 
    7273  
    7374      CASE ( 1 )           ! enhance vertical eddy viscosity and diffusivity (if rn2<-1.e-12) 
     75 
     76         zavm_evd(:,:,:) = avm(:,:,:)           ! set avm prior to evd application 
     77         ! 
    7478         DO jk = 1, jpkm1  
    7579#if defined key_vectopt_loop 
     
    99103         CALL lbc_lnk( avmu, 'U', 1. )   ;   CALL lbc_lnk( avmv, 'V', 1. ) 
    100104         ! 
     105         zavm_evd(:,:,:) = avm(:,:,:) - zavm_evd(:,:,:)   ! change in avm due to evd 
     106         CALL iom_put( "avm_evd", zavm_evd )              ! output this change 
     107         ! 
    101108      CASE DEFAULT         ! enhance vertical eddy diffusivity only (if rn2<-1.e-12)  
    102109         DO jk = 1, jpkm1 
     
    122129      END SELECT  
    123130 
    124       zav_evd(:,:,:) = avt(:,:,:) - zav_evd(:,:,:)      ! change in avt due to evd 
    125       CALL iom_put( "avt_evd", zav_evd )                ! output this change 
     131      zavt_evd(:,:,:) = avt(:,:,:) - zavt_evd(:,:,:)   ! change in avt due to evd 
     132      CALL iom_put( "avt_evd", zavt_evd )              ! output this change 
    126133      ! 
    127134   END SUBROUTINE zdf_evd 
Note: See TracChangeset for help on using the changeset viewer.