Changeset 10963 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/trcstp.F90
- Timestamp:
- 2019-05-10T16:15:43+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/trcstp.F90
r10922 r10963 77 77 IF( .NOT.ln_linssh ) THEN ! update ocean volume due to ssh temporal evolution 78 78 DO jk = 1, jpk 79 cvol(:,:,jk) = e1e2t(:,:) * e3t _n(:,:,jk) * tmask(:,:,jk)79 cvol(:,:,jk) = e1e2t(:,:) * e3t(:,:,jk,Kmm) * tmask(:,:,jk) 80 80 END DO 81 81 IF ( ln_ctl .OR. kt == nitrst .OR. ( ln_check_mass .AND. kt == nitend ) & … … 87 87 IF( l_trcdm2dc ) CALL trc_mean_qsr( kt ) 88 88 ! 89 IF( nn_dttrc /= 1 ) CALL trc_sub_stp( kt, K mm) ! averaging physical variables for sub-stepping89 IF( nn_dttrc /= 1 ) CALL trc_sub_stp( kt, Kbb, Kmm, Krhs ) ! averaging physical variables for sub-stepping 90 90 ! 91 91 IF( MOD( kt , nn_dttrc ) == 0 ) THEN ! only every nn_dttrc time step … … 96 96 ENDIF 97 97 ! 98 tr a(:,:,:,:) = 0.e098 tr(:,:,:,:,Krhs) = 0.e0 99 99 ! 100 100 CALL trc_rst_opn ( kt ) ! Open tracer restart file 101 101 IF( lrst_trc ) CALL trc_rst_cal ( kt, 'WRITE' ) ! calendar 102 CALL trc_wri ( kt ) ! output of passive tracers with iom I/O manager103 CALL trc_sms ( kt ) ! tracers: sinks and sources104 CALL trc_trp ( kt, Kbb, Kmm, Krhs, Kaa ) 102 CALL trc_wri ( kt ) ! output of passive tracers with iom I/O manager 103 CALL trc_sms ( kt, Kmm, Krhs ) ! tracers: sinks and sources 104 CALL trc_trp ( kt, Kbb, Kmm, Krhs, Kaa ) ! transport of passive tracers 105 105 IF( kt == nittrc000 ) THEN 106 106 CALL iom_close( numrtr ) ! close input tracer restart file 107 107 IF(lwm) CALL FLUSH( numont ) ! flush namelist output 108 108 ENDIF 109 IF( lrst_trc ) CALL trc_rst_wri ( kt ) ! write tracer restart file110 IF( lk_trdmxl_trc ) CALL trd_mxl_trc ( kt ) ! trends: Mixed-layer111 ! 112 IF( nn_dttrc /= 1 ) CALL trc_sub_reset( kt ) ! resetting physical variables when sub-stepping109 IF( lrst_trc ) CALL trc_rst_wri ( kt, Kbb, Kmm, Krhs ) ! write tracer restart file 110 IF( lk_trdmxl_trc ) CALL trd_mxl_trc ( kt ) ! trends: Mixed-layer 111 ! 112 IF( nn_dttrc /= 1 ) CALL trc_sub_reset( kt, Kbb, Kmm, Krhs ) ! resetting physical variables when sub-stepping 113 113 ! 114 114 ENDIF … … 117 117 ztrai = 0._wp ! content of all tracers 118 118 DO jn = 1, jptra 119 ztrai = ztrai + glob_sum( 'trcstp', tr n(:,:,:,jn) * cvol(:,:,:) )119 ztrai = ztrai + glob_sum( 'trcstp', tr(:,:,:,jn,Kmm) * cvol(:,:,:) ) 120 120 END DO 121 121 IF( lwm ) WRITE(numstr,9300) kt, ztrai / areatot
Note: See TracChangeset
for help on using the changeset viewer.