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

Changeset 2904


Ignore:
Timestamp:
2011-10-12T10:00:57+02:00 (13 years ago)
Author:
kpedwards
Message:

Added code for e3?_temp - needed for averaging.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/TOP_SRC/trcsub.F90

    r2892 r2904  
    122122         END DO 
    123123#elif defined key_traldf_c2d 
    124             DO jj = 1, jpj 
    125                DO ji = 1, jpi 
    126                   ahtt_tm(ji,jj)= ahtt_tm(ji,jj)+ahtt(ji,jj) 
    127                   ahtu_tm(ji,jj)= ahtu_tm(ji,jj)+ahtu(ji,jj) 
    128                   ahtv_tm(ji,jj)= ahtv_tm(ji,jj)+ahtv(ji,jj) 
    129                   ahtw_tm(ji,jj)= ahtw_tm(ji,jj)+ahtw(ji,jj) 
    130                END DO 
    131             END DO 
     124         DO jj = 1, jpj 
     125            DO ji = 1, jpi 
     126               ahtt_tm(ji,jj)= ahtt_tm(ji,jj)+ahtt(ji,jj) 
     127               ahtu_tm(ji,jj)= ahtu_tm(ji,jj)+ahtu(ji,jj) 
     128               ahtv_tm(ji,jj)= ahtv_tm(ji,jj)+ahtv(ji,jj) 
     129               ahtw_tm(ji,jj)= ahtw_tm(ji,jj)+ahtw(ji,jj) 
     130            END DO 
     131         END DO 
    132132#elif defined key_traldf_c1d 
    133           DO jk = 1, jpk 
    134              ahtt_tm(jk)= ahtt_tm(jk)+ahtt(jk) 
    135              ahtu_tm(jk)= ahtu_tm(jk)+ahtu(jk) 
    136              ahtv_tm(jk)= ahtv_tm(jk)+ahtv(jk) 
    137              ahtw_tm(jk)= ahtw_tm(jk)+ahtw(jk) 
    138           END DO 
     133         DO jk = 1, jpk 
     134            ahtt_tm(jk)= ahtt_tm(jk)+ahtt(jk) 
     135            ahtu_tm(jk)= ahtu_tm(jk)+ahtu(jk) 
     136            ahtv_tm(jk)= ahtv_tm(jk)+ahtv(jk) 
     137            ahtw_tm(jk)= ahtw_tm(jk)+ahtw(jk) 
     138         END DO 
    139139#else 
    140              ahtt_tm= ahtt_tm+ahtt 
    141              ahtu_tm= ahtu_tm+ahtu 
    142              ahtv_tm= ahtv_tm+ahtv 
    143              ahtw_tm= ahtw_tm+ahtw 
     140         ahtt_tm= ahtt_tm+ahtt 
     141         ahtu_tm= ahtu_tm+ahtu 
     142         ahtv_tm= ahtv_tm+ahtv 
     143         ahtw_tm= ahtw_tm+ahtw 
    144144#endif 
    145145#if defined key_traldf_eiv 
    … …  
    155155         END DO 
    156156# elif defined key_traldf_c2d 
    157             DO jj = 1, jpj 
    158                DO ji = 1, jpi 
    159                   aeiu_tm(ji,jj)= aeiu_tm(ji,jj)+aeiu(ji,jj) 
    160                   aeiv_tm(ji,jj)= aeiv_tm(ji,jj)+aeiv(ji,jj) 
    161                   aeiw_tm(ji,jj)= aeiw_tm(ji,jj)+aeiw(ji,jj) 
    162                END DO 
    163             END DO 
     157         DO jj = 1, jpj 
     158            DO ji = 1, jpi 
     159               aeiu_tm(ji,jj)= aeiu_tm(ji,jj)+aeiu(ji,jj) 
     160               aeiv_tm(ji,jj)= aeiv_tm(ji,jj)+aeiv(ji,jj) 
     161               aeiw_tm(ji,jj)= aeiw_tm(ji,jj)+aeiw(ji,jj) 
     162            END DO 
     163         END DO 
    164164# elif defined key_traldf_c1d 
    165           DO jk = 1, jpk 
    166              aeiu_tm(jk)= aeiu_tm(jk)+aeiu(jk) 
    167              aeiv_tm(jk)= aeiv_tm(jk)+aeiv(jk) 
    168              aeiw_tm(jk)= aeiw_tm(jk)+aeiw(jk) 
    169           END DO 
     165         DO jk = 1, jpk 
     166            aeiu_tm(jk)= aeiu_tm(jk)+aeiu(jk) 
     167            aeiv_tm(jk)= aeiv_tm(jk)+aeiv(jk) 
     168            aeiw_tm(jk)= aeiw_tm(jk)+aeiw(jk) 
     169         END DO 
    170170# else 
    171              aeiu_tm= aeiu_tm+aeiu 
    172              aeiv_tm= aeiv_tm+aeiv 
    173              aeiw_tm= aeiw_tm+aeiw 
     171         aeiu_tm= aeiu_tm+aeiu 
     172         aeiv_tm= aeiv_tm+aeiv 
     173         aeiw_tm= aeiw_tm+aeiw 
    174174# endif 
    175175# endif 
    … …  
    178178!   it is time to substep 
    179179!   1. set temporary arrays to hold physics variables 
    180             sshn_temp(:,:)=sshn(:,:) 
    181             sshu_n_temp(:,:)=sshu_n(:,:) 
    182             sshv_n_temp(:,:)=sshv_n(:,:) 
    183             sshf_n_temp(:,:)=sshf_n(:,:) 
    184             sshb_temp(:,:)=sshb(:,:) 
    185             sshu_b_temp(:,:)=sshu_b(:,:) 
    186             sshv_b_temp(:,:)=sshv_b(:,:) 
    187             ssha_temp(:,:)=ssha(:,:) 
    188             sshu_a_temp(:,:)=sshu_a(:,:) 
    189             sshv_a_temp(:,:)=sshv_a(:,:) 
    190             rnf_temp(:,:)=rnf(:,:) 
    191             h_rnf_temp(:,:)=h_rnf(:,:) 
    192             rotn_temp(:,:,:)=rotn(:,:,:) 
    193             hdivn_temp(:,:,:)=hdivn(:,:,:) 
    194             rotb_temp(:,:,:)=rotb(:,:,:) 
    195             hdivb_temp(:,:,:)=hdivb(:,:,:) 
    196             un_temp(:,:,:)=un(:,:,:) 
    197             vn_temp(:,:,:)=vn(:,:,:) 
    198             wn_temp(:,:,:)=wn(:,:,:) 
    199             tsn_temp(:,:,:,:)=tsn(:,:,:,:) 
    200             avt_temp(:,:,:)=avt(:,:,:) 
    201             hmld_temp(:,:)=hmld(:,:) 
    202             qsr_temp(:,:)=qsr(:,:) 
    203             wndm_temp(:,:)=wndm(:,:) 
    204             emp_temp(:,:)=emp(:,:) 
    205             emp_b_temp(:,:)=emp_b(:,:) 
    206             emps_temp(:,:)=emps(:,:) 
     180         sshn_temp(:,:)=sshn(:,:) 
     181         sshu_n_temp(:,:)=sshu_n(:,:) 
     182         sshv_n_temp(:,:)=sshv_n(:,:) 
     183         sshf_n_temp(:,:)=sshf_n(:,:) 
     184         sshb_temp(:,:)=sshb(:,:) 
     185         sshu_b_temp(:,:)=sshu_b(:,:) 
     186         sshv_b_temp(:,:)=sshv_b(:,:) 
     187         ssha_temp(:,:)=ssha(:,:) 
     188         sshu_a_temp(:,:)=sshu_a(:,:) 
     189         sshv_a_temp(:,:)=sshv_a(:,:) 
     190         rnf_temp(:,:)=rnf(:,:) 
     191         h_rnf_temp(:,:)=h_rnf(:,:) 
     192         rotn_temp(:,:,:)=rotn(:,:,:) 
     193         hdivn_temp(:,:,:)=hdivn(:,:,:) 
     194         rotb_temp(:,:,:)=rotb(:,:,:) 
     195         hdivb_temp(:,:,:)=hdivb(:,:,:) 
     196         un_temp(:,:,:)=un(:,:,:) 
     197         vn_temp(:,:,:)=vn(:,:,:) 
     198         wn_temp(:,:,:)=wn(:,:,:) 
     199         tsn_temp(:,:,:,:)=tsn(:,:,:,:) 
     200         avt_temp(:,:,:)=avt(:,:,:) 
     201         hmld_temp(:,:)=hmld(:,:) 
     202         qsr_temp(:,:)=qsr(:,:) 
     203         wndm_temp(:,:)=wndm(:,:) 
     204         emp_temp(:,:)=emp(:,:) 
     205         emp_b_temp(:,:)=emp_b(:,:) 
     206         emps_temp(:,:)=emps(:,:) 
    207207#if defined key_traldf_c3d 
    208             ahtu_temp(:,:,:)=ahtu(:,:,:) 
    209             ahtv_temp(:,:,:)=ahtv(:,:,:) 
    210             ahtw_temp(:,:,:)=ahtw(:,:,:) 
    211             ahtt_temp(:,:,:)=ahtt(:,:,:) 
     208         ahtu_temp(:,:,:)=ahtu(:,:,:) 
     209         ahtv_temp(:,:,:)=ahtv(:,:,:) 
     210         ahtw_temp(:,:,:)=ahtw(:,:,:) 
     211         ahtt_temp(:,:,:)=ahtt(:,:,:) 
    212212#elif defined key_traldf_c2d 
    213             ahtu_temp(:,:)=ahtu(:,:) 
    214             ahtv_temp(:,:)=ahtv(:,:) 
    215             ahtw_temp(:,:)=ahtw(:,:) 
    216             ahtt_temp(:,:)=ahtt(:,:) 
     213         ahtu_temp(:,:)=ahtu(:,:) 
     214         ahtv_temp(:,:)=ahtv(:,:) 
     215         ahtw_temp(:,:)=ahtw(:,:) 
     216         ahtt_temp(:,:)=ahtt(:,:) 
    217217#elif defined key_traldf_c1d 
    218             ahtu_temp(:)=ahtu(:) 
    219             ahtv_temp(:)=ahtv(:) 
    220             ahtw_temp(:)=ahtw(:) 
    221             ahtt_temp(:)=ahtt(:) 
     218         ahtu_temp(:)=ahtu(:) 
     219         ahtv_temp(:)=ahtv(:) 
     220         ahtw_temp(:)=ahtw(:) 
     221         ahtt_temp(:)=ahtt(:) 
    222222#else 
    223             ahtu_temp=ahtu 
    224             ahtv_temp=ahtv 
    225             ahtw_temp=ahtw 
    226             ahtt_temp=ahtt 
     223         ahtu_temp=ahtu 
     224         ahtv_temp=ahtv 
     225         ahtw_temp=ahtw 
     226         ahtt_temp=ahtt 
    227227#endif 
    228228# if defined key_zdfddm 
    229             avs_temp(:,:,:)=avs(:,:,:) 
     229         avs_temp(:,:,:)=avs(:,:,:) 
    230230# endif 
    231231#if defined key_traldf_eiv 
    232232#if defined key_traldf_c3d 
    233             aeiu_temp(:,:,:)=aeiu(:,:,:) 
    234             aeiv_temp(:,:,:)=aeiv(:,:,:) 
    235             aeiw_temp(:,:,:)=aeiw(:,:,:) 
     233         aeiu_temp(:,:,:)=aeiu(:,:,:) 
     234         aeiv_temp(:,:,:)=aeiv(:,:,:) 
     235         aeiw_temp(:,:,:)=aeiw(:,:,:) 
    236236#elif defined key_traldf_c2d 
    237             aeiu_temp(:,:)=aeiu(:,:) 
    238             aeiv_temp(:,:)=aeiv(:,:) 
    239             aeiw_temp(:,:)=aeiw(:,:) 
     237         aeiu_temp(:,:)=aeiu(:,:) 
     238         aeiv_temp(:,:)=aeiv(:,:) 
     239         aeiw_temp(:,:)=aeiw(:,:) 
    240240#elif defined key_traldf_c1d 
    241             aeiu_temp(:)=aeiu(:) 
    242             aeiv_temp(:)=aeiv(:) 
    243             aeiw_temp(:)=aeiw(:) 
     241         aeiu_temp(:)=aeiu(:) 
     242         aeiv_temp(:)=aeiv(:) 
     243         aeiw_temp(:)=aeiw(:) 
    244244#else 
    245             aeiu_temp=aeiu 
    246             aeiv_temp=aeiv 
    247             aeiw_temp=aeiw 
     245         aeiu_temp=aeiu 
     246         aeiv_temp=aeiv 
     247         aeiw_temp=aeiw 
    248248#endif 
    249249#endif 
    250250#if defined key_ldfslp 
    251             wslpi_temp(:,:,:)=wslpi(:,:,:) 
    252             wslpj_temp(:,:,:)=wslpj(:,:,:) 
    253             uslp_temp(:,:,:)=uslp(:,:,:) 
    254             vslp_temp(:,:,:)=vslp(:,:,:) 
    255 #endif 
    256             hu_temp(:,:)=hu(:,:) 
    257             hv_temp(:,:)=hv(:,:) 
    258             hur_temp(:,:)=hur(:,:) 
    259             hvr_temp(:,:)=hvr(:,:) 
     251         wslpi_temp(:,:,:)=wslpi(:,:,:) 
     252         wslpj_temp(:,:,:)=wslpj(:,:,:) 
     253         uslp_temp(:,:,:)=uslp(:,:,:) 
     254         vslp_temp(:,:,:)=vslp(:,:,:) 
     255#endif 
     256         hu_temp(:,:)=hu(:,:) 
     257         hv_temp(:,:)=hv(:,:) 
     258         hur_temp(:,:)=hur(:,:) 
     259         hvr_temp(:,:)=hvr(:,:) 
    260260 
    261261!   2. create averages and reassign variables 
    262             DO jj = 1, jpj 
    263                DO ji = 1, jpi 
    264                   sshn_tm(ji,jj)= sshn_tm(ji,jj)+sshn(ji,jj) 
    265                   sshn(ji,jj)=sshn_tm(ji,jj)/(nn_dttrc+1) 
     262         DO jj = 1, jpj 
     263            DO ji = 1, jpi 
     264               sshn_tm(ji,jj)= sshn_tm(ji,jj)+sshn(ji,jj) 
     265               sshn(ji,jj)=sshn_tm(ji,jj)/(nn_dttrc+1) 
    266266         
    267                   sshu_n_tm(ji,jj)= sshu_n_tm(ji,jj)+sshu_n(ji,jj) 
    268                   sshu_n(ji,jj)=sshu_n_tm(ji,jj)/(nn_dttrc+1) 
    269          
    270                   sshv_n_tm(ji,jj)= sshv_n_tm(ji,jj)+sshv_n(ji,jj) 
    271                   sshv_n(ji,jj)=sshv_n_tm(ji,jj)/(nn_dttrc+1) 
    272  
    273         sshb(ji,jj)=sshb_hold(ji,jj) 
    274  
    275         sshu_b(ji,jj)=sshu_b_hold(ji,jj) 
    276  
    277         sshv_b(ji,jj)=sshv_b_hold(ji,jj) 
    278  
    279                   rnf_tm(ji,jj)= rnf_tm(ji,jj)+rnf(ji,jj) 
    280                   rnf(ji,jj)=rnf_tm(ji,jj)/(nn_dttrc+1) 
     267               sshu_n_tm(ji,jj)= sshu_n_tm(ji,jj)+sshu_n(ji,jj) 
     268               sshu_n(ji,jj)=sshu_n_tm(ji,jj)/(nn_dttrc+1) 
     269 
     270               sshv_n_tm(ji,jj)= sshv_n_tm(ji,jj)+sshv_n(ji,jj) 
     271               sshv_n(ji,jj)=sshv_n_tm(ji,jj)/(nn_dttrc+1) 
     272 
     273               sshb(ji,jj)=sshb_hold(ji,jj) 
     274 
     275               sshu_b(ji,jj)=sshu_b_hold(ji,jj) 
     276 
     277          sshv_b(ji,jj)=sshv_b_hold(ji,jj) 
     278 
     279               rnf_tm(ji,jj)= rnf_tm(ji,jj)+rnf(ji,jj) 
     280               rnf(ji,jj)=rnf_tm(ji,jj)/(nn_dttrc+1) 
    281281  
    282                   h_rnf_tm(ji,jj)= h_rnf_tm(ji,jj)+h_rnf(ji,jj) 
    283                   h_rnf(ji,jj)=h_rnf_tm(ji,jj)/(nn_dttrc+1) 
    284  
    285                   hmld_tm(ji,jj)= hmld_tm(ji,jj)+hmld(ji,jj) 
    286                   hmld(ji,jj)=hmld_tm(ji,jj)/(nn_dttrc+1) 
     282               h_rnf_tm(ji,jj)= h_rnf_tm(ji,jj)+h_rnf(ji,jj) 
     283               h_rnf(ji,jj)=h_rnf_tm(ji,jj)/(nn_dttrc+1) 
     284 
     285               hmld_tm(ji,jj)= hmld_tm(ji,jj)+hmld(ji,jj) 
     286               hmld(ji,jj)=hmld_tm(ji,jj)/(nn_dttrc+1) 
    287287  
    288                   fr_i_tm(ji,jj)= fr_i_tm(ji,jj)+fr_i(ji,jj) 
    289                   fr_i(ji,jj)=fr_i_tm(ji,jj)/(nn_dttrc+1) 
     288               fr_i_tm(ji,jj)= fr_i_tm(ji,jj)+fr_i(ji,jj) 
     289               fr_i(ji,jj)=fr_i_tm(ji,jj)/(nn_dttrc+1) 
    290290 
    291291!  variables that are initialized after averages initialized 
    292                   wndm_tm(ji,jj)= wndm_tm(ji,jj)+wndm(ji,jj) 
    293                   IF(kt == nittrc000)THEN 
    294                       wndm(ji,jj)=wndm_tm(ji,jj)/(nn_dttrc) 
    295                   ELSE 
    296                       wndm(ji,jj)=wndm_tm(ji,jj)/(nn_dttrc+1) 
    297                   ENDIF  
    298  
    299                   qsr_tm(ji,jj)= qsr_tm(ji,jj)+qsr(ji,jj) 
    300                   IF(kt == nittrc000)THEN 
    301                       qsr(ji,jj)=qsr_tm(ji,jj)/(nn_dttrc) 
    302                   ELSE 
    303                       qsr(ji,jj)=qsr_tm(ji,jj)/(nn_dttrc+1) 
    304                   ENDIF  
    305  
    306         emp_b(ji,jj)=emp_b_hold(ji,jj) 
    307  
    308                   emp_tm(ji,jj)= emp_tm(ji,jj)+emp(ji,jj) 
    309                   IF(kt == nittrc000)THEN 
    310                      emp(ji,jj)=emp_tm(ji,jj)/(nn_dttrc) 
    311                   ELSE 
    312                      emp(ji,jj)=emp_tm(ji,jj)/(nn_dttrc+1) 
    313                   ENDIF  
    314  
    315                   emps_tm(ji,jj)= emps_tm(ji,jj)+emps(ji,jj) 
    316                   IF(kt == nittrc000)THEN 
    317                      emps(ji,jj)=emps_tm(ji,jj)/(nn_dttrc) 
    318                   ELSE 
    319                      emps(ji,jj)=emps_tm(ji,jj)/(nn_dttrc+1) 
    320                   ENDIF 
    321  
    322                END DO 
    323             END DO 
    324          DO jk = 1, jpk 
    325              IF( lk_vvl ) THEN                           !  Update Now Vertical coord.  !   (only in vvl case) 
    326              !                                           !------------------------------! 
     292               wndm_tm(ji,jj)= wndm_tm(ji,jj)+wndm(ji,jj) 
     293               IF(kt == nittrc000)THEN 
     294                   wndm(ji,jj)=wndm_tm(ji,jj)/(nn_dttrc) 
     295               ELSE 
     296                   wndm(ji,jj)=wndm_tm(ji,jj)/(nn_dttrc+1) 
     297               ENDIF  
     298 
     299               qsr_tm(ji,jj)= qsr_tm(ji,jj)+qsr(ji,jj) 
     300               IF(kt == nittrc000)THEN 
     301                   qsr(ji,jj)=qsr_tm(ji,jj)/(nn_dttrc) 
     302               ELSE 
     303                   qsr(ji,jj)=qsr_tm(ji,jj)/(nn_dttrc+1) 
     304               ENDIF  
     305 
     306          emp_b(ji,jj)=emp_b_hold(ji,jj) 
     307 
     308               emp_tm(ji,jj)= emp_tm(ji,jj)+emp(ji,jj) 
     309               IF(kt == nittrc000)THEN 
     310                  emp(ji,jj)=emp_tm(ji,jj)/(nn_dttrc) 
     311               ELSE 
     312                  emp(ji,jj)=emp_tm(ji,jj)/(nn_dttrc+1) 
     313               ENDIF  
     314 
     315               emps_tm(ji,jj)= emps_tm(ji,jj)+emps(ji,jj) 
     316               IF(kt == nittrc000)THEN 
     317                  emps(ji,jj)=emps_tm(ji,jj)/(nn_dttrc) 
     318               ELSE 
     319                  emps(ji,jj)=emps_tm(ji,jj)/(nn_dttrc+1) 
     320               ENDIF 
     321 
     322            END DO 
     323         END DO 
     324         DO jk = 1, jpk 
     325            e3t_temp(:,:,jk)= fse3t(:,:,jk)  
     326            e3u_temp(:,:,jk)= fse3u(:,:,jk)  
     327            e3v_temp(:,:,jk)= fse3v(:,:,jk) 
     328            e3w_temp(:,:,jk)= fse3w(:,:,jk)   
     329         END DO 
     330         IF( lk_vvl ) THEN                           !  Update Now Vertical coord.  !   (only in vvl case) 
     331         !                                           !------------------------------! 
     332            DO jk = 1, jpk 
    327333                 fse3t (:,:,jk) = fse3t_n (:,:,jk)          ! vertical scale factors stored in fse3. arrays 
    328334                 fse3u (:,:,jk) = fse3u_n (:,:,jk) 
    329335                 fse3v (:,:,jk) = fse3v_n (:,:,jk) 
    330                  fse3f (:,:,jk) = fse3f_n (:,:,jk) 
    331336                 fse3w (:,:,jk) = fse3w_n (:,:,jk) 
    332                  fse3uw(:,:,jk) = fse3uw_n(:,:,jk) 
    333                  fse3vw(:,:,jk) = fse3vw_n(:,:,jk) 
    334              ENDIF 
    335          END DO 
     337            END DO 
     338         ENDIF 
    336339         DO jk = 1, jpk 
    337340            DO jj = 1, jpj 
    … …  
    382385            DO jj = 1, jpj 
    383386               DO ji = 1, jpi 
    384                   ahtt_tm(ji,jj,jk)= ahtt_tm(ji,jj,jk)+ahtt(ji,jj,jk)*fse3t(ji,jj,jk) 
     387                  ahtt_tm(ji,jj,jk)= ahtt_tm(ji,jj,jk)+ahtt(ji,jj,jk)*e3t_temp(ji,jj,jk) 
    385388                  ahtt(ji,jj,jk)=(ahtt_tm(ji,jj,jk)/(nn_dttrc+1))/fse3t(ji,jj,jk)    
    386389 
    387                   ahtu_tm(ji,jj,jk)= ahtu_tm(ji,jj,jk)+ahtu(ji,jj,jk)*fse3u(ji,jj,jk) 
     390                  ahtu_tm(ji,jj,jk)= ahtu_tm(ji,jj,jk)+ahtu(ji,jj,jk)*e3u_temp(ji,jj,jk) 
    388391                  ahtu(ji,jj,jk)=(ahtu_tm(ji,jj,jk)/(nn_dttrc+1))/fse3u(ji,jj,jk)    
    389392 
    390                   ahtv_tm(ji,jj,jk)= ahtv_tm(ji,jj,jk)+ahtv(ji,jj,jk)*fse3v(ji,jj,jk) 
     393                  ahtv_tm(ji,jj,jk)= ahtv_tm(ji,jj,jk)+ahtv(ji,jj,jk)*e3v_temp(ji,jj,jk) 
    391394                  ahtv(ji,jj,jk)=(ahtv_tm(ji,jj,jk)/(nn_dttrc+1))/fse3v(ji,jj,jk)    
    392395 
    393                   ahtw_tm(ji,jj,jk)= ahtw_tm(ji,jj,jk)+ahtw(ji,jj,jk)*fse3w(ji,jj,jk) 
     396                  ahtw_tm(ji,jj,jk)= ahtw_tm(ji,jj,jk)+ahtw(ji,jj,jk)*e3w_temp(ji,jj,jk) 
    394397                  ahtw(ji,jj,jk)=(ahtw_tm(ji,jj,jk)/(nn_dttrc+1))/fse3w(ji,jj,jk)    
    395398 
    396               END DO 
     399               END DO 
    397400            END DO 
    398401         END DO 
    399402#elif defined key_traldf_c2d 
     403         DO jj = 1, jpj 
     404            DO ji = 1, jpi 
     405               ahtt_tm(ji,jj)= ahtt_tm(ji,jj)+ahtt(ji,jj) 
     406               ahtt(ji,jj)=ahtt_tm(ji,jj)/(nn_dttrc+1) 
     407 
     408               ahtu_tm(ji,jj)= ahtu_tm(ji,jj)+ahtu(ji,jj) 
     409               ahtu(ji,jj)=ahtu_tm(ji,jj)/(nn_dttrc+1) 
     410 
     411               ahtv_tm(ji,jj)= ahtv_tm(ji,jj)+ahtv(ji,jj) 
     412               ahtv(ji,jj)=ahtv_tm(ji,jj)/(nn_dttrc+1) 
     413 
     414               ahtw_tm(ji,jj)= ahtw_tm(ji,jj)+ahtw(ji,jj) 
     415               ahtw(ji,jj)=ahtw_tm(ji,jj)/(nn_dttrc+1) 
     416            END DO 
     417         END DO 
     418#elif defined key_traldf_c1d 
     419         DO jk = 1, jpk 
     420            ahtt_tm(jk)= ahtt_tm(jk)+ahtt(jk) 
     421            ahtt(jk)=ahtt_tm(jk)/(nn_dttrc+1) 
     422 
     423            ahtu_tm(jk)= ahtu_tm(jk)+ahtu(jk) 
     424            ahtu(jk)=ahtu_tm(jk)/(nn_dttrc+1) 
     425 
     426            ahtv_tm(jk)= ahtv_tm(jk)+ahtv(jk) 
     427            ahtv(jk)=ahtv_tm(jk)/(nn_dttrc+1) 
     428 
     429            ahtw_tm(jk)= ahtw_tm(jk)+ahtw(jk) 
     430            ahtw(jk)=ahtw_tm(jk)/(nn_dttrc+1) 
     431         END DO 
     432#else 
     433         ahtt_tm= ahtt_tm+ahtt 
     434         ahtt=ahtt_tm/(nn_dttrc+1) 
     435 
     436         ahtu_tm= ahtu_tm+ahtu 
     437         ahtu=ahtu_tm/(nn_dttrc+1) 
     438 
     439         ahtv_tm= ahtv_tm+ahtv 
     440         ahtv=ahtv_tm/(nn_dttrc+1) 
     441 
     442         ahtw_tm= ahtw_tm+ahtw 
     443         ahtw=ahtw_tm/(nn_dttrc+1) 
     444 
     445#endif 
     446#if defined key_traldf_eiv 
     447# if defined key_traldf_c3d 
     448         DO jk = 1, jpk 
    400449            DO jj = 1, jpj 
    401450               DO ji = 1, jpi 
    402                   ahtt_tm(ji,jj)= ahtt_tm(ji,jj)+ahtt(ji,jj) 
    403                   ahtt(ji,jj)=ahtt_tm(ji,jj)/(nn_dttrc+1) 
    404  
    405                   ahtu_tm(ji,jj)= ahtu_tm(ji,jj)+ahtu(ji,jj) 
    406                   ahtu(ji,jj)=ahtu_tm(ji,jj)/(nn_dttrc+1) 
    407  
    408                   ahtv_tm(ji,jj)= ahtv_tm(ji,jj)+ahtv(ji,jj) 
    409                   ahtv(ji,jj)=ahtv_tm(ji,jj)/(nn_dttrc+1) 
    410  
    411                   ahtw_tm(ji,jj)= ahtw_tm(ji,jj)+ahtw(ji,jj) 
    412                   ahtw(ji,jj)=ahtw_tm(ji,jj)/(nn_dttrc+1) 
    413  
    414              END DO 
    415             END DO 
    416 #elif defined key_traldf_c1d 
    417           DO jk = 1, jpk 
    418              ahtt_tm(jk)= ahtt_tm(jk)+ahtt(jk) 
    419              ahtt(jk)=ahtt_tm(jk)/(nn_dttrc+1) 
    420  
    421              ahtu_tm(jk)= ahtu_tm(jk)+ahtu(jk) 
    422              ahtu(jk)=ahtu_tm(jk)/(nn_dttrc+1) 
    423  
    424              ahtv_tm(jk)= ahtv_tm(jk)+ahtv(jk) 
    425              ahtv(jk)=ahtv_tm(jk)/(nn_dttrc+1) 
    426  
    427              ahtw_tm(jk)= ahtw_tm(jk)+ahtw(jk) 
    428              ahtw(jk)=ahtw_tm(jk)/(nn_dttrc+1) 
    429           END DO 
    430 #else 
    431              ahtt_tm= ahtt_tm+ahtt 
    432              ahtt=ahtt_tm/(nn_dttrc+1) 
    433  
    434              ahtu_tm= ahtu_tm+ahtu 
    435              ahtu=ahtu_tm/(nn_dttrc+1) 
    436  
    437              ahtv_tm= ahtv_tm+ahtv 
    438              ahtv=ahtv_tm/(nn_dttrc+1) 
    439  
    440              ahtw_tm= ahtw_tm+ahtw 
    441              ahtw=ahtw_tm/(nn_dttrc+1) 
    442  
    443 #endif 
    444 #if defined key_traldf_eiv 
    445 # if defined key_traldf_c3d 
    446          DO jk = 1, jpk 
    447             DO jj = 1, jpj 
    448                DO ji = 1, jpi 
    449                   aeiu_tm(ji,jj,jk)= aeiu_tm(ji,jj,jk)+aeiu(ji,jj,jk)*fse3u(ji,jj,jk) 
     451                  aeiu_tm(ji,jj,jk)= aeiu_tm(ji,jj,jk)+aeiu(ji,jj,jk)*e3u_temp(ji,jj,jk) 
    450452                  aeiu(ji,jj,jk)=(aeiu_tm(ji,jj,jk)/(nn_dttrc+1))/fse3u(ji,jj,jk)    
    451453 
    452                   aeiv_tm(ji,jj,jk)= aeiv_tm(ji,jj,jk)+aeiv(ji,jj,jk)*fse3v(ji,jj,jk) 
     454                  aeiv_tm(ji,jj,jk)= aeiv_tm(ji,jj,jk)+aeiv(ji,jj,jk)*e3v_temp(ji,jj,jk) 
    453455                  aeiv(ji,jj,jk)=(aeiv_tm(ji,jj,jk)/(nn_dttrc+1))/fse3v(ji,jj,jk)    
    454456 
    455                   aeiw_tm(ji,jj,jk)= aeiw_tm(ji,jj,jk)+aeiw(ji,jj,jk)*fse3w(ji,jj,jk) 
     457                  aeiw_tm(ji,jj,jk)= aeiw_tm(ji,jj,jk)+aeiw(ji,jj,jk)*e3w_temp(ji,jj,jk) 
    456458                  aeiw(ji,jj,jk)=(aeiw_tm(ji,jj,jk)/(nn_dttrc+1))/fse3w(ji,jj,jk)    
    457459 
    458              END DO 
     460               END DO 
    459461            END DO 
    460462         END DO 
    461463# elif defined key_traldf_c2d 
    462             DO jj = 1, jpj 
    463                DO ji = 1, jpi 
    464                   aeiu_tm(ji,jj)= aeiu_tm(ji,jj)+aeiu(ji,jj) 
    465                   aeiu(ji,jj)=aeiu_tm(ji,jj)/(nn_dttrc+1) 
    466  
    467                   aeiv_tm(ji,jj)= aeiv_tm(ji,jj)+aeiv(ji,jj) 
    468                   aeiv(ji,jj)=aeiv_tm(ji,jj)/(nn_dttrc+1) 
    469  
    470                   aeiw_tm(ji,jj)= aeiw_tm(ji,jj)+aeiw(ji,jj) 
    471                   aeiw(ji,jj)=aeiw_tm(ji,jj)/(nn_dttrc+1) 
    472              END DO 
    473             END DO 
     464         DO jj = 1, jpj 
     465            DO ji = 1, jpi 
     466               aeiu_tm(ji,jj)= aeiu_tm(ji,jj)+aeiu(ji,jj) 
     467               aeiu(ji,jj)=aeiu_tm(ji,jj)/(nn_dttrc+1) 
     468 
     469               aeiv_tm(ji,jj)= aeiv_tm(ji,jj)+aeiv(ji,jj) 
     470               aeiv(ji,jj)=aeiv_tm(ji,jj)/(nn_dttrc+1) 
     471 
     472               aeiw_tm(ji,jj)= aeiw_tm(ji,jj)+aeiw(ji,jj) 
     473               aeiw(ji,jj)=aeiw_tm(ji,jj)/(nn_dttrc+1) 
     474            END DO 
     475         END DO 
    474476# elif defined key_traldf_c1d 
    475           DO jk = 1, jpk 
    476              aeiu_tm(jk)= aeiu_tm(jk)+aeiu(jk) 
    477              aeiu(jk)=aeiu_tm(jk)/(nn_dttrc+1) 
    478  
    479              aeiv_tm(jk)= aeiv_tm(jk)+aeiv(jk) 
    480              aeiv(jk)=aeiv_tm(jk)/(nn_dttrc+1) 
    481  
    482              aeiw_tm(jk)= aeiw_tm(jk)+aeiw(jk) 
    483              aeiw(jk)=aeiw_tm(jk)/(nn_dttrc+1) 
    484           END DO 
     477         DO jk = 1, jpk 
     478            aeiu_tm(jk)= aeiu_tm(jk)+aeiu(jk) 
     479            aeiu(jk)=aeiu_tm(jk)/(nn_dttrc+1) 
     480 
     481            aeiv_tm(jk)= aeiv_tm(jk)+aeiv(jk) 
     482            aeiv(jk)=aeiv_tm(jk)/(nn_dttrc+1) 
     483 
     484            aeiw_tm(jk)= aeiw_tm(jk)+aeiw(jk) 
     485            aeiw(jk)=aeiw_tm(jk)/(nn_dttrc+1) 
     486         END DO 
    485487# else 
    486              aeiu_tm= aeiu_tm+aeiu 
    487              aeiu=aeiu_tm/(nn_dttrc+1) 
    488  
    489              aeiv_tm= aeiv_tm+aeiv 
    490              aeiv=aeiv_tm/(nn_dttrc+1) 
    491  
    492              aeiw_tm= aeiw_tm+aeiw 
    493              aeiw=aeiw_tm/(nn_dttrc+1) 
     488         aeiu_tm= aeiu_tm+aeiu 
     489         aeiu=aeiu_tm/(nn_dttrc+1) 
     490 
     491         aeiv_tm= aeiv_tm+aeiv 
     492         aeiv=aeiv_tm/(nn_dttrc+1) 
     493 
     494         aeiw_tm= aeiw_tm+aeiw 
     495         aeiw=aeiw_tm/(nn_dttrc+1) 
    494496# endif 
    495497# endif 
    … …  
    535537#endif 
    536538# if defined key_zdfddm 
    537           CALL lbc_lnk( avs, 'W', 1. )  
     539         CALL lbc_lnk( avs, 'W', 1. )  
    538540# endif 
    539541#if defined key_traldf_eiv 
    … …  
    748750       emp_b(:,:)=emp_b_temp(:,:) 
    749751#if defined key_traldf_c3d 
    750             ahtu_temp(:,:,:)=ahtu(:,:,:) 
    751             ahtv_temp(:,:,:)=ahtv(:,:,:) 
    752             ahtw_temp(:,:,:)=ahtw(:,:,:) 
    753             ahtt_temp(:,:,:)=ahtt(:,:,:) 
     752       ahtu_temp(:,:,:)=ahtu(:,:,:) 
     753       ahtv_temp(:,:,:)=ahtv(:,:,:) 
     754       ahtw_temp(:,:,:)=ahtw(:,:,:) 
     755       ahtt_temp(:,:,:)=ahtt(:,:,:) 
    754756#elif defined key_traldf_c2d 
    755             ahtu_temp(:,:)=ahtu(:,:) 
    756             ahtv_temp(:,:)=ahtv(:,:) 
    757             ahtw_temp(:,:)=ahtw(:,:) 
    758             ahtt_temp(:,:)=ahtt(:,:) 
     757       ahtu_temp(:,:)=ahtu(:,:) 
     758       ahtv_temp(:,:)=ahtv(:,:) 
     759       ahtw_temp(:,:)=ahtw(:,:) 
     760       ahtt_temp(:,:)=ahtt(:,:) 
    759761#elif defined key_traldf_c1d 
    760             ahtu_temp(:)=ahtu(:) 
    761             ahtv_temp(:)=ahtv(:) 
    762             ahtw_temp(:)=ahtw(:) 
    763             ahtt_temp(:)=ahtt(:) 
     762       ahtu_temp(:)=ahtu(:) 
     763       ahtv_temp(:)=ahtv(:) 
     764       ahtw_temp(:)=ahtw(:) 
     765       ahtt_temp(:)=ahtt(:) 
    764766#else 
    765             ahtu_temp=ahtu 
    766             ahtv_temp=ahtv 
    767             ahtw_temp=ahtw 
    768             ahtt_temp=ahtt 
     767       ahtu_temp=ahtu 
     768       ahtv_temp=ahtv 
     769       ahtw_temp=ahtw 
     770       ahtt_temp=ahtt 
    769771#endif 
    770772#if defined key_traldf_eiv 
    771773#if defined key_traldf_c3d 
    772             aeiu_temp(:,:,:)=aeiu(:,:,:) 
    773             aeiv_temp(:,:,:)=aeiv(:,:,:) 
    774             aeiw_temp(:,:,:)=aeiw(:,:,:) 
     774       aeiu_temp(:,:,:)=aeiu(:,:,:) 
     775       aeiv_temp(:,:,:)=aeiv(:,:,:) 
     776       aeiw_temp(:,:,:)=aeiw(:,:,:) 
    775777#elif defined key_traldf_c2d 
    776             aeiu_temp(:,:)=aeiu(:,:) 
    777             aeiv_temp(:,:)=aeiv(:,:) 
    778             aeiw_temp(:,:)=aeiw(:,:) 
     778       aeiu_temp(:,:)=aeiu(:,:) 
     779       aeiv_temp(:,:)=aeiv(:,:) 
     780       aeiw_temp(:,:)=aeiw(:,:) 
    779781#elif defined key_traldf_c1d 
    780             aeiu_temp(:)=aeiu(:) 
    781             aeiv_temp(:)=aeiv(:) 
    782             aeiw_temp(:)=aeiw(:) 
     782       aeiu_temp(:)=aeiu(:) 
     783       aeiv_temp(:)=aeiv(:) 
     784       aeiw_temp(:)=aeiw(:) 
    783785#else 
    784             aeiu_temp=aeiu 
    785             aeiv_temp=aeiv 
    786             aeiw_temp=aeiw 
     786       aeiu_temp=aeiu 
     787       aeiv_temp=aeiv 
     788       aeiw_temp=aeiw 
    787789#endif 
    788790#endif         
    … …  
    790792 
    791793       !                                           !------------------------------! 
     794       DO jk = 1, jpk 
     795          fse3t(:,:,jk)= e3t_temp(:,:,jk) !0.  
     796          fse3u(:,:,jk)= e3u_temp(:,:,jk) !0.  
     797          fse3v(:,:,jk)= e3v_temp(:,:,jk)!0.  
     798          fse3w(:,:,jk)= e3w_temp(:,:,jk) !0.  
     799       END DO 
    792800       IF( lk_vvl ) THEN                           !  Update Now Vertical coord.  !   (only in vvl case) 
    793801       !                                           !------------------------------! 
    … …  
    808816 
    809817!   start new averages 
    810          DO jk = 1, jpk 
    811             DO jj = 1, jpj 
    812                DO ji = 1, jpi 
    813                   un_tm(ji,jj,jk)=un(ji,jj,jk)*fse3u(ji,jj,jk)  !0.0 
    814                   vn_tm(ji,jj,jk)=vn(ji,jj,jk)*fse3v(ji,jj,jk)  !0.0 
    815                   wn_tm(ji,jj,jk)=wn(ji,jj,jk)*fse3w(ji,jj,jk)  !0.0 
    816                   avt_tm(ji,jj,jk)=avt(ji,jj,jk)*fse3t(ji,jj,jk)  !0.0 
    817                   tsn_tm(ji,jj,jk,jp_tem)=tsn(ji,jj,jk,jp_tem)*fse3t(ji,jj,jk)  !0.0 
    818                   tsn_tm(ji,jj,jk,jp_sal)=tsn(ji,jj,jk,jp_sal)*fse3t(ji,jj,jk)  !0.0 
     818      DO jk = 1, jpk 
     819         DO jj = 1, jpj 
     820            DO ji = 1, jpi 
     821               un_tm(ji,jj,jk)=un(ji,jj,jk)*fse3u(ji,jj,jk)  !0.0 
     822               vn_tm(ji,jj,jk)=vn(ji,jj,jk)*fse3v(ji,jj,jk)  !0.0 
     823               wn_tm(ji,jj,jk)=wn(ji,jj,jk)*fse3w(ji,jj,jk)  !0.0 
     824               avt_tm(ji,jj,jk)=avt(ji,jj,jk)*fse3t(ji,jj,jk)  !0.0 
     825               tsn_tm(ji,jj,jk,jp_tem)=tsn(ji,jj,jk,jp_tem)*fse3t(ji,jj,jk)  !0.0 
     826               tsn_tm(ji,jj,jk,jp_sal)=tsn(ji,jj,jk,jp_sal)*fse3t(ji,jj,jk)  !0.0 
    819827# if defined key_zdfddm 
    820                   avs_tm==avs(ji,jj,jk)*fse3w(ji,jj,jk)  !0.0 
     828               avs_tm==avs(ji,jj,jk)*fse3w(ji,jj,jk)  !0.0 
    821829# endif 
    822830#if defined key_ldfslp 
    823                   wslpi_tm(ji,jj,jk)= wslpi(ji,jj,jk)*fse3w(ji,jj,jk) !0.  
    824                   wslpj_tm(ji,jj,jk)=wslpj(ji,jj,jk)*fse3w(ji,jj,jk) !0.  
    825                   uslp_tm(ji,jj,jk)=uslp(ji,jj,jk)*fse3u(ji,jj,jk) !0.  
    826                   vslp_tm(ji,jj,jk)=vslp(ji,jj,jk)*fse3v(ji,jj,jk) !0.  
    827 #endif 
    828           END DO 
    829             END DO 
    830          END DO 
     831               wslpi_tm(ji,jj,jk)= wslpi(ji,jj,jk)*fse3w(ji,jj,jk) !0.  
     832               wslpj_tm(ji,jj,jk)=wslpj(ji,jj,jk)*fse3w(ji,jj,jk) !0.  
     833               uslp_tm(ji,jj,jk)=uslp(ji,jj,jk)*fse3u(ji,jj,jk) !0.  
     834               vslp_tm(ji,jj,jk)=vslp(ji,jj,jk)*fse3v(ji,jj,jk) !0.  
     835#endif 
     836            END DO 
     837         END DO 
     838      END DO 
    831839      sshn_tm(:,:) = sshn(:,:) !0. 
    832840      sshu_n_tm(:,:) = sshu_n(:,:) !0. 
    … …  
    845853      wndm_tm(:,:)=wndm(:,:) 
    846854#if defined key_traldf_c3d 
    847          DO jk = 1, jpk 
    848             DO jj = 1, jpj 
    849                DO ji = 1, jpi 
    850                   ahtt_tm(ji,jj,jk)= ahtt(ji,jj,jk)*fse3t(ji,jj,jk) 
    851                   ahtu_tm(ji,jj,jk)= ahtu(ji,jj,jk)*fse3u(ji,jj,jk) 
    852                   ahtv_tm(ji,jj,jk)= ahtv(ji,jj,jk)*fse3v(ji,jj,jk) 
    853                   ahtw_tm(ji,jj,jk)= ahtw(ji,jj,jk)*fse3w(ji,jj,jk) 
    854              END DO 
    855             END DO 
    856          END DO 
     855      DO jk = 1, jpk 
     856         DO jj = 1, jpj 
     857            DO ji = 1, jpi 
     858               ahtt_tm(ji,jj,jk)= ahtt(ji,jj,jk)*fse3t(ji,jj,jk) 
     859               ahtu_tm(ji,jj,jk)= ahtu(ji,jj,jk)*fse3u(ji,jj,jk) 
     860               ahtv_tm(ji,jj,jk)= ahtv(ji,jj,jk)*fse3v(ji,jj,jk) 
     861               ahtw_tm(ji,jj,jk)= ahtw(ji,jj,jk)*fse3w(ji,jj,jk) 
     862            END DO 
     863         END DO 
     864      END DO 
    857865#elif defined key_traldf_c2d 
    858             DO jj = 1, jpj 
    859                DO ji = 1, jpi 
    860                   ahtt_tm(ji,jj)= ahtt(ji,jj) 
    861                   ahtu_tm(ji,jj)= ahtu(ji,jj) 
    862                   ahtv_tm(ji,jj)= ahtv(ji,jj) 
    863                   ahtw_tm(ji,jj)= ahtw(ji,jj) 
    864              END DO 
    865             END DO 
     866      DO jj = 1, jpj 
     867         DO ji = 1, jpi 
     868            ahtt_tm(ji,jj)= ahtt(ji,jj) 
     869            ahtu_tm(ji,jj)= ahtu(ji,jj) 
     870            ahtv_tm(ji,jj)= ahtv(ji,jj) 
     871            ahtw_tm(ji,jj)= ahtw(ji,jj) 
     872         END DO 
     873      END DO 
    866874#elif defined key_traldf_c1d 
    867           DO jk = 1, jpk 
    868              ahtt_tm(jk)= ahtt(jk) 
    869              ahtu_tm(jk)= ahtu(jk) 
    870              ahtv_tm(jk)= ahtv(jk) 
    871              ahtw_tm(jk)= ahtw(jk) 
    872           END DO 
     875      DO jk = 1, jpk 
     876         ahtt_tm(jk)= ahtt(jk) 
     877         ahtu_tm(jk)= ahtu(jk) 
     878         ahtv_tm(jk)= ahtv(jk) 
     879         ahtw_tm(jk)= ahtw(jk) 
     880      END DO 
    873881#else 
    874              ahtt_tm= ahtt 
    875              ahtu_tm= ahtu 
    876              ahtv_tm= ahtv 
    877              ahtw_tm= ahtw 
     882      ahtt_tm= ahtt 
     883      ahtu_tm= ahtu 
     884      ahtv_tm= ahtv 
     885      ahtw_tm= ahtw 
    878886#endif 
    879887#if defined key_traldf_eiv 
    880888# if defined key_traldf_c3d 
    881          DO jk = 1, jpk 
    882             DO jj = 1, jpj 
    883                DO ji = 1, jpi 
    884                   aeiu_tm(ji,jj,jk)= aeiu(ji,jj,jk)*fse3u(ji,jj,jk) 
    885                   aeiv_tm(ji,jj,jk)= aeiv(ji,jj,jk)*fse3v(ji,jj,jk) 
    886                   aeiw_tm(ji,jj,jk)= aeiw(ji,jj,jk)*fse3w(ji,jj,jk) 
    887              END DO 
    888             END DO 
    889          END DO 
     889      DO jk = 1, jpk 
     890         DO jj = 1, jpj 
     891            DO ji = 1, jpi 
     892               aeiu_tm(ji,jj,jk)= aeiu(ji,jj,jk)*fse3u(ji,jj,jk) 
     893               aeiv_tm(ji,jj,jk)= aeiv(ji,jj,jk)*fse3v(ji,jj,jk) 
     894               aeiw_tm(ji,jj,jk)= aeiw(ji,jj,jk)*fse3w(ji,jj,jk) 
     895            END DO 
     896         END DO 
     897      END DO 
    890898# elif defined key_traldf_c2d 
    891             DO jj = 1, jpj 
    892                DO ji = 1, jpi 
    893                   aeiu_tm(ji,jj)= aeiu(ji,jj) 
    894                   aeiv_tm(ji,jj)= aeiv(ji,jj) 
    895                   aeiw_tm(ji,jj)= aeiw(ji,jj) 
    896              END DO 
    897             END DO 
     899      DO jj = 1, jpj 
     900         DO ji = 1, jpi 
     901            aeiu_tm(ji,jj)= aeiu(ji,jj) 
     902            aeiv_tm(ji,jj)= aeiv(ji,jj) 
     903            aeiw_tm(ji,jj)= aeiw(ji,jj) 
     904         END DO 
     905      END DO 
    898906# elif defined key_traldf_c1d 
    899           DO jk = 1, jpk 
    900              aeiu_tm(jk)= aeiu(jk) 
    901              aeiv_tm(jk)= aeiv(jk) 
    902              aeiw_tm(jk)= aeiw(jk) 
    903           END DO 
     907      DO jk = 1, jpk 
     908         aeiu_tm(jk)= aeiu(jk) 
     909         aeiv_tm(jk)= aeiv(jk) 
     910         aeiw_tm(jk)= aeiw(jk) 
     911      END DO 
    904912# else 
    905              aeiu_tm= aeiu 
    906              aeiv_tm= aeiv 
    907              aeiw_tm= aeiw 
     913      aeiu_tm= aeiu 
     914      aeiv_tm= aeiv 
     915      aeiw_tm= aeiw 
    908916# endif 
    909917# endif 
    … …  
    10661074         &      sshv_a_temp(jpi,jpj)        ,  sshf_a_temp(jpi,jpj),     & 
    10671075         &      sshf_n_temp(jpi,jpj)        ,  sshf_b_temp(jpi,jpj),     & 
     1076         &      e3t_temp(jpi,jpj,jpk)       ,  e3u_temp(jpi,jpj,jpk),    & 
     1077         &      e3v_temp(jpi,jpj,jpk)       ,  e3w_temp(jpi,jpj,jpk),    & 
    10681078#if defined key_ldfslp 
    10691079         &      wslpi_temp(jpi,jpj,jpk)     ,  wslpj_temp(jpi,jpj,jpk),  & 
Note: See TracChangeset for help on using the changeset viewer.