- Timestamp:
- 2020-09-29T12:41:06+02:00 (3 years ago)
- Location:
- NEMO/branches/2020/r12377_ticket2386
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r12377_ticket2386
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@13507 sette
-
- Property svn:externals
-
NEMO/branches/2020/r12377_ticket2386/src/TOP/TRP/trcsbc.F90
r12511 r13540 18 18 USE oce_trc ! ocean dynamics and active tracers variables 19 19 USE trc ! ocean passive tracers variables 20 USE prtctl _trc! Print control for debbuging20 USE prtctl ! Print control for debbuging 21 21 USE iom 22 22 USE trd_oce … … 30 30 !! * Substitutions 31 31 # include "do_loop_substitute.h90" 32 # include "domzgr_substitute.h90" 32 33 !!---------------------------------------------------------------------- 33 34 !! NEMO/TOP 4.0 , NEMO Consortium (2018) … … 49 50 !! The surface freshwater flux modify the ocean volume 50 51 !! and thus the concentration of a tracer as : 51 !! tr(Krhs) = tr(Krhs) + emp * tr(Kmm) / e3t for k=152 !! tr(Krhs) = tr(Krhs) + emp * tr(Kmm) / e3t_ for k=1 52 53 !! where emp, the surface freshwater budget (evaporation minus 53 54 !! precipitation ) given in kg/m2/s is divided … … 87 88 zfact = 0.5_wp 88 89 DO jn = 1, jptra 89 CALL iom_get( numrtr, jpdom_auto glo, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn) ) ! before tracer content sbc90 CALL iom_get( numrtr, jpdom_auto, 'sbc_'//TRIM(ctrcnm(jn))//'_b', sbc_trc_b(:,:,jn) ) ! before tracer content sbc 90 91 END DO 91 92 ELSE ! No restart or restart not found: Euler forward time stepping … … 120 121 ! 121 122 DO jn = 1, jptra 122 DO_2D _01_00123 DO_2D( 0, 1, 0, 0 ) 123 124 sbc_trc(ji,jj,jn) = zsfx(ji,jj) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 124 125 END_2D … … 128 129 ! 129 130 DO jn = 1, jptra 130 DO_2D _01_00131 DO_2D( 0, 1, 0, 0 ) 131 132 sbc_trc(ji,jj,jn) = ( zsfx(ji,jj) + fmmflx(ji,jj) ) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 132 133 END_2D … … 136 137 ! 137 138 DO jn = 1, jptra 138 DO_2D _01_00139 DO_2D( 0, 1, 0, 0 ) 139 140 zse3t = 1. / e3t(ji,jj,1,Kmm) 140 141 ! tracer flux at the ice/ocean interface (tracer/m2/s) … … 154 155 END SELECT 155 156 ! 156 CALL lbc_lnk( 'trcsbc', sbc_trc(:,:,:), 'T', 1. )157 CALL lbc_lnk( 'trcsbc', sbc_trc(:,:,:), 'T', 1.0_wp ) 157 158 ! Concentration dilution effect on tracers due to evaporation & precipitation 158 159 DO jn = 1, jptra … … 160 161 IF( l_trdtrc ) ztrtrd(:,:,:) = ptr(:,:,:,jn,Krhs) ! save trends 161 162 ! 162 DO_2D _01_00163 DO_2D( 0, 1, 0, 0 ) 163 164 zse3t = zfact / e3t(ji,jj,1,Kmm) 164 165 ptr(ji,jj,1,jn,Krhs) = ptr(ji,jj,1,jn,Krhs) + ( sbc_trc_b(ji,jj,jn) + sbc_trc(ji,jj,jn) ) * zse3t … … 186 187 ! 187 188 IF( sn_cfctl%l_prttrc ) THEN 188 WRITE(charout, FMT="('sbc ')") ; CALL prt_ctl_ trc_info(charout)189 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' ) 190 191 ENDIF 191 192 IF( l_trdtrc ) DEALLOCATE( ztrtrd )
Note: See TracChangeset
for help on using the changeset viewer.