- Timestamp:
- 2015-09-24T08:31:40+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90
r5120 r5758 48 48 CONTAINS 49 49 50 SUBROUTINE trc_trp( k stp)50 SUBROUTINE trc_trp( kt ) 51 51 !!---------------------------------------------------------------------- 52 52 !! *** ROUTINE trc_trp *** … … 57 57 !! - Update the passive tracers 58 58 !!---------------------------------------------------------------------- 59 INTEGER, INTENT( in ) :: k stp! ocean time-step index59 INTEGER, INTENT( in ) :: kt ! ocean time-step index 60 60 !! --------------------------------------------------------------------- 61 61 ! … … 64 64 IF( .NOT. lk_c1d ) THEN 65 65 ! 66 CALL trc_sbc( kstp ) ! surface boundary condition 67 IF( lk_trabbl ) CALL trc_bbl( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme 68 IF( ln_trcdmp ) CALL trc_dmp( kstp ) ! internal damping trends 69 IF( ln_trcdmp_clo ) CALL trc_dmp_clo( kstp ) ! internal damping trends on closed seas only 70 CALL trc_adv( kstp ) ! horizontal & vertical advection 71 CALL trc_ldf( kstp ) ! lateral mixing 72 IF( .NOT. lk_offline .AND. lk_zdfkpp ) & 73 & CALL trc_kpp( kstp ) ! KPP non-local tracer fluxes 66 CALL trc_sbc ( kt ) ! surface boundary condition 67 IF( lk_trabbl ) CALL trc_bbl ( kt ) ! advective (and/or diffusive) bottom boundary layer scheme 68 IF( ln_trcdmp ) CALL trc_dmp ( kt ) ! internal damping trends 69 IF( ln_trcdmp_clo ) CALL trc_dmp_clo( kt ) ! internal damping trends on closed seas only 70 CALL trc_adv ( kt ) ! horizontal & vertical advection 71 CALL trc_ldf ( kt ) ! lateral mixing 74 72 #if defined key_agrif 75 IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_trc 73 IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_trc ! tracers sponge 76 74 #endif 77 CALL trc_zdf ( kstp )! vertical mixing and after tracer fields78 CALL trc_nxt ( kstp )! tracer fields at next time step79 IF( ln_trcrad ) CALL trc_rad ( kstp )! Correct artificial negative concentrations75 CALL trc_zdf ( kt ) ! vertical mixing and after tracer fields 76 CALL trc_nxt ( kt ) ! tracer fields at next time step 77 IF( ln_trcrad ) CALL trc_rad ( kt ) ! Correct artificial negative concentrations 80 78 81 79 #if defined key_agrif 82 IF( .NOT. Agrif_Root()) CALL Agrif_Update_Trc( kstp )! Update tracer at AGRIF zoom boundaries : children only80 IF( .NOT.Agrif_Root()) CALL Agrif_Update_Trc( kt ) ! Update tracer at AGRIF zoom boundaries : children only 83 81 #endif 84 82 85 IF( ln_zps .AND. .NOT. ln_isfcav) & 86 & CALL zps_hde ( kstp, jptra, trn, gtru, gtrv ) ! Partial steps: now horizontal gradient of passive 87 IF( ln_zps .AND. ln_isfcav) & 88 & CALL zps_hde_isf( kstp, jptra, trn, pgtu=gtru, pgtv=gtrv, pgtui=gtrui, pgtvi=gtrvi ) ! Partial steps: now horizontal gradient of passive 89 ! tracers at the bottom ocean level 83 ! ! Partial top/bottom cell: GRADh( trn ) 84 IF( ln_isfcav .AND. ln_zps ) THEN ; CALL zps_hde_isf( kt, jptra, trn, gtru, gtrv, gtrui, gtrvi ) ! both top & bottom 85 ELSEIF( ln_zps ) THEN ; CALL zps_hde ( kt, jptra, trn, gtru, gtrv ) ! only bottom 86 ENDIF 87 !!gm IF( ln_zps ) THEN 88 ! & CALL zps_hde ( kt, jptra, trn, gtru, gtrv ) ! Partial steps: now horizontal gradient of passive 89 ! IF( ln_isfcav) & 90 ! & CALL zps_hde_isf( kt, jptra, trn, pgtu=gtru, pgtv=gtrv, pgtui=gtrui, pgtvi=gtrvi ) ! Partial steps: now horizontal gradient of passive 91 !!gm ENDIF 90 92 ! 91 93 ELSE ! 1D vertical configuration 92 CALL trc_sbc( k stp) ! surface boundary condition94 CALL trc_sbc( kt ) ! surface boundary condition 93 95 IF( .NOT. lk_offline .AND. lk_zdfkpp ) & 94 & CALL trc_kpp( k stp) ! KPP non-local tracer fluxes95 CALL trc_zdf( k stp) ! vertical mixing and after tracer fields96 CALL trc_nxt( k stp) ! tracer fields at next time step97 IF( ln_trcrad ) CALL trc_rad( k stp) ! Correct artificial negative concentrations96 & CALL trc_kpp( kt ) ! KPP non-local tracer fluxes 97 CALL trc_zdf( kt ) ! vertical mixing and after tracer fields 98 CALL trc_nxt( kt ) ! tracer fields at next time step 99 IF( ln_trcrad ) CALL trc_rad( kt ) ! Correct artificial negative concentrations 98 100 ! 99 101 END IF … … 108 110 !!---------------------------------------------------------------------- 109 111 CONTAINS 110 SUBROUTINE trc_trp( k stp) ! Empty routine111 INTEGER, INTENT(in) :: k stp112 WRITE(*,*) 'trc_trp: You should not have seen this print! error?', k stp112 SUBROUTINE trc_trp( kt ) ! Empty routine 113 INTEGER, INTENT(in) :: kt 114 WRITE(*,*) 'trc_trp: You should not have seen this print! error?', kt 113 115 END SUBROUTINE trc_trp 114 116 #endif
Note: See TracChangeset
for help on using the changeset viewer.