Ignore:
Timestamp:
2021-02-08T21:55:49+01:00 (6 months ago)
Author:
techene
Message:

#2605 correct the truncature change of ww induced by diawri when ln_zad_Aimp=T along a time step

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DIA/diawri.F90

    r14239 r14418  
    250250         CALL iom_put( "sbv", z2d )                ! bottom j-current 
    251251      ENDIF 
    252  
    253       IF( ln_zad_Aimp ) ww = ww + wi               ! Recombine explicit and implicit parts of vertical velocity for diagnostic output 
    254       CALL iom_put( "woce", ww )                   ! vertical velocity 
    255  
     252       
     253      IF( ln_zad_Aimp ) THEN                       ! Recombine explicit and implicit parts of vertical velocity for diagnostic output 
     254         CALL iom_put( "woce", ww + wi )           ! vertical velocity 
     255      ELSE 
     256         CALL iom_put( "woce", ww )                ! vertical velocity 
     257      ENDIF 
     258       
    256259      IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN   ! vertical mass transport & its square value 
    257260         ! Caution: in the VVL case, it only correponds to the baroclinic mass transport. 
    258261         z2d(:,:) = rho0 * e1e2t(:,:) 
    259          DO jk = 1, jpk 
    260             z3d(:,:,jk) = ww(:,:,jk) * z2d(:,:) 
    261          END DO 
     262         IF( ln_zad_Aimp ) THEN 
     263            DO jk = 1, jpk 
     264               z3d(:,:,jk) = ( ww(:,:,jk) + wi(:,:,jk) ) * z2d(:,:) 
     265            END DO 
     266         ELSE 
     267            DO jk = 1, jpk 
     268               z3d(:,:,jk) = ww(:,:,jk) * z2d(:,:) 
     269            END DO 
     270         ENDIF 
    262271         CALL iom_put( "w_masstr" , z3d )   
    263272         IF( iom_use('w_masstr2') )   CALL iom_put( "w_masstr2", z3d(:,:,:) * z3d(:,:,:) ) 
    264273      ENDIF 
    265       ! 
    266       IF( ln_zad_Aimp ) ww = ww - wi               ! Remove implicit part of vertical velocity that was added for diagnostic output 
    267274 
    268275      CALL iom_put( "avt" , avt )                  ! T vert. eddy diff. coef. 
Note: See TracChangeset for help on using the changeset viewer.