Changeset 6234 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2016-01-13T14:36:18+01:00 (9 years ago)
- Location:
- branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90
r5729 r6234 65 65 ! 66 66 CALL trc_sbc( kstp ) ! surface boundary condition 67 # if defined key_debug_medusa 68 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_sbc at kt =', kt 69 CALL trc_rst_tra_stat 70 CALL flush(numout) 71 # endif 67 72 IF( lk_trabbl ) CALL trc_bbl( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme 68 73 IF( ln_trcdmp ) CALL trc_dmp( kstp ) ! internal damping trends 69 74 IF( ln_trcdmp_clo ) CALL trc_dmp_clo( kstp ) ! internal damping trends on closed seas only 70 75 CALL trc_adv( kstp ) ! horizontal & vertical advection 76 # if defined key_debug_medusa 77 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_adv at kt =', kt 78 CALL trc_rst_tra_stat 79 CALL flush(numout) 80 # endif 71 81 CALL trc_ldf( kstp ) ! lateral mixing 72 82 IF( .NOT. lk_offline .AND. lk_zdfkpp ) & … … 76 86 #endif 77 87 CALL trc_zdf( kstp ) ! vertical mixing and after tracer fields 88 # if defined key_debug_medusa 89 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_zdf at kt =', kt 90 CALL trc_rst_tra_stat 91 CALL flush(numout) 92 # endif 78 93 CALL trc_nxt( kstp ) ! tracer fields at next time step 94 # if defined key_debug_medusa 95 IF(lwp) WRITE(numout,*) ' MEDUSA trc_trp after trc_nxt at kt =', kt 96 CALL trc_rst_tra_stat 97 CALL flush(numout) 98 # endif 79 99 IF( ln_trcrad ) CALL trc_rad( kstp ) ! Correct artificial negative concentrations 80 100 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r6233 r6234 49 49 PUBLIC trc_rst_stat 50 50 PUBLIC trc_rst_dia_stat 51 PUBLIC trc_rst_tra_stat 51 52 52 53 !! * Substitutions … … 569 570 570 571 572 SUBROUTINE trc_rst_tra_stat 573 !!---------------------------------------------------------------------- 574 !! *** trc_rst_tra_stat *** 575 !! 576 !! ** purpose : Compute tracers statistics - check where crazy values appears 577 !!---------------------------------------------------------------------- 578 INTEGER :: jk, jn 579 REAL(wp) :: ztraf, zmin, zmax, zmean, zdrift 580 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zvol 581 !!---------------------------------------------------------------------- 582 583 IF( lwp ) THEN 584 WRITE(numout,*) 585 WRITE(numout,*) ' ----SURFACE TRA STAT---- ' 586 WRITE(numout,*) 587 ENDIF 588 ! 589 DO jk = 1, jpk 590 zvol(:,:,jk) = e1e2t(:,:) * fse3t_a(:,:,jk) * tmask(:,:,jk) 591 END DO 592 ! 593 DO jn = 1, jptra 594 ztraf = glob_sum( tra(:,:,1,jn) * zvol(:,:,:) ) 595 zmin = MINVAL( tra(:,:,1,jn), mask= ((tmask(:,:,1).NE.0.)) ) 596 zmax = MAXVAL( tra(:,:,1,jn), mask= ((tmask(:,:,1).NE.0.)) ) 597 IF( lk_mpp ) THEN 598 CALL mpp_min( zmin ) ! min over the global domain 599 CALL mpp_max( zmax ) ! max over the global domain 600 END IF 601 zmean = ztraf / areatot 602 IF(lwp) WRITE(numout,9001) jn, TRIM( ctrcnm(jn) ), zmean, zmin, zmax 603 END DO 604 IF(lwp) WRITE(numout,*) 605 9001 FORMAT(' tracer nb :',i2,' name :',a10,' mean :',e18.10,' min :',e18.10, & 606 & ' max :',e18.10) 607 ! 608 END SUBROUTINE trc_rst_tra_stat 609 610 611 571 612 SUBROUTINE trc_rst_dia_stat( dgtr, names) 572 613 !!---------------------------------------------------------------------- … … 595 636 END IF 596 637 zmean = ztraf / areatot 597 IF(lwp) WRITE(numout,900 1) TRIM( names ), zmean, zmin, zmax638 IF(lwp) WRITE(numout,9002) TRIM( names ), zmean, zmin, zmax 598 639 ! 599 640 IF(lwp) WRITE(numout,*) 600 900 1FORMAT(' tracer name :',a10,' mean :',e18.10,' min :',e18.10, &641 9002 FORMAT(' tracer name :',a10,' mean :',e18.10,' min :',e18.10, & 601 642 & ' max :',e18.10 ) 602 643 ! -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcstp.F90
r6223 r6234 94 94 # if defined key_debug_medusa 95 95 CALL trc_rst_stat 96 CALL trc_rst_tra_stat 96 97 # endif 97 98 IF( lrst_trc ) CALL trc_rst_cal ( kt, 'WRITE' ) ! calendar … … 103 104 IF(lwp) WRITE(numout,*) ' MEDUSA trc_stp SMS complete at kt =', kt 104 105 CALL trc_rst_stat 106 CALL trc_rst_tra_stat 105 107 CALL flush(numout) 106 108 # endif … … 109 111 IF(lwp) WRITE(numout,*) ' MEDUSA trc_stp transport complete at kt =', kt 110 112 CALL trc_rst_stat 113 CALL trc_rst_tra_stat 111 114 CALL flush(numout) 112 115 # endif
Note: See TracChangeset
for help on using the changeset viewer.