Changeset 13150


Ignore:
Timestamp:
2020-06-24T14:21:34+02:00 (6 weeks ago)
Author:
jpalmier
Message:

update MEDUSA-external revision

Location:
NEMO/branches/NERC/dev_r4.0.2_NERC_Externals
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/NERC/dev_r4.0.2_NERC_Externals

    • Property svn:externals
      •  

        old new  
        66^/vendors/FCM@10134           ext/FCM 
        77^/vendors/IOIPSL@9655         ext/IOIPSL 
        8 https://code.metoffice.gov.uk/svn/medusa/main/trunk/src/MEDUSA@25 src/TOP/MEDUSA 
        9 https://code.metoffice.gov.uk/svn/medusa/main/trunk/cfgs/ORCA2_MEDUSA@25 cfgs/ORCA2_MEDUSA 
         8https://code.metoffice.gov.uk/svn/medusa/main/trunk/src/MEDUSA@64 src/TOP/MEDUSA 
         9https://code.metoffice.gov.uk/svn/medusa/main/trunk/cfgs/ORCA2_MEDUSA@64 cfgs/ORCA2_MEDUSA 
  • NEMO/branches/NERC/dev_r4.0.2_NERC_Externals/src/TOP/trcrst.F90

    r11536 r13150  
    2424   USE daymod 
    2525   USE lib_mpp 
     26   USE sms_medusa 
    2627    
    2728   IMPLICIT NONE 
     
    3233   PUBLIC   trc_rst_wri       ! called by ??? 
    3334   PUBLIC   trc_rst_cal 
     35   PUBLIC   trc_rst_stat 
     36   PUBLIC   trc_rst_conserve  ! Conservation Checks  
     37 
    3438 
    3539   !!---------------------------------------------------------------------- 
     
    149153      IF( kt == nitrst ) THEN 
    150154          CALL trc_rst_stat            ! statistics 
     155          CALL trc_rst_conserve        ! Conservation Checks  
    151156          CALL iom_close( numrtw )     ! close the restart file (only at last time step) 
    152157#if ! defined key_trdmxl_trc 
     
    338343   END SUBROUTINE trc_rst_stat 
    339344 
     345   SUBROUTINE trc_rst_conserve 
     346      !!---------------------------------------------------------------------- 
     347      !!                    ***  trc_rst_conserve  *** 
     348      !! 
     349      !! ** purpose  :   Compute tracers conservation statistics 
     350      !! 
     351      !! AXY (17/11/2017) 
     352      !! This routine calculates the "now" inventories of the elemental  
     353      !! cycles of MEDUSA and compares them to those calculate when the 
     354      !! model was initialised / restarted; the cycles calculated are: 
     355      !!    nitrogen, silicon, iron, carbon, alkalinity and oxygen 
     356      !!---------------------------------------------------------------------- 
     357      INTEGER  :: ji, jj, jk, jn 
     358      REAL(wp) :: zsum3d, zsum2d, zinvt, zdelta, zratio 
     359      REAL(wp), DIMENSION(jpi,jpj,jpk) :: z3d, zvol 
     360      REAL(wp), DIMENSION(jpi,jpj)     :: z2d, zarea 
     361      REAL(wp), DIMENSION(6)           :: loc_cycletot3, loc_cycletot2 
     362      !!---------------------------------------------------------------------- 
     363      ! 
     364      IF( lwp ) THEN 
     365         WRITE(numout,*)  
     366         WRITE(numout,*) '           ----TRACER CONSERVATION----             ' 
     367         WRITE(numout,*)  
     368      ENDIF 
     369      ! 
     370      ! ocean volume 
     371      DO jk = 1, jpk 
     372         zvol(:,:,jk) = e1e2t(:,:) * e3t_a(:,:,jk) * tmask(:,:,jk) 
     373      END DO 
     374      ! 
     375      ! ocean area (for sediments) 
     376      zarea(:,:)      = e1e2t(:,:) * tmask(:,:,1) 
     377      ! 
     378      !---------------------------------------------------------------------- 
     379      ! nitrogen 
     380      z3d(:,:,:) = trn(:,:,:,jpphn) + trn(:,:,:,jpphd) + trn(:,:,:,jpzmi) + & 
     381                   trn(:,:,:,jpzme) + trn(:,:,:,jpdet) + trn(:,:,:,jpdin) 
     382      z2d(:,:)   = zn_sed_n(:,:) 
     383      zsum3d     = glob_sum( z3d(:,:,:) * zvol(:,:,:) ) 
     384      zsum2d     = glob_sum( z2d(:,:) * zarea(:,:) ) 
     385      ! total tracer, and delta 
     386      zinvt      = zsum3d + zsum2d 
     387      zdelta     = zinvt - cycletot(1) 
     388      zratio     = 1.0e2 * zdelta / cycletot(1) 
     389      ! 
     390      IF ( lwp ) WRITE(numout,9010) 'nitrogen', zsum3d, zsum2d, zinvt,   & 
     391         cycletot(1), zdelta, zratio 
     392      IF ( lwp ) WRITE(numout,*)  
     393      ! 
     394      !---------------------------------------------------------------------- 
     395      ! silicon 
     396      z3d(:,:,:) = trn(:,:,:,jppds) + trn(:,:,:,jpsil) 
     397      z2d(:,:)   = zn_sed_si(:,:) 
     398      zsum3d     = glob_sum( z3d(:,:,:) * zvol(:,:,:) ) 
     399      zsum2d     = glob_sum( z2d(:,:) * zarea(:,:) ) 
     400      ! total tracer, and delta 
     401      zinvt      = zsum3d + zsum2d 
     402      zdelta     = zinvt - cycletot(2) 
     403      zratio     = 1.0e2 * zdelta / cycletot(2) 
     404      ! 
     405      IF ( lwp ) WRITE(numout,9010) 'silicon', zsum3d, zsum2d, zinvt,    & 
     406         cycletot(2), zdelta, zratio 
     407      IF ( lwp ) WRITE(numout,*)  
     408      ! 
     409      !---------------------------------------------------------------------- 
     410      ! iron 
     411      z3d(:,:,:) = ((trn(:,:,:,jpphn) + trn(:,:,:,jpphd) + trn(:,:,:,jpzmi) + & 
     412            trn(:,:,:,jpzme) + trn(:,:,:,jpdet)) * xrfn) + trn(:,:,:,jpfer) 
     413      z2d(:,:)   = zn_sed_fe(:,:) 
     414      zsum3d     = glob_sum( z3d(:,:,:) * zvol(:,:,:) ) 
     415      zsum2d     = glob_sum( z2d(:,:) * zarea(:,:) ) 
     416      ! total tracer, and delta 
     417      zinvt      = zsum3d + zsum2d 
     418      zdelta     = zinvt - cycletot(3) 
     419      zratio     = 1.0e2 * zdelta / cycletot(3) 
     420      ! 
     421      IF ( lwp ) WRITE(numout,9010) 'iron', zsum3d, zsum2d, zinvt,       & 
     422         cycletot(3), zdelta, zratio 
     423      IF ( lwp ) WRITE(numout,*)  
     424      ! 
     425      !---------------------------------------------------------------------- 
     426      ! carbon 
     427      z3d(:,:,:) = (trn(:,:,:,jpphn) * xthetapn)  + (trn(:,:,:,jpphd) * xthetapd)  + & 
     428                   (trn(:,:,:,jpzmi) * xthetazmi) + (trn(:,:,:,jpzme) * xthetazme) + & 
     429                   trn(:,:,:,jpdtc) + trn(:,:,:,jpdic) 
     430      z2d(:,:)   = zn_sed_c(:,:) + zn_sed_ca(:,:) 
     431      zsum3d     = glob_sum( z3d(:,:,:) * zvol(:,:,:) ) 
     432      zsum2d     = glob_sum( z2d(:,:) * zarea(:,:) ) 
     433      ! total tracer, and delta 
     434      zinvt      = zsum3d + zsum2d 
     435      zdelta     = zinvt - cycletot(4) 
     436      zratio     = 1.0e2 * zdelta / cycletot(4) 
     437      ! 
     438      IF ( lwp ) WRITE(numout,9010) 'carbon', zsum3d, zsum2d, zinvt,     & 
     439         cycletot(4), zdelta, zratio 
     440      IF ( lwp ) WRITE(numout,*)  
     441      ! 
     442      !---------------------------------------------------------------------- 
     443      ! alkalinity 
     444      z3d(:,:,:) = trn(:,:,:,jpalk) 
     445      z2d(:,:)   = zn_sed_ca(:,:) * 2.0 
     446      zsum3d     = glob_sum( z3d(:,:,:) * zvol(:,:,:) ) 
     447      zsum2d     = glob_sum( z2d(:,:) * zarea(:,:) ) 
     448      ! total tracer, and delta 
     449      zinvt      = zsum3d + zsum2d 
     450      zdelta     = zinvt - cycletot(5) 
     451      zratio     = 1.0e2 * zdelta / cycletot(5) 
     452      ! 
     453      IF ( lwp ) WRITE(numout,9010) 'alkalinity', zsum3d, zsum2d, zinvt, & 
     454         cycletot(5), zdelta, zratio 
     455      IF ( lwp ) WRITE(numout,*)  
     456      ! 
     457      !---------------------------------------------------------------------- 
     458      ! oxygen 
     459      z3d(:,:,:) = trn(:,:,:,jpoxy) 
     460      z2d(:,:)   = 0.0 
     461      zsum3d     = glob_sum( z3d(:,:,:) * zvol(:,:,:) ) 
     462      zsum2d     = glob_sum( z2d(:,:) * zarea(:,:) ) 
     463      ! total tracer, and delta 
     464      zinvt      = zsum3d + zsum2d 
     465      zdelta     = zinvt - cycletot(6) 
     466      zratio     = 1.0e2 * zdelta / cycletot(6) 
     467      ! 
     468      IF ( lwp ) WRITE(numout,9010) 'oxygen', zsum3d, zsum2d, zinvt,     & 
     469         cycletot(6), zdelta, zratio 
     470      ! 
     471      !---------------------------------------------------------------------- 
     472      ! Check  
     473      zsum3d        = glob_sum( zvol(:,:,:) ) 
     474      zsum2d        = glob_sum( zarea(:,:) ) 
     475      IF ( lwp ) THEN  
     476         WRITE(numout,*) 
     477         WRITE(numout,*) ' check : cvol    : ', zsum3d 
     478         WRITE(numout,*) ' check : carea   : ', zsum2d 
     479         WRITE(numout,*) 
     480      ENDIF 
     481      ! 
     4829010  FORMAT(' element:',a10,                     & 
     483             ' 3d sum:',e18.10,' 2d sum:',e18.10, & 
     484             ' total:',e18.10,' initial:',e18.10, & 
     485             ' delta:',e18.10,' %:',e18.10) 
     486      ! 
     487   END SUBROUTINE trc_rst_conserve  
     488 
    340489#else 
    341490   !!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.