Ignore:
Timestamp:
2019-05-15T21:19:35+02:00 (18 months ago)
Author:
acc
Message:

2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps : interface changes to tra and trc routines for design compliance and consistency. Fully SETTE tested (non-AGRIF, only)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/TRA/traqsr.F90

    r10954 r10985  
    7575CONTAINS 
    7676 
    77    SUBROUTINE tra_qsr( kt, Kmm, Krhs ) 
     77   SUBROUTINE tra_qsr( kt, Kmm, pts, Krhs ) 
    7878      !!---------------------------------------------------------------------- 
    7979      !!                  ***  ROUTINE tra_qsr  *** 
     
    101101      !!              Morel, A. et Berthon, JF, 1989, Limnol Oceanogr 34(8), 1545-1562 
    102102      !!---------------------------------------------------------------------- 
    103       INTEGER, INTENT(in) ::   kt     ! ocean time-step 
    104       INTEGER, INTENT(in) ::   Kmm, Krhs     ! time level indices 
     103      INTEGER,                                   INTENT(in   ) :: kt            ! ocean time-step 
     104      INTEGER,                                   INTENT(in   ) :: Kmm, Krhs     ! time level indices 
     105      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts,jpt), INTENT(inout) :: pts           ! active tracers and RHS of tracer equation 
    105106      ! 
    106107      INTEGER  ::   ji, jj, jk               ! dummy loop indices 
     
    127128      IF( l_trdtra ) THEN      ! trends diagnostic: save the input temperature trend 
    128129         ALLOCATE( ztrdt(jpi,jpj,jpk) )  
    129          ztrdt(:,:,:) = ts(:,:,:,jp_tem,Krhs) 
     130         ztrdt(:,:,:) = pts(:,:,:,jp_tem,Krhs) 
    130131      ENDIF 
    131132      ! 
     
    262263         DO jj = 2, jpjm1        !-----------------------------! 
    263264            DO ji = fs_2, fs_jpim1   ! vector opt. 
    264                ts(ji,jj,jk,jp_tem,Krhs) = ts(ji,jj,jk,jp_tem,Krhs)   & 
    265                   &                 + z1_2 * ( qsr_hc_b(ji,jj,jk) + qsr_hc(ji,jj,jk) ) / e3t(ji,jj,jk,Kmm) 
     265               pts(ji,jj,jk,jp_tem,Krhs) = pts(ji,jj,jk,jp_tem,Krhs)   & 
     266                  &                      + z1_2 * ( qsr_hc_b(ji,jj,jk) + qsr_hc(ji,jj,jk) ) / e3t(ji,jj,jk,Kmm) 
    266267            END DO 
    267268         END DO 
     
    296297      ! 
    297298      IF( l_trdtra ) THEN     ! qsr tracers trends saved for diagnostics 
    298          ztrdt(:,:,:) = ts(:,:,:,jp_tem,Krhs) - ztrdt(:,:,:) 
     299         ztrdt(:,:,:) = pts(:,:,:,jp_tem,Krhs) - ztrdt(:,:,:) 
    299300         CALL trd_tra( kt, Kmm, Krhs, 'TRA', jp_tem, jptra_qsr, ztrdt ) 
    300301         DEALLOCATE( ztrdt )  
    301302      ENDIF 
    302303      !                       ! print mean trends (used for debugging) 
    303       IF(ln_ctl)   CALL prt_ctl( tab3d_1=ts(:,:,:,jp_tem,Krhs), clinfo1=' qsr  - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 
     304      IF(ln_ctl)   CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' qsr  - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 
    304305      ! 
    305306      IF( ln_timing )   CALL timing_stop('tra_qsr') 
Note: See TracChangeset for help on using the changeset viewer.