Changeset 6223 for branches/NERC
- Timestamp:
- 2016-01-07T16:48:36+01:00 (9 years ago)
- Location:
- branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r6221 r6223 47 47 PUBLIC trc_rst_wri ! called by ??? 48 48 PUBLIC trc_rst_cal 49 PUBLIC trc_rst_stat 49 50 50 51 !! * Substitutions … … 296 297 !! 297 298 !! calculate stats on these fields 299 call trc_rst_dia_stat(qint_idtra(:,:,1), 'qint_IDTRA') 298 300 IF(lwp) WRITE(numout,*) ' IDTRA averaged properties stats (min, max, sum) ...' 299 301 fq0 = MINVAL(qint_idtra(:,:,1)) … … 319 321 !! 320 322 !! calculate stats on these fields 323 call trc_rst_dia_stat(qint_cfc(:,:,jl), 'qint_'//ctrcnm(jn)) 321 324 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...' 322 325 fq0 = MINVAL(qint_cfc(:,:,jl)) … … 457 460 !! calculate stats on these fields 458 461 IF(lwp) WRITE(numout,*) ' IDTRA averaged properties stats (min, max, sum) ...' 462 call trc_rst_dia_stat(qint_idtra(:,:,1), 'qint_IDTRA') 459 463 fq0 = MINVAL(qint_idtra(:,:,1)) 460 464 fq1 = MAXVAL(qint_idtra(:,:,1)) … … 470 474 !! 471 475 !! calculate stats on these fields 476 call trc_rst_dia_stat(qint_cfc(:,:,jl), 'qint_'//ctrcnm(jn)) 472 477 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...' 473 478 fq0 = MINVAL(qint_cfc(:,:,jl)) … … 640 645 IF(lwp) WRITE(numout,9000) jn, TRIM( ctrcnm(jn) ), zmean, zmin, zmax, zdrift 641 646 END DO 642 WRITE(numout,*)647 IF(lwp) WRITE(numout,*) 643 648 9000 FORMAT(' tracer nb :',i2,' name :',a10,' mean :',e18.10,' min :',e18.10, & 644 649 & ' max :',e18.10,' drift :',e18.10, ' %') 645 650 ! 646 651 END SUBROUTINE trc_rst_stat 652 653 654 SUBROUTINE trc_rst_dia_stat 655 !!---------------------------------------------------------------------- 656 !! *** trc_rst_dia_stat *** 657 !! 658 !! ** purpose : Compute tracers statistics 659 !!---------------------------------------------------------------------- 660 REAL(wp), DIMENSION(jpi,jpj) , INTENT(in) :: dgtr ! ocean time-step 661 CHARACTER(len=*) , INTENT(in) :: names ! "READ"/"WRITE" flag 662 !!--------------------------------------------------------------------- 663 INTEGER :: jk, jn 664 REAL(wp) :: ztraf, zmin, zmax, zmean, areasf 665 REAL(wp), DIMENSION(jpi,jpj) :: zvol 666 !!---------------------------------------------------------------------- 667 668 IF( lwp ) WRITE(numout,*) 'STAT- ', names 669 ! 670 zvol(:,:) = e1e2t(:,:) * fse3t_a(:,:,1) * tmask(:,:,1) 671 ztraf = glob_sum( dgtr(:,:) * zvol(:,:) ) 672 areasf = glob_sum(e1e2t(:,:) * tmask(:,:,1) ) 673 zmin = MINVAL( dgtr(:,:), mask= ((tmask(:,:,1)*SPREAD(tmask_i,DIM=2,NCOPIES=1).NE.0.)) ) 674 zmax = MAXVAL( dgtr(:,:), mask= ((tmask(:,:,1)*SPREAD(tmask_i,DIM=2,NCOPIES=1).NE.0.)) ) 675 IF( lk_mpp ) THEN 676 CALL mpp_min( zmin ) ! min over the global domain 677 CALL mpp_max( zmax ) ! max over the global domain 678 END IF 679 zmean = ztraf / areatot 680 IF(lwp) WRITE(numout,9001) TRIM( names ), zmean, zmin, zmax 681 ! 682 IF(lwp) WRITE(numout,*) 683 9001 FORMAT(' tracer name :',a10,' mean :',e18.10,' min :',e18.10, & 684 & ' max :',e18.10 ) 685 ! 686 END SUBROUTINE trc_rst_dia_stat 687 647 688 648 689 #else -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcstp.F90
r5726 r6223 92 92 # endif 93 93 CALL trc_rst_opn ( kt ) ! Open tracer restart file 94 # if defined key_debug_medusa 95 CALL trc_rst_stat 96 # endif 94 97 IF( lrst_trc ) CALL trc_rst_cal ( kt, 'WRITE' ) ! calendar 95 98 IF( lk_iomput ) THEN ; CALL trc_wri ( kt ) ! output of passive tracers with iom I/O manager … … 99 102 # if defined key_debug_medusa 100 103 IF(lwp) WRITE(numout,*) ' MEDUSA trc_stp SMS complete at kt =', kt 104 CALL trc_rst_stat 101 105 CALL flush(numout) 102 106 # endif … … 104 108 # if defined key_debug_medusa 105 109 IF(lwp) WRITE(numout,*) ' MEDUSA trc_stp transport complete at kt =', kt 110 CALL trc_rst_stat 106 111 CALL flush(numout) 107 112 # endif
Note: See TracChangeset
for help on using the changeset viewer.