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 13899 for NEMO/branches/2020/tickets_icb_1900/src/SWE/dynldf_lap_blp.F90 – NEMO

Ignore:
Timestamp:
2020-11-27T17:26:33+01:00 (4 years ago)
Author:
mathiot
Message:

ticket #1900: update branch to trunk and add ICB test case

Location:
NEMO/branches/2020/tickets_icb_1900
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/tickets_icb_1900

    • Property svn:externals
      •  

        old new  
        22^/utils/build/makenemo@HEAD   makenemo 
        33^/utils/build/mk@HEAD         mk 
        4 ^/utils/tools/@HEAD           tools 
         4^/utils/tools@HEAD            tools 
        55^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
         
        88 
        99# SETTE 
        10 ^/utils/CI/sette@12931        sette 
         10^/utils/CI/sette@13559        sette 
  • NEMO/branches/2020/tickets_icb_1900/src/SWE/dynldf_lap_blp.F90

    r12983 r13899  
    9797            DO jk = 1, jpkm1                                 ! Horizontal slab 
    9898               !                                              
    99                DO_2D_01_01 
     99               DO_2D( 0, 1, 0, 1 ) 
    100100               !                                      ! ahm * e3 * curl  (computed from 1 to jpim1/jpjm1) 
    101101!!gm open question here : e3f  at before or now ?    probably now...  
     
    112112               END_2D 
    113113               ! 
    114                DO_2D_00_00 
     114               DO_2D( 0, 0, 0, 0 ) 
    115115                  pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * (                                             & 
    116116                     &              - ( zcur(ji  ,jj) - zcur(ji,jj-1) ) * r1_e2u(ji,jj) / e3u(ji,jj,jk,Kmm)   & 
     
    128128            DO jk = 1, jpkm1                                 ! Horizontal slab 
    129129               ! 
    130                DO_2D_01_01 
     130               DO_2D( 0, 1, 0, 1 ) 
    131131                  !                                      ! shearing stress component (F-point)   NB : ahmf has already been multiplied by fmask 
    132132                  zshe(ji-1,jj-1) = ahmf(ji-1,jj-1,jk)                                                              & 
     
    143143               END_2D 
    144144               ! 
    145                DO_2D_00_00 
     145               DO_2D( 0, 0, 0, 0 ) 
    146146                  pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * r1_e1e2u(ji,jj) / e3u(ji,jj,jk,Kmm)                               & 
    147147                     &    * (   (   zten(ji+1,jj  ) * e2t(ji+1,jj  )*e2t(ji+1,jj  ) * e3t(ji+1,jj  ,jk,Kmm)                       & 
     
    164164            DO jk = 1, jpkm1                                 ! Horizontal slab 
    165165               ! 
    166                DO_2D_01_01 
     166               DO_2D( 0, 1, 0, 1 ) 
    167167                  !                                      ! shearing stress component (F-point)   NB : ahmf has already been multiplied by fmask 
    168168                  zshe(ji-1,jj-1) = ahmf(ji-1,jj-1,jk)                                           & 
     
    175175               END_2D 
    176176               ! 
    177                DO_2D_00_00 
     177               DO_2D( 0, 0, 0, 0 ) 
    178178                  pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * r1_e1e2u(ji,jj) / e3u(ji,jj,jk,Kmm)   & 
    179179                     &    * (   zten(ji+1,jj  ) * e2t(ji+1,jj  ) * e3t(ji+1,jj  ,jk,Kmm)              & 
Note: See TracChangeset for help on using the changeset viewer.