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/TOP/TRP/trcsbc.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/TOP/TRP/trcsbc.F90

    r13237 r13899  
    1818   USE oce_trc         ! ocean dynamics and active tracers variables 
    1919   USE trc             ! ocean  passive tracers variables 
    20    USE prtctl_trc      ! Print control for debbuging 
     20   USE prtctl          ! Print control for debbuging 
    2121   USE iom 
    2222   USE trd_oce 
     
    8888            zfact = 0.5_wp 
    8989            DO jn = 1, jptra 
    90                CALL iom_get( numrtr, jpdom_autoglo, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn) )   ! before tracer content sbc 
     90               CALL iom_get( numrtr, jpdom_auto, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn) )   ! before tracer content sbc 
    9191            END DO 
    9292         ELSE                                         ! No restart or restart not found: Euler forward time stepping 
     
    121121         ! 
    122122         DO jn = 1, jptra 
    123             DO_2D_01_00 
     123            DO_2D( 0, 1, 0, 0 ) 
    124124               sbc_trc(ji,jj,jn) = zsfx(ji,jj) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 
    125125            END_2D 
     
    129129         ! 
    130130         DO jn = 1, jptra 
    131             DO_2D_01_00 
     131            DO_2D( 0, 1, 0, 0 ) 
    132132               sbc_trc(ji,jj,jn) = ( zsfx(ji,jj) + fmmflx(ji,jj) ) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 
    133133            END_2D 
     
    137137         ! 
    138138         DO jn = 1, jptra 
    139             DO_2D_01_00 
     139            DO_2D( 0, 1, 0, 0 ) 
    140140               zse3t = 1. / e3t(ji,jj,1,Kmm) 
    141141               ! tracer flux at the ice/ocean interface (tracer/m2/s) 
     
    161161         IF( l_trdtrc )   ztrtrd(:,:,:) = ptr(:,:,:,jn,Krhs)  ! save trends 
    162162         ! 
    163          DO_2D_01_00 
     163         DO_2D( 0, 1, 0, 0 ) 
    164164            zse3t = zfact / e3t(ji,jj,1,Kmm) 
    165165            ptr(ji,jj,1,jn,Krhs) = ptr(ji,jj,1,jn,Krhs) + ( sbc_trc_b(ji,jj,jn) + sbc_trc(ji,jj,jn) ) * zse3t 
     
    187187      ! 
    188188      IF( sn_cfctl%l_prttrc )   THEN 
    189          WRITE(charout, FMT="('sbc ')") ;  CALL prt_ctl_trc_info(charout) 
    190                                            CALL prt_ctl_trc( tab4d=ptr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
     189         WRITE(charout, FMT="('sbc ')") ;  CALL prt_ctl_info( charout, cdcomp = 'top' ) 
     190                                           CALL prt_ctl( tab4d_1=ptr(:,:,:,:,Krhs), mask1=tmask, clinfo=ctrcnm, clinfo3='trd' ) 
    191191      ENDIF 
    192192      IF( l_trdtrc )  DEALLOCATE( ztrtrd ) 
Note: See TracChangeset for help on using the changeset viewer.