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 11463 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST – NEMO

Ignore:
Timestamp:
2019-08-20T14:14:56+02:00 (5 years ago)
Author:
acc
Message:

Branch: dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap. Minor bugfix in step.F90 to enable AGRIF SETTE tests to run. Also merged prettification changes to NST routines from the dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps branch. Neither AGRIF_DEMO nor VORTEX restart perfectly (drifting after 8 and 121 timesteps, respectively).

Location:
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_oce_interp.F90

    r11099 r11463  
    110110            DO jk = 1, jpkm1 
    111111               DO jj = 1, jpj 
    112                   uu_b(ibdy1:ibdy2,jj,Krhs_a) = uu_b(ibdy1:ibdy2,jj,Krhs_a) &  
    113                       & + e3u(ibdy1:ibdy2,jj,jk,Krhs_a) * uu(ibdy1:ibdy2,jj,jk,Krhs_a) * umask(ibdy1:ibdy2,jj,jk) 
     112                  uu_b(ibdy1:ibdy2,jj,Krhs_a) = uu_b(ibdy1:ibdy2,jj,   Krhs_a) &  
     113                      &                        + e3u(ibdy1:ibdy2,jj,jk,Krhs_a) & 
     114                      &                        *  uu(ibdy1:ibdy2,jj,jk,Krhs_a) * umask(ibdy1:ibdy2,jj,jk) 
    114115               END DO 
    115116            END DO 
    … …  
    122123            DO jk=1,jpkm1              ! Smooth 
    123124               DO jj=j1,j2 
    124                   uu(ibdy2,jj,jk,Krhs_a) = 0.25_wp*(uu(ibdy2-1,jj,jk,Krhs_a)+2._wp*uu(ibdy2,jj,jk,Krhs_a)+uu(ibdy2+1,jj,jk,Krhs_a)) 
     125                  uu(ibdy2,jj,jk,Krhs_a) = 0.25_wp*( uu(ibdy2-1,jj,jk,Krhs_a)+2._wp*uu(ibdy2,jj,jk,Krhs_a) & 
     126                      &                            + uu(ibdy2+1,jj,jk,Krhs_a) ) 
    125127               END DO 
    126128            END DO 
    … …  
    130132         DO jk = 1, jpkm1 
    131133            DO jj = 1, jpj 
    132                zub(ibdy1:ibdy2,jj) = zub(ibdy1:ibdy2,jj) &  
    133                   & + e3u(ibdy1:ibdy2,jj,jk,Krhs_a)  * uu(ibdy1:ibdy2,jj,jk,Krhs_a)*umask(ibdy1:ibdy2,jj,jk) 
     134               zub(ibdy1:ibdy2,jj) = zub(ibdy1:ibdy2,jj)           +   e3u(ibdy1:ibdy2,jj,jk,Krhs_a)  &  
     135                 &                  * uu(ibdy1:ibdy2,jj,jk,Krhs_a) * umask(ibdy1:ibdy2,jj,jk) 
    134136            END DO 
    135137         END DO 
    … …  
    140142         DO jk = 1, jpkm1 
    141143            DO jj = 1, jpj 
    142                uu(ibdy1:ibdy2,jj,jk,Krhs_a) = ( uu(ibdy1:ibdy2,jj,jk,Krhs_a) & 
    143                  & + uu_b(ibdy1:ibdy2,jj,Krhs_a)-zub(ibdy1:ibdy2,jj)) * umask(ibdy1:ibdy2,jj,jk) 
     144               uu(ibdy1:ibdy2,jj,jk,Krhs_a) = (    uu(ibdy1:ibdy2,jj,jk,Krhs_a)   & 
     145                 &                             + uu_b(ibdy1:ibdy2,jj,   Krhs_a)   & 
     146                 &                             -  zub(ibdy1:ibdy2,jj)           ) & 
     147                 &                            * umask(ibdy1:ibdy2,jj,jk) 
    144148            END DO 
    145149         END DO 
    … …  
    149153            DO jk = 1, jpkm1 
    150154               DO jj = 1, jpj 
    151                   zvb(ibdy1:ibdy2,jj) = zvb(ibdy1:ibdy2,jj) &  
    152                      & + e3v(ibdy1:ibdy2,jj,jk,Krhs_a) * vv(ibdy1:ibdy2,jj,jk,Krhs_a) * vmask(ibdy1:ibdy2,jj,jk) 
     155                  zvb(ibdy1:ibdy2,jj) = zvb(ibdy1:ibdy2,jj)           +   e3v(ibdy1:ibdy2,jj,jk,Krhs_a)  & 
     156                    &                  * vv(ibdy1:ibdy2,jj,jk,Krhs_a) * vmask(ibdy1:ibdy2,jj,jk) 
    153157               END DO 
    154158            END DO 
    … …  
    158162            DO jk = 1, jpkm1 
    159163               DO jj = 1, jpj 
    160                   vv(ibdy1:ibdy2,jj,jk,Krhs_a) = ( vv(ibdy1:ibdy2,jj,jk,Krhs_a) &  
    161                     & + vv_b(ibdy1:ibdy2,jj,Krhs_a)-zvb(ibdy1:ibdy2,jj))*vmask(ibdy1:ibdy2,jj,jk) 
     164                  vv(ibdy1:ibdy2,jj,jk,Krhs_a) = (    vv(ibdy1:ibdy2,jj,jk,Krhs_a)   &  
     165                    &                             + vv_b(ibdy1:ibdy2,jj,   Krhs_a)   & 
     166                    &                             -  zvb(ibdy1:ibdy2,jj)           ) & 
     167                    &                            * vmask(ibdy1:ibdy2,jj,jk) 
    162168               END DO 
    163169            END DO 
    … …  
    181187            DO jk = 1, jpkm1 
    182188               DO jj = 1, jpj 
    183                   uu_b(ibdy1:ibdy2,jj,Krhs_a) = uu_b(ibdy1:ibdy2,jj,Krhs_a) &  
    184                       & + e3u(ibdy1:ibdy2,jj,jk,Krhs_a) * uu(ibdy1:ibdy2,jj,jk,Krhs_a) * umask(ibdy1:ibdy2,jj,jk) 
     189                  uu_b(ibdy1:ibdy2,jj,Krhs_a) =   uu_b(ibdy1:ibdy2,jj,   Krhs_a) &  
     190                      &                        +   e3u(ibdy1:ibdy2,jj,jk,Krhs_a) & 
     191                      &                        *    uu(ibdy1:ibdy2,jj,jk,Krhs_a) & 
     192                      &                        * umask(ibdy1:ibdy2,jj,jk) 
    185193               END DO 
    186194            END DO 
    … …  
    193201            DO jk=1,jpkm1              ! Smooth 
    194202               DO jj=j1,j2 
    195                   uu(ibdy1,jj,jk,Krhs_a) = 0.25_wp*(uu(ibdy1-1,jj,jk,Krhs_a)+2._wp*uu(ibdy1,jj,jk,Krhs_a)+uu(ibdy1+1,jj,jk,Krhs_a)) 
     203                  uu(ibdy1,jj,jk,Krhs_a) = 0.25_wp*(        uu(ibdy1-1,jj,jk,Krhs_a)  & 
     204                    &                               + 2._wp*uu(ibdy1  ,jj,jk,Krhs_a)  & 
     205                    &                                     + uu(ibdy1+1,jj,jk,Krhs_a) ) 
    196206               END DO 
    197207            END DO 
    … …  
    201211         DO jk = 1, jpkm1 
    202212            DO jj = 1, jpj 
    203                zub(ibdy1:ibdy2,jj) = zub(ibdy1:ibdy2,jj) &  
    204                   & + e3u(ibdy1:ibdy2,jj,jk,Krhs_a)  * uu(ibdy1:ibdy2,jj,jk,Krhs_a) * umask(ibdy1:ibdy2,jj,jk) 
     213               zub(ibdy1:ibdy2,jj) =  zub(ibdy1:ibdy2,jj)                                      &  
     214                  &                 + e3u(ibdy1:ibdy2,jj,jk,Krhs_a)                            & 
     215                  &                 *  uu(ibdy1:ibdy2,jj,jk,Krhs_a) * umask(ibdy1:ibdy2,jj,jk) 
    205216            END DO 
    206217         END DO 
    … …  
    211222         DO jk = 1, jpkm1 
    212223            DO jj = 1, jpj 
    213                uu(ibdy1:ibdy2,jj,jk,Krhs_a) = ( uu(ibdy1:ibdy2,jj,jk,Krhs_a) &  
    214                  & + uu_b(ibdy1:ibdy2,jj,Krhs_a)-zub(ibdy1:ibdy2,jj))*umask(ibdy1:ibdy2,jj,jk) 
     224               uu(ibdy1:ibdy2,jj,jk,Krhs_a) = (      uu(ibdy1:ibdy2,jj,jk,Krhs_a) &  
     225                 &                             +   uu_b(ibdy1:ibdy2,jj,   Krhs_a) & 
     226                 &                             -    zub(ibdy1:ibdy2,jj)           & 
     227                 &                            ) * umask(ibdy1:ibdy2,jj,jk) 
    215228            END DO 
    216229         END DO 
    … …  
    222235            DO jk = 1, jpkm1 
    223236               DO jj = 1, jpj 
    224                   zvb(ibdy1:ibdy2,jj) = zvb(ibdy1:ibdy2,jj) & 
    225                      & + e3v(ibdy1:ibdy2,jj,jk,Krhs_a) * vv(ibdy1:ibdy2,jj,jk,Krhs_a) * vmask(ibdy1:ibdy2,jj,jk) 
     237                  zvb(ibdy1:ibdy2,jj) =    zvb(ibdy1:ibdy2,jj)                     & 
     238                     &                 +   e3v(ibdy1:ibdy2,jj,jk,Krhs_a)           & 
     239                     &                 *    vv(ibdy1:ibdy2,jj,jk,Krhs_a)           & 
     240                     &                 * vmask(ibdy1:ibdy2,jj,jk) 
    226241               END DO 
    227242            END DO 
    … …  
    231246            DO jk = 1, jpkm1 
    232247               DO jj = 1, jpj 
    233                   vv(ibdy1:ibdy2,jj,jk,Krhs_a) = ( vv(ibdy1:ibdy2,jj,jk,Krhs_a) &  
    234                       & + vv_b(ibdy1:ibdy2,jj,Krhs_a)-zvb(ibdy1:ibdy2,jj)) * vmask(ibdy1:ibdy2,jj,jk) 
     248                  vv(ibdy1:ibdy2,jj,jk,Krhs_a) = (      vv(ibdy1:ibdy2,jj,jk,Krhs_a) &  
     249                      &                           +   vv_b(ibdy1:ibdy2,jj,   Krhs_a) & 
     250                      &                           -    zvb(ibdy1:ibdy2,jj)           & 
     251                      &                          ) * vmask(ibdy1:ibdy2,jj,jk) 
    235252               END DO 
    236253            END DO 
    … …  
    254271            DO jk = 1, jpkm1 
    255272               DO ji = 1, jpi 
    256                   vv_b(ji,jbdy1:jbdy2,Krhs_a) = vv_b(ji,jbdy1:jbdy2,Krhs_a) &  
    257                       & + e3v(ji,jbdy1:jbdy2,jk,Krhs_a) * vv(ji,jbdy1:jbdy2,jk,Krhs_a) * vmask(ji,jbdy1:jbdy2,jk) 
     273                  vv_b(ji,jbdy1:jbdy2,Krhs_a) =   vv_b(ji,jbdy1:jbdy2,   Krhs_a) &  
     274                      &                        +   e3v(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     275                      &                        *    vv(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     276                      &                        * vmask(ji,jbdy1:jbdy2,jk) 
    258277               END DO 
    259278            END DO 
    … …  
    266285            DO jk = 1, jpkm1           ! Smooth 
    267286               DO ji = i1, i2 
    268                   vv(ji,jbdy2,jk,Krhs_a) = 0.25_wp*(vv(ji,jbdy2-1,jk,Krhs_a)+2._wp*vv(ji,jbdy2,jk,Krhs_a)+vv(ji,jbdy2+1,jk,Krhs_a)) 
     287                  vv(ji,jbdy2,jk,Krhs_a) = 0.25_wp*(        vv(ji,jbdy2-1,jk,Krhs_a)  & 
     288                    &                               + 2._wp*vv(ji,jbdy2  ,jk,Krhs_a)  & 
     289                    &                               +       vv(ji,jbdy2+1,jk,Krhs_a) ) 
    269290               END DO 
    270291            END DO 
    … …  
    274295         DO jk=1,jpkm1 
    275296            DO ji=1,jpi 
    276                zvb(ji,jbdy1:jbdy2) = zvb(ji,jbdy1:jbdy2) &  
    277                   & + e3v(ji,jbdy1:jbdy2,jk,Krhs_a) * vv(ji,jbdy1:jbdy2,jk,Krhs_a) * vmask(ji,jbdy1:jbdy2,jk) 
     297               zvb(ji,jbdy1:jbdy2) =    zvb(ji,jbdy1:jbdy2)           &  
     298                  &                 +   e3v(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     299                  &                 *    vv(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     300                  &                 * vmask(ji,jbdy1:jbdy2,jk) 
    278301            END DO 
    279302         END DO 
    … …  
    284307         DO jk = 1, jpkm1 
    285308            DO ji = 1, jpi 
    286                vv(ji,jbdy1:jbdy2,jk,Krhs_a) = ( vv(ji,jbdy1:jbdy2,jk,Krhs_a) &  
    287                  & + vv_b(ji,jbdy1:jbdy2,Krhs_a) - zvb(ji,jbdy1:jbdy2) ) * vmask(ji,jbdy1:jbdy2,jk) 
     309               vv(ji,jbdy1:jbdy2,jk,Krhs_a) = (      vv(ji,jbdy1:jbdy2,jk,Krhs_a) &  
     310                 &                             +   vv_b(ji,jbdy1:jbdy2,   Krhs_a) & 
     311                 &                             -    zvb(ji,jbdy1:jbdy2)           & 
     312                 &                            ) * vmask(ji,jbdy1:jbdy2,jk) 
    288313            END DO 
    289314         END DO 
    … …  
    293318            DO jk = 1, jpkm1 
    294319               DO ji = 1, jpi 
    295                   zub(ji,jbdy1:jbdy2) = zub(ji,jbdy1:jbdy2) &  
    296                      & + e3u(ji,jbdy1:jbdy2,jk,Krhs_a) * uu(ji,jbdy1:jbdy2,jk,Krhs_a) * umask(ji,jbdy1:jbdy2,jk) 
     320                  zub(ji,jbdy1:jbdy2) =    zub(ji,jbdy1:jbdy2)           &  
     321                     &                 +   e3u(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     322                     &                 *    uu(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     323                     &                 * umask(ji,jbdy1:jbdy2,jk) 
    297324               END DO 
    298325            END DO 
    … …  
    303330            DO jk = 1, jpkm1 
    304331               DO ji = 1, jpi 
    305                   uu(ji,jbdy1:jbdy2,jk,Krhs_a) = ( uu(ji,jbdy1:jbdy2,jk,Krhs_a) &  
    306                     & + uu_b(ji,jbdy1:jbdy2,Krhs_a) - zub(ji,jbdy1:jbdy2) ) * umask(ji,jbdy1:jbdy2,jk) 
     332                  uu(ji,jbdy1:jbdy2,jk,Krhs_a) = (      uu(ji,jbdy1:jbdy2,jk,Krhs_a) &  
     333                    &                             +   uu_b(ji,jbdy1:jbdy2,   Krhs_a) & 
     334                    &                             -    zub(ji,jbdy1:jbdy2)           & 
     335                    &                            ) * umask(ji,jbdy1:jbdy2,jk) 
    307336               END DO 
    308337            END DO 
    … …  
    326355            DO jk = 1, jpkm1 
    327356               DO ji = 1, jpi 
    328                   vv_b(ji,jbdy1:jbdy2,Krhs_a) = vv_b(ji,jbdy1:jbdy2,Krhs_a) &  
    329                       & + e3v(ji,jbdy1:jbdy2,jk,Krhs_a) * vv(ji,jbdy1:jbdy2,jk,Krhs_a) * vmask(ji,jbdy1:jbdy2,jk) 
     357                  vv_b(ji,jbdy1:jbdy2,Krhs_a) =   vv_b(ji,jbdy1:jbdy2,   Krhs_a) &  
     358                      &                        +   e3v(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     359                      &                        *    vv(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     360                      &                        * vmask(ji,jbdy1:jbdy2,jk) 
    330361               END DO 
    331362            END DO 
    … …  
    338369            DO jk = 1, jpkm1           ! Smooth 
    339370               DO ji = i1, i2 
    340                   vv(ji,jbdy1,jk,Krhs_a) = 0.25_wp*(vv(ji,jbdy1-1,jk,Krhs_a)+2._wp*vv(ji,jbdy1,jk,Krhs_a)+vv(ji,jbdy1+1,jk,Krhs_a)) 
     371                  vv(ji,jbdy1,jk,Krhs_a) = 0.25_wp*(        vv(ji,jbdy1-1,jk,Krhs_a)  & 
     372                    &                               + 2._wp*vv(ji,jbdy1  ,jk,Krhs_a)  & 
     373                    &                               +       vv(ji,jbdy1+1,jk,Krhs_a) ) 
    341374               END DO 
    342375            END DO 
    … …  
    346379         DO jk=1,jpkm1 
    347380            DO ji=1,jpi 
    348                zvb(ji,jbdy1:jbdy2) = zvb(ji,jbdy1:jbdy2) &  
    349                   & + e3v(ji,jbdy1:jbdy2,jk,Krhs_a) * vv(ji,jbdy1:jbdy2,jk,Krhs_a) * vmask(ji,jbdy1:jbdy2,jk) 
     381               zvb(ji,jbdy1:jbdy2) =    zvb(ji,jbdy1:jbdy2)           &  
     382                  &                 +   e3v(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     383                  &                 *    vv(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     384                  &                 * vmask(ji,jbdy1:jbdy2,jk) 
    350385            END DO 
    351386         END DO 
    … …  
    356391         DO jk = 1, jpkm1 
    357392            DO ji = 1, jpi 
    358                vv(ji,jbdy1:jbdy2,jk,Krhs_a) = ( vv(ji,jbdy1:jbdy2,jk,Krhs_a) &  
    359                  & + vv_b(ji,jbdy1:jbdy2,Krhs_a) - zvb(ji,jbdy1:jbdy2) ) * vmask(ji,jbdy1:jbdy2,jk) 
     393               vv(ji,jbdy1:jbdy2,jk,Krhs_a) = (      vv(ji,jbdy1:jbdy2,jk,Krhs_a) &  
     394                 &                             +   vv_b(ji,jbdy1:jbdy2,   Krhs_a) & 
     395                 &                             -    zvb(ji,jbdy1:jbdy2)           & 
     396                 &                            ) * vmask(ji,jbdy1:jbdy2,jk) 
    360397            END DO 
    361398         END DO 
    … …  
    367404            DO jk = 1, jpkm1 
    368405               DO ji = 1, jpi 
    369                   zub(ji,jbdy1:jbdy2) = zub(ji,jbdy1:jbdy2) &  
    370                      & + e3u(ji,jbdy1:jbdy2,jk,Krhs_a) * uu(ji,jbdy1:jbdy2,jk,Krhs_a) * umask(ji,jbdy1:jbdy2,jk) 
     406                  zub(ji,jbdy1:jbdy2) =    zub(ji,jbdy1:jbdy2)           &  
     407                     &                 +   e3u(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     408                     &                 *    uu(ji,jbdy1:jbdy2,jk,Krhs_a) & 
     409                     &                 * umask(ji,jbdy1:jbdy2,jk) 
    371410               END DO 
    372411            END DO 
    … …  
    377416            DO jk = 1, jpkm1 
    378417               DO ji = 1, jpi 
    379                   uu(ji,jbdy1:jbdy2,jk,Krhs_a) = ( uu(ji,jbdy1:jbdy2,jk,Krhs_a) &  
    380                     & + uu_b(ji,jbdy1:jbdy2,Krhs_a) - zub(ji,jbdy1:jbdy2) ) * umask(ji,jbdy1:jbdy2,jk) 
     418                  uu(ji,jbdy1:jbdy2,jk,Krhs_a) = (      uu(ji,jbdy1:jbdy2,jk,Krhs_a) &  
     419                    &                             +   uu_b(ji,jbdy1:jbdy2,   Krhs_a) & 
     420                    &                             -    zub(ji,jbdy1:jbdy2)           & 
     421                    &                            ) * umask(ji,jbdy1:jbdy2,jk) 
    381422               END DO 
    382423            END DO 
    … …  
    684725               iref = ji 
    685726               jref = jj 
    686                if(western_side) iref=MAX(2,ji) 
    687                if(eastern_side) iref=MIN(nlci-1,ji) 
    688                if(southern_side) jref=MAX(2,jj) 
     727               if(western_side)  iref=MAX(2     ,ji) 
     728               if(eastern_side)  iref=MIN(nlci-1,ji) 
     729               if(southern_side) jref=MAX(2     ,jj) 
    689730               if(northern_side) jref=MIN(nlcj-1,jj) 
    690731               N_in = 0 
    … …  
    693734                  N_in = N_in + 1 
    694735                  tabin(jk,:) = ptab(ji,jj,jk,n1:n2-1) 
    695                   h_in(N_in) = ptab(ji,jj,jk,n2) 
     736                  h_in(N_in)  = ptab(ji,jj,jk,n2) 
    696737               END DO 
    697738               N_out = 0 
    … …  
    722763            !  
    723764            ! Remove CORNERS 
    724             IF((nbondj == -1).OR.(nbondj == 2)) jmin = 2 + nbghostcells 
     765            IF((nbondj == -1).OR.(nbondj == 2)) jmin = 2    + nbghostcells 
    725766            IF((nbondj == +1).OR.(nbondj == 2)) jmax = nlcj - nbghostcells - 1 
    726             IF((nbondi == -1).OR.(nbondi == 2)) imin = 2 + nbghostcells 
     767            IF((nbondi == -1).OR.(nbondi == 2)) imin = 2    + nbghostcells 
    727768            IF((nbondi == +1).OR.(nbondi == 2)) imax = nlci - nbghostcells - 1       
    728769            ! 
    … …  
    737778               ibdy = nlci-nbghostcells 
    738779               DO jn = 1, jpts 
    739                   ts(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn) + z2 * ptab_child(ibdy,jmin:jmax,1:jpkm1,jn) 
     780                  ts(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn) & 
     781                    &                                     + z2 * ptab_child(ibdy  ,jmin:jmax,1:jpkm1,jn) 
    740782                  DO jk = 1, jpkm1 
    741783                     DO jj = jmin,jmax 
    … …  
    743785                           ts(ibdy,jj,jk,jn,Krhs_a) = ts(ibdy+1,jj,jk,jn,Krhs_a) * tmask(ibdy,jj,jk) 
    744786                        ELSE 
    745                            ts(ibdy,jj,jk,jn,Krhs_a)=(z4*ts(ibdy+1,jj,jk,jn,Krhs_a)+z3*ts(ibdy-1,jj,jk,jn,Krhs_a))*tmask(ibdy,jj,jk) 
     787                           ts(ibdy,jj,jk,jn,Krhs_a) = (   z4 * ts(ibdy+1,jj,jk,jn,Krhs_a) & 
     788                    &                                  +  z3 * ts(ibdy-1,jj,jk,jn,Krhs_a) & 
     789                    &                                 ) *   tmask(ibdy  ,jj,jk) 
    746790                           IF( uu(ibdy-1,jj,jk,Kmm_a) > 0._wp ) THEN 
    747                               ts(ibdy,jj,jk,jn,Krhs_a)=( z6*ts(ibdy-1,jj,jk,jn,Krhs_a)+z5*ts(ibdy+1,jj,jk,jn,Krhs_a) &  
    748                                                  + z7*ts(ibdy-2,jj,jk,jn,Krhs_a) ) * tmask(ibdy,jj,jk) 
     791                              ts(ibdy,jj,jk,jn,Krhs_a) = (  z6 * ts(ibdy-1,jj,jk,jn,Krhs_a) & 
     792                    &                                     + z5 * ts(ibdy+1,jj,jk,jn,Krhs_a) &  
     793                    &                                     + z7 * ts(ibdy-2,jj,jk,jn,Krhs_a) & 
     794                    &                                    ) *  tmask(ibdy  ,jj,jk) 
    749795                           ENDIF 
    750796                        ENDIF 
    … …  
    752798                  END DO 
    753799                  ! Restore ghost points: 
    754                   ts(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn) * tmask(ibdy+1,jmin:jmax,1:jpkm1) 
     800                  ts(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn) & 
     801                    &                                     *     tmask(ibdy+1,jmin:jmax,1:jpkm1) 
    755802               END DO 
    756803            ENDIF 
    … …  
    766813               jbdy = nlcj-nbghostcells          
    767814               DO jn = 1, jpts 
    768                   ts(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(imin:imax,jbdy+1,1:jpkm1,jn) + z2 * ptab_child(imin:imax,jbdy,1:jpkm1,jn) 
     815                  ts(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(imin:imax,jbdy+1,1:jpkm1,jn) & 
     816                    &                                     + z2 * ptab_child(imin:imax,jbdy  ,1:jpkm1,jn) 
    769817                  DO jk = 1, jpkm1 
    770818                     DO ji = imin,imax 
    … …  
    772820                           ts(ji,jbdy,jk,jn,Krhs_a) = ts(ji,jbdy+1,jk,jn,Krhs_a) * tmask(ji,jbdy,jk) 
    773821                        ELSE 
    774                            ts(ji,jbdy,jk,jn,Krhs_a)=(z4*ts(ji,jbdy+1,jk,jn,Krhs_a)+z3*ts(ji,jbdy-1,jk,jn,Krhs_a))*tmask(ji,jbdy,jk)         
     822                           ts(ji,jbdy,jk,jn,Krhs_a)=(  z4 * ts(ji,jbdy+1,jk,jn,Krhs_a)   & 
     823                    &                                + z3 * ts(ji,jbdy-1,jk,jn,Krhs_a)   & 
     824                    &                               ) *  tmask(ji,jbdy  ,jk)         
    775825                           IF (vv(ji,jbdy-1,jk,Kmm_a) > 0._wp ) THEN 
    776                               ts(ji,jbdy,jk,jn,Krhs_a)=( z6*ts(ji,jbdy-1,jk,jn,Krhs_a)+z5*ts(ji,jbdy+1,jk,jn,Krhs_a)  & 
    777                                                  + z7*ts(ji,jbdy-2,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk) 
     826                              ts(ji,jbdy,jk,jn,Krhs_a)=(  z6 * ts(ji,jbdy-1,jk,jn,Krhs_a) & 
     827                    &                                   + z5 * ts(ji,jbdy+1,jk,jn,Krhs_a) & 
     828                    &                                   + z7 * ts(ji,jbdy-2,jk,jn,Krhs_a) & 
     829                    &                                  ) *  tmask(ji,jbdy  ,jk) 
    778830                           ENDIF 
    779831                        ENDIF 
    … …  
    781833                  END DO 
    782834                  ! Restore ghost points: 
    783                   ts(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy+1,1:jpkm1,jn) * tmask(imin:imax,jbdy+1,1:jpkm1) 
     835                  ts(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy+1,1:jpkm1,jn) & 
     836                    &                                     *     tmask(imin:imax,jbdy+1,1:jpkm1) 
    784837               END DO 
    785838            ENDIF 
    … …  
    795848               ibdy = 1+nbghostcells        
    796849               DO jn = 1, jpts 
    797                   ts(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn) + z2 * ptab_child(ibdy,jmin:jmax,1:jpkm1,jn) 
     850                  ts(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn) & 
     851                    &                                     + z2 * ptab_child(ibdy  ,jmin:jmax,1:jpkm1,jn) 
    798852                  DO jk = 1, jpkm1 
    799853                     DO jj = jmin,jmax 
    … …  
    801855                           ts(ibdy,jj,jk,jn,Krhs_a) = ts(ibdy-1,jj,jk,jn,Krhs_a) * tmask(ibdy,jj,jk) 
    802856                        ELSE 
    803                            ts(ibdy,jj,jk,jn,Krhs_a)=(z4*ts(ibdy-1,jj,jk,jn,Krhs_a)+z3*ts(ibdy+1,jj,jk,jn,Krhs_a))*tmask(ibdy,jj,jk)         
     857                           ts(ibdy,jj,jk,jn,Krhs_a) = (  z4 * ts(ibdy-1,jj,jk,jn,Krhs_a) & 
     858                    &                                  + z3 * ts(ibdy+1,jj,jk,jn,Krhs_a) & 
     859                    &                                 ) *  tmask(ibdy  ,jj,jk)         
    804860                           IF( uu(ibdy,jj,jk,Kmm_a) < 0._wp ) THEN 
    805                               ts(ibdy,jj,jk,jn,Krhs_a)=( z6*ts(ibdy+1,jj,jk,jn,Krhs_a)+z5*ts(ibdy-1,jj,jk,jn,Krhs_a) & 
    806                                                  + z7*ts(ibdy+2,jj,jk,jn,Krhs_a) ) * tmask(ibdy,jj,jk) 
     861                              ts(ibdy,jj,jk,jn,Krhs_a) = (  z6 * ts(ibdy+1,jj,jk,jn,Krhs_a) & 
     862                    &                                     + z5 * ts(ibdy-1,jj,jk,jn,Krhs_a) & 
     863                    &                                     + z7 * ts(ibdy+2,jj,jk,jn,Krhs_a) & 
     864                    &                                    ) *  tmask(ibdy,jj,jk) 
    807865                           ENDIF 
    808866                        ENDIF 
    … …  
    810868                  END DO 
    811869                  ! Restore ghost points: 
    812                   ts(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn) * tmask(ibdy-1,jmin:jmax,1:jpkm1) 
     870                  ts(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn) & 
     871                    &                                     *     tmask(ibdy-1,jmin:jmax,1:jpkm1) 
    813872               END DO 
    814873            ENDIF 
    … …  
    824883               jbdy=1+nbghostcells         
    825884               DO jn = 1, jpts 
    826                   ts(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(imin:imax,jbdy-1,1:jpkm1,jn) + z2 * ptab_child(imin:imax,jbdy,1:jpkm1,jn) 
     885                  ts(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(imin:imax,jbdy-1,1:jpkm1,jn) & 
     886                    &                                     + z2 * ptab_child(imin:imax,jbdy  ,1:jpkm1,jn) 
    827887                  DO jk = 1, jpkm1       
    828888                     DO ji = imin,imax 
    829889                        IF( vmask(ji,jbdy,jk) == 0._wp ) THEN 
    830                            ts(ji,jbdy,jk,jn,Krhs_a)=ts(ji,jbdy-1,jk,jn,Krhs_a) * tmask(ji,jbdy,jk) 
     890                           ts(ji,jbdy,jk,jn,Krhs_a) = ts(ji,jbdy-1,jk,jn,Krhs_a) * tmask(ji,jbdy,jk) 
    831891                        ELSE 
    832                            ts(ji,jbdy,jk,jn,Krhs_a)=(z4*ts(ji,jbdy-1,jk,jn,Krhs_a)+z3*ts(ji,jbdy+1,jk,jn,Krhs_a))*tmask(ji,jbdy,jk) 
     892                           ts(ji,jbdy,jk,jn,Krhs_a) = (  z4 * ts(ji,jbdy-1,jk,jn,Krhs_a) & 
     893                    &                                  + z3 * ts(ji,jbdy+1,jk,jn,Krhs_a) & 
     894                    &                                 ) *  tmask(ji,jbdy  ,jk) 
    833895                           IF( vv(ji,jbdy,jk,Kmm_a) < 0._wp ) THEN 
    834                               ts(ji,jbdy,jk,jn,Krhs_a)=( z6*ts(ji,jbdy+1,jk,jn,Krhs_a)+z5*ts(ji,jbdy-1,jk,jn,Krhs_a) &  
    835                                                  + z7*ts(ji,jbdy+2,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk) 
     896                              ts(ji,jbdy,jk,jn,Krhs_a) = (  z6 * ts(ji,jbdy+1,jk,jn,Krhs_a) & 
     897                    &                                     + z5 * ts(ji,jbdy-1,jk,jn,Krhs_a) &  
     898                    &                                     + z7 * ts(ji,jbdy+2,jk,jn,Krhs_a) & 
     899                    &                                    ) *  tmask(ji,jbdy  ,jk) 
    836900                           ENDIF 
    837901                        ENDIF 
    … …  
    839903                  END DO 
    840904                  ! Restore ghost points: 
    841                   ts(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy-1,1:jpkm1,jn) * tmask(imin:imax,jbdy-1,1:jpkm1) 
     905                  ts(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy-1,1:jpkm1,jn) & 
     906                    &                                     *     tmask(imin:imax,jbdy-1,1:jpkm1) 
    842907               END DO 
    843908            ENDIF 
    … …  
    900965            DO jj=j1,j2 
    901966               DO ji=i1,i2 
    902                   ptab(ji,jj,jk,1) = (e2u(ji,jj) * e3u(ji,jj,jk,Kmm_a) * uu(ji,jj,jk,Kmm_a)*umask(ji,jj,jk))  
     967                  ptab(ji,jj,jk,1) = (  e2u(ji,jj)          *   e3u(ji,jj,jk,Kmm_a)  & 
     968                    &                 *  uu(ji,jj,jk,Kmm_a) * umask(ji,jj,jk)      )  
    903969# if defined key_vertical 
    904970                  ptab(ji,jj,jk,2) = (umask(ji,jj,jk) * e2u(ji,jj) * e3u(ji,jj,jk,Kmm_a)) 
    … …  
    923989                  IF (ptab(ji,jj,jk,2) == 0) EXIT 
    924990                  N_in = N_in + 1 
    925                   tabin(jk) = ptab(ji,jj,jk,1)/ptab(ji,jj,jk,2) 
     991                  tabin(jk)  = ptab(ji,jj,jk,1)/ptab(ji,jj,jk,2) 
    926992                  h_in(N_in) = ptab(ji,jj,jk,2)/(e2u(ji,jj)*zrhoy)  
    927993              ENDDO 
    … …  
    9451011          
    9461012              IF (N_in * N_out > 0) THEN 
    947                  h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
     1013                 h_diff = SUM( h_out(1:N_out) ) - SUM( h_in(1:N_in) ) 
    9481014! Should be able to remove the next IF/ELSEIF statement once scale factors are dealt with properly 
    9491015                 if (h_diff < -1.e4) then 
    950                     print *,'CHECK YOUR BATHY ...', h_diff, sum(h_out(1:N_out)), sum(h_in(1:N_in)) 
     1016                    print *,'CHECK YOUR BATHY ...', h_diff, SUM( h_out(1:N_out) ), SUM( h_in(1:N_in) ) 
    9511017!                    stop 
    9521018                 endif 
    9531019              ENDIF 
    954               call reconstructandremap(tabin(1:N_in),h_in(1:N_in),uu(ji,jj,1:N_out,Krhs_a),h_out(1:N_out),N_in,N_out) 
     1020              call reconstructandremap( tabin(1:N_in) , h_in(1:N_in), uu(ji,jj,1:N_out,Krhs_a), & 
     1021                    &                   h_out(1:N_out), N_in        , N_out                    ) 
    9551022            ENDDO 
    9561023         ENDDO 
    … …  
    9921059            DO jj=j1,j2 
    9931060               DO ji=i1,i2 
    994                   ptab(ji,jj,jk,1) = (e1v(ji,jj) * e3v(ji,jj,jk,Kmm_a) * vv(ji,jj,jk,Kmm_a)*vmask(ji,jj,jk)) 
     1061                  ptab(ji,jj,jk,1) = (  e1v(ji,jj)          *   e3v(ji,jj,jk,Kmm_a) & 
     1062                    &                 *  vv(ji,jj,jk,Kmm_a) * vmask(ji,jj,jk)      ) 
    9951063# if defined key_vertical 
    9961064                  ptab(ji,jj,jk,2) = vmask(ji,jj,jk) * e1v(ji,jj) * e3v(ji,jj,jk,Kmm_a) 
    … …  
    10151083                  if (ptab(ji,jj,jk,2) == 0) EXIT 
    10161084                  N_in = N_in + 1 
    1017                   tabin(jk) = ptab(ji,jj,jk,1)/ptab(ji,jj,jk,2) 
     1085                  tabin(jk)  = ptab(ji,jj,jk,1)/ptab(ji,jj,jk,2) 
    10181086                  h_in(N_in) = ptab(ji,jj,jk,2)/(e1v(ji,jj)*zrhox) 
    10191087               END DO 
    … …  
    10331101                 CYCLE 
    10341102               ENDIF 
    1035                call reconstructandremap(tabin(1:N_in),h_in(1:N_in),vv(ji,jj,1:N_out,Krhs_a),h_out(1:N_out),N_in,N_out) 
     1103               call reconstructandremap( tabin(1:N_in) , h_in(1:N_in), vv(ji,jj,1:N_out,Krhs_a), & 
     1104                    &                    h_out(1:N_out), N_in        , N_out                    ) 
    10361105            END DO 
    10371106         END DO 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_oce_sponge.F90

    r11053 r11463  
    393393          
    394394              IF (N_in * N_out > 0) THEN 
    395                  h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
     395                 h_diff = SUM( h_out(1:N_out) ) - SUM( h_in(1:N_in) ) 
    396396                 if (h_diff < -1.e4) then 
    397                     print *,'CHECK YOUR BATHY ...', h_diff, sum(h_out(1:N_out)), sum(h_in(1:N_in)) 
     397                    print *,'CHECK YOUR BATHY ...', h_diff, SUM( h_out(1:N_out) ), SUM( h_in(1:N_in) ) 
    398398                 endif 
    399399              ENDIF 
    … …  
    403403         ENDDO 
    404404 
    405          ubdiff(i1:i2,j1:j2,:) = (uu(i1:i2,j1:j2,:,Kbb_a) - tabres_child(i1:i2,j1:j2,:))*umask(i1:i2,j1:j2,:) 
     405         ubdiff(i1:i2,j1:j2,:) = ( uu(i1:i2,j1:j2,:,Kbb_a) - tabres_child(i1:i2,j1:j2,:  ) )*umask(i1:i2,j1:j2,:) 
    406406#else 
    407          ubdiff(i1:i2,j1:j2,:) = (uu(i1:i2,j1:j2,:,Kbb_a) - tabres(i1:i2,j1:j2,:,1))*umask(i1:i2,j1:j2,:) 
     407         ubdiff(i1:i2,j1:j2,:) = ( uu(i1:i2,j1:j2,:,Kbb_a) -       tabres(i1:i2,j1:j2,:,1) )*umask(i1:i2,j1:j2,:) 
    408408#endif 
    409409         ! 
    … …  
    540540          
    541541              IF (N_in * N_out > 0) THEN 
    542                  h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
     542                 h_diff = SUM( h_out(1:N_out) ) - SUM( h_in(1:N_in) ) 
    543543                 if (h_diff < -1.e4) then 
    544                     print *,'CHECK YOUR BATHY ...', h_diff, sum(h_out(1:N_out)), sum(h_in(1:N_in)) 
     544                    print *,'CHECK YOUR BATHY ...', h_diff, SUM( h_out(1:N_out) ), SUM( h_in(1:N_in) ) 
    545545                 endif 
    546546              ENDIF 
    … …  
    549549         ENDDO 
    550550 
    551          vbdiff(i1:i2,j1:j2,:) = (vv(i1:i2,j1:j2,:,Kbb_a) - tabres_child(i1:i2,j1:j2,:))*vmask(i1:i2,j1:j2,:)   
     551         vbdiff(i1:i2,j1:j2,:) = ( vv(i1:i2,j1:j2,:,Kbb_a) - tabres_child(i1:i2,j1:j2,:  ) )*vmask(i1:i2,j1:j2,:)   
    552552# else 
    553          vbdiff(i1:i2,j1:j2,:) = (vv(i1:i2,j1:j2,:,Kbb_a) - tabres(i1:i2,j1:j2,:,1))*vmask(i1:i2,j1:j2,:) 
     553         vbdiff(i1:i2,j1:j2,:) = ( vv(i1:i2,j1:j2,:,Kbb_a) -       tabres(i1:i2,j1:j2,:,1) )*vmask(i1:i2,j1:j2,:) 
    554554# endif 
    555555         ! 
    … …  
    570570               DO ji = i1,i2-1   ! vector opt. 
    571571                  zbtr = r1_e1e2f(ji,jj) * e3f(ji,jj,jk) * fsahm_spf(ji,jj) 
    572                   rotdiff(ji,jj,jk) = ( e2v(ji+1,jj) * vbdiff(ji+1,jj,jk) &  
    573                                     &  -e2v(ji  ,jj) * vbdiff(ji  ,jj,jk)  ) * fmask(ji,jj,jk) * zbtr 
     572                  rotdiff(ji,jj,jk) = (  e2v(ji+1,jj) * vbdiff(ji+1,jj,jk) &  
     573                                    &  - e2v(ji  ,jj) * vbdiff(ji  ,jj,jk)  ) * fmask(ji,jj,jk) * zbtr 
    574574               END DO 
    575575            END DO 
    … …  
    586586               IF( .NOT. tabspongedone_u(ji,jj) ) THEN 
    587587                  DO jk = 1, jpkm1 
    588                      uu(ji,jj,jk,Krhs_a) = uu(ji,jj,jk,Krhs_a)                                                               & 
    589                         & - ( rotdiff (ji  ,jj,jk) - rotdiff (ji,jj-1,jk)) / ( e2u(ji,jj) * e3u(ji,jj,jk,Kmm_a) )  & 
    590                         & + ( hdivdiff(ji+1,jj,jk) - hdivdiff(ji,jj  ,jk)) * r1_e1u(ji,jj) 
     588                     uu(ji,jj,jk,Krhs_a) =          uu(ji  ,jj,jk,Krhs_a)                                   & 
     589                        &                 - (  rotdiff(ji  ,jj,jk) -  rotdiff(ji,jj-1,jk) )                 & 
     590                        &                   / (    e2u(ji  ,jj)    *      e3u(ji,jj  ,jk,Kmm_a) )           & 
     591                        &                 + ( hdivdiff(ji+1,jj,jk) - hdivdiff(ji,jj  ,jk) ) * r1_e1u(ji,jj) 
    591592                  END DO 
    592593               ENDIF 
    … …  
    600601               IF( .NOT. tabspongedone_v(ji,jj) ) THEN 
    601602                  DO jk = 1, jpkm1 
    602                      vv(ji,jj,jk,Krhs_a) = vv(ji,jj,jk,Krhs_a)                                                                  & 
    603                         &  + ( rotdiff (ji,jj  ,jk) - rotdiff (ji-1,jj,jk) ) / ( e1v(ji,jj) * e3v(ji,jj,jk,Kmm_a) )   & 
    604                         &  + ( hdivdiff(ji,jj+1,jk) - hdivdiff(ji  ,jj,jk) ) * r1_e2v(ji,jj) 
     603                     vv(ji,jj,jk,Krhs_a) =          vv(ji,jj  ,jk,Krhs_a)                                   & 
     604                        &                 + (  rotdiff(ji,jj  ,jk) -  rotdiff(ji-1,jj,jk) )                 & 
     605                        &                   / (    e1v(ji,jj  ) *         e3v(ji  ,jj,jk,Kmm_a) )           & 
     606                        &                 + ( hdivdiff(ji,jj+1,jk) - hdivdiff(ji  ,jj,jk) ) * r1_e2v(ji,jj) 
    605607                  END DO 
    606608               ENDIF 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_oce_update.F90

    r11099 r11463  
    315315               DO jj=j1,j2 
    316316                  DO ji=i1,i2 
    317                      tabres(ji,jj,jk,jn) = (ts(ji,jj,jk,jn,Kmm_a) * e3t(ji,jj,jk,Kmm_a) ) & 
    318                                            * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 
     317                     tabres(ji,jj,jk,jn) = (    ts(ji,jj,jk,jn,Kmm_a) * e3t(ji,jj,jk,Kmm_a) )          & 
     318                       &                   * tmask(ji,jj,jk)       + (tmask(ji,jj,jk) - 1._wp)*999._wp 
    319319                  END DO 
    320320               END DO 
    … …  
    325325               DO ji=i1,i2 
    326326                  tabres(ji,jj,jk,n2) = tmask(ji,jj,jk) * e3t(ji,jj,jk,Kmm_a) & 
    327                                            + (tmask(ji,jj,jk)-1)*999._wp 
     327                    &                      + (tmask(ji,jj,jk) - 1._wp)*999._wp 
    328328               END DO 
    329329            END DO 
    330330         END DO 
    331331      ELSE 
    332          tabres_child(:,:,:,:) = 0. 
     332         tabres_child(:,:,:,:) = 0._wp 
    333333         AGRIF_SpecialValue = 0._wp 
    334334         DO jj=j1,j2 
    … …  
    356356                  ENDIF 
    357357                  DO jn=n1,n2-1 
    358                      CALL reconstructandremap(tabin(1:N_in,jn),h_in(1:N_in),tabres_child(ji,jj,1:N_out,jn),h_out(1:N_out),N_in,N_out) 
     358                     CALL reconstructandremap( tabin(1:N_in,jn), h_in(1:N_in), tabres_child(ji,jj,1:N_out,jn),  & 
     359                       &                       h_out(1:N_out)  , N_in        , N_out                           ) 
    359360                  ENDDO 
    360361               ENDIF 
    … …  
    369370                     DO ji=i1,i2 
    370371                        IF( tabres_child(ji,jj,jk,jn) .NE. 0. ) THEN 
    371                            ts(ji,jj,jk,jn,Kbb_a) = ts(ji,jj,jk,jn,Kbb_a) &  
    372                                  & + atfp * ( tabres_child(ji,jj,jk,jn) & 
    373                                  &          - ts(ji,jj,jk,jn,Kmm_a) ) * tmask(ji,jj,jk) 
     372                           ts(ji,jj,jk,jn,Kbb_a) =                     ts(ji,jj,jk,jn,Kbb_a) &  
     373                                 &                + atfp * ( tabres_child(ji,jj,jk,jn)       & 
     374                                 &                                   - ts(ji,jj,jk,jn,Kmm_a) & 
     375                                 &                         ) * tmask(ji,jj,jk) 
    374376                        ENDIF 
    375377                     ENDDO 
    … …  
    414416!> jc tmp 
    415417                     tabres(ji,jj,jk,jn) = ts(ji,jj,jk,jn,Kmm_a)  * e3t(ji,jj,jk,Kmm_a) / e3t_0(ji,jj,jk) 
    416 !                     tabres(ji,jj,jk,jn) = ts(ji,jj,jk,jn,Kmm_a)  * e3t(ji,jj,jk,Kmm_a) 
     418!                    tabres(ji,jj,jk,jn) = ts(ji,jj,jk,jn,Kmm_a)  * e3t(ji,jj,jk,Kmm_a) 
    417419!< jc tmp 
    418420                  END DO 
    … …  
    438440                           ztno = ts(ji,jj,jk,jn,Kmm_a) * e3t(ji,jj,jk,Krhs_a) 
    439441                           ts(ji,jj,jk,jn,Kbb_a) = ( ztb + atfp * ( ztnu - ztno) )  &  
    440                                      &        * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb_a) 
     442                                     &            * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb_a) 
    441443                        ENDIF 
    442444                     END DO 
    … …  
    496498               DO ji=i1,i2 
    497499                  tabres(ji,jj,jk,1) = zrhoy * e2u(ji,jj) * e3u(ji,jj,jk,Kmm_a) * umask(ji,jj,jk) * uu(ji,jj,jk,Kmm_a)  & 
    498                                        + (umask(ji,jj,jk)-1)*999._wp 
     500                    &                  + (umask(ji,jj,jk)-1._wp)*999._wp 
    499501                  tabres(ji,jj,jk,2) = zrhoy * umask(ji,jj,jk) * e2u(ji,jj) * e3u(ji,jj,jk,Kmm_a)  & 
    500                                        + (umask(ji,jj,jk)-1)*999._wp 
     502                    &                  + (umask(ji,jj,jk)-1._wp)*999._wp 
    501503               END DO 
    502504            END DO 
    … …  
    513515                  IF( tabres(ji,jj,jk,2) < -900) EXIT 
    514516                  N_in = N_in + 1 
    515                   tabin(jk) = tabres(ji,jj,jk,1)/tabres(ji,jj,jk,2) 
     517                  tabin(jk)  = tabres(ji,jj,jk,1)/tabres(ji,jj,jk,2) 
    516518                  h_in(N_in) = tabres(ji,jj,jk,2)/e2u(ji,jj) 
    517519               ENDDO 
    … …  
    523525               ENDDO 
    524526               IF (N_in * N_out > 0) THEN 
    525                   h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
     527                  h_diff = SUM( h_out(1:N_out) ) - SUM( h_in(1:N_in) ) 
    526528                  IF (h_diff < -1.e-4) THEN 
    527529!Even if bathy at T points match it's possible for the U points to be deeper in the child grid.  
    … …  
    540542                     ENDDO 
    541543                  ENDIF 
    542                   CALL reconstructandremap(tabin(1:N_in),h_in(1:N_in),tabres_child(ji,jj,1:N_out),h_out(1:N_out),N_in,N_out) 
    543                   tabres_child(ji,jj,N_out) = tabres_child(ji,jj,N_out) + excess/(e2u(ji,jj)*h_out(N_out)) 
     544                  CALL reconstructandremap( tabin(1:N_in) , h_in(1:N_in), tabres_child(ji,jj,1:N_out),       & 
     545                    &                       h_out(1:N_out), N_in        , N_out                            ) 
     546                  tabres_child(ji,jj,N_out) = tabres_child(ji,jj,N_out) + excess/( e2u(ji,jj)*h_out(N_out) ) 
    544547               ENDIF 
    545548            ENDDO 
    … …  
    550553               DO ji=i1,i2 
    551554                  IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) THEN ! Add asselin part 
    552                      uu(ji,jj,jk,Kbb_a) = uu(ji,jj,jk,Kbb_a) &  
    553                            & + atfp * ( tabres_child(ji,jj,jk) - uu(ji,jj,jk,Kmm_a) ) * umask(ji,jj,jk) 
     555                     uu(ji,jj,jk,Kbb_a) =                     uu(ji,jj,jk,Kbb_a)                                  &  
     556                           &             + atfp * ( tabres_child(ji,jj,jk) - uu(ji,jj,jk,Kmm_a) ) * umask(ji,jj,jk) 
    554557                  ENDIF 
    555558                  ! 
    … …  
    579582         zrhoy = Agrif_Rhoy() 
    580583         DO jk = k1, k2 
    581             tabres(i1:i2,j1:j2,jk,1) = zrhoy * e2u(i1:i2,j1:j2) * e3u(i1:i2,j1:j2,jk,Kmm_a) * uu(i1:i2,j1:j2,jk,Kmm_a) 
     584            tabres(i1:i2,j1:j2,jk,1) = zrhoy * e2u(i1:i2,j1:j2) * e3u(i1:i2,j1:j2,jk,Kmm_a) & 
     585              &                                                 *  uu(i1:i2,j1:j2,jk,Kmm_a) 
    582586         END DO 
    583587      ELSE 
    … …  
    591595                     zuno = uu(ji,jj,jk,Kmm_a) * e3u(ji,jj,jk,Krhs_a) 
    592596                     zunu = tabres(ji,jj,jk,1) 
    593                      uu(ji,jj,jk,Kbb_a) = ( zub + atfp * ( zunu - zuno) ) &       
    594                                     & * umask(ji,jj,jk) / e3u(ji,jj,jk,Kbb_a) 
     597                     uu(ji,jj,jk,Kbb_a) = ( zub + atfp * ( zunu - zuno ) )      &       
     598                       &                  * umask(ji,jj,jk) / e3u(ji,jj,jk,Kbb_a) 
    595599                  ENDIF 
    596600                  ! 
    … …  
    682686            DO jj=j1,j2 
    683687               DO ji=i1,i2 
    684                   tabres(ji,jj,jk,1) = zrhox * e1v(ji,jj) * e3v(ji,jj,jk,Kmm_a) * vmask(ji,jj,jk) * vv(ji,jj,jk,Kmm_a) & 
    685                                        + (vmask(ji,jj,jk)-1)*999._wp 
     688                  tabres(ji,jj,jk,1) = zrhox * e1v(ji,jj) * e3v(ji,jj,jk,Kmm_a)                   & 
     689                    &                  *  vmask(ji,jj,jk) *  vv(ji,jj,jk,Kmm_a)                   & 
     690                    &                  + (vmask(ji,jj,jk)-1)*999._wp 
    686691                  tabres(ji,jj,jk,2) = vmask(ji,jj,jk) * zrhox * e1v(ji,jj) * e3v(ji,jj,jk,Kmm_a) & 
    687692                                       + (vmask(ji,jj,jk)-1)*999._wp 
    … …  
    708713               ENDDO 
    709714               IF (N_in * N_out > 0) THEN 
    710                   h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
     715                  h_diff = SUM( h_out(1:N_out) ) - SUM( h_in(1:N_in) ) 
    711716                  IF (h_diff < -1.e-4) then 
    712717!Even if bathy at T points match it's possible for the U points to be deeper in the child grid.  
    … …  
    714719                     excess = 0._wp 
    715720                     DO jk=N_in,1,-1 
    716                         thick = MIN(-1*h_diff, h_in(jk)) 
     721                        thick = MIN( -1._wp*h_diff, h_in(jk) ) 
    717722                        excess = excess + tabin(jk)*thick*e2u(ji,jj) 
    718                         tabin(jk) = tabin(jk)*(1. - thick/h_in(jk)) 
     723                        tabin(jk) = tabin(jk)*(1._wp - thick/h_in(jk)) 
    719724                        h_diff = h_diff + thick 
    720725                        IF ( h_diff == 0) THEN 
    … …  
    725730                     ENDDO 
    726731                  ENDIF 
    727                   CALL reconstructandremap(tabin(1:N_in),h_in(1:N_in),tabres_child(ji,jj,1:N_out),h_out(1:N_out),N_in,N_out) 
     732                  CALL reconstructandremap( tabin(1:N_in) , h_in(1:N_in), tabres_child(ji,jj,1:N_out),     & 
     733                    &                       h_out(1:N_out), N_in        , N_out                          ) 
    728734                  tabres_child(ji,jj,N_out) = tabres_child(ji,jj,N_out) + excess/(e1v(ji,jj)*h_out(N_out)) 
    729735               ENDIF 
    … …  
    736742                  ! 
    737743                  IF( .NOT.(lk_agrif_fstep.AND.(neuler==0)) ) THEN ! Add asselin part 
    738                      vv(ji,jj,jk,Kbb_a) = vv(ji,jj,jk,Kbb_a) &  
    739                            & + atfp * ( tabres_child(ji,jj,jk) - vv(ji,jj,jk,Kmm_a) ) * vmask(ji,jj,jk) 
     744                     vv(ji,jj,jk,Kbb_a) =                     vv(ji,jj,jk,Kbb_a)                  &  
     745                           &             + atfp * ( tabres_child(ji,jj,jk) - vv(ji,jj,jk,Kmm_a) ) & 
     746                           &                      *        vmask(ji,jj,jk) 
    740747                  ENDIF 
    741748                  ! 
    … …  
    781788                     zvno = vv(ji,jj,jk,Kmm_a) * e3v(ji,jj,jk,Krhs_a) 
    782789                     zvnu = tabres(ji,jj,jk,1) 
    783                      vv(ji,jj,jk,Kbb_a) = ( zvb + atfp * ( zvnu - zvno) ) &       
    784                                     & * vmask(ji,jj,jk) / e3v(ji,jj,jk,Kbb_a) 
     790                     vv(ji,jj,jk,Kbb_a) = ( zvb + atfp * ( zvnu - zvno) )        &       
     791                       &                  * vmask(ji,jj,jk) / e3v(ji,jj,jk,Kbb_a) 
    785792                  ENDIF 
    786793                  ! 
    … …  
    883890               ! Update barotropic velocities: 
    884891               IF ( .NOT.ln_dynspg_ts .OR. (ln_dynspg_ts.AND.(.NOT.ln_bt_fw)) ) THEN 
    885                   IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) THEN ! Add asselin part 
     892                  IF ( .NOT.( lk_agrif_fstep .AND. (neuler==0) ) ) THEN                          ! Add asselin part 
    886893                     zcorr = (tabres(ji,jj) - uu_b(ji,jj,Kmm_a) * hu(ji,jj,Krhs_a)) * r1_hu(ji,jj,Kbb_a) 
    887894                     uu_b(ji,jj,Kbb_a) = uu_b(ji,jj,Kbb_a) + atfp * zcorr * umask(ji,jj,1) 
    … …  
    948955               ! 
    949956               ! Update barotropic velocities: 
    950                IF ( .NOT.ln_dynspg_ts .OR. (ln_dynspg_ts.AND.(.NOT.ln_bt_fw)) ) THEN 
    951                   IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) THEN ! Add asselin part 
     957               IF ( .NOT.ln_dynspg_ts .OR. ( ln_dynspg_ts .AND. ( .NOT.ln_bt_fw ) ) ) THEN 
     958                  IF ( .NOT. ( lk_agrif_fstep .AND. ( neuler==0 ) ) ) THEN                      ! Add asselin part 
    952959                     zcorr = (tabres(ji,jj) - vv_b(ji,jj,Kmm_a) * hv(ji,jj,Krhs_a)) * r1_hv(ji,jj,Kbb_a) 
    953                      vv_b(ji,jj,Kbb_a) = vv_b(ji,jj,Kbb_a) + atfp * zcorr * vmask(ji,jj,1) 
     960                     vv_b(ji,jj,Kbb_a) =       vv_b(ji,jj,Kbb_a) + atfp * zcorr  *  vmask(ji,jj,1) 
    954961                  END IF 
    955962               ENDIF               
    … …  
    10001007            DO jj=j1,j2 
    10011008               DO ji=i1,i2 
    1002                   ssh(ji,jj,Kbb_a) =   ssh(ji,jj,Kbb_a) & 
    1003                         & + atfp * ( tabres(ji,jj) - ssh(ji,jj,Kmm_a) ) * tmask(ji,jj,1) 
     1009                  ssh(ji,jj,Kbb_a) =              ssh(ji,jj,Kbb_a)                & 
     1010                    &               + atfp * ( tabres(ji,jj) - ssh(ji,jj,Kmm_a) ) & 
     1011                    &                        *  tmask(ji,jj,1) 
    10041012               END DO 
    10051013            END DO 
    … …  
    10951103               zcor = rdt * r1_e1e2t(i1  ,jj) * e2u(i1,jj) * (ub2_b(i1,jj)-tabres(i1,jj))  
    10961104               ssh(i1  ,jj,Kmm_a) = ssh(i1  ,jj,Kmm_a) + zcor 
    1097                IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) ssh(i1  ,jj,Kbb_a) = ssh(i1  ,jj,Kbb_a) + atfp * zcor 
     1105               IF ( .NOT. ( lk_agrif_fstep .AND. ( neuler==0 ) ) )  & 
     1106                 &                  ssh(i1  ,jj,Kbb_a) = ssh(i1  ,jj,Kbb_a) + atfp * zcor 
    10981107            END DO 
    10991108         ENDIF 
    … …  
    11021111               zcor = - rdt * r1_e1e2t(i2+1,jj) * e2u(i2,jj) * (ub2_b(i2,jj)-tabres(i2,jj)) 
    11031112               ssh(i2+1,jj,Kmm_a) = ssh(i2+1,jj,Kmm_a) + zcor 
    1104                IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) ssh(i2+1,jj,Kbb_a) = ssh(i2+1,jj,Kbb_a) + atfp * zcor 
     1113               IF ( .NOT. ( lk_agrif_fstep .AND. ( neuler==0 ) ) )   & 
     1114                 &                  ssh(i2+1,jj,Kbb_a) = ssh(i2+1,jj,Kbb_a) + atfp * zcor 
    11051115            END DO 
    11061116         ENDIF 
    … …  
    11811191         IF (southern_side) THEN 
    11821192            DO ji=i1,i2 
    1183                zcor = rdt * r1_e1e2t(ji,j1  ) * e1v(ji,j1  ) * (vb2_b(ji,j1)-tabres(ji,j1)) 
     1193               zcor = rdt * r1_e1e2t(ji,j1  ) * e1v(ji,j1  ) * ( vb2_b(ji,j1)-tabres(ji,j1) ) 
    11841194               ssh(ji,j1  ,Kmm_a) = ssh(ji,j1  ,Kmm_a) + zcor 
    1185                IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) ssh(ji,j1  ,Kbb_a) = ssh(ji,j1,Kbb_a) + atfp * zcor 
     1195               IF ( .NOT. ( lk_agrif_fstep .AND. ( neuler==0 ) ) )   & 
     1196                 &                 ssh(ji,j1  ,Kbb_a) = ssh(ji,j1,Kbb_a) + atfp * zcor 
    11861197            END DO 
    11871198         ENDIF 
    11881199         IF (northern_side) THEN                
    11891200            DO ji=i1,i2 
    1190                zcor = - rdt * r1_e1e2t(ji,j2+1) * e1v(ji,j2  ) * (vb2_b(ji,j2)-tabres(ji,j2)) 
     1201               zcor = - rdt * r1_e1e2t(ji,j2+1) * e1v(ji,j2  ) * ( vb2_b(ji,j2)-tabres(ji,j2) ) 
    11911202               ssh(ji,j2+1,Kmm_a) = ssh(ji,j2+1,Kmm_a) + zcor 
    1192                IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) ssh(ji,j2+1,Kbb_a) = ssh(ji,j2+1,Kbb_a) + atfp * zcor 
     1203               IF ( .NOT. ( lk_agrif_fstep .AND. ( neuler==0 ) ) )   & 
     1204                 &                  ssh(ji,j2+1,Kbb_a) = ssh(ji,j2+1,Kbb_a) + atfp * zcor 
    11931205            END DO 
    11941206         ENDIF 
    … …  
    12231235            END DO 
    12241236         END DO 
    1225          tabres(:,:,:,1)=tabres(:,:,:,1)*Agrif_Rhox()*Agrif_Rhoy() 
    1226          tabres(:,:,:,2)=tabres(:,:,:,2)*Agrif_Rhox() 
    1227          tabres(:,:,:,3)=tabres(:,:,:,3)*Agrif_Rhoy() 
     1237         tabres(:,:,:,1) = tabres(:,:,:,1)*Agrif_Rhox()*Agrif_Rhoy() 
     1238         tabres(:,:,:,2) = tabres(:,:,:,2)*Agrif_Rhox() 
     1239         tabres(:,:,:,3) = tabres(:,:,:,3)*Agrif_Rhoy() 
    12281240      ELSE 
    12291241         DO jk=k1,k2 
    … …  
    12341246                     print *,'VAL2 = ',ji,jj,jk,tabres(ji,jj,jk,2),e1t(ji,jj)*tmask(ji,jj,jk) 
    12351247                     print *,'VAL3 = ',ji,jj,jk,tabres(ji,jj,jk,3),e2t(ji,jj)*tmask(ji,jj,jk) 
    1236                      ztemp = sqrt(tabres(ji,jj,jk,1)/(tabres(ji,jj,jk,2)*tabres(ji,jj,jk,3))) 
     1248                     ztemp = SQRT( tabres(ji,jj,jk,1)/(tabres(ji,jj,jk,2)*tabres(ji,jj,jk,3)) ) 
    12371249                     print *,'CORR = ',ztemp-1. 
    12381250                     print *,'NEW VALS = ',tabres(ji,jj,jk,2)*ztemp,tabres(ji,jj,jk,3)*ztemp, & 
    12391251                           tabres(ji,jj,jk,2)*ztemp*tabres(ji,jj,jk,3)*ztemp 
    1240                      e1t(ji,jj) = tabres(ji,jj,jk,2)*ztemp 
    1241                      e2t(ji,jj) = tabres(ji,jj,jk,3)*ztemp 
     1252                     e1t(ji,jj) = tabres(ji,jj,jk,2) * ztemp 
     1253                     e2t(ji,jj) = tabres(ji,jj,jk,3) * ztemp 
    12421254                  END IF 
    12431255               END DO 
    … …  
    13191331            DO jj=j1,j2 
    13201332               DO ji=i1,i2 
    1321                   ptab(ji,jj,jk) = e3t_0(ji,jj,jk) * (1._wp + ssh(ji,jj,Kmm_a) & 
    1322                                      & *ssmask(ji,jj)/(ht_0(ji,jj)-1._wp + ssmask(ji,jj))) 
     1333                  ptab(ji,jj,jk) =   e3t_0(ji,jj,jk)                                    & 
     1334                    &                * ( 1._wp + ssh(ji,jj,Kmm_a)                       & 
     1335                    &                          * ssmask(ji,jj)                          & 
     1336                    &                          / ( ht_0(ji,jj)-1._wp + ssmask(ji,jj) ) ) 
    13231337               END DO 
    13241338            END DO 
    … …  
    13391353               DO jj=j1,j2 
    13401354                  DO ji=i1,i2 
    1341                      e3t(ji,jj,jk,Kbb_a) =  e3t(ji,jj,jk,Kbb_a) & 
    1342                            & + atfp * ( ptab(ji,jj,jk) - e3t(ji,jj,jk,Kmm_a) ) 
     1355                     e3t(ji,jj,jk,Kbb_a) =             e3t(ji,jj,jk,Kbb_a)                  & 
     1356                       &                  + atfp * ( ptab(ji,jj,jk) - e3t(ji,jj,jk,Kmm_a) ) 
    13431357                  END DO 
    13441358               END DO 
    … …  
    13531367                  DO ji = i1,i2             
    13541368                     zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk)) 
    1355                      e3w(ji,jj,jk,Kbb_a)  = e3w_0(ji,jj,jk) + ( 1.0_wp - 0.5_wp * tmask(ji,jj,jk) ) *        &  
    1356                      &                                        ( e3t(ji,jj,jk-1,Kbb_a) - e3t_0(ji,jj,jk-1) )  & 
    1357                      &                                  +            0.5_wp * tmask(ji,jj,jk)   *        & 
    1358                      &                                        ( e3t(ji,jj,jk  ,Kbb_a) - e3t_0(ji,jj,jk  ) ) 
    1359                      gdepw(ji,jj,jk,Kbb_a) = gdepw(ji,jj,jk-1,Kbb_a) + e3t(ji,jj,jk-1,Kbb_a) 
    1360                      gdept(ji,jj,jk,Kbb_a) =      zcoef  * ( gdepw(ji,jj,jk  ,Kbb_a) + 0.5 * e3w(ji,jj,jk,Kbb_a))  & 
    1361                          &               + (1-zcoef) * ( gdept(ji,jj,jk-1,Kbb_a) +       e3w(ji,jj,jk,Kbb_a))  
     1369                     e3w(ji,jj,jk,Kbb_a)  = e3w_0(ji,jj,jk) + ( 1.0_wp - 0.5_wp * tmask(ji,jj,jk) )          &  
     1370                       &                                    * ( e3t(ji,jj,jk-1,Kbb_a) - e3t_0(ji,jj,jk-1) )  & 
     1371                       &                                    + 0.5_wp * tmask(ji,jj,jk)                       & 
     1372                       &                                    * ( e3t(ji,jj,jk  ,Kbb_a) - e3t_0(ji,jj,jk  ) ) 
     1373                     gdepw(ji,jj,jk,Kbb_a) =                       gdepw(ji,jj,jk-1,Kbb_a) +      e3t(ji,jj,jk-1,Kbb_a) 
     1374                     gdept(ji,jj,jk,Kbb_a) =            zcoef  * ( gdepw(ji,jj,jk  ,Kbb_a) + 0.5 * e3w(ji,jj,jk  ,Kbb_a) )  & 
     1375                       &                     + (1._wp - zcoef) * ( gdept(ji,jj,jk-1,Kbb_a) +       e3w(ji,jj,jk  ,Kbb_a) )  
    13621376                  END DO 
    13631377               END DO 
    … …  
    13791393         ! 
    13801394         ! Update vertical scale factor at W-points and depths: 
    1381          e3w (i1:i2,j1:j2,1,Kmm_a) = e3w_0(i1:i2,j1:j2,1) + e3t(i1:i2,j1:j2,1,Kmm_a) - e3t_0(i1:i2,j1:j2,1) 
     1395           e3w(i1:i2,j1:j2,1,Kmm_a) = e3w_0(i1:i2,j1:j2,1) + e3t(i1:i2,j1:j2,1,Kmm_a) - e3t_0(i1:i2,j1:j2,1) 
    13821396         gdept(i1:i2,j1:j2,1,Kmm_a) = 0.5_wp * e3w(i1:i2,j1:j2,1,Kmm_a) 
    13831397         gdepw(i1:i2,j1:j2,1,Kmm_a) = 0.0_wp 
    1384          gde3w(i1:i2,j1:j2,1) = gdept(i1:i2,j1:j2,1,Kmm_a) - (ht(i1:i2,j1:j2)-ht_0(i1:i2,j1:j2)) ! Last term in the rhs is ssh 
     1398         gde3w(i1:i2,j1:j2,1      ) = gdept(i1:i2,j1:j2,1,Kmm_a) - ( ht(i1:i2,j1:j2) - ht_0(i1:i2,j1:j2) ) ! Last term in the rhs is ssh 
    13851399         ! 
    13861400         DO jk = 2, jpk 
    … …  
    13881402               DO ji = i1,i2             
    13891403               zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk)) 
    1390                e3w(ji,jj,jk,Kmm_a)  = e3w_0(ji,jj,jk) + ( 1.0_wp - 0.5_wp * tmask(ji,jj,jk) ) * ( e3t(ji,jj,jk-1,Kmm_a) - e3t_0(ji,jj,jk-1) )   & 
    1391                &                                  +            0.5_wp * tmask(ji,jj,jk)   * ( e3t(ji,jj,jk  ,Kmm_a) - e3t_0(ji,jj,jk  ) ) 
    1392                gdepw(ji,jj,jk,Kmm_a) = gdepw(ji,jj,jk-1,Kmm_a) + e3t(ji,jj,jk-1,Kmm_a) 
    1393                gdept(ji,jj,jk,Kmm_a) =      zcoef  * ( gdepw(ji,jj,jk  ,Kmm_a) + 0.5 * e3w(ji,jj,jk,Kmm_a))  & 
    1394                    &               + (1-zcoef) * ( gdept(ji,jj,jk-1,Kmm_a) +       e3w(ji,jj,jk,Kmm_a))  
    1395                gde3w(ji,jj,jk) = gdept(ji,jj,jk,Kmm_a) - (ht(ji,jj)-ht_0(ji,jj)) ! Last term in the rhs is ssh 
     1404                e3w(ji,jj,jk,Kmm_a)  = e3w_0(ji,jj,jk) + ( 1.0_wp - 0.5_wp * tmask(ji,jj,jk) )          & 
     1405                 &                                     * ( e3t(ji,jj,jk-1,Kmm_a) - e3t_0(ji,jj,jk-1) )  & 
     1406                 &                                     + 0.5_wp * tmask(ji,jj,jk)                       & 
     1407                 &                                     * ( e3t(ji,jj,jk  ,Kmm_a) - e3t_0(ji,jj,jk  ) ) 
     1408               gdepw(ji,jj,jk,Kmm_a) =                         gdepw(ji,jj,jk-1,Kmm_a) +       e3t(ji,jj,jk-1,Kmm_a) 
     1409               gdept(ji,jj,jk,Kmm_a) =             zcoef   * ( gdepw(ji,jj,jk  ,Kmm_a) + 0.5 * e3w(ji,jj,jk  ,Kmm_a) )  & 
     1410                 &                     + ( 1._wp - zcoef ) * ( gdept(ji,jj,jk-1,Kmm_a) +       e3w(ji,jj,jk  ,Kmm_a) )  
     1411               gde3w(ji,jj,jk      ) =                         gdept(ji,jj,jk  ,Kmm_a) - ( ht(ji,jj)-ht_0(ji,jj) )    ! Last term in the rhs is ssh 
    13961412               END DO 
    13971413            END DO 
    … …  
    13991415         ! 
    14001416         IF  ((neuler==0).AND.(Agrif_Nb_Step()==0) ) THEN 
    1401             e3t (i1:i2,j1:j2,1:jpk,Kbb_a)  = e3t (i1:i2,j1:j2,1:jpk,Kmm_a) 
    1402             e3w (i1:i2,j1:j2,1:jpk,Kbb_a)  = e3w (i1:i2,j1:j2,1:jpk,Kmm_a) 
     1417              e3t(i1:i2,j1:j2,1:jpk,Kbb_a) =   e3t(i1:i2,j1:j2,1:jpk,Kmm_a) 
     1418              e3w(i1:i2,j1:j2,1:jpk,Kbb_a) =   e3w(i1:i2,j1:j2,1:jpk,Kmm_a) 
    14031419            gdepw(i1:i2,j1:j2,1:jpk,Kbb_a) = gdepw(i1:i2,j1:j2,1:jpk,Kmm_a) 
    14041420            gdept(i1:i2,j1:j2,1:jpk,Kbb_a) = gdept(i1:i2,j1:j2,1:jpk,Kmm_a) 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_top_interp.F90

    r11053 r11463  
    151151               ibdy = nlci-nbghostcells 
    152152               DO jn = 1, jptra 
    153                   tr(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn) + z2 * ptab_child(ibdy,jmin:jmax,1:jpkm1,jn) 
     153                  tr(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn)  & 
     154                    &                                     + z2 * ptab_child(ibdy  ,jmin:jmax,1:jpkm1,jn) 
    154155                  DO jk = 1, jpkm1 
    155156                     DO jj = jmin,jmax 
    … …  
    157158                           tr(ibdy,jj,jk,jn,Krhs_a) = tr(ibdy+1,jj,jk,jn,Krhs_a) * tmask(ibdy,jj,jk) 
    158159                        ELSE 
    159                            tr(ibdy,jj,jk,jn,Krhs_a)=(z4*tr(ibdy+1,jj,jk,jn,Krhs_a)+z3*tr(ibdy-1,jj,jk,jn,Krhs_a))*tmask(ibdy,jj,jk) 
     160                           tr(ibdy,jj,jk,jn,Krhs_a) = (  z4 * tr(ibdy+1,jj,jk,jn,Krhs_a)     & 
     161                             &                         + z3 * tr(ibdy-1,jj,jk,jn,Krhs_a)     & 
     162                             &                        ) *  tmask(ibdy  ,jj,jk) 
    160163                           IF( uu(ibdy-1,jj,jk,Kmm_a) > 0._wp ) THEN 
    161                               tr(ibdy,jj,jk,jn,Krhs_a)=( z6*tr(ibdy-1,jj,jk,jn,Krhs_a)+z5*tr(ibdy+1,jj,jk,jn,Krhs_a) &  
    162                                                  + z7*tr(ibdy-2,jj,jk,jn,Krhs_a) ) * tmask(ibdy,jj,jk) 
    163                            ENDIF 
    164                         ENDIF 
    165                      END DO 
    166                   END DO 
    167                   ! Restore ghost points: 
    168                   tr(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn) * tmask(ibdy+1,jmin:jmax,1:jpkm1) 
     164                              tr(ibdy,jj,jk,jn,Krhs_a) = (  z6 * tr(ibdy-1,jj,jk,jn,Krhs_a)  & 
     165                                &                         + z5 * tr(ibdy+1,jj,jk,jn,Krhs_a)  &  
     166                                &                         + z7 * tr(ibdy-2,jj,jk,jn,Krhs_a)  & 
     167                                &                         ) * tmask(ibdy  ,jj,jk) 
     168                           ENDIF 
     169                        ENDIF 
     170                     END DO 
     171                  END DO 
     172                  ! Restore ghost points: 
     173                  tr(ibdy+1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy+1,jmin:jmax,1:jpkm1,jn) & 
     174                    &                                     *     tmask(ibdy+1,jmin:jmax,1:jpkm1) 
    169175               END DO 
    170176            ENDIF 
    … …  
    180186               jbdy = nlcj-nbghostcells          
    181187               DO jn = 1, jptra 
    182                   tr(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(imin:imax,jbdy+1,1:jpkm1,jn) + z2 * ptab_child(imin:imax,jbdy,1:jpkm1,jn) 
     188                  tr(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(imin:imax,jbdy+1,1:jpkm1,jn) & 
     189                    &                                     + z2 * ptab_child(imin:imax,jbdy  ,1:jpkm1,jn) 
    183190                  DO jk = 1, jpkm1 
    184191                     DO ji = imin,imax 
    … …  
    186193                           tr(ji,jbdy,jk,jn,Krhs_a) = tr(ji,jbdy+1,jk,jn,Krhs_a) * tmask(ji,jbdy,jk) 
    187194                        ELSE 
    188                            tr(ji,jbdy,jk,jn,Krhs_a)=(z4*tr(ji,jbdy+1,jk,jn,Krhs_a)+z3*tr(ji,jbdy-1,jk,jn,Krhs_a))*tmask(ji,jbdy,jk)         
     195                           tr(ji,jbdy,jk,jn,Krhs_a) = ( z4 * tr(ji,jbdy+1,jk,jn,Krhs_a)  
     196                             &                        + z3 * tr(ji,jbdy-1,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk)         
    189197                           IF (vv(ji,jbdy-1,jk,Kmm_a) > 0._wp ) THEN 
    190                               tr(ji,jbdy,jk,jn,Krhs_a)=( z6*tr(ji,jbdy-1,jk,jn,Krhs_a)+z5*tr(ji,jbdy+1,jk,jn,Krhs_a)  & 
    191                                                  + z7*tr(ji,jbdy-2,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk) 
    192                            ENDIF 
    193                         ENDIF 
    194                      END DO 
    195                   END DO 
    196                   ! Restore ghost points: 
    197                   tr(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy+1,1:jpkm1,jn) * tmask(imin:imax,jbdy+1,1:jpkm1) 
     198                              tr(ji,jbdy,jk,jn,Krhs_a) = (  z6 * tr(ji,jbdy-1,jk,jn,Krhs_a)                       & 
     199                                &                         + z5 * tr(ji,jbdy+1,jk,jn,Krhs_a)                       & 
     200                                &                         + z7 * tr(ji,jbdy-2,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk) 
     201                           ENDIF 
     202                        ENDIF 
     203                     END DO 
     204                  END DO 
     205                  ! Restore ghost points: 
     206                  tr(imin:imax,jbdy+1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy+1,1:jpkm1,jn)  & 
     207                    &                                     *     tmask(imin:imax,jbdy+1,1:jpkm1) 
    198208               END DO 
    199209            ENDIF 
    … …  
    209219               ibdy = 1+nbghostcells        
    210220               DO jn = 1, jptra 
    211                   tr(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn) + z2 * ptab_child(ibdy,jmin:jmax,1:jpkm1,jn) 
     221                  tr(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn)         & 
     222                    &                                     + z2 * ptab_child(ibdy  ,jmin:jmax,1:jpkm1,jn) 
    212223                  DO jk = 1, jpkm1 
    213224                     DO jj = jmin,jmax 
    … …  
    215226                           tr(ibdy,jj,jk,jn,Krhs_a) = tr(ibdy-1,jj,jk,jn,Krhs_a) * tmask(ibdy,jj,jk) 
    216227                        ELSE 
    217                            tr(ibdy,jj,jk,jn,Krhs_a)=(z4*tr(ibdy-1,jj,jk,jn,Krhs_a)+z3*tr(ibdy+1,jj,jk,jn,Krhs_a))*tmask(ibdy,jj,jk)         
     228                           tr(ibdy,jj,jk,jn,Krhs_a) = (  z4 * tr(ibdy-1,jj,jk,jn,Krhs_a)                          & 
     229                             &                         + z3 * tr(ibdy+1,jj,jk,jn,Krhs_a) ) * tmask(ibdy,jj,jk)         
    218230                           IF( uu(ibdy,jj,jk,Kmm_a) < 0._wp ) THEN 
    219                               tr(ibdy,jj,jk,jn,Krhs_a)=( z6*tr(ibdy+1,jj,jk,jn,Krhs_a)+z5*tr(ibdy-1,jj,jk,jn,Krhs_a) & 
    220                                                  + z7*tr(ibdy+2,jj,jk,jn,Krhs_a) ) * tmask(ibdy,jj,jk) 
    221                            ENDIF 
    222                         ENDIF 
    223                      END DO 
    224                   END DO 
    225                   ! Restore ghost points: 
    226                   tr(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn) * tmask(ibdy-1,jmin:jmax,1:jpkm1) 
     231                              tr(ibdy,jj,jk,jn,Krhs_a) = (  z6 * tr(ibdy+1,jj,jk,jn,Krhs_a)                       & 
     232                                &                         + z5 * tr(ibdy-1,jj,jk,jn,Krhs_a)                       & 
     233                                &                         + z7 * tr(ibdy+2,jj,jk,jn,Krhs_a) ) * tmask(ibdy,jj,jk) 
     234                           ENDIF 
     235                        ENDIF 
     236                     END DO 
     237                  END DO 
     238                  ! Restore ghost points: 
     239                  tr(ibdy-1,jmin:jmax,1:jpkm1,jn,Krhs_a) = ptab_child(ibdy-1,jmin:jmax,1:jpkm1,jn)                & 
     240                    &                                     *     tmask(ibdy-1,jmin:jmax,1:jpkm1) 
    227241               END DO 
    228242            ENDIF 
    … …  
    238252               jbdy=1+nbghostcells         
    239253               DO jn = 1, jptra 
    240                   tr(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) = z1 * ptab_child(imin:imax,jbdy-1,1:jpkm1,jn) + z2 * ptab_child(imin:imax,jbdy,1:jpkm1,jn) 
     254                  tr(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) =  z1 * ptab_child(imin:imax,jbdy-1,1:jpkm1,jn)          & 
     255                    &                                     + z2 * ptab_child(imin:imax,jbdy  ,1:jpkm1,jn) 
    241256                  DO jk = 1, jpkm1       
    242257                     DO ji = imin,imax 
    243258                        IF( vmask(ji,jbdy,jk) == 0._wp ) THEN 
    244                            tr(ji,jbdy,jk,jn,Krhs_a)=tr(ji,jbdy-1,jk,jn,Krhs_a) * tmask(ji,jbdy,jk) 
    245                         ELSE 
    246                            tr(ji,jbdy,jk,jn,Krhs_a)=(z4*tr(ji,jbdy-1,jk,jn,Krhs_a)+z3*tr(ji,jbdy+1,jk,jn,Krhs_a))*tmask(ji,jbdy,jk) 
     259                           tr(ji,jbdy,jk,jn,Krhs_a) = tr(ji,jbdy-1,jk,jn,Krhs_a) * tmask(ji,jbdy,jk) 
     260                        ELSE 
     261                           tr(ji,jbdy,jk,jn,Krhs_a) = (  z4 * tr(ji,jbdy-1,jk,jn,Krhs_a)                          & 
     262                             &                         + z3 * tr(ji,jbdy+1,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk) 
    247263                           IF( vv(ji,jbdy,jk,Kmm_a) < 0._wp ) THEN 
    248                               tr(ji,jbdy,jk,jn,Krhs_a)=( z6*tr(ji,jbdy+1,jk,jn,Krhs_a)+z5*tr(ji,jbdy-1,jk,jn,Krhs_a) &  
    249                                                  + z7*tr(ji,jbdy+2,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk) 
    250                            ENDIF 
    251                         ENDIF 
    252                      END DO 
    253                   END DO 
    254                   ! Restore ghost points: 
    255                   tr(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy-1,1:jpkm1,jn) * tmask(imin:imax,jbdy-1,1:jpkm1) 
     264                              tr(ji,jbdy,jk,jn,Krhs_a) = (  z6 * tr(ji,jbdy+1,jk,jn,Krhs_a)                       & 
     265                                &                         + z5 * tr(ji,jbdy-1,jk,jn,Krhs_a)                       &  
     266                                &                         + z7 * tr(ji,jbdy+2,jk,jn,Krhs_a) ) * tmask(ji,jbdy,jk) 
     267                           ENDIF 
     268                        ENDIF 
     269                     END DO 
     270                  END DO 
     271                  ! Restore ghost points: 
     272                  tr(imin:imax,jbdy-1,1:jpkm1,jn,Krhs_a) = ptab_child(imin:imax,jbdy-1,1:jpkm1,jn)                & 
     273                    &                                     *     tmask(imin:imax,jbdy-1,1:jpkm1) 
    256274               END DO 
    257275            ENDIF 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_top_sponge.F90

    r10989 r11463  
    108108                  N_in = N_in + 1 
    109109                  tabin(jk,:) = tabres(ji,jj,jk,n1:n2-1) 
    110                   h_in(N_in) = tabres(ji,jj,jk,n2) 
     110                  h_in(N_in)  = tabres(ji,jj,jk,n2) 
    111111               END DO 
    112112               N_out = 0 
    113113               DO jk=1,jpk ! jpk of child grid 
    114114                  IF (tmask(ji,jj,jk) == 0) EXIT  
    115                   N_out = N_out + 1 
     115                  N_out     = N_out + 1 
    116116                  h_out(jk) = e3t(ji,jj,jk,Kmm) !Child grid scale factors. Could multiply by e1e2t here instead of division above 
    117117               ENDDO 
    118118               IF (N_in > 0) THEN 
    119                   h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
     119                  h_diff = SUM( h_out(1:N_out) ) - SUM( h_in(1:N_in) ) 
    120120                  tabres(ji,jj,k2,:) = tabres(ji,jj,k2-1,:) !what is this line for????? 
    121121                  DO jn=1,jptra 
    … …  
    133133                  trbdiff(ji,jj,jk,1:jptra) = tr(ji,jj,jk,1:jptra,Kbb) - tabres_child(ji,jj,jk,1:jptra) 
    134134# else 
    135                   trbdiff(ji,jj,jk,1:jptra) = tr(ji,jj,jk,1:jptra,Kbb) - tabres(ji,jj,jk,1:jptra) 
     135                  trbdiff(ji,jj,jk,1:jptra) = tr(ji,jj,jk,1:jptra,Kbb) -       tabres(ji,jj,jk,1:jptra) 
    136136# endif 
    137137               ENDDO 
    … …  
    153153                  DO ji = i1+1,i2-1 
    154154                     IF( .NOT. tabspongedone_trn(ji,jj) ) THEN  
    155                         tr(ji,jj,jk,jn,Krhs) = tr(ji,jj,jk,jn,Krhs) + (  ztu(ji,jj) - ztu(ji-1,jj  )     & 
    156                            &                                   + ztv(ji,jj) - ztv(ji  ,jj-1)  )  & 
    157                            &                                * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 
     155                        tr(ji,jj,jk,jn,Krhs) = tr(ji,jj,jk,jn,Krhs) + (  ztu(ji,jj) - ztu(ji-1,jj  )       & 
     156                           &                                           + ztv(ji,jj) - ztv(ji  ,jj-1)  )    & 
     157                           &                                          * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 
    158158                     ENDIF 
    159159                  END DO 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_top_update.F90

    r10989 r11463  
    8484               DO jj=j1,j2 
    8585                  DO ji=i1,i2 
    86                      tabres(ji,jj,jk,jn) = (tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) ) & 
    87                                            * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 
     86                     tabres(ji,jj,jk,jn) = (    tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) ) & 
     87                       &                   * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 
    8888                  END DO 
    8989               END DO 
    … …  
    9393            DO jj=j1,j2 
    9494               DO ji=i1,i2 
    95                   tabres(ji,jj,jk,n2) = tmask(ji,jj,jk) * e3t(ji,jj,jk,Kmm) & 
    96                                            + (tmask(ji,jj,jk)-1)*999._wp 
     95                  tabres(ji,jj,jk,n2) =   tmask(ji,jj,jk) * e3t(ji,jj,jk,Kmm) & 
     96                    &                  + (tmask(ji,jj,jk)-1)*999._wp 
    9797               END DO 
    9898            END DO 
    … …  
    138138                     DO ji=i1,i2 
    139139                        IF( tabres_child(ji,jj,jk,jn) .NE. 0. ) THEN 
    140                            tr(ji,jj,jk,jn,Kbb) = ts(ji,jj,jk,jn,Kbb) &  
    141                                  & + atfp * ( tabres_child(ji,jj,jk,jn) & 
    142                                  &          - tr(ji,jj,jk,jn,Kmm) ) * tmask(ji,jj,jk) 
     140                           tr(ji,jj,jk,jn,Kbb) = ts(ji,jj,jk,jn,Kbb)                     &  
     141                             &                  + atfp * ( tabres_child(ji,jj,jk,jn)     & 
     142                             &                                     - tr(ji,jj,jk,jn,Kmm) & 
     143                             &                           ) *      tmask(ji,jj,jk) 
    143144                        ENDIF 
    144145                     ENDDO 
    … …  
    183184                  DO ji=i1,i2 
    184185!> jc tmp 
    185                      tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm)  * e3t(ji,jj,jk,Kmm) / e3t_0(ji,jj,jk) 
    186 !                     tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) 
     186                     tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) / e3t_0(ji,jj,jk) 
     187!                    tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) 
    187188!< jc tmp 
    188189                  END DO 
    … …  
    194195         DO jn = n1,n2 
    195196            tabres(i1:i2,j1:j2,k1:k2,jn) =  tabres(i1:i2,j1:j2,k1:k2,jn) * e3t_0(i1:i2,j1:j2,k1:k2) & 
    196                                          & * tmask(i1:i2,j1:j2,k1:k2) 
     197              &                            * tmask(i1:i2,j1:j2,k1:k2) 
    197198         ENDDO 
    198199!< jc tmp 
    … …  
    204205                     DO ji=i1,i2 
    205206                        IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN 
    206                            ztb  = tr(ji,jj,jk,jn,Kbb) * e3t(ji,jj,jk,Kbb) ! fse3t_b prior update should be used 
     207                           ztb  =     tr(ji,jj,jk,jn,Kbb) * e3t(ji,jj,jk,Kbb) ! fse3t_b prior update should be used 
    207208                           ztnu = tabres(ji,jj,jk,jn) 
    208                            ztno = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Krhs) 
    209                            tr(ji,jj,jk,jn,Kbb) = ( ztb + atfp * ( ztnu - ztno) )  &  
    210                                      &        * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb) 
     209                           ztno =     tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Krhs) 
     210                           tr(ji,jj,jk,jn,Kbb) = ( ztb + atfp * ( ztnu - ztno) ) * tmask(ji,jj,jk)     & 
     211                             &                                                   /  e3t(ji,jj,jk,Kbb) 
    211212                        ENDIF 
    212213                     ENDDO 
Note: See TracChangeset for help on using the changeset viewer.