- Timestamp:
- 2019-11-22T15:29:17+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Property svn:mergeinfo deleted
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/TRA/trabbc.F90
r11536 r11949 51 51 CONTAINS 52 52 53 SUBROUTINE tra_bbc( kt )53 SUBROUTINE tra_bbc( kt, Kmm, pts, Krhs ) 54 54 !!---------------------------------------------------------------------- 55 55 !! *** ROUTINE tra_bbc *** … … 73 73 !! Emile-Geay and Madec, 2009, Ocean Science. 74 74 !!---------------------------------------------------------------------- 75 INTEGER, INTENT(in) :: kt ! ocean time-step index 75 INTEGER, INTENT(in ) :: kt ! ocean time-step index 76 INTEGER, INTENT(in ) :: Kmm, Krhs ! time level indices 77 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts,jpt), INTENT(inout) :: pts ! active tracers and RHS of tracer equation 76 78 ! 77 79 INTEGER :: ji, jj ! dummy loop indices … … 83 85 IF( l_trdtra ) THEN ! Save the input temperature trend 84 86 ALLOCATE( ztrdt(jpi,jpj,jpk) ) 85 ztrdt(:,:,:) = tsa(:,:,:,jp_tem)87 ztrdt(:,:,:) = pts(:,:,:,jp_tem,Krhs) 86 88 ENDIF 87 89 ! ! Add the geothermal trend on temperature 88 90 DO jj = 2, jpjm1 89 91 DO ji = 2, jpim1 90 tsa(ji,jj,mbkt(ji,jj),jp_tem) = tsa(ji,jj,mbkt(ji,jj),jp_tem) + qgh_trd0(ji,jj) / e3t_n(ji,jj,mbkt(ji,jj))92 pts(ji,jj,mbkt(ji,jj),jp_tem,Krhs) = pts(ji,jj,mbkt(ji,jj),jp_tem,Krhs) + qgh_trd0(ji,jj) / e3t(ji,jj,mbkt(ji,jj),Kmm) 91 93 END DO 92 94 END DO 93 95 ! 94 CALL lbc_lnk( 'trabbc', tsa(:,:,:,jp_tem) , 'T', 1. )96 CALL lbc_lnk( 'trabbc', pts(:,:,:,jp_tem,Krhs) , 'T', 1. ) 95 97 ! 96 98 IF( l_trdtra ) THEN ! Send the trend for diagnostics 97 ztrdt(:,:,:) = tsa(:,:,:,jp_tem) - ztrdt(:,:,:)98 CALL trd_tra( kt, 'TRA', jp_tem, jptra_bbc, ztrdt )99 ztrdt(:,:,:) = pts(:,:,:,jp_tem,Krhs) - ztrdt(:,:,:) 100 CALL trd_tra( kt, Kmm, Krhs, 'TRA', jp_tem, jptra_bbc, ztrdt ) 99 101 DEALLOCATE( ztrdt ) 100 102 ENDIF 101 103 ! 102 IF(ln_ctl) CALL prt_ctl( tab3d_1= tsa(:,:,:,jp_tem), clinfo1=' bbc - Ta: ', mask1=tmask, clinfo3='tra-ta' )104 IF(ln_ctl) CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' bbc - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 103 105 ! 104 106 IF( ln_timing ) CALL timing_stop('tra_bbc')
Note: See TracChangeset
for help on using the changeset viewer.