Changeset 10954 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/TRA/trabbc.F90
- Timestamp:
- 2019-05-09T18:12:29+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/TRA/trabbc.F90
r10946 r10954 84 84 IF( l_trdtra ) THEN ! Save the input temperature trend 85 85 ALLOCATE( ztrdt(jpi,jpj,jpk) ) 86 ztrdt(:,:,:) = ts a(:,:,:,jp_tem)86 ztrdt(:,:,:) = ts(:,:,:,jp_tem,Krhs) 87 87 ENDIF 88 88 ! ! Add the geothermal trend on temperature 89 89 DO jj = 2, jpjm1 90 90 DO ji = 2, jpim1 91 ts a(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))91 ts(ji,jj,mbkt(ji,jj),jp_tem,Krhs) = ts(ji,jj,mbkt(ji,jj),jp_tem,Krhs) + qgh_trd0(ji,jj) / e3t(ji,jj,mbkt(ji,jj),Kmm) 92 92 END DO 93 93 END DO 94 94 ! 95 CALL lbc_lnk( 'trabbc', ts a(:,:,:,jp_tem) , 'T', 1. )95 CALL lbc_lnk( 'trabbc', ts(:,:,:,jp_tem,Krhs) , 'T', 1. ) 96 96 ! 97 97 IF( l_trdtra ) THEN ! Send the trend for diagnostics 98 ztrdt(:,:,:) = ts a(:,:,:,jp_tem) - ztrdt(:,:,:)98 ztrdt(:,:,:) = ts(:,:,:,jp_tem,Krhs) - ztrdt(:,:,:) 99 99 CALL trd_tra( kt, Kmm, Krhs, 'TRA', jp_tem, jptra_bbc, ztrdt ) 100 100 DEALLOCATE( ztrdt ) 101 101 ENDIF 102 102 ! 103 IF(ln_ctl) CALL prt_ctl( tab3d_1=ts a(:,:,:,jp_tem), clinfo1=' bbc - Ta: ', mask1=tmask, clinfo3='tra-ta' )103 IF(ln_ctl) CALL prt_ctl( tab3d_1=ts(:,:,:,jp_tem,Krhs), clinfo1=' bbc - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 104 104 ! 105 105 IF( ln_timing ) CALL timing_stop('tra_bbc')
Note: See TracChangeset
for help on using the changeset viewer.