- Timestamp:
- 2018-06-21T11:58:42+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_nemo2cice_prints/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90
r9816 r9817 27 27 USE trcsbc ! surface boundary condition (trc_sbc routine) 28 28 USE zpshde ! partial step: hor. derivative (zps_hde routine) 29 # if defined key_debug_medusa 30 USE trcrst 31 # endif 32 29 33 30 34 #if defined key_agrif … … 65 69 ! 66 70 CALL trc_sbc( kstp ) ! surface boundary condition 71 # if defined key_debug_medusa 72 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_sbc at kt =', kstp 73 CALL trc_rst_tra_stat 74 CALL flush(numout) 75 # endif 67 76 IF( lk_trabbl ) CALL trc_bbl( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme 68 77 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 only70 78 CALL trc_adv( kstp ) ! horizontal & vertical advection 79 # if defined key_debug_medusa 80 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_adv at kt =', kstp 81 CALL trc_rst_tra_stat 82 CALL flush(numout) 83 # endif 84 85 IF( ln_zps ) THEN 86 IF( ln_isfcav ) THEN ; CALL zps_hde_isf( kstp, jptra, trb, pgtu=gtru, pgtv=gtrv, pgtui=gtrui, pgtvi=gtrvi ) ! both top & bottom 87 ELSE ; CALL zps_hde ( kstp, jptra, trb, gtru, gtrv ) ! only bottom 88 ENDIF 89 ENDIF 71 90 CALL trc_ldf( kstp ) ! lateral mixing 72 91 IF( .NOT. lk_offline .AND. lk_zdfkpp ) & … … 76 95 #endif 77 96 CALL trc_zdf( kstp ) ! vertical mixing and after tracer fields 97 # if defined key_debug_medusa 98 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_zdf at kt =', kstp 99 CALL trc_rst_tra_stat 100 CALL flush(numout) 101 # endif 78 102 CALL trc_nxt( kstp ) ! tracer fields at next time step 103 # if defined key_debug_medusa 104 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_nxt at kt =', kstp 105 CALL trc_rst_tra_stat 106 CALL flush(numout) 107 # endif 79 108 IF( ln_trcrad ) CALL trc_rad( kstp ) ! Correct artificial negative concentrations 109 IF( ln_trcdmp_clo ) CALL trc_dmp_clo( kstp ) ! internal damping trends on closed seas only 80 110 81 111 #if defined key_agrif … … 83 113 #endif 84 114 85 IF( ln_zps .AND. .NOT. ln_isfcav) &86 & CALL zps_hde ( kstp, jptra, trn, gtru, gtrv ) ! Partial steps: now horizontal gradient of passive87 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 passive89 ! tracers at the bottom ocean level90 !91 115 ELSE ! 1D vertical configuration 92 116 CALL trc_sbc( kstp ) ! surface boundary condition … … 100 124 ! 101 125 IF( nn_timing == 1 ) CALL timing_stop('trc_trp') 126 ! 127 9400 FORMAT(a25,i4,D23.16) 102 128 ! 103 129 END SUBROUTINE trc_trp
Note: See TracChangeset
for help on using the changeset viewer.