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

Changeset 9043


Ignore:
Timestamp:
2017-12-14T12:49:18+01:00 (6 years ago)
Author:
acc
Message:

Branch 2017/dev_merge_2017. revert renamed constants in dynspg_ts.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r9037 r9043  
    8787   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::  ftsw, ftse          ! (only used with een vorticity scheme) 
    8888 
    89    REAL(wp) ::   z1_12 = 1._wp / 12._wp   ! local ratios 
    90    REAL(wp) ::   z1_8  = 0.125_wp         ! 
    91    REAL(wp) ::   z1_4  = 0.25_wp          ! 
    92    REAL(wp) ::   z1_2  = 0.5_wp           ! 
     89   REAL(wp) ::   r1_12 = 1._wp / 12._wp   ! local ratios 
     90   REAL(wp) ::   r1_8  = 0.125_wp         ! 
     91   REAL(wp) ::   r1_4  = 0.25_wp          ! 
     92   REAL(wp) ::   r1_2  = 0.5_wp           ! 
    9393 
    9494   !! * Substitutions 
     
    157157      INTEGER  ::   ikbu, iktu, noffset   ! local integers 
    158158      INTEGER  ::   ikbv, iktv            !   -      - 
    159       REAL(wp) ::   z1_2dt_b, z2dt_bf        ! local scalars 
     159      REAL(wp) ::   r1_2dt_b, z2dt_bf        ! local scalars 
    160160      REAL(wp) ::   zx1, zx2, zu_spg, zhura  !   -      - 
    161161      REAL(wp) ::   zy1, zy2, zv_spg, zhvra  !   -      - 
     
    193193      ELSE                                        ;   z2dt_bf = 2.0_wp * rdt 
    194194      ENDIF 
    195       z1_2dt_b = 1.0_wp / z2dt_bf  
     195      r1_2dt_b = 1.0_wp / z2dt_bf  
    196196      ! 
    197197      ll_init     = ln_bt_av                    ! if no time averaging, then no specific restart  
     
    390390               zx2 = ( zwx(ji  ,jj) + zwx(ji  ,jj+1) ) * r1_e2v(ji,jj) 
    391391               ! energy conserving formulation for planetary vorticity term 
    392                zu_trd(ji,jj) =   z1_4 * ( zwz(ji  ,jj-1) * zy1 + zwz(ji,jj) * zy2 ) 
    393                zv_trd(ji,jj) = - z1_4 * ( zwz(ji-1,jj  ) * zx1 + zwz(ji,jj) * zx2 ) 
     392               zu_trd(ji,jj) =   r1_4 * ( zwz(ji  ,jj-1) * zy1 + zwz(ji,jj) * zy2 ) 
     393               zv_trd(ji,jj) = - r1_4 * ( zwz(ji-1,jj  ) * zx1 + zwz(ji,jj) * zx2 ) 
    394394            END DO 
    395395         END DO 
     
    398398         DO jj = 2, jpjm1 
    399399            DO ji = fs_2, fs_jpim1   ! vector opt. 
    400                zy1 =   z1_8 * ( zwy(ji  ,jj-1) + zwy(ji+1,jj-1) & 
     400               zy1 =   r1_8 * ( zwy(ji  ,jj-1) + zwy(ji+1,jj-1) & 
    401401                 &            + zwy(ji  ,jj  ) + zwy(ji+1,jj  ) ) * r1_e1u(ji,jj) 
    402                zx1 = - z1_8 * ( zwx(ji-1,jj  ) + zwx(ji-1,jj+1) & 
     402               zx1 = - r1_8 * ( zwx(ji-1,jj  ) + zwx(ji-1,jj+1) & 
    403403                 &            + zwx(ji  ,jj  ) + zwx(ji  ,jj+1) ) * r1_e2v(ji,jj) 
    404404               zu_trd(ji,jj)  = zy1 * ( zwz(ji  ,jj-1) + zwz(ji,jj) ) 
     
    410410         DO jj = 2, jpjm1 
    411411            DO ji = fs_2, fs_jpim1   ! vector opt. 
    412                zu_trd(ji,jj) = + z1_12 * r1_e1u(ji,jj) * (  ftne(ji,jj  ) * zwy(ji  ,jj  ) & 
     412               zu_trd(ji,jj) = + r1_12 * r1_e1u(ji,jj) * (  ftne(ji,jj  ) * zwy(ji  ,jj  ) & 
    413413                &                                         + ftnw(ji+1,jj) * zwy(ji+1,jj  ) & 
    414414                &                                         + ftse(ji,jj  ) * zwy(ji  ,jj-1) & 
    415415                &                                         + ftsw(ji+1,jj) * zwy(ji+1,jj-1) ) 
    416                zv_trd(ji,jj) = - z1_12 * r1_e2v(ji,jj) * (  ftsw(ji,jj+1) * zwx(ji-1,jj+1) & 
     416               zv_trd(ji,jj) = - r1_12 * r1_e2v(ji,jj) * (  ftsw(ji,jj+1) * zwx(ji-1,jj+1) & 
    417417                &                                         + ftse(ji,jj+1) * zwx(ji  ,jj+1) & 
    418418                &                                         + ftnw(ji,jj  ) * zwx(ji-1,jj  ) & 
     
    585585         END DO 
    586586      ELSE 
    587          zztmp = r1_rau0 * z1_2 
     587         zztmp = r1_rau0 * r1_2 
    588588         DO jj = 2, jpjm1 
    589589            DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    605605            END DO 
    606606         ELSE 
    607             zztmp = grav * z1_2 
     607            zztmp = grav * r1_2 
    608608            DO jj = 2, jpjm1               
    609609               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    624624         zssh_frc(:,:) = r1_rau0 * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) 
    625625      ELSE 
    626          zztmp = r1_rau0 * z1_2 
     626         zztmp = r1_rau0 * r1_2 
    627627         zssh_frc(:,:) = zztmp * (  emp(:,:) + emp_b(:,:) - rnf(:,:) - rnf_b(:,:)   & 
    628628                &                 + fwfisf(:,:) + fwfisf_b(:,:)                     ) 
     
    762762            DO jj = 2, jpjm1                                    ! Sea Surface Height at u- & v-points 
    763763               DO ji = 2, fs_jpim1   ! Vector opt. 
    764                   zwx(ji,jj) = z1_2 * ssumask(ji,jj)  * r1_e1e2u(ji,jj)     & 
     764                  zwx(ji,jj) = r1_2 * ssumask(ji,jj)  * r1_e1e2u(ji,jj)     & 
    765765                     &              * ( e1e2t(ji  ,jj) * zsshp2_e(ji  ,jj)  & 
    766766                     &              +   e1e2t(ji+1,jj) * zsshp2_e(ji+1,jj) ) 
    767                   zwy(ji,jj) = z1_2 * ssvmask(ji,jj)  * r1_e1e2v(ji,jj)     & 
     767                  zwy(ji,jj) = r1_2 * ssvmask(ji,jj)  * r1_e1e2v(ji,jj)     & 
    768768                     &              * ( e1e2t(ji,jj  ) * zsshp2_e(ji,jj  )  & 
    769769                     &              +   e1e2t(ji,jj+1) * zsshp2_e(ji,jj+1) ) 
     
    873873            DO jj = 2, jpjm1 
    874874               DO ji = 2, jpim1      ! NO Vector Opt. 
    875                   zsshu_a(ji,jj) = z1_2 * ssumask(ji,jj) * r1_e1e2u(ji,jj)    & 
     875                  zsshu_a(ji,jj) = r1_2 * ssumask(ji,jj) * r1_e1e2u(ji,jj)    & 
    876876                     &              * ( e1e2t(ji  ,jj  )  * ssha_e(ji  ,jj  ) & 
    877877                     &              +   e1e2t(ji+1,jj  )  * ssha_e(ji+1,jj  ) ) 
    878                   zsshv_a(ji,jj) = z1_2 * ssvmask(ji,jj) * r1_e1e2v(ji,jj)    & 
     878                  zsshv_a(ji,jj) = r1_2 * ssvmask(ji,jj) * r1_e1e2v(ji,jj)    & 
    879879                     &              * ( e1e2t(ji  ,jj  )  * ssha_e(ji  ,jj  ) & 
    880880                     &              +   e1e2t(ji  ,jj+1)  * ssha_e(ji  ,jj+1) ) 
     
    961961            DO jj = 2, jpjm1                             
    962962               DO ji = 2, jpim1 
    963                   zx1 = z1_2 * ssumask(ji  ,jj) *  r1_e1e2u(ji  ,jj)    & 
     963                  zx1 = r1_2 * ssumask(ji  ,jj) *  r1_e1e2u(ji  ,jj)    & 
    964964                     &      * ( e1e2t(ji  ,jj  ) * zsshp2_e(ji  ,jj)    & 
    965965                     &      +   e1e2t(ji+1,jj  ) * zsshp2_e(ji+1,jj  ) ) 
    966                   zy1 = z1_2 * ssvmask(ji  ,jj) *  r1_e1e2v(ji  ,jj  )  & 
     966                  zy1 = r1_2 * ssvmask(ji  ,jj) *  r1_e1e2v(ji  ,jj  )  & 
    967967                     &       * ( e1e2t(ji ,jj  ) * zsshp2_e(ji  ,jj  )  & 
    968968                     &       +   e1e2t(ji ,jj+1) * zsshp2_e(ji  ,jj+1) ) 
     
    988988                  zx1 = ( zwx(ji-1,jj  ) + zwx(ji-1,jj+1) ) * r1_e2v(ji,jj) 
    989989                  zx2 = ( zwx(ji  ,jj  ) + zwx(ji  ,jj+1) ) * r1_e2v(ji,jj) 
    990                   zu_trd(ji,jj) = z1_4 * ( zwz(ji  ,jj-1) * zy1 + zwz(ji,jj) * zy2 ) 
    991                   zv_trd(ji,jj) =-z1_4 * ( zwz(ji-1,jj  ) * zx1 + zwz(ji,jj) * zx2 ) 
     990                  zu_trd(ji,jj) = r1_4 * ( zwz(ji  ,jj-1) * zy1 + zwz(ji,jj) * zy2 ) 
     991                  zv_trd(ji,jj) =-r1_4 * ( zwz(ji-1,jj  ) * zx1 + zwz(ji,jj) * zx2 ) 
    992992               END DO 
    993993            END DO 
     
    996996            DO jj = 2, jpjm1 
    997997               DO ji = fs_2, fs_jpim1   ! vector opt. 
    998                   zy1 =   z1_8 * ( zwy(ji  ,jj-1) + zwy(ji+1,jj-1) & 
     998                  zy1 =   r1_8 * ( zwy(ji  ,jj-1) + zwy(ji+1,jj-1) & 
    999999                   &             + zwy(ji  ,jj  ) + zwy(ji+1,jj  ) ) * r1_e1u(ji,jj) 
    1000                   zx1 = - z1_8 * ( zwx(ji-1,jj  ) + zwx(ji-1,jj+1) & 
     1000                  zx1 = - r1_8 * ( zwx(ji-1,jj  ) + zwx(ji-1,jj+1) & 
    10011001                   &             + zwx(ji  ,jj  ) + zwx(ji  ,jj+1) ) * r1_e2v(ji,jj) 
    10021002                  zu_trd(ji,jj)  = zy1 * ( zwz(ji  ,jj-1) + zwz(ji,jj) ) 
     
    10081008            DO jj = 2, jpjm1 
    10091009               DO ji = fs_2, fs_jpim1   ! vector opt. 
    1010                   zu_trd(ji,jj) = + z1_12 * r1_e1u(ji,jj) * (  ftne(ji,jj  ) * zwy(ji  ,jj  ) & 
     1010                  zu_trd(ji,jj) = + r1_12 * r1_e1u(ji,jj) * (  ftne(ji,jj  ) * zwy(ji  ,jj  ) & 
    10111011                     &                                       + ftnw(ji+1,jj) * zwy(ji+1,jj  ) & 
    10121012                     &                                       + ftse(ji,jj  ) * zwy(ji  ,jj-1) &  
    10131013                     &                                       + ftsw(ji+1,jj) * zwy(ji+1,jj-1) ) 
    1014                   zv_trd(ji,jj) = - z1_12 * r1_e2v(ji,jj) * (  ftsw(ji,jj+1) * zwx(ji-1,jj+1) &  
     1014                  zv_trd(ji,jj) = - r1_12 * r1_e2v(ji,jj) * (  ftsw(ji,jj+1) * zwx(ji-1,jj+1) &  
    10151015                     &                                       + ftse(ji,jj+1) * zwx(ji  ,jj+1) & 
    10161016                     &                                       + ftnw(ji,jj  ) * zwx(ji-1,jj  ) &  
     
    11741174         zwy(:,:) = vn_adv(:,:) 
    11751175         IF( .NOT.( kt == nit000 .AND. neuler==0 ) ) THEN 
    1176             un_adv(:,:) = z1_2 * ( ub2_b(:,:) + zwx(:,:) - atfp * un_bf(:,:) ) 
    1177             vn_adv(:,:) = z1_2 * ( vb2_b(:,:) + zwy(:,:) - atfp * vn_bf(:,:) ) 
     1176            un_adv(:,:) = r1_2 * ( ub2_b(:,:) + zwx(:,:) - atfp * un_bf(:,:) ) 
     1177            vn_adv(:,:) = r1_2 * ( vb2_b(:,:) + zwy(:,:) - atfp * vn_bf(:,:) ) 
    11781178            ! 
    11791179            ! Update corrective fluxes for next time step: 
     
    11941194      IF( ln_dynadv_vec .OR. ln_linssh ) THEN 
    11951195         DO jk=1,jpkm1 
    1196             ua(:,:,jk) = ua(:,:,jk) + ( ua_b(:,:) - ub_b(:,:) ) * z1_2dt_b 
    1197             va(:,:,jk) = va(:,:,jk) + ( va_b(:,:) - vb_b(:,:) ) * z1_2dt_b 
     1196            ua(:,:,jk) = ua(:,:,jk) + ( ua_b(:,:) - ub_b(:,:) ) * r1_2dt_b 
     1197            va(:,:,jk) = va(:,:,jk) + ( va_b(:,:) - vb_b(:,:) ) * r1_2dt_b 
    11981198         END DO 
    11991199      ELSE 
     
    12011201         DO jj = 1, jpjm1 
    12021202            DO ji = 1, jpim1      ! NO Vector Opt. 
    1203                zsshu_a(ji,jj) = z1_2 * umask(ji,jj,1)  * r1_e1e2u(ji,jj) & 
     1203               zsshu_a(ji,jj) = r1_2 * umask(ji,jj,1)  * r1_e1e2u(ji,jj) & 
    12041204                  &              * ( e1e2t(ji  ,jj) * ssha(ji  ,jj)    & 
    12051205                  &              +   e1e2t(ji+1,jj) * ssha(ji+1,jj) ) 
    1206                zsshv_a(ji,jj) = z1_2 * vmask(ji,jj,1)  * r1_e1e2v(ji,jj) & 
     1206               zsshv_a(ji,jj) = r1_2 * vmask(ji,jj,1)  * r1_e1e2v(ji,jj) & 
    12071207                  &              * ( e1e2t(ji,jj  ) * ssha(ji,jj  )    & 
    12081208                  &              +   e1e2t(ji,jj+1) * ssha(ji,jj+1) ) 
     
    12121212         ! 
    12131213         DO jk=1,jpkm1 
    1214             ua(:,:,jk) = ua(:,:,jk) + r1_hu_n(:,:) * ( ua_b(:,:) - ub_b(:,:) * hu_b(:,:) ) * z1_2dt_b 
    1215             va(:,:,jk) = va(:,:,jk) + r1_hv_n(:,:) * ( va_b(:,:) - vb_b(:,:) * hv_b(:,:) ) * z1_2dt_b 
     1214            ua(:,:,jk) = ua(:,:,jk) + r1_hu_n(:,:) * ( ua_b(:,:) - ub_b(:,:) * hu_b(:,:) ) * r1_2dt_b 
     1215            va(:,:,jk) = va(:,:,jk) + r1_hv_n(:,:) * ( va_b(:,:) - vb_b(:,:) * hv_b(:,:) ) * r1_2dt_b 
    12161216         END DO 
    12171217         ! Save barotropic velocities not transport: 
Note: See TracChangeset for help on using the changeset viewer.