Changeset 14822 for NEMO/branches/2021/dev_r14122_HPC-08_Mueller_OSMOSIS_streamlining/src/OCE/TRA/trasbc.F90
- Timestamp:
- 2021-05-10T13:32:16+02:00 (3 years ago)
- Location:
- NEMO/branches/2021/dev_r14122_HPC-08_Mueller_OSMOSIS_streamlining
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14122_HPC-08_Mueller_OSMOSIS_streamlining
- Property svn:externals
-
old new 1 ^/utils/build/arch@ 14122arch2 ^/utils/build/makenemo@ 14122makenemo3 ^/utils/build/mk@ 14122mk4 ^/utils/tools@ 14122tools5 ^/vendors/AGRIF/dev@ 14122ext/AGRIF6 ^/vendors/FCM@ 14122ext/FCM7 ^/vendors/IOIPSL@ 14122ext/IOIPSL8 ^/vendors/PPR@ 14122ext/PPR1 ^/utils/build/arch@HEAD arch 2 ^/utils/build/makenemo@HEAD makenemo 3 ^/utils/build/mk@HEAD mk 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev@HEAD ext/AGRIF 6 ^/vendors/FCM@HEAD ext/FCM 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 ^/vendors/PPR@HEAD ext/PPR 9 9 10 10 # SETTE 11 ^/utils/CI/sette _wave@13990sette11 ^/utils/CI/sette@14244 sette
-
- Property svn:externals
-
NEMO/branches/2021/dev_r14122_HPC-08_Mueller_OSMOSIS_streamlining/src/OCE/TRA/trasbc.F90
r14072 r14822 98 98 ENDIF 99 99 ! 100 ! NOTE: [tiling-comms-merge] Many DO loop bounds changed to avoid changing results when using tiling. Some bounds were also adjusted to account for those changed in tra_atf101 100 IF( ntsi == Nis0 ) THEN ; isi = nn_hls ; ELSE ; isi = 0 ; ENDIF ! Avoid double-counting when using tiling 102 101 IF( ntsj == Njs0 ) THEN ; isj = nn_hls ; ELSE ; isj = 0 ; ENDIF … … 106 105 !!gm This should be moved into sbcmod.F90 module ? (especially now that ln_traqsr is read in namsbc namelist) 107 106 IF( .NOT.ln_traqsr ) THEN ! no solar radiation penetration 108 DO_2D( is j, iej, isi, iei)107 DO_2D( isi, iei, isj, iej ) 109 108 qns(ji,jj) = qns(ji,jj) + qsr(ji,jj) ! total heat flux in qns 110 109 qsr(ji,jj) = 0._wp ! qsr set to zero … … 127 126 ELSE ! No restart or restart not found: Euler forward time stepping 128 127 zfact = 1._wp 129 DO_2D( is j, iej, isi, iei)128 DO_2D( isi, iei, isj, iej ) 130 129 sbc_tsc(ji,jj,:) = 0._wp 131 130 sbc_tsc_b(ji,jj,:) = 0._wp … … 134 133 ELSE !* other time-steps: swap of forcing fields 135 134 zfact = 0.5_wp 136 DO_2D( is j, iej, isi, iei)135 DO_2D( isi, iei, isj, iej ) 137 136 sbc_tsc_b(ji,jj,:) = sbc_tsc(ji,jj,:) 138 137 END_2D 139 138 ENDIF 140 139 ! !== Now sbc tracer content fields ==! 141 DO_2D( is j, iej, isi, iei)140 DO_2D( isi, iei, isj, iej ) 142 141 sbc_tsc(ji,jj,jp_tem) = r1_rho0_rcp * qns(ji,jj) ! non solar heat flux 143 142 sbc_tsc(ji,jj,jp_sal) = r1_rho0 * sfx(ji,jj) ! salt flux due to freezing/melting 144 143 END_2D 145 144 IF( ln_linssh ) THEN !* linear free surface 146 DO_2D( is j, iej, isi, iei) !==>> add concentration/dilution effect due to constant volume cell145 DO_2D( isi, iei, isj, iej ) !==>> add concentration/dilution effect due to constant volume cell 147 146 sbc_tsc(ji,jj,jp_tem) = sbc_tsc(ji,jj,jp_tem) + r1_rho0 * emp(ji,jj) * pts(ji,jj,1,jp_tem,Kmm) 148 147 sbc_tsc(ji,jj,jp_sal) = sbc_tsc(ji,jj,jp_sal) + r1_rho0 * emp(ji,jj) * pts(ji,jj,1,jp_sal,Kmm) … … 155 154 ! 156 155 DO jn = 1, jpts !== update tracer trend ==! 157 ! NOTE: [tiling-comms-merge] This looped over nn_hls, which changes the results when using tiling158 156 DO_2D( 0, 0, 0, 0 ) 159 157 pts(ji,jj,1,jn,Krhs) = pts(ji,jj,1,jn,Krhs) + zfact * ( sbc_tsc_b(ji,jj,jn) + sbc_tsc(ji,jj,jn) ) &
Note: See TracChangeset
for help on using the changeset viewer.