Changeset 12377 for NEMO/trunk/src/OCE/TRA/trabbc.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/OCE/TRA/trabbc.F90
r12276 r12377 44 44 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_qgh ! structure of input qgh (file informations, fields read) 45 45 46 !! * Substitutions 47 # include "do_loop_substitute.h90" 46 48 !!---------------------------------------------------------------------- 47 49 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 51 53 CONTAINS 52 54 53 SUBROUTINE tra_bbc( kt )55 SUBROUTINE tra_bbc( kt, Kmm, pts, Krhs ) 54 56 !!---------------------------------------------------------------------- 55 57 !! *** ROUTINE tra_bbc *** … … 73 75 !! Emile-Geay and Madec, 2009, Ocean Science. 74 76 !!---------------------------------------------------------------------- 75 INTEGER, INTENT(in) :: kt ! ocean time-step index 77 INTEGER, INTENT(in ) :: kt ! ocean time-step index 78 INTEGER, INTENT(in ) :: Kmm, Krhs ! time level indices 79 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts,jpt), INTENT(inout) :: pts ! active tracers and RHS of tracer equation 76 80 ! 77 81 INTEGER :: ji, jj ! dummy loop indices … … 83 87 IF( l_trdtra ) THEN ! Save the input temperature trend 84 88 ALLOCATE( ztrdt(jpi,jpj,jpk) ) 85 ztrdt(:,:,:) = tsa(:,:,:,jp_tem)89 ztrdt(:,:,:) = pts(:,:,:,jp_tem,Krhs) 86 90 ENDIF 87 91 ! ! Add the geothermal trend on temperature 88 DO jj = 2, jpjm1 89 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)) 91 END DO 92 END DO 92 DO_2D_00_00 93 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) 94 END_2D 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 104 CALL iom_put ( "hfgeou" , rau0_rcp * qgh_trd0(:,:) ) 103 ! 104 IF(ln_ctl) CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' bbc - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 105 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' bbc - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 105 106 ! 106 107 IF( ln_timing ) CALL timing_stop('tra_bbc') … … 135 136 !!---------------------------------------------------------------------- 136 137 ! 137 REWIND( numnam_ref )138 138 READ ( numnam_ref, nambbc, IOSTAT = ios, ERR = 901) 139 139 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambbc in reference namelist' ) 140 140 ! 141 REWIND( numnam_cfg )142 141 READ ( numnam_cfg, nambbc, IOSTAT = ios, ERR = 902 ) 143 142 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nambbc in configuration namelist' )
Note: See TracChangeset
for help on using the changeset viewer.