Changeset 14644 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/TRA/trasbc.F90
- Timestamp:
- 2021-03-26T15:33:49+01:00 (3 years ago)
- 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 9 9 10 10 # SETTE 11 ^/utils/CI/sette _wave@13990sette11 ^/utils/CI/sette@14244 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/TRA/trasbc.F90
r14219 r14644 106 106 !!gm This should be moved into sbcmod.F90 module ? (especially now that ln_traqsr is read in namsbc namelist) 107 107 IF( .NOT.ln_traqsr ) THEN ! no solar radiation penetration 108 DO_2D( is j, iej, isi, iei)108 DO_2D( isi, iei, isj, iej ) 109 109 qns(ji,jj) = qns(ji,jj) + qsr(ji,jj) ! total heat flux in qns 110 110 qsr(ji,jj) = 0._wp ! qsr set to zero … … 127 127 ELSE ! No restart or restart not found: Euler forward time stepping 128 128 zfact = 1._wp 129 DO_2D( is j, iej, isi, iei)129 DO_2D( isi, iei, isj, iej ) 130 130 sbc_tsc(ji,jj,:) = 0._wp 131 131 sbc_tsc_b(ji,jj,:) = 0._wp … … 134 134 ELSE !* other time-steps: swap of forcing fields 135 135 zfact = 0.5_wp 136 DO_2D( is j, iej, isi, iei)136 DO_2D( isi, iei, isj, iej ) 137 137 sbc_tsc_b(ji,jj,:) = sbc_tsc(ji,jj,:) 138 138 END_2D 139 139 ENDIF 140 140 ! !== Now sbc tracer content fields ==! 141 DO_2D( is j, iej, isi, iei)141 DO_2D( isi, iei, isj, iej ) 142 142 sbc_tsc(ji,jj,jp_tem) = r1_rho0_rcp * qns(ji,jj) ! non solar heat flux 143 143 sbc_tsc(ji,jj,jp_sal) = r1_rho0 * sfx(ji,jj) ! salt flux due to freezing/melting 144 144 END_2D 145 145 IF( ln_linssh ) THEN !* linear free surface 146 DO_2D( is j, iej, isi, iei) !==>> add concentration/dilution effect due to constant volume cell146 DO_2D( isi, iei, isj, iej ) !==>> add concentration/dilution effect due to constant volume cell 147 147 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 148 sbc_tsc(ji,jj,jp_sal) = sbc_tsc(ji,jj,jp_sal) + r1_rho0 * emp(ji,jj) * pts(ji,jj,1,jp_sal,Kmm)
Note: See TracChangeset
for help on using the changeset viewer.