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 14644 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/TRA/traqsr.F90 – NEMO

Ignore:
Timestamp:
2021-03-26T15:33:49+01:00 (3 years ago)
Author:
sparonuz
Message:

Merge trunk -r14642:HEAD

Location:
NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final

    • Property svn:externals

      •  

        old new  
        99 
        1010# SETTE 
        11 ^/utils/CI/sette_wave@13990         sette 
         11^/utils/CI/sette@14244        sette 

  • NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/TRA/traqsr.F90

    r14219 r14644  
    152152         ELSE                                           ! No restart or Euler forward at 1st time step 
    153153            z1_2 = 1._wp 
    154             DO_3D( isj, iej, isi, iei, 1, jpk ) 
     154            DO_3D( isi, iei, isj, iej, 1, jpk ) 
    155155               qsr_hc_b(ji,jj,jk) = 0._wp 
    156156            END_3D 
    … …  
    158158      ELSE                             !==  Swap of qsr heat content  ==! 
    159159         z1_2 = 0.5_wp 
    160          DO_3D( isj, iej, isi, iei, 1, jpk ) 
     160         DO_3D( isi, iei, isj, iej, 1, jpk ) 
    161161            qsr_hc_b(ji,jj,jk) = qsr_hc(ji,jj,jk) 
    162162         END_3D 
    … …  
    169169      CASE( np_BIO )                   !==  bio-model fluxes  ==! 
    170170         ! 
    171          DO_3D( isj, iej, isi, iei, 1, nksr ) 
     171         DO_3D( isi, iei, isj, iej, 1, nksr ) 
    172172            qsr_hc(ji,jj,jk) = r1_rho0_rcp * ( etot3(ji,jj,jk) - etot3(ji,jj,jk+1) ) 
    173173         END_3D 
    … …  
    191191            ! most expensive calculations) 
    192192            ! 
    193             DO_2D( isj, iej, isi, iei ) 
     193            DO_2D( isi, iei, isj, iej ) 
    194194                       ! zlogc = log(zchl) 
    195195               zlogc = LOG ( MIN( 10. , MAX( 0.03, sf_chl(1)%fnow(ji,jj,1) ) ) ) 
    … …  
    210210 
    211211! 
    212             DO_3D( isj, iej, isi, iei, 1, nksr + 1 ) 
     212            DO_3D( isi, iei, isj, iej, 1, nksr + 1 ) 
    213213               ! zchl    = ALOG( ze0(ji,jj) ) 
    214214               zlogc = ze0(ji,jj) 
    … …  
    240240         ! 
    241241         zcoef  = ( 1. - rn_abs ) / 3._wp    !* surface equi-partition in R-G-B 
    242          DO_2D( isj, iej, isi, iei ) 
     242         DO_2D( isi, iei, isj, iej ) 
    243243            ze0(ji,jj) = rn_abs * qsr(ji,jj) 
    244244            ze1(ji,jj) = zcoef  * qsr(ji,jj) 
    … …  
    251251         ! 
    252252         !                                    !* interior equi-partition in R-G-B depending on vertical profile of Chl 
    253          DO_3D( isj, iej, isi, iei, 2, nksr + 1 ) 
     253         DO_3D( isi, iei, isj, iej, 2, nksr + 1 ) 
    254254            ze3t = e3t(ji,jj,jk-1,Kmm) 
    255255            irgb = NINT( ztmp3d(ji,jj,jk) ) 
    … …  
    265265         END_3D 
    266266         ! 
    267          DO_3D( isj, iej, isi, iei, 1, nksr )          !* now qsr induced heat content 
     267         DO_3D( isi, iei, isj, iej, 1, nksr )          !* now qsr induced heat content 
    268268            qsr_hc(ji,jj,jk) = r1_rho0_rcp * ( ztmp3d(ji,jj,jk) - ztmp3d(ji,jj,jk+1) ) 
    269269         END_3D 
    … …  
    275275         zz0 =        rn_abs   * r1_rho0_rcp      ! surface equi-partition in 2-bands 
    276276         zz1 = ( 1. - rn_abs ) * r1_rho0_rcp 
    277          DO_3D( isj, iej, isi, iei, 1, nksr )          !* now qsr induced heat content 
     277         DO_3D( isi, iei, isj, iej, 1, nksr )          !* now qsr induced heat content 
    278278            zc0 = zz0 * EXP( -gdepw(ji,jj,jk  ,Kmm)*xsi0r ) + zz1 * EXP( -gdepw(ji,jj,jk  ,Kmm)*xsi1r ) 
    279279            zc1 = zz0 * EXP( -gdepw(ji,jj,jk+1,Kmm)*xsi0r ) + zz1 * EXP( -gdepw(ji,jj,jk+1,Kmm)*xsi1r ) 
    … …  
    293293      ! 
    294294      ! sea-ice: store the 1st ocean level attenuation coefficient 
    295       DO_2D( isj, iej, isi, iei ) 
     295      DO_2D( isi, iei, isj, iej ) 
    296296         IF( qsr(ji,jj) /= 0._wp ) THEN   ;   fraqsr_1lev(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rho0_rcp * qsr(ji,jj) ) 
    297297         ELSE                             ;   fraqsr_1lev(ji,jj) = 1._wp 
Note: See TracChangeset for help on using the changeset viewer.