Ignore:
Timestamp:
2020-04-09T21:06:01+02:00 (11 months ago)
Author:
techene
Message:

some cleaning and proper module/routine name, mini bug introduced and corrected in sbcice_cice

File:
1 moved

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/traatf_qco.F90

    r12731 r12732  
    5353 
    5454   PUBLIC   tra_atf_qco       ! routine called by step.F90 
    55    PUBLIC   tra_atf_fix_lf   ! to be used in trcnxt !!st WARNING discrepancy here interpol is used 
    56    PUBLIC   tra_atf_qco_lf    ! to be used in trcnxt !!st WARNING discrepancy here interpol is used 
     55   PUBLIC   tra_atf_fix_lf    ! to be used in trcnxt !!st WARNING discrepancy here interpol is used by PISCES 
     56   PUBLIC   tra_atf_qco_lf    ! to be used in trcnxt !!st WARNING discrepancy here interpol is used by PISCES 
    5757 
    5858   !! * Substitutions 
     
    103103         IF(lwp) WRITE(numout,*) '~~~~~~~' 
    104104      ENDIF 
    105  
    106 !       ! Update after tracer on domain lateral boundaries 
    107 !       ! 
    108 ! #if defined key_agrif 
    109 !       CALL Agrif_tra                     ! AGRIF zoom boundaries 
    110 ! #endif 
    111 !       !                                              ! local domain boundaries  (T-point, unchanged sign) 
    112 !       CALL lbc_lnk_multi( 'traatfLF', pts(:,:,:,jp_tem,Kaa), 'T', 1., pts(:,:,:,jp_sal,Kaa), 'T', 1. ) 
    113 !       ! 
    114 !       IF( ln_bdy )   CALL bdy_tra( kt, Kbb, pts, Kaa )  ! BDY open boundaries 
     105!!st  Update after tracer on domain lateral boundaries as been removed outside  
    115106 
    116107      ! trends computation initialisation 
     
    127118         ! G Nurser 23 Mar 2017. Recalculate trend as Delta(e3t*T)/e3tn; e3tn cancel from pts(Kmm) terms 
    128119         DO jk = 1, jpkm1 
    129             ! ztrdt(:,:,jk) = ( pts(:,:,jk,jp_tem,Kaa)*e3t(:,:,jk,Kaa) / e3t(:,:,jk,Kmm) - pts(:,:,jk,jp_tem,Kmm)) * zfact 
    130             ! ztrds(:,:,jk) = ( pts(:,:,jk,jp_sal,Kaa)*e3t(:,:,jk,Kaa) / e3t(:,:,jk,Kmm) - pts(:,:,jk,jp_sal,Kmm)) * zfact 
    131120            ztrdt(:,:,jk) = ( pts(:,:,jk,jp_tem,Kaa) * (1._wp + r3t(:,:,Kaa) * tmask(:,:,jk))/(1._wp + r3t(:,:,Kmm) * tmask(:,:,jk))  & 
    132121               &            - pts(:,:,jk,jp_tem,Kmm) ) * zfact 
     
    251240      INTEGER  ::   ji, jj, jk, jn              ! dummy loop indices 
    252241      REAL(wp) ::   zfact, zfact1, ztc_a , ztc_n , ztc_b , ztc_f , ztc_d    ! local scalar 
    253       REAL(wp) ::   zfact2, ze3t_b, ze3t_n, ze3t_a, ze3t_f, ze3t_d  !   -      - 
     242      REAL(wp) ::   zfact2, ze3t_b, ze3t_n, ze3t_a, ze3t_f                  !   -      - 
    254243      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) ::   ztrd_atf 
    255244      !!---------------------------------------------------------------------- 
     
    288277            ztc_a  = pt(ji,jj,jk,jn,Kaa) * ze3t_a 
    289278            ! 
    290             ze3t_d = ze3t_a - 2. * ze3t_n + ze3t_b 
    291279            ztc_d  = ztc_a  - 2. * ztc_n  + ztc_b 
    292280            ! 
Note: See TracChangeset for help on using the changeset viewer.