- Timestamp:
- 2017-06-25T12:26:32+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/TOP_SRC/trcstp.F90
r7812 r8215 8 8 #if defined key_top 9 9 !!---------------------------------------------------------------------- 10 !! trc_stp : passive tracer system time-stepping11 !!---------------------------------------------------------------------- 12 USE oce_trc 10 !! trc_stp : passive tracer system time-stepping 11 !!---------------------------------------------------------------------- 12 USE oce_trc ! ocean dynamics and active tracers variables 13 13 USE sbc_oce 14 14 USE trc 15 USE trctrp 16 USE trcsms 15 USE trctrp ! passive tracers transport 16 USE trcsms ! passive tracers sources and sinks 17 17 USE trcwri 18 18 USE trcrst 19 USE trcsub ! 19 20 USE trdtrc_oce 20 21 USE trdmxl_trc 21 USE prtctl_trc ! Print control for debbuging22 USE iom23 USE i n_out_manager24 USE trcsub22 ! 23 USE prtctl_trc ! Print control for debbuging 24 USE iom ! 25 USE in_out_manager ! 25 26 26 27 IMPLICIT NONE … … 29 30 PUBLIC trc_stp ! called by step 30 31 31 REAL(wp), DIMENSION(:,:,:), SAVE, ALLOCATABLE :: qsr_arr ! save qsr during TOP time-step32 REAL(wp) :: rdt_sampl33 INTEGER :: nb_rec_per_day, ktdcy34 REAL(wp) :: rsecfst, rseclast35 LOGICAL :: llnew32 LOGICAL :: llnew ! ??? 33 REAL(wp) :: rdt_sampl ! ??? 34 INTEGER :: nb_rec_per_day, ktdcy ! ??? 35 REAL(wp) :: rsecfst, rseclast ! ??? 36 REAL(wp), DIMENSION(:,:,:), SAVE, ALLOCATABLE :: qsr_arr ! save qsr during TOP time-step 36 37 37 38 !!---------------------------------------------------------------------- … … 46 47 !! *** ROUTINE trc_stp *** 47 48 !! 48 !! ** Purpose : Time loop of opa for passive tracer49 !! ** Purpose : Time loop of opa for passive tracer 49 50 !! 50 !! ** Method : 51 !! Compute the passive tracers trends 52 !! Update the passive tracers 51 !! ** Method : Compute the passive tracers trends 52 !! Update the passive tracers 53 53 !!------------------------------------------------------------------- 54 INTEGER, INTENT( in ) :: kt! ocean time-step index55 INTEGER :: jk, jn ! dummy loop indices56 REAL(wp) :: ztrai57 CHARACTER (len=25) :: charout58 54 INTEGER, INTENT( in ) :: kt ! ocean time-step index 55 ! 56 INTEGER :: jk, jn ! dummy loop indices 57 REAL(wp):: ztrai ! local scalar 58 CHARACTER (len=25) :: charout ! 59 59 !!------------------------------------------------------------------- 60 60 ! … … 115 115 ! 116 116 END SUBROUTINE trc_stp 117 117 118 118 119 SUBROUTINE trc_mean_qsr( kt ) … … 128 129 !! In coupled mode, the sampling is done at every coupling frequency 129 130 !!---------------------------------------------------------------------- 130 INTEGER, INTENT(in) :: kt 131 INTEGER :: jn 132 REAL(wp) :: zkt, zrec 133 CHARACTER(len=1) :: cl1 ! 1 character 134 CHARACTER(len=2) :: cl2 ! 2 characters 135 131 INTEGER, INTENT( in ) :: kt ! ocean time-step index 132 ! 133 INTEGER :: jn ! dummy loop indices 134 REAL(wp) :: zkt, zrec ! local scalars 135 CHARACTER(len=1) :: cl1 ! 1 character 136 CHARACTER(len=2) :: cl2 ! 2 characters 137 !!---------------------------------------------------------------------- 138 ! 136 139 IF( kt == nittrc000 ) THEN 137 140 IF( ln_cpl ) THEN … … 143 146 ENDIF 144 147 ! 145 IF( lwp) THEN148 IF(lwp) THEN 146 149 WRITE(numout,*) 147 150 WRITE(numout,*) ' Sampling frequency dt = ', rdt_sampl, 's',' Number of sampling per day nrec = ', nb_rec_per_day … … 171 174 CALL iom_get( numrtr, jpdom_autoglo, 'qsr_arr_'//cl2, qsr_arr(:,:,jn) ) ! A mean of qsr 172 175 ENDIF 173 END DO176 END DO 174 177 ELSE 175 178 DO jn = 1, nb_rec_per_day … … 184 187 DO jn = 1, nb_rec_per_day 185 188 qsr_arr(:,:,jn) = qsr_mean(:,:) 186 END DO189 END DO 187 190 ENDIF 188 191 ! … … 220 223 CALL iom_rstput( kt, nitrst, numrtw, 'qsr_arr_'//cl2, qsr_arr(:,:,jn) ) 221 224 ENDIF 222 END DO225 END DO 223 226 CALL iom_rstput( kt, nitrst, numrtw, 'qsr_mean', qsr_mean(:,:) ) 224 227 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.