Changeset 9987 for branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90
- Timestamp:
- 2018-07-23T11:33:03+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90
r7960 r9987 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( ln_zps ) THEN 80 IF( ln_isfcav ) THEN ; CALL zps_hde_isf( kstp, jptra, trb, pgtu=gtru, pgtv=gtrv, pgtui=gtrui, pgtvi=gtrvi ) ! both top & bottom 81 ELSE ; CALL zps_hde ( kstp, jptra, trb, gtru, gtrv ) ! only bottom 82 ENDIF 83 ENDIF 71 84 CALL trc_ldf( kstp ) ! lateral mixing 72 85 IF( .NOT. lk_offline .AND. lk_zdfkpp ) & … … 77 90 CALL trc_zdf( kstp ) ! vertical mixing and after tracer fields 78 91 CALL trc_nxt( kstp ) ! tracer fields at next time step 92 # if defined key_debug_medusa 93 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_nxt at kt =', kstp 94 CALL trc_rst_tra_stat 95 CALL flush(numout) 96 # endif 79 97 IF( ln_trcrad ) CALL trc_rad( kstp ) ! Correct artificial negative concentrations 98 IF( ln_trcdmp_clo ) CALL trc_dmp_clo( kstp ) ! internal damping trends on closed seas only 80 99 81 100 #if defined key_agrif … … 83 102 #endif 84 103 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 104 ELSE ! 1D vertical configuration 92 105 CALL trc_sbc( kstp ) ! surface boundary condition … … 100 113 ! 101 114 IF( nn_timing == 1 ) CALL timing_stop('trc_trp') 115 ! 116 9400 FORMAT(a25,i4,D23.16) 102 117 ! 103 118 END SUBROUTINE trc_trp
Note: See TracChangeset
for help on using the changeset viewer.