Changeset 10966 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/TRP/trcrad.F90
- Timestamp:
- 2019-05-10T18:43:09+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/TOP/TRP/trcrad.F90
r10946 r10966 37 37 CONTAINS 38 38 39 SUBROUTINE trc_rad( kt, K mm, Krhs )39 SUBROUTINE trc_rad( kt, Kbb, Kmm, Krhs ) 40 40 !!---------------------------------------------------------------------- 41 41 !! *** ROUTINE trc_rad *** … … 52 52 !! (the total CFC content is not strictly preserved) 53 53 !!---------------------------------------------------------------------- 54 INTEGER, INTENT(in) :: kt ! ocean time-step index55 INTEGER, INTENT(in) :: K mm, Krhs ! time level indices54 INTEGER, INTENT(in) :: kt ! ocean time-step index 55 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 56 56 ! 57 57 CHARACTER (len=22) :: charout … … 60 60 IF( ln_timing ) CALL timing_start('trc_rad') 61 61 ! 62 IF( ln_age ) CALL trc_rad_sms( kt, Kmm, Krhs, tr b, trn, jp_age , jp_age ) ! AGE63 IF( ll_cfc ) CALL trc_rad_sms( kt, Kmm, Krhs, tr b, trn, jp_cfc0, jp_cfc1 ) ! CFC model64 IF( ln_c14 ) CALL trc_rad_sms( kt, Kmm, Krhs, tr b, trn, jp_c14 , jp_c14 ) ! C1465 IF( ln_pisces ) CALL trc_rad_sms( kt, Kmm, Krhs, tr b, trn, jp_pcs0, jp_pcs1, cpreserv='Y' ) ! PISCES model66 IF( ln_my_trc ) CALL trc_rad_sms( kt, Kmm, Krhs, tr b, trn, jp_myt0, jp_myt1 ) ! MY_TRC model62 IF( ln_age ) CALL trc_rad_sms( kt, Kmm, Krhs, tr(:,:,:,:,Kbb), tr(:,:,:,:,Kmm), jp_age , jp_age ) ! AGE 63 IF( ll_cfc ) CALL trc_rad_sms( kt, Kmm, Krhs, tr(:,:,:,:,Kbb), tr(:,:,:,:,Kmm), jp_cfc0, jp_cfc1 ) ! CFC model 64 IF( ln_c14 ) CALL trc_rad_sms( kt, Kmm, Krhs, tr(:,:,:,:,Kbb), tr(:,:,:,:,Kmm), jp_c14 , jp_c14 ) ! C14 65 IF( ln_pisces ) CALL trc_rad_sms( kt, Kmm, Krhs, tr(:,:,:,:,Kbb), tr(:,:,:,:,Kmm), jp_pcs0, jp_pcs1, cpreserv='Y' ) ! PISCES model 66 IF( ln_my_trc ) CALL trc_rad_sms( kt, Kmm, Krhs, tr(:,:,:,:,Kbb), tr(:,:,:,:,Kmm), jp_myt0, jp_myt1 ) ! MY_TRC model 67 67 ! 68 68 IF(ln_ctl) THEN ! print mean trends (used for debugging) 69 69 WRITE(charout, FMT="('rad')") 70 70 CALL prt_ctl_trc_info( charout ) 71 CALL prt_ctl_trc( tab4d=tr n, mask=tmask, clinfo=ctrcnm )71 CALL prt_ctl_trc( tab4d=tr(:,:,:,:,Kmm), mask=tmask, clinfo=ctrcnm ) 72 72 ENDIF 73 73 ! … … 160 160 DO jn = jp_sms0, jp_sms1 161 161 ! 162 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrb(:,:,:,jn) ! save input tr bfor trend computation162 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrb(:,:,:,jn) ! save input tr(:,:,:,:,Kbb) for trend computation 163 163 ! 164 164 DO jk = 1, jpkm1 … … 210 210 DO jn = jp_sms0, jp_sms1 211 211 ! 212 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrn(:,:,:,jn) ! save input tr bfor trend computation212 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrn(:,:,:,jn) ! save input tr for trend computation 213 213 ! 214 214 DO jk = 1, jpkm1 … … 257 257 DO jn = jp_sms0, jp_sms1 258 258 ! 259 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrb(:,:,:,jn) ! save input tr bfor trend computation259 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrb(:,:,:,jn) ! save input tr for trend computation 260 260 ! 261 261 WHERE( ptrb(:,:,:,jn) < 0. ) ptrb(:,:,:,jn) = 0. … … 266 266 ENDIF 267 267 ! 268 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrn(:,:,:,jn) ! save input tr nfor trend computation268 IF( l_trdtrc ) ztrtrd(:,:,:) = ptrn(:,:,:,jn) ! save input tr for trend computation 269 269 ! 270 270 WHERE( ptrn(:,:,:,jn) < 0. ) ptrn(:,:,:,jn) = 0. … … 288 288 !!---------------------------------------------------------------------- 289 289 CONTAINS 290 SUBROUTINE trc_rad( kt ) ! Empty routine290 SUBROUTINE trc_rad( kt, Kbb, Kmm, Krhs ) ! Empty routine 291 291 INTEGER, INTENT(in) :: kt 292 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 292 293 WRITE(*,*) 'trc_rad: You should not have seen this print! error?', kt 293 294 END SUBROUTINE trc_rad
Note: See TracChangeset
for help on using the changeset viewer.