New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7894 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 – NEMO

Ignore:
Timestamp:
2017-04-11T13:58:38+02:00 (7 years ago)
Author:
jpalmier
Message:

JPALM -- 11-04-2017 -- MEDUSA spring tidy-up refreshning session

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

    r7766 r7894  
    1212   !!              -   !  2015-11  (J. Palmieri) added iom_use for diags 
    1313   !!              -   !  2016-11  (A. Yool) updated diags for CMIP6 
     14   !!              -   !  2017-03  (A. Yool) minor tidy 
    1415   !!---------------------------------------------------------------------- 
    1516#if defined key_medusa 
     
    6768      ! 
    6869      NAMELIST/natbio/ xxi,xaln,xald,jphy,xvpn,xvpd,          & 
    69       &    xsin0,xnsi0,xuif,jliebig, jq10,                    & 
     70      &    xsin0,xnsi0,xuif,jliebig, xq10,                    & 
    7071      &    xthetam,xthetamd,xnln,xnld,xsld,xfln,xfld,         & 
    7172      &  xgmi,xgme,xkmi,xkme,xphi,xbetan,xbetac,xkc,          & 
     
    136137# endif 
    137138 
    138      IF( ( .NOT.lk_iomput .AND. ln_diatrc ) .OR. ( ln_diatrc .AND. lk_medusa ) ) THEN 
     139      !! AXY (06/04/17): can this next section be deleted? I think it can 
     140      IF( ( .NOT.lk_iomput .AND. ln_diatrc ) .OR. ( ln_diatrc .AND. lk_medusa ) ) THEN 
    139141         ! 
    140142         ! Namelist nammeddia 
     
    221223      ! ------------------------------------------- 
    222224       
    223       xxi         = 0. 
    224       xaln        = 0. 
    225       xald        = 0. 
    226       jphy        = 0 
    227       xvpn        = 0. 
    228       xvpd        = 0. 
    229       xthetam     = 0. 
    230       xthetamd    = 0. 
    231 !! 
    232       xsin0       = 0. 
    233       xnsi0       = 0. 
    234       xuif        = 0. 
    235 !! 
    236       jliebig     = 0 
    237       jq10        = 0. 
    238       xnln        = 0. 
    239       xnld        = 0. 
    240       xsld        = 0. 
    241       xfln        = 0. 
    242       xfld        = 0. 
    243 !! 
    244       xgmi        = 0. 
    245       xgme        = 0. 
    246       xkmi        = 0. 
    247       xkme        = 0. 
    248       xphi    = 0. 
    249       xbetan      = 0. 
    250       xbetac      = 0. 
    251       xkc         = 0. 
    252       xpmipn      = 0. 
    253       xpmid       = 0. 
    254       xpmepn      = 0. 
    255       xpmepd      = 0. 
    256       xpmezmi     = 0. 
    257       xpmed       = 0. 
    258 !! 
    259       xmetapn     = 0. 
    260       xmetapd     = 0. 
    261       xmetazmi    = 0. 
    262       xmetazme    = 0. 
    263 !! 
    264       jmpn        = 0 
    265       xmpn        = 0. 
    266       xkphn       = 0. 
    267       jmpd        = 0 
    268       xmpd        = 0. 
    269       xkphd       = 0. 
    270       jmzmi       = 0 
    271       xmzmi       = 0. 
    272       xkzmi       = 0. 
    273       jmzme       = 0 
    274       xmzme       = 0. 
    275       xkzme       = 0. 
    276 !! 
    277       jmd         = 0 
    278       jsfd        = 0 
    279       xmd         = 0. 
    280       xmdc        = 0. 
    281 !! 
    282       xthetapn    = 0. 
    283       xthetapd    = 0. 
    284       xthetazmi   = 0. 
    285       xthetazme   = 0. 
    286       xthetad     = 0. 
    287       xrfn        = 0. 
    288       xrsn        = 0.  !: (NOT USED HERE; RETAINED FOR LOBSTER) 
    289 !! 
    290       jiron       = 0 
    291       xfe_mass    = 0. 
    292       xfe_sol     = 0. 
    293       xfe_sed     = 0. 
    294       xLgT        = 0. 
    295       xk_FeL     = 0. 
    296       xk_sc_Fe    = 0. 
    297 !! 
    298       jexport     = 0 
    299       jfdfate     = 0 
    300       jrratio     = 0 
    301       jocalccd    = 0 
    302       xridg_r0    = 0. 
    303       xfdfrac1   = 0. 
    304       xfdfrac2   = 0. 
    305       xfdfrac3   = 0. 
    306       xcaco3a    = 0. 
    307       xcaco3b    = 0. 
    308       xmassc     = 0. 
    309       xmassca    = 0. 
    310       xmasssi    = 0. 
    311       xprotca    = 0. 
    312       xprotsi    = 0. 
    313       xfastc     = 0. 
    314       xfastca    = 0. 
    315       xfastsi    = 0. 
    316 !! 
    317       jorgben     = 0 
    318       jinorgben   = 0 
    319       xsedn       = 0. 
    320       xsedfe      = 0. 
    321       xsedsi      = 0. 
    322       xsedc       = 0. 
    323       xsedca      = 0. 
    324       xburial     = 0. 
    325 !! 
    326       jriver_n    = 0 
    327       jriver_si   = 0 
    328       jriver_c    = 0 
    329       jriver_alk  = 0 
    330       jriver_dep  = 1 
    331 !! 
    332       xsdiss     = 0. 
    333 !! 
    334       vsed        = 0. 
    335       xhr         = 0. 
    336 !! 
    337       sedlam     = 0. 
    338       sedlostpoc  = 0. 
    339       jpkb    = 0. 
    340       jdms        = 0 
    341       jdms_input  = 0 
    342       jdms_model  = 0 
     225!! Primary production and chl related quantities 
     226      xxi         = 0.  !! conversion factor from gC to mmolN 
     227      xaln        = 0.  !! Chl-a specific initial slope of P-I curve for non-diatoms 
     228      xald        = 0.  !! Chl-a specific initial slope of P-I curve for diatoms 
     229      jphy        = 0   !! phytoplankton T-dependent growth switch 
     230      xvpn        = 0.  !! maximum growth rate for non-diatoms 
     231      xvpd        = 0.  !! maximum growth rate for diatoms 
     232      xthetam     = 0.  !! maximum Chl to C ratio for non-diatoms 
     233      xthetamd    = 0.  !! maximum Chl to C ratio for diatoms 
     234!! 
     235!! Diatom silicon parameters 
     236      xsin0       = 0.  !! minimum diatom Si:N ratio 
     237      xnsi0       = 0.  !! minimum diatom N:Si ratio 
     238      xuif        = 0.  !! hypothetical growth ratio at infinite Si:N ratio 
     239!! 
     240!! Nutrient limitation 
     241      jliebig     = 0   !! Liebig nutrient uptake switch 
     242      xq10        = 0.  !! specific Q10 value (if jphy = 2)  
     243      xnln        = 0.  !! half-sat constant for DIN uptake by non-diatoms 
     244      xnld        = 0.  !! half-sat constant for DIN uptake by diatoms 
     245      xsld        = 0.  !! half-sat constant for  Si uptake by diatoms 
     246      xfln        = 0.  !! half-sat constant for  Fe uptake by non-datoms 
     247      xfld        = 0.  !! half-sat constant for  Fe uptake by diatoms 
     248!! 
     249!! Grazing 
     250      xgmi        = 0.  !! microzoo maximum growth rate 
     251      xgme        = 0.  !! mesozoo maximum growth rate 
     252      xkmi        = 0.  !! microzoo grazing half-sat parameter 
     253      xkme        = 0.  !! mesozoo grazing half-sat parameter 
     254      xphi    = 0.  !! micro/mesozoo grazing inefficiency 
     255      xbetan      = 0.  !! micro/mesozoo N assimilation efficiency 
     256      xbetac      = 0.  !! micro/mesozoo C assimilation efficiency 
     257      xkc         = 0.  !! micro/mesozoo net C growth efficiency 
     258      xpmipn      = 0.  !! grazing preference of microzoo for non-diatoms 
     259      xpmid       = 0.  !! grazing preference of microzoo for diatoms 
     260      xpmepn      = 0.  !! grazing preference of mesozoo for non-diatoms 
     261      xpmepd      = 0.  !! grazing preference of mesozoo for diatoms 
     262      xpmezmi     = 0.  !! grazing preference of mesozoo for microzoo 
     263      xpmed       = 0.  !! grazing preference of mesozoo for detritus 
     264!! 
     265!! Metabolic losses 
     266      xmetapn     = 0.  !! non-diatom metabolic loss rate 
     267      xmetapd     = 0.  !! diatom     metabolic loss rate 
     268      xmetazmi    = 0.  !! microzoo   metabolic loss rate 
     269      xmetazme    = 0.  !! mesozoo    metabolic loss rate 
     270!! 
     271!! Mortality/Remineralisation 
     272      jmpn        = 0   !! non-diatom mortality functional form 
     273      xmpn        = 0.  !! non-diatom mortality rate 
     274      xkphn       = 0.  !! non-diatom mortality half-sat constant 
     275      jmpd        = 0   !! diatom     mortality functional form 
     276      xmpd        = 0.  !! diatom     mortality rate 
     277      xkphd       = 0.  !! diatom     mortality half-sat constant 
     278      jmzmi       = 0   !! microzoo   mortality functional form 
     279      xmzmi       = 0.  !! microzoo   mortality rate 
     280      xkzmi       = 0.  !! microzoo   mortality half-sat constant 
     281      jmzme       = 0   !! mesozoo    mortality functional form 
     282      xmzme       = 0.  !! mesozoo    mortality rate 
     283      xkzme       = 0.  !! mesozoo    mortality half-sat constant 
     284!! 
     285!! Remineralisation 
     286      jmd         = 0   !! detritus T-dependent remineralisation switch 
     287      jsfd        = 0   !! accelerate seafloor detritus remin. switch 
     288      xmd         = 0.  !! detrital nitrogen remineralisation rate 
     289      xmdc        = 0.  !! detrital carbon remineralisation rate 
     290!! 
     291!! Stochiometric ratios 
     292      xthetapn    = 0.  !! non-diatom C:N ratio 
     293      xthetapd    = 0.  !! diatom     C:N ratio 
     294      xthetazmi   = 0.  !! microzoo   C:N ratio 
     295      xthetazme   = 0.  !! mesozoo    C:N ratio 
     296      xthetad     = 0.  !! detritus   C:N ratio (MEDUSA-1 only) 
     297      xrfn        = 0.  !! phytoplankton Fe:N ratio 
     298      xrsn        = 0.  !! diatom Si:N ratio (NOT USED HERE; LOBSTER ARTIFACT) 
     299!! 
     300!! Iron parameters 
     301      jiron       = 0   !! iron scavenging submodel switch 
     302      xfe_mass    = 0.  !! iron atomic mass 
     303      xfe_sol     = 0.  !! aeolian iron solubility 
     304      xfe_sed     = 0.  !! sediment iron input 
     305      xLgT        = 0.  !! total ligand concentration (umol/m3) 
     306      xk_FeL     = 0.  !! dissociation constant for (Fe + L) 
     307      xk_sc_Fe    = 0.  !! scavenging rate of "free" iron 
     308!! 
     309!! Fast-sinking detritus parameters 
     310      jexport     = 0   !! fast detritus remineralisation switch 
     311      jfdfate     = 0   !! fate of fast detritus at seafloor switch 
     312      jrratio     = 0   !! rain ratio switch 
     313      jocalccd    = 0   !! CCD switch 
     314      xridg_r0    = 0.  !! Ridgwell rain ratio coefficient 
     315      xfdfrac1   = 0.  !! fast-sinking fraction of diatom nat. mort. losses 
     316      xfdfrac2   = 0.  !! fast-sinking fraction of meszooplankton mort. losses 
     317      xfdfrac3   = 0.  !! fast-sinking fraction of diatom silicon grazing losses 
     318      xcaco3a    = 0.  !! polar (high latitude) CaCO3 fraction 
     319      xcaco3b    = 0.  !! equatorial (low latitude) CaCO3 fraction 
     320      xmassc     = 0.  !! organic C mass:mole ratio, C106 H175 O40 N16 P1 
     321      xmassca    = 0.  !! calcium carbonate mass:mole ratio, CaCO3 
     322      xmasssi    = 0.  !! biogenic silicon mass:mole ratio, (H2SiO3)n 
     323      xprotca    = 0.  !! calcium carbonate protection ratio 
     324      xprotsi    = 0.  !! biogenic silicon protection ratio 
     325      xfastc     = 0.  !! organic C remineralisation length scale 
     326      xfastca    = 0.  !! calcium carbonate dissolution length scale 
     327      xfastsi    = 0.  !! biogenic silicon dissolution length scale 
     328!! 
     329!! Benthic 
     330      jorgben     = 0   !! does   organic detritus go to the benthos? 
     331      jinorgben   = 0   !! does inorganic detritus go to the benthos? 
     332      xsedn       = 0.  !! organic   nitrogen sediment remineralisation rate 
     333      xsedfe      = 0.  !! organic   iron     sediment remineralisation rate 
     334      xsedsi      = 0.  !! inorganic silicon  sediment dissolution      rate 
     335      xsedc       = 0.  !! organic   carbon   sediment remineralisation rate 
     336      xsedca      = 0.  !! inorganic carbon   sediment dissolution      rate 
     337      xburial     = 0.  !! burial rate of seafloor detritus 
     338!! 
     339!! Riverine inputs 
     340      jriver_n    = 0   !! riverine N          input? 
     341      jriver_si   = 0   !! riverine Si         input? 
     342      jriver_c    = 0   !! riverine C          input? 
     343      jriver_alk  = 0   !! riverine alkalinity input? 
     344      jriver_dep  = 1   !! depth of riverine   input? 
     345!! 
     346!! Miscellaneous 
     347      xsdiss     = 0.  !! diatom frustule dissolution rate 
     348!! 
     349!! Gravitational sinking 
     350      vsed        = 0.  !! detritus gravitational sinking rate 
     351      xhr         = 0.  !! coeff for Martin's remineralisation profile (LOBSTER) 
     352!! 
     353!! Additional parameters 
     354      sedlam     = 0.  !! time coeff of POC in sediments (LOBSTER) 
     355      sedlostpoc  = 0.  !! sediment geol loss for POC (LOBSTER) 
     356      jpkb    = 0.  !! vertical layer for diagnostic of the vertical flux (LOBSTER) 
     357!! 
     358!! UKESM1 - new DMS parameters  !! Jpalm 
     359      jdms        = 0   !! include dms diagnostics 
     360      jdms_input  = 0   !! use instant (0) or diel-avg (1) inputs 
     361      jdms_model  = 0   !! choice of DMS model passed to atmosphere 
     362                        !! 1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 
    343363             
    344364      !REWIND(numnatm) 
    345365      !READ(numnatm,natbio) 
    346          ! Namelist natbio 
    347          ! ------------------- 
    348          REWIND( numnatp_ref )              ! Namelist natbio in reference namelist : MEDUSA diagnostics 
    349          READ  ( numnatp_ref, natbio, IOSTAT = ios, ERR = 903) 
    350 903      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in reference namelist', lwp ) 
    351  
    352          REWIND( numnatp_cfg )              ! Namelist natbio in configuration namelist : MEDUSA diagnostics 
    353          READ  ( numnatp_cfg, natbio, IOSTAT = ios, ERR = 904 ) 
    354 904      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in configuration namelist', lwp ) 
    355          IF(lwm) WRITE ( numonp, natbio ) 
    356  
    357 !! Primary production and chl related quantities 
    358 !!       xxi         :  conversion factor from gC to mmolN  
    359 !!       xaln        :  Chl-a specific initial slope of P-I curve for non-diatoms  
    360 !!       xald        :  Chl-a specific initial slope of P-I curve for diatoms 
    361 !!       jphy        :  phytoplankton T-dependent growth switch 
    362 !!       xvpn        :  maximum growth rate for non-diatoms 
    363 !!       xvpd        :  maximum growth rate for diatoms 
    364 !!       xthetam     :  maximum Chl to C ratio for non-diatoms       
    365 !!       xthetamd    :  maximum Chl to C ratio for diatoms       
    366 !! 
    367 !! Diatom silicon parameters 
    368 !!       xsin0       :  minimum diatom Si:N ratio 
    369 !!       xnsi0       :  minimum diatom N:Si ratio 
    370 !!       xuif        :  hypothetical growth ratio at infinite Si:N ratio 
    371 !! 
    372 !! Nutrient limitation 
    373 !!       jliebig     :  Liebig nutrient uptake switch 
    374 !!       xnln        :  half-sat constant for DIN uptake by non-diatoms  
    375 !!       xnld        :  half-sat constant for DIN uptake by diatoms  
    376 !!       xsl         :  half-sat constant for Si uptake by diatoms  
    377 !!       xfld        :  half-sat constant for Fe uptake by diatoms   
    378 !!       xfln        :  half-sat constant for Fe uptake by non-datoms  
    379 !! 
    380 !! Grazing 
    381 !!       xgmi        :  microzoo maximum growth rate  
    382 !!       xgme        :  mesozoo maximum growth rate  
    383 !!       xkmi        :  microzoo grazing half-sat parameter 
    384 !!       xkme        :  mesozoo grazing half-sat parameter 
    385 !!       xphi        :  micro/mesozoo grazing inefficiency 
    386 !!       xbetan      :  micro/mesozoo N assimilation efficiency 
    387 !!       xbetac      :  micro/mesozoo C assimilation efficiency 
    388 !!       xkc         :  micro/mesozoo net C growth efficiency 
    389 !!       xpmipn      :  grazing preference of microzoo for non-diatoms 
    390 !!       xpmid       :  grazing preference of microzoo for diatoms 
    391 !!       xpmepn      :  grazing preference of mesozoo for non-diatoms  
    392 !!       xpmepd      :  grazing preference of mesozoo for diatoms 
    393 !!       xpmezmi     :  grazing preference of mesozoo for microzoo 
    394 !!       xpmed       :  grazing preference of mesozoo for detritus 
    395 !! 
    396 !! Metabolic losses 
    397 !!       xmetapn     :  non-diatom metabolic loss rate 
    398 !!       xmetapd     :  diatom     metabolic loss rate 
    399 !!       xmetazmi    :  microzoo   metabolic loss rate 
    400 !!       xmetazme    :  mesozoo    metabolic loss rate 
    401 !! 
    402 !! Mortality/Remineralisation 
    403 !!       jmpn        :  non-diatom mortality functional form 
    404 !!       xmpn        :  non-diatom mortality rate 
    405 !!       xkphn       :  non-diatom mortality half-sat constant 
    406 !!       jmpd        :  diatom     mortality functional form 
    407 !!       xmpd        :  diatom mortality rate 
    408 !!       xkphd       :  diatom mortality half-sat constant 
    409 !!       jmzmi       :  microzoo   mortality functional form 
    410 !!       xmzmi       :  microzoo mortality rate 
    411 !!       xkzmi       :  microzoo mortality half-sat constant 
    412 !!       jmzme       :  mesozoo    mortality functional form 
    413 !!       xmzme       :  mesozoo mortality rate 
    414 !!       xkzme       :  mesozoo mortality half-sat constant 
    415 !! 
    416 !! Remineralisation 
    417 !!       jmd         :  detritus T-dependent remineralisation switch 
    418 !!       jsfd        :  accelerate seafloor detritus remin. switch 
    419 !!       xmd         :  detrital nitrogen remineralisation rate  
    420 !!       xmdc        :  detrital carbon remineralisation rate  
    421 !! 
    422 !! Stochiometric ratios 
    423 !!       xthetapn    :  non-diatom C:N ratio 
    424 !!       xthetapd    :  diatom C:N ratio 
    425 !!       xthetazmi   :  microzoo C:N ratio 
    426 !!       xthetazme   :  mesozoo C:N ratio 
    427 !!       xthetad     :  detritus C:N ratio 
    428 !!       xrfn        :  phytoplankton Fe:N ratio 
    429 !!  xrsn        :  diatom Si:N ratio (*NOT* used) 
    430 !! 
    431 !! Iron parameters 
    432 !!       jiron       :  iron scavenging submodel switch 
    433 !!       xfe_mass    :  iron atomic mass 
    434 !!  xfe_sol     :  aeolian iron solubility 
    435 !!  xfe_sed     :  sediment iron input 
    436 !!  xLgT      :  total ligand concentration (umol/m3) 
    437 !!  xk_FeL       :  dissociation constant for (Fe + L) 
    438 !!  xk_sc_Fe    :  scavenging rate of "free" iron 
    439 !!   
    440 !! Fast-sinking detritus parameters 
    441 !!       jexport     :  fast detritus remineralisation switch 
    442 !!       jfdfate     :  fate of fast detritus at seafloor switch 
    443 !!       jrratio     :  rain ratio switch 
    444 !!       jocalccd    :  CCD switch 
    445 !!       xridg_r0    :  Ridgwell rain ratio coefficient 
    446 !!       xfdfrac1    :  fast-sinking fraction of diatom nat. mort. losses 
    447 !!       xfdfrac2    :  fast-sinking fraction of meszooplankton mort. losses 
    448 !!       xfdfrac3    :  fast-sinking fraction of diatom silicon grazing losses 
    449 !!       xcaco3a     :  polar (high latitude) CaCO3 fraction 
    450 !!       xcaco3b     :  equatorial (low latitude) CaCO3 fraction 
    451 !!       xmassc      :  organic C mass:mole ratio, C106 H175 O40 N16 P1 
    452 !!       xmassca     :  calcium carbonate mass:mole ratio, CaCO3 
    453 !!       xmasssi     :  biogenic silicon mass:mole ratio, (H2SiO3)n 
    454 !!       xprotca     :  calcium carbonate protection ratio 
    455 !!       xprotsi     :  biogenic silicon protection ratio 
    456 !!       xfastc      :  organic C remineralisation length scale 
    457 !!       xfastca     :  calcium carbonate dissolution length scale 
    458 !!       xfastsi     :  biogenic silicon dissolution length scale 
    459 !! 
    460 !! Benthic  
    461 !!       jorgben     :  does   organic detritus go to the benthos? 
    462 !!       jinorgben   :  does inorganic detritus go to the benthos? 
    463 !!       xsedn       :  organic   nitrogen sediment remineralisation rate  
    464 !!       xsedfe      :  organic   iron     sediment remineralisation rate  
    465 !!       xsedsi      :  inorganic silicon  sediment dissolution      rate  
    466 !!       xsedc       :  organic   carbon   sediment remineralisation rate  
    467 !!       xsedca      :  inorganic carbon   sediment dissolution      rate  
    468 !!       xburial     :  burial rate of seafloor detritus 
    469 !! 
    470 !! Riverine inputs 
    471 !!       jriver_n    :  riverine N          input? 
    472 !!       jriver_si   :  riverine Si         input? 
    473 !!       jriver_c    :  riverine C          input? 
    474 !!       jriver_alk  :  riverine alkalinity input? 
    475 !!       jriver_dep  :  depth of riverine   input? 
    476 !! 
    477 !! Miscellaneous 
    478 !!       xsdiss      :  diatom frustule dissolution rate 
    479 !! 
    480 !! Gravitational sinking       
    481 !!       vsed        :  detritus gravitational sinking rate  
    482 !!       xhr         :  coeff for Martin's remineralisation profile 
    483 !! 
    484 !! Additional parameters 
    485 !!       sedlam      :  time coeff of POC in sediments 
    486 !!      sedlostpoc   :  sediment geol loss for POC 
    487 !!       jpkb        :  vertical layer for diagnostic of the vertical flux  
    488 !!                      NOTE that in LOBSTER it is a first vertical layers where  
    489 !!                      biology is active   
    490 !! 
    491 !! UKESM1 - new diagnostics  !! Jpalm 
    492 !!       jdms        :  include dms diagnostics 
    493 !!  jdms_input  :  use instant (0) or diel-avg (1) inputs 
    494 !!       jdms_model  :  choice of DMS model passed to atmosphere 
    495 !!                      1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 
    496 !! 
     366      ! Namelist natbio 
     367      ! ------------------- 
     368      REWIND( numnatp_ref )     ! Namelist natbio in reference namelist : MEDUSA diagnostics 
     369      READ  ( numnatp_ref, natbio, IOSTAT = ios, ERR = 903) 
     370 903  IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in reference namelist', lwp ) 
     371 
     372      REWIND( numnatp_cfg )     ! Namelist natbio in configuration namelist : MEDUSA diagnostics 
     373      READ  ( numnatp_cfg, natbio, IOSTAT = ios, ERR = 904 ) 
     374 904  IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in configuration namelist', lwp ) 
     375      IF(lwm) WRITE ( numonp, natbio ) 
     376 
    497377      IF(lwp) THEN 
    498 !! 
     378 
    499379!! AXY (08/11/13): compilation key notification 
    500380         WRITE(numout,*) '=== Compilation keys' 
     
    512392         WRITE(numout,*)     & 
    513393         &   ' key_axy_carbchem                                                       = INACTIVE' 
    514 #endif         
    515 #if defined key_mocsy 
    516          WRITE(numout,*)     & 
    517          &   ' key_mocsy                                                              = ACTIVE' 
    518 #else 
    519          WRITE(numout,*)     & 
    520          &   ' key_mocsy                                                              = INACTIVE' 
    521394#endif         
    522395#if defined key_avgqsr_medusa 
     
    601474         &   ' maximum Chl to C ratio for diatoms,                        xthetamd    = ', xthetamd 
    602475         WRITE(numout,*)     & 
    603          &   ' specific Q10 value (jphy==2),                                  jq10    = ', jq10 
     476         &   ' specific Q10 value (jphy==2),                                  xq10    = ', xq10 
    604477!! 
    605478!! Diatom silicon parameters 
     
    791664            WRITE(numout,*)     & 
    792665            &   ' Dutkiewicz et al. (2005) iron scavenging                   jiron       = ', jiron 
    793          elseif (jiron.eq.2) then 
    794             WRITE(numout,*)     & 
    795             &   ' Moore et al. (2004) iron scavenging                        jiron       = ', jiron 
    796          elseif (jiron.eq.3) then 
    797             WRITE(numout,*)     & 
    798             &   ' Moore et al. (2008) iron scavenging                        jiron       = ', jiron 
    799          elseif (jiron.eq.4) then 
    800             WRITE(numout,*)     & 
    801             &   ' Galbraith et al. (2010) iron scavenging                    jiron       = ', jiron 
    802666         else 
    803667            WRITE(numout,*)     & 
     
    1004868!! 
    1005869!! UKESM1 - new diagnostics  !! Jpalm; AXY (08/07/15) 
    1006          WRITE(numout,*) '=== UKESM1-related parameters' 
     870         WRITE(numout,*) '=== UKESM1 DMS parameters' 
    1007871         WRITE(numout,*)     & 
    1008872         &   ' include DMS diagnostic?,                                   jdms        = ', jdms 
     
    1084948 
    1085949      !READ(numnatm,natroam) 
    1086          ! Namelist natroam 
    1087          ! ------------------- 
    1088          REWIND( numnatp_ref )              ! Namelist natroam in reference namelist : MEDUSA diagnostics 
    1089          READ  ( numnatp_ref, natroam, IOSTAT = ios, ERR = 905) 
    1090 905      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in reference namelist', lwp ) 
    1091  
    1092          REWIND( numnatp_cfg )              ! Namelist natroam in configuration namelist : MEDUSA diagnostics 
    1093          READ  ( numnatp_cfg, natroam, IOSTAT = ios, ERR = 906 ) 
    1094 906      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in configuration namelist', lwp ) 
    1095          IF(lwm) WRITE ( numonp, natroam ) 
     950      ! Namelist natroam 
     951      ! ------------------- 
     952      REWIND( numnatp_ref )     ! Namelist natroam in reference namelist : MEDUSA diagnostics 
     953      READ  ( numnatp_ref, natroam, IOSTAT = ios, ERR = 905) 
     954 905  IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in reference namelist', lwp ) 
     955 
     956      REWIND( numnatp_cfg )     ! Namelist natroam in configuration namelist : MEDUSA diagnostics 
     957      READ  ( numnatp_cfg, natroam, IOSTAT = ios, ERR = 906 ) 
     958 906  IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in configuration namelist', lwp ) 
     959      IF(lwm) WRITE ( numonp, natroam ) 
    1096960 
    1097961!! ROAM carbon, alkalinity and oxygen cycle parameters 
     
    1101965!!       xthetarem :  oxygen consumption by carbon remineralisation 
    1102966!!       xo2min    :  oxygen minimum concentration 
    1103  
    1104967      IF(lwp) THEN 
    1105968          WRITE(numout,*) 'natroam' 
    1106969          WRITE(numout,*) ' ' 
    1107 !! 
    1108 !! ROAM carbon, alkalinity and oxygen cycle parameters 
    1109970          WRITE(numout,*) '=== ROAM carbon, alkalinity and oxygen cycle parameters' 
    1110971          WRITE(numout,*)     & 
     
    1136997 
    1137998      !READ(numnatm,natopt) 
    1138          ! Namelist natopt 
    1139          ! ------------------- 
    1140          REWIND( numnatp_ref )              ! Namelist natopt in reference namelist : MEDUSA diagnostics 
    1141          READ  ( numnatp_ref, natopt, IOSTAT = ios, ERR = 907) 
    1142 907      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in reference namelist', lwp ) 
    1143  
    1144          REWIND( numnatp_cfg )              ! Namelist natopt in configuration namelist : MEDUSA diagnostics 
    1145          READ  ( numnatp_cfg, natopt, IOSTAT = ios, ERR = 908 ) 
    1146 908      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in configuration namelist', lwp ) 
    1147          IF(lwm) WRITE ( numonp, natopt ) 
     999      ! Namelist natopt 
     1000      ! ------------------- 
     1001      REWIND( numnatp_ref )     ! Namelist natopt in reference namelist : MEDUSA diagnostics 
     1002      READ  ( numnatp_ref, natopt, IOSTAT = ios, ERR = 907) 
     1003 907  IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in reference namelist', lwp ) 
     1004       
     1005      REWIND( numnatp_cfg )     ! Namelist natopt in configuration namelist : MEDUSA diagnostics 
     1006      READ  ( numnatp_cfg, natopt, IOSTAT = ios, ERR = 908 ) 
     1007 908  IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in configuration namelist', lwp ) 
     1008      IF(lwm) WRITE ( numonp, natopt ) 
    11481009 
    11491010      IF(lwp) THEN 
     
    11581019         WRITE(numout,*) ' chla/chla+phea ratio            rpig  = ',rpig 
    11591020         WRITE(numout,*) ' ' 
    1160  
    11611021      ENDIF 
    11621022 
     
    11741034         endif 
    11751035         WRITE(numout,*) ' ' 
    1176        ENDIF 
     1036      ENDIF 
    11771037 
    11781038   END SUBROUTINE trc_nam_medusa 
     
    11901050      !!--------------------------------------------------------------------- 
    11911051      !! 
    1192       !! 
    11931052      !!----------------------------------------------------------------------             
    11941053      !! Variable conventions 
     
    11961055      !! 
    11971056      IF (iom_use("INVTN")) THEN  
    1198           med_diag%INVTN%dgsave = .TRUE. 
    1199       ELSE  
    1200           med_diag%INVTN%dgsave = .FALSE. 
    1201       ENDIF 
    1202       IF  (iom_use("INVTSI")) THEN  
    1203           med_diag%INVTSI%dgsave = .TRUE. 
    1204       ELSE  
    1205           med_diag%INVTSI%dgsave = .FALSE. 
    1206       ENDIF 
    1207       IF  (iom_use("INVTFE")) THEN  
    1208           med_diag%INVTFE%dgsave = .TRUE. 
    1209       ELSE  
    1210           med_diag%INVTFE%dgsave = .FALSE. 
    1211       ENDIF 
    1212       IF  (iom_use("PRN")) THEN  
    1213           med_diag%PRN%dgsave = .TRUE. 
    1214       ELSE  
    1215           med_diag%PRN%dgsave = .FALSE. 
    1216       ENDIF 
    1217       IF  (iom_use("MPN")) THEN  
    1218           med_diag%MPN%dgsave = .TRUE. 
    1219       ELSE  
    1220           med_diag%MPN%dgsave = .FALSE. 
    1221       ENDIF 
    1222       IF  (iom_use("PRD")) THEN  
    1223           med_diag%PRD%dgsave = .TRUE. 
    1224       ELSE  
    1225           med_diag%PRD%dgsave = .FALSE. 
    1226       ENDIF 
    1227       IF  (iom_use("MPD")) THEN  
    1228           med_diag%MPD%dgsave = .TRUE. 
    1229       ELSE  
    1230           med_diag%MPD%dgsave = .FALSE. 
    1231       ENDIF 
    1232       IF  (iom_use("DSED")) THEN  
    1233           med_diag%DSED%dgsave = .TRUE. 
    1234       ELSE  
    1235           med_diag%DSED%dgsave = .FALSE. 
    1236       ENDIF 
    1237       IF  (iom_use("OPAL")) THEN  
    1238           med_diag%OPAL%dgsave = .TRUE. 
    1239       ELSE  
    1240           med_diag%OPAL%dgsave = .FALSE. 
    1241       ENDIF 
    1242       IF  (iom_use("OPALDISS")) THEN  
    1243           med_diag%OPALDISS%dgsave = .TRUE. 
    1244       ELSE  
    1245           med_diag%OPALDISS%dgsave = .FALSE. 
    1246       ENDIF 
    1247       IF  (iom_use("GMIPn")) THEN  
    1248           med_diag%GMIPn%dgsave = .TRUE. 
    1249       ELSE  
    1250           med_diag%GMIPn%dgsave = .FALSE. 
    1251       ENDIF 
    1252       IF  (iom_use("GMID")) THEN  
    1253           med_diag%GMID%dgsave = .TRUE. 
    1254       ELSE  
    1255           med_diag%GMID%dgsave = .FALSE. 
    1256       ENDIF 
    1257       IF  (iom_use("MZMI")) THEN  
    1258           med_diag%MZMI%dgsave = .TRUE. 
    1259       ELSE  
    1260           med_diag%MZMI%dgsave = .FALSE. 
    1261       ENDIF 
    1262       IF  (iom_use("GMEPN")) THEN  
    1263           med_diag%GMEPN%dgsave = .TRUE. 
    1264       ELSE  
    1265           med_diag%GMEPN%dgsave = .FALSE. 
    1266       ENDIF 
    1267       IF  (iom_use("GMEPD")) THEN  
    1268           med_diag%GMEPD%dgsave = .TRUE. 
    1269       ELSE  
    1270           med_diag%GMEPD%dgsave = .FALSE. 
    1271       ENDIF 
    1272       IF  (iom_use("GMEZMI")) THEN  
    1273           med_diag%GMEZMI%dgsave = .TRUE. 
    1274       ELSE  
    1275           med_diag%GMEZMI%dgsave = .FALSE. 
    1276       ENDIF 
    1277       IF  (iom_use("GMED")) THEN  
    1278           med_diag%GMED%dgsave = .TRUE. 
    1279       ELSE  
    1280           med_diag%GMED%dgsave = .FALSE. 
    1281       ENDIF 
    1282       IF  (iom_use("MZME")) THEN  
    1283           med_diag%MZME%dgsave = .TRUE. 
    1284       ELSE  
    1285           med_diag%MZME%dgsave = .FALSE. 
    1286       ENDIF 
    1287       IF  (iom_use("DEXP")) THEN  
    1288           med_diag%DEXP%dgsave = .TRUE. 
    1289       ELSE  
    1290           med_diag%DEXP%dgsave = .FALSE. 
    1291       ENDIF 
    1292       IF  (iom_use("DETN")) THEN  
    1293           med_diag%DETN%dgsave = .TRUE. 
    1294       ELSE  
    1295           med_diag%DETN%dgsave = .FALSE. 
    1296       ENDIF 
    1297       IF  (iom_use("MDET")) THEN  
    1298           med_diag%MDET%dgsave = .TRUE. 
    1299       ELSE  
    1300           med_diag%MDET%dgsave = .FALSE. 
    1301       ENDIF 
    1302       IF  (iom_use("AEOLIAN")) THEN  
    1303           med_diag%AEOLIAN%dgsave = .TRUE. 
    1304       ELSE  
    1305           med_diag%AEOLIAN%dgsave = .FALSE. 
    1306       ENDIF 
    1307       IF  (iom_use("BENTHIC")) THEN  
    1308           med_diag%BENTHIC%dgsave = .TRUE. 
    1309       ELSE  
    1310           med_diag%BENTHIC%dgsave = .FALSE. 
    1311       ENDIF 
    1312       IF  (iom_use("SCAVENGE")) THEN  
    1313           med_diag%SCAVENGE%dgsave = .TRUE. 
    1314       ELSE  
    1315           med_diag%SCAVENGE%dgsave = .FALSE. 
    1316       ENDIF 
    1317       IF  (iom_use("PN_JLIM")) THEN  
    1318           med_diag%PN_JLIM%dgsave = .TRUE. 
    1319       ELSE  
    1320           med_diag%PN_JLIM%dgsave = .FALSE. 
    1321       ENDIF 
    1322       IF  (iom_use("PN_NLIM")) THEN  
    1323           med_diag%PN_NLIM%dgsave = .TRUE. 
    1324       ELSE  
    1325           med_diag%PN_NLIM%dgsave = .FALSE. 
    1326       ENDIF 
    1327       IF  (iom_use("PN_FELIM")) THEN  
    1328           med_diag%PN_FELIM%dgsave = .TRUE. 
    1329       ELSE  
    1330           med_diag%PN_FELIM%dgsave = .FALSE. 
    1331       ENDIF 
    1332       IF  (iom_use("PD_JLIM")) THEN  
    1333           med_diag%PD_JLIM%dgsave = .TRUE. 
    1334       ELSE  
    1335           med_diag%PD_JLIM%dgsave = .FALSE. 
    1336       ENDIF 
    1337       IF  (iom_use("PD_NLIM")) THEN  
    1338           med_diag%PD_NLIM%dgsave = .TRUE. 
    1339       ELSE  
    1340           med_diag%PD_NLIM%dgsave = .FALSE. 
    1341       ENDIF 
    1342       IF  (iom_use("PD_FELIM")) THEN  
    1343           med_diag%PD_FELIM%dgsave = .TRUE. 
    1344       ELSE  
    1345           med_diag%PD_FELIM%dgsave = .FALSE. 
    1346       ENDIF 
    1347       IF  (iom_use("PD_SILIM")) THEN  
    1348           med_diag%PD_SILIM%dgsave = .TRUE. 
    1349       ELSE  
    1350           med_diag%PD_SILIM%dgsave = .FALSE. 
    1351       ENDIF 
    1352       IF  (iom_use("PDSILIM2")) THEN  
    1353           med_diag%PDSILIM2%dgsave = .TRUE. 
    1354       ELSE  
    1355           med_diag%PDSILIM2%dgsave = .FALSE. 
    1356       ENDIF 
    1357       IF  (iom_use("SDT__100")) THEN  
    1358           med_diag%SDT__100%dgsave = .TRUE. 
    1359       ELSE  
    1360           med_diag%SDT__100%dgsave = .FALSE. 
    1361       ENDIF 
    1362       IF  (iom_use("SDT__200")) THEN  
    1363           med_diag%SDT__200%dgsave = .TRUE. 
    1364       ELSE  
    1365           med_diag%SDT__200%dgsave = .FALSE. 
    1366       ENDIF 
    1367       IF  (iom_use("SDT__500")) THEN  
    1368           med_diag%SDT__500%dgsave = .TRUE. 
    1369       ELSE  
    1370           med_diag%SDT__500%dgsave = .FALSE. 
    1371       ENDIF 
    1372       IF  (iom_use("SDT_1000")) THEN  
    1373           med_diag%SDT_1000%dgsave = .TRUE. 
    1374       ELSE  
    1375           med_diag%SDT_1000%dgsave = .FALSE. 
    1376       ENDIF 
    1377       IF  (iom_use("TOTREG_N")) THEN  
    1378           med_diag%TOTREG_N%dgsave = .TRUE. 
    1379       ELSE  
    1380           med_diag%TOTREG_N%dgsave = .FALSE. 
    1381       ENDIF 
    1382       IF  (iom_use("TOTRG_SI")) THEN  
    1383           med_diag%TOTRG_SI%dgsave = .TRUE. 
    1384       ELSE  
    1385           med_diag%TOTRG_SI%dgsave = .FALSE. 
    1386       ENDIF 
    1387       IF  (iom_use("REG__100")) THEN  
    1388           med_diag%REG__100%dgsave = .TRUE. 
    1389       ELSE  
    1390           med_diag%REG__100%dgsave = .FALSE. 
    1391       ENDIF 
    1392       IF  (iom_use("REG__200")) THEN  
    1393           med_diag%REG__200%dgsave = .TRUE. 
    1394       ELSE  
    1395           med_diag%REG__200%dgsave = .FALSE. 
    1396       ENDIF 
    1397       IF  (iom_use("REG__500")) THEN  
    1398           med_diag%REG__500%dgsave = .TRUE. 
    1399       ELSE  
    1400           med_diag%REG__500%dgsave = .FALSE. 
    1401       ENDIF 
    1402       IF  (iom_use("REG_1000")) THEN  
    1403           med_diag%REG_1000%dgsave = .TRUE. 
    1404       ELSE  
    1405           med_diag%REG_1000%dgsave = .FALSE. 
    1406       ENDIF 
    1407       IF  (iom_use("FASTN")) THEN  
    1408           med_diag%FASTN%dgsave = .TRUE. 
    1409       ELSE  
    1410           med_diag%FASTN%dgsave = .FALSE. 
    1411       ENDIF 
    1412       IF  (iom_use("FASTSI")) THEN  
    1413           med_diag%FASTSI%dgsave = .TRUE. 
    1414       ELSE  
    1415           med_diag%FASTSI%dgsave = .FALSE. 
    1416       ENDIF 
    1417       IF  (iom_use("FASTFE")) THEN  
    1418           med_diag%FASTFE%dgsave = .TRUE. 
    1419       ELSE  
    1420           med_diag%FASTFE%dgsave = .FALSE. 
    1421       ENDIF 
    1422       IF  (iom_use("FASTC")) THEN  
    1423           med_diag%FASTC%dgsave = .TRUE. 
    1424       ELSE  
    1425           med_diag%FASTC%dgsave = .FALSE. 
    1426       ENDIF 
    1427       IF  (iom_use("FASTCA")) THEN  
    1428           med_diag%FASTCA%dgsave = .TRUE. 
    1429       ELSE  
    1430           med_diag%FASTCA%dgsave = .FALSE. 
    1431       ENDIF 
    1432       IF  (iom_use("FDT__100")) THEN  
    1433           med_diag%FDT__100%dgsave = .TRUE. 
    1434       ELSE  
    1435           med_diag%FDT__100%dgsave = .FALSE. 
    1436       ENDIF 
    1437       IF  (iom_use("FDT__200")) THEN  
    1438           med_diag%FDT__200%dgsave = .TRUE. 
    1439       ELSE  
    1440           med_diag%FDT__200%dgsave = .FALSE. 
    1441       ENDIF 
    1442       IF  (iom_use("FDT__500")) THEN  
    1443           med_diag%FDT__500%dgsave = .TRUE. 
    1444       ELSE  
    1445           med_diag%FDT__500%dgsave = .FALSE. 
    1446       ENDIF 
    1447       IF  (iom_use("FDT_1000")) THEN  
    1448           med_diag%FDT_1000%dgsave = .TRUE. 
    1449       ELSE  
    1450           med_diag%FDT_1000%dgsave = .FALSE. 
    1451       ENDIF 
    1452       IF  (iom_use("RG__100F")) THEN  
    1453           med_diag%RG__100F%dgsave = .TRUE. 
    1454       ELSE  
    1455           med_diag%RG__100F%dgsave = .FALSE. 
    1456       ENDIF 
    1457       IF  (iom_use("RG__200F")) THEN  
    1458           med_diag%RG__200F%dgsave = .TRUE. 
    1459       ELSE  
    1460           med_diag%RG__200F%dgsave = .FALSE. 
    1461       ENDIF 
    1462       IF  (iom_use("RG__500F")) THEN  
    1463           med_diag%RG__500F%dgsave = .TRUE. 
    1464       ELSE  
    1465           med_diag%RG__500F%dgsave = .FALSE. 
    1466       ENDIF 
    1467       IF  (iom_use("RG_1000F")) THEN  
    1468           med_diag%RG_1000F%dgsave = .TRUE. 
    1469       ELSE  
    1470           med_diag%RG_1000F%dgsave = .FALSE. 
    1471       ENDIF 
    1472       IF  (iom_use("FDS__100")) THEN  
    1473           med_diag%FDS__100%dgsave = .TRUE. 
    1474       ELSE  
    1475           med_diag%FDS__100%dgsave = .FALSE. 
    1476       ENDIF 
    1477       IF  (iom_use("FDS__200")) THEN  
    1478           med_diag%FDS__200%dgsave = .TRUE. 
    1479       ELSE  
    1480           med_diag%FDS__200%dgsave = .FALSE. 
    1481       ENDIF 
    1482       IF  (iom_use("FDS__500")) THEN  
    1483           med_diag%FDS__500%dgsave = .TRUE. 
    1484       ELSE  
    1485           med_diag%FDS__500%dgsave = .FALSE. 
    1486       ENDIF 
    1487       IF  (iom_use("FDS_1000")) THEN  
    1488           med_diag%FDS_1000%dgsave = .TRUE. 
    1489       ELSE  
    1490           med_diag%FDS_1000%dgsave = .FALSE. 
    1491       ENDIF 
    1492       IF  (iom_use("RGS_100F")) THEN  
    1493           med_diag%RGS_100F%dgsave = .TRUE. 
    1494       ELSE  
    1495           med_diag%RGS_100F%dgsave = .FALSE. 
    1496       ENDIF 
    1497       IF  (iom_use("RGS_200F")) THEN  
    1498           med_diag%RGS_200F%dgsave = .TRUE. 
    1499       ELSE  
    1500           med_diag%RGS_200F%dgsave = .FALSE. 
    1501       ENDIF 
    1502       IF  (iom_use("RGS_500F")) THEN  
    1503           med_diag%RGS_500F%dgsave = .TRUE. 
    1504       ELSE  
    1505           med_diag%RGS_500F%dgsave = .FALSE. 
    1506       ENDIF 
    1507       IF  (iom_use("RGS1000F")) THEN  
    1508           med_diag%RGS1000F%dgsave = .TRUE. 
    1509       ELSE  
    1510           med_diag%RGS1000F%dgsave = .FALSE. 
    1511       ENDIF 
    1512       IF  (iom_use("REMINN")) THEN  
    1513           med_diag%REMINN%dgsave = .TRUE. 
    1514       ELSE  
    1515           med_diag%REMINN%dgsave = .FALSE. 
    1516       ENDIF 
    1517       IF  (iom_use("REMINSI")) THEN  
    1518           med_diag%REMINSI%dgsave = .TRUE. 
    1519       ELSE  
    1520           med_diag%REMINSI%dgsave = .FALSE. 
    1521       ENDIF 
    1522       IF  (iom_use("REMINFE")) THEN  
    1523           med_diag%REMINFE%dgsave = .TRUE. 
    1524       ELSE  
    1525           med_diag%REMINFE%dgsave = .FALSE. 
    1526       ENDIF 
    1527       IF  (iom_use("REMINC")) THEN  
    1528           med_diag%REMINC%dgsave = .TRUE. 
    1529       ELSE  
    1530           med_diag%REMINC%dgsave = .FALSE. 
    1531       ENDIF 
    1532       IF  (iom_use("REMINCA")) THEN  
    1533           med_diag%REMINCA%dgsave = .TRUE. 
    1534       ELSE  
    1535           med_diag%REMINCA%dgsave = .FALSE. 
    1536       ENDIF 
    1537       IF  (iom_use("SEAFLRN")) THEN  
    1538           med_diag%SEAFLRN%dgsave = .TRUE. 
    1539       ELSE  
    1540           med_diag%SEAFLRN%dgsave = .FALSE. 
    1541       ENDIF 
    1542       IF  (iom_use("SEAFLRSI")) THEN  
    1543           med_diag%SEAFLRSI%dgsave = .TRUE. 
    1544       ELSE  
    1545           med_diag%SEAFLRSI%dgsave = .FALSE. 
    1546       ENDIF 
    1547       IF  (iom_use("SEAFLRFE")) THEN  
    1548           med_diag%SEAFLRFE%dgsave = .TRUE. 
    1549       ELSE  
    1550           med_diag%SEAFLRFE%dgsave = .FALSE. 
    1551       ENDIF 
    1552       IF  (iom_use("SEAFLRC")) THEN  
    1553           med_diag%SEAFLRC%dgsave = .TRUE. 
    1554       ELSE  
    1555           med_diag%SEAFLRC%dgsave = .FALSE. 
    1556       ENDIF 
    1557       IF  (iom_use("SEAFLRCA")) THEN  
    1558           med_diag%SEAFLRCA%dgsave = .TRUE. 
    1559       ELSE  
    1560           med_diag%SEAFLRCA%dgsave = .FALSE. 
    1561       ENDIF 
    1562       IF  (iom_use("MED_QSR")) THEN  
    1563           med_diag%MED_QSR%dgsave = .TRUE. 
    1564       ELSE  
    1565           med_diag%MED_QSR%dgsave = .FALSE. 
    1566       ENDIF 
    1567       IF  (iom_use("MED_XPAR")) THEN  
    1568           med_diag%MED_XPAR%dgsave = .TRUE. 
    1569       ELSE  
    1570           med_diag%MED_XPAR%dgsave = .FALSE. 
    1571       ENDIF 
    1572       IF  (iom_use("INTFLX_N")) THEN  
    1573           med_diag%INTFLX_N%dgsave = .TRUE. 
    1574       ELSE  
    1575           med_diag%INTFLX_N%dgsave = .FALSE. 
    1576       ENDIF 
    1577       IF  (iom_use("INTFLX_SI")) THEN  
    1578           med_diag%INTFLX_SI%dgsave = .TRUE. 
    1579       ELSE  
    1580           med_diag%INTFLX_SI%dgsave = .FALSE. 
    1581       ENDIF 
    1582       IF  (iom_use("INTFLX_FE")) THEN  
    1583           med_diag%INTFLX_FE%dgsave = .TRUE. 
    1584       ELSE  
    1585           med_diag%INTFLX_FE%dgsave = .FALSE. 
    1586       ENDIF 
    1587       IF  (iom_use("INT_PN")) THEN  
    1588           med_diag%INT_PN%dgsave = .TRUE. 
    1589       ELSE  
    1590           med_diag%INT_PN%dgsave = .FALSE. 
    1591       ENDIF 
    1592       IF  (iom_use("INT_PD")) THEN  
    1593           med_diag%INT_PD%dgsave = .TRUE. 
    1594       ELSE  
    1595           med_diag%INT_PD%dgsave = .FALSE. 
    1596       ENDIF 
    1597       IF  (iom_use("ML_PRN")) THEN  
    1598           med_diag%ML_PRN%dgsave = .TRUE. 
    1599       ELSE  
    1600           med_diag%ML_PRN%dgsave = .FALSE. 
    1601       ENDIF 
    1602       IF  (iom_use("ML_PRD")) THEN  
    1603           med_diag%ML_PRD%dgsave = .TRUE. 
    1604       ELSE  
    1605           med_diag%ML_PRD%dgsave = .FALSE. 
    1606       ENDIF 
    1607       IF  (iom_use("OCAL_CCD")) THEN  
    1608           med_diag%OCAL_CCD%dgsave = .TRUE. 
    1609       ELSE  
    1610           med_diag%OCAL_CCD%dgsave = .FALSE. 
    1611       ENDIF 
    1612       IF  (iom_use("OCAL_LVL")) THEN  
    1613           med_diag%OCAL_LVL%dgsave = .TRUE. 
    1614       ELSE  
    1615           med_diag%OCAL_LVL%dgsave = .FALSE. 
    1616       ENDIF 
    1617       IF  (iom_use("FE_0000")) THEN  
    1618           med_diag%FE_0000%dgsave = .TRUE. 
    1619       ELSE  
    1620           med_diag%FE_0000%dgsave = .FALSE. 
    1621       ENDIF 
    1622       IF  (iom_use("FE_0100")) THEN  
    1623           med_diag%FE_0100%dgsave = .TRUE. 
    1624       ELSE  
    1625           med_diag%FE_0100%dgsave = .FALSE. 
    1626       ENDIF 
    1627       IF  (iom_use("FE_0200")) THEN  
    1628           med_diag%FE_0200%dgsave = .TRUE. 
    1629       ELSE  
    1630           med_diag%FE_0200%dgsave = .FALSE. 
    1631       ENDIF 
    1632       IF  (iom_use("FE_0500")) THEN  
    1633           med_diag%FE_0500%dgsave = .TRUE. 
    1634       ELSE  
    1635           med_diag%FE_0500%dgsave = .FALSE. 
    1636       ENDIF 
    1637       IF  (iom_use("FE_1000")) THEN  
    1638           med_diag%FE_1000%dgsave = .TRUE. 
    1639       ELSE  
    1640           med_diag%FE_1000%dgsave = .FALSE. 
    1641       ENDIF 
    1642       IF  (iom_use("MED_XZE")) THEN  
    1643           med_diag%MED_XZE%dgsave = .TRUE. 
    1644       ELSE  
    1645           med_diag%MED_XZE%dgsave = .FALSE. 
    1646       ENDIF 
    1647       IF  (iom_use("WIND")) THEN  
    1648           med_diag%WIND%dgsave = .TRUE. 
    1649       ELSE  
    1650           med_diag%WIND%dgsave = .FALSE. 
    1651       ENDIF 
    1652       IF  (iom_use("ATM_PCO2")) THEN  
    1653           med_diag%ATM_PCO2%dgsave = .TRUE. 
    1654       ELSE  
    1655           med_diag%ATM_PCO2%dgsave = .FALSE. 
    1656       ENDIF 
    1657       IF  (iom_use("OCN_PH")) THEN  
    1658           med_diag%OCN_PH%dgsave = .TRUE. 
    1659       ELSE  
    1660           med_diag%OCN_PH%dgsave = .FALSE. 
    1661       ENDIF 
    1662       IF  (iom_use("OCN_PCO2")) THEN  
    1663           med_diag%OCN_PCO2%dgsave = .TRUE. 
    1664       ELSE  
    1665           med_diag%OCN_PCO2%dgsave = .FALSE. 
    1666       ENDIF 
    1667       IF  (iom_use("OCNH2CO3")) THEN  
    1668           med_diag%OCNH2CO3%dgsave = .TRUE. 
    1669       ELSE  
    1670           med_diag%OCNH2CO3%dgsave = .FALSE. 
    1671       ENDIF 
    1672       IF  (iom_use("OCN_HCO3")) THEN  
    1673           med_diag%OCN_HCO3%dgsave = .TRUE. 
    1674       ELSE  
    1675           med_diag%OCN_HCO3%dgsave = .FALSE. 
    1676       ENDIF 
    1677       IF  (iom_use("OCN_CO3")) THEN  
    1678           med_diag%OCN_CO3%dgsave = .TRUE. 
    1679       ELSE  
    1680           med_diag%OCN_CO3%dgsave = .FALSE. 
    1681       ENDIF 
    1682       IF  (iom_use("CO2FLUX")) THEN  
    1683           med_diag%CO2FLUX%dgsave = .TRUE. 
    1684       ELSE  
    1685           med_diag%CO2FLUX%dgsave = .FALSE. 
    1686       ENDIF 
    1687       IF  (iom_use("OM_CAL")) THEN  
    1688           med_diag%OM_CAL%dgsave = .TRUE. 
    1689       ELSE  
    1690           med_diag%OM_CAL%dgsave = .FALSE. 
    1691       ENDIF 
    1692       IF  (iom_use("OM_ARG")) THEN  
    1693           med_diag%OM_ARG%dgsave = .TRUE. 
    1694       ELSE  
    1695           med_diag%OM_ARG%dgsave = .FALSE. 
    1696       ENDIF 
    1697       IF  (iom_use("TCO2")) THEN  
    1698           med_diag%TCO2%dgsave = .TRUE. 
    1699       ELSE  
    1700           med_diag%TCO2%dgsave = .FALSE. 
    1701       ENDIF 
    1702       IF  (iom_use("TALK")) THEN  
    1703           med_diag%TALK%dgsave = .TRUE. 
    1704       ELSE  
    1705           med_diag%TALK%dgsave = .FALSE. 
    1706       ENDIF 
    1707       IF  (iom_use("KW660")) THEN  
    1708           med_diag%KW660%dgsave = .TRUE. 
    1709       ELSE  
    1710           med_diag%KW660%dgsave = .FALSE. 
    1711       ENDIF 
    1712       IF  (iom_use("ATM_PP0")) THEN  
    1713           med_diag%ATM_PP0%dgsave = .TRUE. 
    1714       ELSE  
    1715           med_diag%ATM_PP0%dgsave = .FALSE. 
    1716       ENDIF 
    1717       IF  (iom_use("O2FLUX")) THEN  
    1718           med_diag%O2FLUX%dgsave = .TRUE. 
    1719       ELSE  
    1720           med_diag%O2FLUX%dgsave = .FALSE. 
    1721       ENDIF 
    1722       IF  (iom_use("O2SAT")) THEN  
    1723           med_diag%O2SAT%dgsave = .TRUE. 
    1724       ELSE  
    1725           med_diag%O2SAT%dgsave = .FALSE. 
    1726       ENDIF 
    1727       IF  (iom_use("CAL_CCD")) THEN  
    1728           med_diag%CAL_CCD%dgsave = .TRUE. 
    1729       ELSE  
    1730           med_diag%CAL_CCD%dgsave = .FALSE. 
    1731       ENDIF 
    1732       IF  (iom_use("ARG_CCD")) THEN  
    1733           med_diag%ARG_CCD%dgsave = .TRUE. 
    1734       ELSE  
    1735           med_diag%ARG_CCD%dgsave = .FALSE. 
    1736       ENDIF 
    1737       IF  (iom_use("SFR_OCAL")) THEN  
    1738           med_diag%SFR_OCAL%dgsave = .TRUE. 
    1739       ELSE  
    1740           med_diag%SFR_OCAL%dgsave = .FALSE. 
    1741       ENDIF 
    1742       IF  (iom_use("SFR_OARG")) THEN  
    1743           med_diag%SFR_OARG%dgsave = .TRUE. 
    1744       ELSE  
    1745           med_diag%SFR_OARG%dgsave = .FALSE. 
    1746       ENDIF 
    1747       IF  (iom_use("N_PROD")) THEN  
    1748           med_diag%N_PROD%dgsave = .TRUE. 
    1749       ELSE  
    1750           med_diag%N_PROD%dgsave = .FALSE. 
    1751       ENDIF 
    1752       IF  (iom_use("N_CONS")) THEN  
    1753           med_diag%N_CONS%dgsave = .TRUE. 
    1754       ELSE  
    1755           med_diag%N_CONS%dgsave = .FALSE. 
    1756       ENDIF 
    1757       IF  (iom_use("C_PROD")) THEN  
    1758           med_diag%C_PROD%dgsave = .TRUE. 
    1759       ELSE  
    1760           med_diag%C_PROD%dgsave = .FALSE. 
    1761       ENDIF 
    1762       IF  (iom_use("C_CONS")) THEN  
    1763           med_diag%C_CONS%dgsave = .TRUE. 
    1764       ELSE  
    1765           med_diag%C_CONS%dgsave = .FALSE. 
    1766       ENDIF 
    1767       IF  (iom_use("O2_PROD")) THEN  
    1768           med_diag%O2_PROD%dgsave = .TRUE. 
    1769       ELSE  
    1770           med_diag%O2_PROD%dgsave = .FALSE. 
    1771       ENDIF 
    1772       IF  (iom_use("O2_CONS")) THEN  
    1773           med_diag%O2_CONS%dgsave = .TRUE. 
    1774       ELSE  
    1775           med_diag%O2_CONS%dgsave = .FALSE. 
    1776       ENDIF 
    1777       IF  (iom_use("O2_ANOX")) THEN  
    1778           med_diag%O2_ANOX%dgsave = .TRUE. 
    1779       ELSE  
    1780           med_diag%O2_ANOX%dgsave = .FALSE. 
    1781       ENDIF 
    1782       IF  (iom_use("RR_0100")) THEN  
    1783           med_diag%RR_0100%dgsave = .TRUE. 
    1784       ELSE  
    1785           med_diag%RR_0100%dgsave = .FALSE. 
    1786       ENDIF 
    1787       IF  (iom_use("RR_0500")) THEN  
    1788           med_diag%RR_0500%dgsave = .TRUE. 
    1789       ELSE  
    1790           med_diag%RR_0500%dgsave = .FALSE. 
    1791       ENDIF 
    1792       IF  (iom_use("RR_1000")) THEN  
    1793           med_diag%RR_1000%dgsave = .TRUE. 
    1794       ELSE  
    1795           med_diag%RR_1000%dgsave = .FALSE. 
    1796       ENDIF 
    1797       IF  (iom_use("IBEN_N")) THEN  
    1798           med_diag%IBEN_N%dgsave = .TRUE. 
    1799       ELSE  
    1800           med_diag%IBEN_N%dgsave = .FALSE. 
    1801       ENDIF 
    1802       IF  (iom_use("IBEN_FE")) THEN  
    1803           med_diag%IBEN_FE%dgsave = .TRUE. 
    1804       ELSE  
    1805           med_diag%IBEN_FE%dgsave = .FALSE. 
    1806       ENDIF 
    1807       IF  (iom_use("IBEN_C")) THEN  
    1808           med_diag%IBEN_C%dgsave = .TRUE. 
    1809       ELSE  
    1810           med_diag%IBEN_C%dgsave = .FALSE. 
    1811       ENDIF 
    1812       IF  (iom_use("IBEN_SI")) THEN  
    1813           med_diag%IBEN_SI%dgsave = .TRUE. 
    1814       ELSE  
    1815           med_diag%IBEN_SI%dgsave = .FALSE. 
    1816       ENDIF 
    1817       IF  (iom_use("IBEN_CA")) THEN  
    1818           med_diag%IBEN_CA%dgsave = .TRUE. 
    1819       ELSE  
    1820           med_diag%IBEN_CA%dgsave = .FALSE. 
    1821       ENDIF 
    1822       IF  (iom_use("OBEN_N")) THEN  
    1823           med_diag%OBEN_N%dgsave = .TRUE. 
    1824       ELSE  
    1825           med_diag%OBEN_N%dgsave = .FALSE. 
    1826       ENDIF 
    1827       IF  (iom_use("OBEN_FE")) THEN  
    1828           med_diag%OBEN_FE%dgsave = .TRUE. 
    1829       ELSE  
    1830           med_diag%OBEN_FE%dgsave = .FALSE. 
    1831       ENDIF 
    1832       IF  (iom_use("OBEN_C")) THEN  
    1833           med_diag%OBEN_C%dgsave = .TRUE. 
    1834       ELSE  
    1835           med_diag%OBEN_C%dgsave = .FALSE. 
    1836       ENDIF 
    1837       IF  (iom_use("OBEN_SI")) THEN  
    1838           med_diag%OBEN_SI%dgsave = .TRUE. 
    1839       ELSE  
    1840           med_diag%OBEN_SI%dgsave = .FALSE. 
    1841       ENDIF 
    1842       IF  (iom_use("OBEN_CA")) THEN  
    1843           med_diag%OBEN_CA%dgsave = .TRUE. 
    1844       ELSE  
    1845           med_diag%OBEN_CA%dgsave = .FALSE. 
    1846       ENDIF 
    1847       IF  (iom_use("BEN_N")) THEN  
    1848           med_diag%BEN_N%dgsave = .TRUE. 
    1849       ELSE  
    1850           med_diag%BEN_N%dgsave = .FALSE. 
    1851       ENDIF 
    1852       IF  (iom_use("BEN_FE")) THEN  
    1853           med_diag%BEN_FE%dgsave = .TRUE. 
    1854       ELSE  
    1855           med_diag%BEN_FE%dgsave = .FALSE. 
    1856       ENDIF 
    1857       IF  (iom_use("BEN_C")) THEN  
    1858           med_diag%BEN_C%dgsave = .TRUE. 
    1859       ELSE  
    1860           med_diag%BEN_C%dgsave = .FALSE. 
    1861       ENDIF 
    1862       IF  (iom_use("BEN_SI")) THEN  
    1863           med_diag%BEN_SI%dgsave = .TRUE. 
    1864       ELSE  
    1865           med_diag%BEN_SI%dgsave = .FALSE. 
    1866       ENDIF 
    1867       IF  (iom_use("BEN_CA")) THEN  
    1868           med_diag%BEN_CA%dgsave = .TRUE. 
    1869       ELSE  
    1870           med_diag%BEN_CA%dgsave = .FALSE. 
    1871       ENDIF 
    1872       IF  (iom_use("RUNOFF")) THEN  
    1873           med_diag%RUNOFF%dgsave = .TRUE. 
    1874       ELSE  
    1875           med_diag%RUNOFF%dgsave = .FALSE. 
    1876       ENDIF 
    1877       IF  (iom_use("RIV_N")) THEN  
    1878           med_diag%RIV_N%dgsave = .TRUE. 
    1879       ELSE  
    1880           med_diag%RIV_N%dgsave = .FALSE. 
    1881       ENDIF 
    1882       IF  (iom_use("RIV_SI")) THEN  
    1883           med_diag%RIV_SI%dgsave = .TRUE. 
    1884       ELSE  
    1885           med_diag%RIV_SI%dgsave = .FALSE. 
    1886       ENDIF 
    1887       IF  (iom_use("RIV_C")) THEN  
    1888           med_diag%RIV_C%dgsave = .TRUE. 
    1889       ELSE  
    1890           med_diag%RIV_C%dgsave = .FALSE. 
    1891       ENDIF 
    1892       IF  (iom_use("RIV_ALK")) THEN  
    1893           med_diag%RIV_ALK%dgsave = .TRUE. 
    1894       ELSE  
    1895           med_diag%RIV_ALK%dgsave = .FALSE. 
    1896       ENDIF 
    1897       IF  (iom_use("DETC")) THEN  
    1898           med_diag%DETC%dgsave = .TRUE. 
    1899       ELSE  
    1900           med_diag%DETC%dgsave = .FALSE. 
    1901       ENDIF 
    1902       IF  (iom_use("SDC__100")) THEN  
    1903           med_diag%SDC__100%dgsave = .TRUE. 
    1904       ELSE  
    1905           med_diag%SDC__100%dgsave = .FALSE. 
    1906       ENDIF 
    1907       IF  (iom_use("SDC__200")) THEN  
    1908           med_diag%SDC__200%dgsave = .TRUE. 
    1909       ELSE  
    1910           med_diag%SDC__200%dgsave = .FALSE. 
    1911       ENDIF 
    1912       IF  (iom_use("SDC__500")) THEN  
    1913           med_diag%SDC__500%dgsave = .TRUE. 
    1914       ELSE  
    1915           med_diag%SDC__500%dgsave = .FALSE. 
    1916       ENDIF 
    1917       IF  (iom_use("SDC_1000")) THEN  
    1918           med_diag%SDC_1000%dgsave = .TRUE. 
    1919       ELSE  
    1920           med_diag%SDC_1000%dgsave = .FALSE. 
    1921       ENDIF 
    1922       IF  (iom_use("INVTC")) THEN  
    1923           med_diag%INVTC%dgsave = .TRUE. 
    1924       ELSE  
    1925           med_diag%INVTC%dgsave = .FALSE. 
    1926       ENDIF 
    1927       IF  (iom_use("INVTALK")) THEN  
    1928           med_diag%INVTALK%dgsave = .TRUE. 
    1929       ELSE  
    1930           med_diag%INVTALK%dgsave = .FALSE. 
    1931       ENDIF 
    1932       IF  (iom_use("INVTO2")) THEN  
    1933           med_diag%INVTO2%dgsave = .TRUE. 
    1934       ELSE  
    1935           med_diag%INVTO2%dgsave = .FALSE. 
    1936       ENDIF 
    1937       IF  (iom_use("LYSO_CA")) THEN  
    1938           med_diag%LYSO_CA%dgsave = .TRUE. 
    1939       ELSE  
    1940           med_diag%LYSO_CA%dgsave = .FALSE. 
    1941       ENDIF 
    1942       IF  (iom_use("COM_RESP")) THEN  
    1943           med_diag%COM_RESP%dgsave = .TRUE. 
    1944       ELSE  
    1945           med_diag%COM_RESP%dgsave = .FALSE. 
    1946       ENDIF 
    1947       IF  (iom_use("PN_LLOSS")) THEN  
    1948           med_diag%PN_LLOSS%dgsave = .TRUE. 
    1949       ELSE  
    1950           med_diag%PN_LLOSS%dgsave = .FALSE. 
    1951       ENDIF 
    1952       IF  (iom_use("PD_LLOSS")) THEN  
    1953           med_diag%PD_LLOSS%dgsave = .TRUE. 
    1954       ELSE  
    1955           med_diag%PD_LLOSS%dgsave = .FALSE. 
    1956       ENDIF 
    1957       IF  (iom_use("ZI_LLOSS")) THEN  
    1958           med_diag%ZI_LLOSS%dgsave = .TRUE. 
    1959       ELSE  
    1960           med_diag%ZI_LLOSS%dgsave = .FALSE. 
    1961       ENDIF 
    1962       IF  (iom_use("ZE_LLOSS")) THEN  
    1963           med_diag%ZE_LLOSS%dgsave = .TRUE. 
    1964       ELSE  
    1965           med_diag%ZE_LLOSS%dgsave = .FALSE. 
    1966       ENDIF 
    1967       IF  (iom_use("ZI_MES_N")) THEN  
    1968           med_diag%ZI_MES_N%dgsave = .TRUE. 
    1969       ELSE  
    1970           med_diag%ZI_MES_N%dgsave = .FALSE. 
    1971       ENDIF 
    1972       IF  (iom_use("ZI_MES_D")) THEN  
    1973           med_diag%ZI_MES_D%dgsave = .TRUE. 
    1974       ELSE  
    1975           med_diag%ZI_MES_D%dgsave = .FALSE. 
    1976       ENDIF 
    1977       IF  (iom_use("ZI_MES_C")) THEN  
    1978           med_diag%ZI_MES_C%dgsave = .TRUE. 
    1979       ELSE  
    1980           med_diag%ZI_MES_C%dgsave = .FALSE. 
    1981       ENDIF 
    1982       IF  (iom_use("ZI_MESDC")) THEN  
    1983           med_diag%ZI_MESDC%dgsave = .TRUE. 
    1984       ELSE  
    1985           med_diag%ZI_MESDC%dgsave = .FALSE. 
    1986       ENDIF 
    1987       IF  (iom_use("ZI_EXCR")) THEN  
    1988           med_diag%ZI_EXCR%dgsave = .TRUE. 
    1989       ELSE  
    1990           med_diag%ZI_EXCR%dgsave = .FALSE. 
    1991       ENDIF 
    1992       IF  (iom_use("ZI_RESP")) THEN  
    1993           med_diag%ZI_RESP%dgsave = .TRUE. 
    1994       ELSE  
    1995           med_diag%ZI_RESP%dgsave = .FALSE. 
    1996       ENDIF 
    1997       IF  (iom_use("ZI_GROW")) THEN  
    1998           med_diag%ZI_GROW%dgsave = .TRUE. 
    1999       ELSE  
    2000           med_diag%ZI_GROW%dgsave = .FALSE. 
    2001       ENDIF 
    2002       IF  (iom_use("ZE_MES_N")) THEN  
    2003           med_diag%ZE_MES_N%dgsave = .TRUE. 
    2004       ELSE  
    2005           med_diag%ZE_MES_N%dgsave = .FALSE. 
    2006       ENDIF 
    2007       IF  (iom_use("ZE_MES_D")) THEN  
    2008           med_diag%ZE_MES_D%dgsave = .TRUE. 
    2009       ELSE  
    2010           med_diag%ZE_MES_D%dgsave = .FALSE. 
    2011       ENDIF 
    2012       IF  (iom_use("ZE_MES_C")) THEN  
    2013           med_diag%ZE_MES_C%dgsave = .TRUE. 
    2014       ELSE  
    2015           med_diag%ZE_MES_C%dgsave = .FALSE. 
    2016       ENDIF 
    2017       IF  (iom_use("ZE_MESDC")) THEN  
    2018           med_diag%ZE_MESDC%dgsave = .TRUE. 
    2019       ELSE  
    2020           med_diag%ZE_MESDC%dgsave = .FALSE. 
    2021       ENDIF 
    2022       IF  (iom_use("ZE_EXCR")) THEN  
    2023           med_diag%ZE_EXCR%dgsave = .TRUE. 
    2024       ELSE  
    2025           med_diag%ZE_EXCR%dgsave = .FALSE. 
    2026       ENDIF 
    2027       IF  (iom_use("ZE_RESP")) THEN  
    2028           med_diag%ZE_RESP%dgsave = .TRUE. 
    2029       ELSE  
    2030           med_diag%ZE_RESP%dgsave = .FALSE. 
    2031       ENDIF 
    2032       IF  (iom_use("ZE_GROW")) THEN  
    2033           med_diag%ZE_GROW%dgsave = .TRUE. 
    2034       ELSE  
    2035           med_diag%ZE_GROW%dgsave = .FALSE. 
    2036       ENDIF 
    2037       IF  (iom_use("MDETC")) THEN  
    2038           med_diag%MDETC%dgsave = .TRUE. 
    2039       ELSE  
    2040           med_diag%MDETC%dgsave = .FALSE. 
    2041       ENDIF 
    2042       IF  (iom_use("GMIDC")) THEN  
    2043           med_diag%GMIDC%dgsave = .TRUE. 
    2044       ELSE  
    2045           med_diag%GMIDC%dgsave = .FALSE. 
    2046       ENDIF 
    2047       IF  (iom_use("GMEDC")) THEN  
    2048           med_diag%GMEDC%dgsave = .TRUE. 
    2049       ELSE  
    2050           med_diag%GMEDC%dgsave = .FALSE. 
    2051       ENDIF 
    2052       IF  (iom_use("INT_ZMI")) THEN  
    2053           med_diag%INT_ZMI%dgsave = .TRUE. 
    2054       ELSE  
    2055           med_diag%INT_ZMI%dgsave = .FALSE. 
    2056       ENDIF 
    2057       IF  (iom_use("INT_ZME")) THEN  
    2058           med_diag%INT_ZME%dgsave = .TRUE. 
    2059       ELSE  
    2060           med_diag%INT_ZME%dgsave = .FALSE. 
    2061       ENDIF 
    2062       IF  (iom_use("INT_DET")) THEN  
    2063           med_diag%INT_DET%dgsave = .TRUE. 
    2064       ELSE  
    2065           med_diag%INT_DET%dgsave = .FALSE. 
    2066       ENDIF 
    2067       IF  (iom_use("INT_DTC")) THEN  
    2068           med_diag%INT_DTC%dgsave = .TRUE. 
    2069       ELSE  
    2070           med_diag%INT_DTC%dgsave = .FALSE. 
    2071       ENDIF 
    2072       IF  (iom_use("DMS_SURF")) THEN  
    2073           med_diag%DMS_SURF%dgsave = .TRUE. 
    2074       ELSE  
    2075           med_diag%DMS_SURF%dgsave = .FALSE. 
    2076       ENDIF 
    2077       IF  (iom_use("DMS_ANDR")) THEN  
    2078           med_diag%DMS_ANDR%dgsave = .TRUE. 
    2079       ELSE  
    2080           med_diag%DMS_ANDR%dgsave = .FALSE. 
    2081       ENDIF 
    2082       IF  (iom_use("DMS_SIMO")) THEN  
    2083           med_diag%DMS_SIMO%dgsave = .TRUE. 
    2084       ELSE  
    2085           med_diag%DMS_SIMO%dgsave = .FALSE. 
    2086       ENDIF 
    2087       IF  (iom_use("DMS_ARAN")) THEN  
    2088           med_diag%DMS_ARAN%dgsave = .TRUE. 
    2089       ELSE  
    2090           med_diag%DMS_ARAN%dgsave = .FALSE. 
    2091       ENDIF 
    2092       IF  (iom_use("DMS_HALL")) THEN  
    2093           med_diag%DMS_HALL%dgsave = .TRUE. 
    2094       ELSE  
    2095           med_diag%DMS_HALL%dgsave = .FALSE. 
    2096       ENDIF 
    2097       IF  (iom_use("ATM_XCO2")) THEN  
    2098           med_diag%ATM_XCO2%dgsave = .TRUE. 
    2099       ELSE  
    2100           med_diag%ATM_XCO2%dgsave = .FALSE. 
    2101       ENDIF 
    2102       IF  (iom_use("OCN_FCO2")) THEN  
    2103           med_diag%OCN_FCO2%dgsave = .TRUE. 
    2104       ELSE  
    2105           med_diag%OCN_FCO2%dgsave = .FALSE. 
    2106       ENDIF 
    2107       IF  (iom_use("ATM_FCO2")) THEN  
    2108           med_diag%ATM_FCO2%dgsave = .TRUE. 
    2109       ELSE  
    2110           med_diag%ATM_FCO2%dgsave = .FALSE. 
    2111       ENDIF 
    2112       IF  (iom_use("OCN_RHOSW")) THEN  
    2113           med_diag%OCN_RHOSW%dgsave = .TRUE. 
    2114       ELSE  
    2115           med_diag%OCN_RHOSW%dgsave = .FALSE. 
    2116       ENDIF 
    2117       IF  (iom_use("OCN_SCHCO2")) THEN  
    2118           med_diag%OCN_SCHCO2%dgsave = .TRUE. 
    2119       ELSE  
    2120           med_diag%OCN_SCHCO2%dgsave = .FALSE. 
    2121       ENDIF 
    2122       IF  (iom_use("OCN_KWCO2")) THEN  
    2123           med_diag%OCN_KWCO2%dgsave = .TRUE. 
    2124       ELSE  
    2125           med_diag%OCN_KWCO2%dgsave = .FALSE. 
    2126       ENDIF 
    2127       IF  (iom_use("OCN_K0")) THEN  
    2128           med_diag%OCN_K0%dgsave = .TRUE. 
    2129       ELSE  
    2130           med_diag%OCN_K0%dgsave = .FALSE. 
    2131       ENDIF 
    2132       IF  (iom_use("CO2STARAIR")) THEN  
    2133           med_diag%CO2STARAIR%dgsave = .TRUE. 
    2134       ELSE  
    2135           med_diag%CO2STARAIR%dgsave = .FALSE. 
    2136       ENDIF 
    2137       IF  (iom_use("OCN_DPCO2")) THEN  
    2138           med_diag%OCN_DPCO2%dgsave = .TRUE. 
    2139       ELSE  
    2140           med_diag%OCN_DPCO2%dgsave = .FALSE. 
     1057         med_diag%INVTN%dgsave = .TRUE. 
     1058      ELSE  
     1059         med_diag%INVTN%dgsave = .FALSE. 
     1060      ENDIF 
     1061      IF (iom_use("INVTSI")) THEN  
     1062         med_diag%INVTSI%dgsave = .TRUE. 
     1063      ELSE  
     1064         med_diag%INVTSI%dgsave = .FALSE. 
     1065      ENDIF 
     1066      IF (iom_use("INVTFE")) THEN  
     1067         med_diag%INVTFE%dgsave = .TRUE. 
     1068      ELSE  
     1069         med_diag%INVTFE%dgsave = .FALSE. 
     1070      ENDIF 
     1071      IF (iom_use("PRN")) THEN  
     1072         med_diag%PRN%dgsave = .TRUE. 
     1073      ELSE  
     1074         med_diag%PRN%dgsave = .FALSE. 
     1075      ENDIF 
     1076      IF (iom_use("MPN")) THEN  
     1077         med_diag%MPN%dgsave = .TRUE. 
     1078      ELSE  
     1079         med_diag%MPN%dgsave = .FALSE. 
     1080      ENDIF 
     1081      IF (iom_use("PRD")) THEN  
     1082         med_diag%PRD%dgsave = .TRUE. 
     1083      ELSE  
     1084         med_diag%PRD%dgsave = .FALSE. 
     1085      ENDIF 
     1086      IF (iom_use("MPD")) THEN  
     1087         med_diag%MPD%dgsave = .TRUE. 
     1088      ELSE  
     1089         med_diag%MPD%dgsave = .FALSE. 
     1090      ENDIF 
     1091      IF (iom_use("DSED")) THEN  
     1092         med_diag%DSED%dgsave = .TRUE. 
     1093      ELSE  
     1094         med_diag%DSED%dgsave = .FALSE. 
     1095      ENDIF 
     1096      IF (iom_use("OPAL")) THEN  
     1097         med_diag%OPAL%dgsave = .TRUE. 
     1098      ELSE  
     1099         med_diag%OPAL%dgsave = .FALSE. 
     1100      ENDIF 
     1101      IF (iom_use("OPALDISS")) THEN  
     1102         med_diag%OPALDISS%dgsave = .TRUE. 
     1103      ELSE  
     1104         med_diag%OPALDISS%dgsave = .FALSE. 
     1105      ENDIF 
     1106      IF (iom_use("GMIPn")) THEN  
     1107         med_diag%GMIPn%dgsave = .TRUE. 
     1108      ELSE  
     1109         med_diag%GMIPn%dgsave = .FALSE. 
     1110      ENDIF 
     1111      IF (iom_use("GMID")) THEN  
     1112         med_diag%GMID%dgsave = .TRUE. 
     1113      ELSE  
     1114         med_diag%GMID%dgsave = .FALSE. 
     1115      ENDIF 
     1116      IF (iom_use("MZMI")) THEN  
     1117         med_diag%MZMI%dgsave = .TRUE. 
     1118      ELSE  
     1119         med_diag%MZMI%dgsave = .FALSE. 
     1120      ENDIF 
     1121      IF (iom_use("GMEPN")) THEN  
     1122         med_diag%GMEPN%dgsave = .TRUE. 
     1123      ELSE  
     1124         med_diag%GMEPN%dgsave = .FALSE. 
     1125      ENDIF 
     1126      IF (iom_use("GMEPD")) THEN  
     1127         med_diag%GMEPD%dgsave = .TRUE. 
     1128      ELSE  
     1129         med_diag%GMEPD%dgsave = .FALSE. 
     1130      ENDIF 
     1131      IF (iom_use("GMEZMI")) THEN  
     1132         med_diag%GMEZMI%dgsave = .TRUE. 
     1133      ELSE  
     1134         med_diag%GMEZMI%dgsave = .FALSE. 
     1135      ENDIF 
     1136      IF (iom_use("GMED")) THEN  
     1137         med_diag%GMED%dgsave = .TRUE. 
     1138      ELSE  
     1139         med_diag%GMED%dgsave = .FALSE. 
     1140      ENDIF 
     1141      IF (iom_use("MZME")) THEN  
     1142         med_diag%MZME%dgsave = .TRUE. 
     1143      ELSE  
     1144         med_diag%MZME%dgsave = .FALSE. 
     1145      ENDIF 
     1146      IF (iom_use("DEXP")) THEN  
     1147         med_diag%DEXP%dgsave = .TRUE. 
     1148      ELSE  
     1149         med_diag%DEXP%dgsave = .FALSE. 
     1150      ENDIF 
     1151      IF (iom_use("DETN")) THEN  
     1152         med_diag%DETN%dgsave = .TRUE. 
     1153      ELSE  
     1154         med_diag%DETN%dgsave = .FALSE. 
     1155      ENDIF 
     1156      IF (iom_use("MDET")) THEN  
     1157         med_diag%MDET%dgsave = .TRUE. 
     1158      ELSE  
     1159         med_diag%MDET%dgsave = .FALSE. 
     1160      ENDIF 
     1161      IF (iom_use("AEOLIAN")) THEN  
     1162         med_diag%AEOLIAN%dgsave = .TRUE. 
     1163      ELSE  
     1164         med_diag%AEOLIAN%dgsave = .FALSE. 
     1165      ENDIF 
     1166      IF (iom_use("BENTHIC")) THEN  
     1167         med_diag%BENTHIC%dgsave = .TRUE. 
     1168      ELSE  
     1169         med_diag%BENTHIC%dgsave = .FALSE. 
     1170      ENDIF 
     1171      IF (iom_use("SCAVENGE")) THEN  
     1172         med_diag%SCAVENGE%dgsave = .TRUE. 
     1173      ELSE  
     1174         med_diag%SCAVENGE%dgsave = .FALSE. 
     1175      ENDIF 
     1176      IF (iom_use("PN_JLIM")) THEN  
     1177         med_diag%PN_JLIM%dgsave = .TRUE. 
     1178      ELSE  
     1179         med_diag%PN_JLIM%dgsave = .FALSE. 
     1180      ENDIF 
     1181      IF (iom_use("PN_NLIM")) THEN  
     1182         med_diag%PN_NLIM%dgsave = .TRUE. 
     1183      ELSE  
     1184         med_diag%PN_NLIM%dgsave = .FALSE. 
     1185      ENDIF 
     1186      IF (iom_use("PN_FELIM")) THEN  
     1187         med_diag%PN_FELIM%dgsave = .TRUE. 
     1188      ELSE  
     1189         med_diag%PN_FELIM%dgsave = .FALSE. 
     1190      ENDIF 
     1191      IF (iom_use("PD_JLIM")) THEN  
     1192         med_diag%PD_JLIM%dgsave = .TRUE. 
     1193      ELSE  
     1194         med_diag%PD_JLIM%dgsave = .FALSE. 
     1195      ENDIF 
     1196      IF (iom_use("PD_NLIM")) THEN  
     1197         med_diag%PD_NLIM%dgsave = .TRUE. 
     1198      ELSE  
     1199         med_diag%PD_NLIM%dgsave = .FALSE. 
     1200      ENDIF 
     1201      IF (iom_use("PD_FELIM")) THEN  
     1202         med_diag%PD_FELIM%dgsave = .TRUE. 
     1203      ELSE  
     1204         med_diag%PD_FELIM%dgsave = .FALSE. 
     1205      ENDIF 
     1206      IF (iom_use("PD_SILIM")) THEN  
     1207         med_diag%PD_SILIM%dgsave = .TRUE. 
     1208      ELSE  
     1209         med_diag%PD_SILIM%dgsave = .FALSE. 
     1210      ENDIF 
     1211      IF (iom_use("PDSILIM2")) THEN  
     1212         med_diag%PDSILIM2%dgsave = .TRUE. 
     1213      ELSE  
     1214         med_diag%PDSILIM2%dgsave = .FALSE. 
     1215      ENDIF 
     1216      IF (iom_use("SDT__100")) THEN  
     1217         med_diag%SDT__100%dgsave = .TRUE. 
     1218      ELSE  
     1219         med_diag%SDT__100%dgsave = .FALSE. 
     1220      ENDIF 
     1221      IF (iom_use("SDT__200")) THEN  
     1222         med_diag%SDT__200%dgsave = .TRUE. 
     1223      ELSE  
     1224         med_diag%SDT__200%dgsave = .FALSE. 
     1225      ENDIF 
     1226      IF (iom_use("SDT__500")) THEN  
     1227         med_diag%SDT__500%dgsave = .TRUE. 
     1228      ELSE  
     1229         med_diag%SDT__500%dgsave = .FALSE. 
     1230      ENDIF 
     1231      IF (iom_use("SDT_1000")) THEN  
     1232         med_diag%SDT_1000%dgsave = .TRUE. 
     1233      ELSE  
     1234         med_diag%SDT_1000%dgsave = .FALSE. 
     1235      ENDIF 
     1236      IF (iom_use("TOTREG_N")) THEN  
     1237         med_diag%TOTREG_N%dgsave = .TRUE. 
     1238      ELSE  
     1239         med_diag%TOTREG_N%dgsave = .FALSE. 
     1240      ENDIF 
     1241      IF (iom_use("TOTRG_SI")) THEN  
     1242         med_diag%TOTRG_SI%dgsave = .TRUE. 
     1243      ELSE  
     1244         med_diag%TOTRG_SI%dgsave = .FALSE. 
     1245      ENDIF 
     1246      IF (iom_use("REG__100")) THEN  
     1247         med_diag%REG__100%dgsave = .TRUE. 
     1248      ELSE  
     1249         med_diag%REG__100%dgsave = .FALSE. 
     1250      ENDIF 
     1251      IF (iom_use("REG__200")) THEN  
     1252         med_diag%REG__200%dgsave = .TRUE. 
     1253      ELSE  
     1254         med_diag%REG__200%dgsave = .FALSE. 
     1255      ENDIF 
     1256      IF (iom_use("REG__500")) THEN  
     1257         med_diag%REG__500%dgsave = .TRUE. 
     1258      ELSE  
     1259         med_diag%REG__500%dgsave = .FALSE. 
     1260      ENDIF 
     1261      IF (iom_use("REG_1000")) THEN  
     1262         med_diag%REG_1000%dgsave = .TRUE. 
     1263      ELSE  
     1264         med_diag%REG_1000%dgsave = .FALSE. 
     1265      ENDIF 
     1266      IF (iom_use("FASTN")) THEN  
     1267         med_diag%FASTN%dgsave = .TRUE. 
     1268      ELSE  
     1269         med_diag%FASTN%dgsave = .FALSE. 
     1270      ENDIF 
     1271      IF (iom_use("FASTSI")) THEN  
     1272         med_diag%FASTSI%dgsave = .TRUE. 
     1273      ELSE  
     1274         med_diag%FASTSI%dgsave = .FALSE. 
     1275      ENDIF 
     1276      IF (iom_use("FASTFE")) THEN  
     1277         med_diag%FASTFE%dgsave = .TRUE. 
     1278      ELSE  
     1279         med_diag%FASTFE%dgsave = .FALSE. 
     1280      ENDIF 
     1281      IF (iom_use("FASTC")) THEN  
     1282         med_diag%FASTC%dgsave = .TRUE. 
     1283      ELSE  
     1284         med_diag%FASTC%dgsave = .FALSE. 
     1285      ENDIF 
     1286      IF (iom_use("FASTCA")) THEN  
     1287         med_diag%FASTCA%dgsave = .TRUE. 
     1288      ELSE  
     1289         med_diag%FASTCA%dgsave = .FALSE. 
     1290      ENDIF 
     1291      IF (iom_use("FDT__100")) THEN  
     1292         med_diag%FDT__100%dgsave = .TRUE. 
     1293      ELSE  
     1294         med_diag%FDT__100%dgsave = .FALSE. 
     1295      ENDIF 
     1296      IF (iom_use("FDT__200")) THEN  
     1297         med_diag%FDT__200%dgsave = .TRUE. 
     1298      ELSE  
     1299         med_diag%FDT__200%dgsave = .FALSE. 
     1300      ENDIF 
     1301      IF (iom_use("FDT__500")) THEN  
     1302         med_diag%FDT__500%dgsave = .TRUE. 
     1303      ELSE  
     1304         med_diag%FDT__500%dgsave = .FALSE. 
     1305      ENDIF 
     1306      IF (iom_use("FDT_1000")) THEN  
     1307         med_diag%FDT_1000%dgsave = .TRUE. 
     1308      ELSE  
     1309         med_diag%FDT_1000%dgsave = .FALSE. 
     1310      ENDIF 
     1311      IF (iom_use("RG__100F")) THEN  
     1312         med_diag%RG__100F%dgsave = .TRUE. 
     1313      ELSE  
     1314         med_diag%RG__100F%dgsave = .FALSE. 
     1315      ENDIF 
     1316      IF (iom_use("RG__200F")) THEN  
     1317         med_diag%RG__200F%dgsave = .TRUE. 
     1318      ELSE  
     1319         med_diag%RG__200F%dgsave = .FALSE. 
     1320      ENDIF 
     1321      IF (iom_use("RG__500F")) THEN  
     1322         med_diag%RG__500F%dgsave = .TRUE. 
     1323      ELSE  
     1324         med_diag%RG__500F%dgsave = .FALSE. 
     1325      ENDIF 
     1326      IF (iom_use("RG_1000F")) THEN  
     1327         med_diag%RG_1000F%dgsave = .TRUE. 
     1328      ELSE  
     1329         med_diag%RG_1000F%dgsave = .FALSE. 
     1330      ENDIF 
     1331      IF (iom_use("FDS__100")) THEN  
     1332         med_diag%FDS__100%dgsave = .TRUE. 
     1333      ELSE  
     1334         med_diag%FDS__100%dgsave = .FALSE. 
     1335      ENDIF 
     1336      IF (iom_use("FDS__200")) THEN  
     1337         med_diag%FDS__200%dgsave = .TRUE. 
     1338      ELSE  
     1339         med_diag%FDS__200%dgsave = .FALSE. 
     1340      ENDIF 
     1341      IF (iom_use("FDS__500")) THEN  
     1342         med_diag%FDS__500%dgsave = .TRUE. 
     1343      ELSE  
     1344         med_diag%FDS__500%dgsave = .FALSE. 
     1345      ENDIF 
     1346      IF (iom_use("FDS_1000")) THEN  
     1347         med_diag%FDS_1000%dgsave = .TRUE. 
     1348      ELSE  
     1349         med_diag%FDS_1000%dgsave = .FALSE. 
     1350      ENDIF 
     1351      IF (iom_use("RGS_100F")) THEN  
     1352         med_diag%RGS_100F%dgsave = .TRUE. 
     1353      ELSE  
     1354         med_diag%RGS_100F%dgsave = .FALSE. 
     1355      ENDIF 
     1356      IF (iom_use("RGS_200F")) THEN  
     1357         med_diag%RGS_200F%dgsave = .TRUE. 
     1358      ELSE  
     1359         med_diag%RGS_200F%dgsave = .FALSE. 
     1360      ENDIF 
     1361      IF (iom_use("RGS_500F")) THEN  
     1362         med_diag%RGS_500F%dgsave = .TRUE. 
     1363      ELSE  
     1364         med_diag%RGS_500F%dgsave = .FALSE. 
     1365      ENDIF 
     1366      IF (iom_use("RGS1000F")) THEN  
     1367         med_diag%RGS1000F%dgsave = .TRUE. 
     1368      ELSE  
     1369         med_diag%RGS1000F%dgsave = .FALSE. 
     1370      ENDIF 
     1371      IF (iom_use("REMINN")) THEN  
     1372         med_diag%REMINN%dgsave = .TRUE. 
     1373      ELSE  
     1374         med_diag%REMINN%dgsave = .FALSE. 
     1375      ENDIF 
     1376      IF (iom_use("REMINSI")) THEN  
     1377         med_diag%REMINSI%dgsave = .TRUE. 
     1378      ELSE  
     1379         med_diag%REMINSI%dgsave = .FALSE. 
     1380      ENDIF 
     1381      IF (iom_use("REMINFE")) THEN  
     1382         med_diag%REMINFE%dgsave = .TRUE. 
     1383      ELSE  
     1384         med_diag%REMINFE%dgsave = .FALSE. 
     1385      ENDIF 
     1386      IF (iom_use("REMINC")) THEN  
     1387         med_diag%REMINC%dgsave = .TRUE. 
     1388      ELSE  
     1389         med_diag%REMINC%dgsave = .FALSE. 
     1390      ENDIF 
     1391      IF (iom_use("REMINCA")) THEN  
     1392         med_diag%REMINCA%dgsave = .TRUE. 
     1393      ELSE  
     1394         med_diag%REMINCA%dgsave = .FALSE. 
     1395      ENDIF 
     1396      IF (iom_use("SEAFLRN")) THEN  
     1397         med_diag%SEAFLRN%dgsave = .TRUE. 
     1398      ELSE  
     1399         med_diag%SEAFLRN%dgsave = .FALSE. 
     1400      ENDIF 
     1401      IF (iom_use("SEAFLRSI")) THEN  
     1402         med_diag%SEAFLRSI%dgsave = .TRUE. 
     1403      ELSE  
     1404         med_diag%SEAFLRSI%dgsave = .FALSE. 
     1405      ENDIF 
     1406      IF (iom_use("SEAFLRFE")) THEN  
     1407         med_diag%SEAFLRFE%dgsave = .TRUE. 
     1408      ELSE  
     1409         med_diag%SEAFLRFE%dgsave = .FALSE. 
     1410      ENDIF 
     1411      IF (iom_use("SEAFLRC")) THEN  
     1412         med_diag%SEAFLRC%dgsave = .TRUE. 
     1413      ELSE  
     1414         med_diag%SEAFLRC%dgsave = .FALSE. 
     1415      ENDIF 
     1416      IF (iom_use("SEAFLRCA")) THEN  
     1417         med_diag%SEAFLRCA%dgsave = .TRUE. 
     1418      ELSE  
     1419         med_diag%SEAFLRCA%dgsave = .FALSE. 
     1420      ENDIF 
     1421      IF (iom_use("MED_QSR")) THEN  
     1422         med_diag%MED_QSR%dgsave = .TRUE. 
     1423      ELSE  
     1424         med_diag%MED_QSR%dgsave = .FALSE. 
     1425      ENDIF 
     1426      IF (iom_use("MED_XPAR")) THEN  
     1427         med_diag%MED_XPAR%dgsave = .TRUE. 
     1428      ELSE  
     1429         med_diag%MED_XPAR%dgsave = .FALSE. 
     1430      ENDIF 
     1431      IF (iom_use("INTFLX_N")) THEN  
     1432         med_diag%INTFLX_N%dgsave = .TRUE. 
     1433      ELSE  
     1434         med_diag%INTFLX_N%dgsave = .FALSE. 
     1435      ENDIF 
     1436      IF (iom_use("INTFLX_SI")) THEN  
     1437         med_diag%INTFLX_SI%dgsave = .TRUE. 
     1438      ELSE  
     1439         med_diag%INTFLX_SI%dgsave = .FALSE. 
     1440      ENDIF 
     1441      IF (iom_use("INTFLX_FE")) THEN  
     1442         med_diag%INTFLX_FE%dgsave = .TRUE. 
     1443      ELSE  
     1444         med_diag%INTFLX_FE%dgsave = .FALSE. 
     1445      ENDIF 
     1446      IF (iom_use("INT_PN")) THEN  
     1447         med_diag%INT_PN%dgsave = .TRUE. 
     1448      ELSE  
     1449         med_diag%INT_PN%dgsave = .FALSE. 
     1450      ENDIF 
     1451      IF (iom_use("INT_PD")) THEN  
     1452         med_diag%INT_PD%dgsave = .TRUE. 
     1453      ELSE  
     1454         med_diag%INT_PD%dgsave = .FALSE. 
     1455      ENDIF 
     1456      IF (iom_use("ML_PRN")) THEN  
     1457         med_diag%ML_PRN%dgsave = .TRUE. 
     1458      ELSE  
     1459         med_diag%ML_PRN%dgsave = .FALSE. 
     1460      ENDIF 
     1461      IF (iom_use("ML_PRD")) THEN  
     1462         med_diag%ML_PRD%dgsave = .TRUE. 
     1463      ELSE  
     1464         med_diag%ML_PRD%dgsave = .FALSE. 
     1465      ENDIF 
     1466      IF (iom_use("OCAL_CCD")) THEN  
     1467         med_diag%OCAL_CCD%dgsave = .TRUE. 
     1468      ELSE  
     1469         med_diag%OCAL_CCD%dgsave = .FALSE. 
     1470      ENDIF 
     1471      IF (iom_use("OCAL_LVL")) THEN  
     1472         med_diag%OCAL_LVL%dgsave = .TRUE. 
     1473      ELSE  
     1474         med_diag%OCAL_LVL%dgsave = .FALSE. 
     1475      ENDIF 
     1476      IF (iom_use("FE_0000")) THEN  
     1477         med_diag%FE_0000%dgsave = .TRUE. 
     1478      ELSE  
     1479         med_diag%FE_0000%dgsave = .FALSE. 
     1480      ENDIF 
     1481      IF (iom_use("FE_0100")) THEN  
     1482         med_diag%FE_0100%dgsave = .TRUE. 
     1483      ELSE  
     1484         med_diag%FE_0100%dgsave = .FALSE. 
     1485      ENDIF 
     1486      IF (iom_use("FE_0200")) THEN  
     1487         med_diag%FE_0200%dgsave = .TRUE. 
     1488      ELSE  
     1489         med_diag%FE_0200%dgsave = .FALSE. 
     1490      ENDIF 
     1491      IF (iom_use("FE_0500")) THEN  
     1492         med_diag%FE_0500%dgsave = .TRUE. 
     1493      ELSE  
     1494         med_diag%FE_0500%dgsave = .FALSE. 
     1495      ENDIF 
     1496      IF (iom_use("FE_1000")) THEN  
     1497         med_diag%FE_1000%dgsave = .TRUE. 
     1498      ELSE  
     1499         med_diag%FE_1000%dgsave = .FALSE. 
     1500      ENDIF 
     1501      IF (iom_use("MED_XZE")) THEN  
     1502         med_diag%MED_XZE%dgsave = .TRUE. 
     1503      ELSE  
     1504         med_diag%MED_XZE%dgsave = .FALSE. 
     1505      ENDIF 
     1506      IF (iom_use("WIND")) THEN  
     1507         med_diag%WIND%dgsave = .TRUE. 
     1508      ELSE  
     1509         med_diag%WIND%dgsave = .FALSE. 
     1510      ENDIF 
     1511      IF (iom_use("ATM_PCO2")) THEN  
     1512         med_diag%ATM_PCO2%dgsave = .TRUE. 
     1513      ELSE  
     1514         med_diag%ATM_PCO2%dgsave = .FALSE. 
     1515      ENDIF 
     1516      IF (iom_use("OCN_PH")) THEN  
     1517         med_diag%OCN_PH%dgsave = .TRUE. 
     1518      ELSE  
     1519         med_diag%OCN_PH%dgsave = .FALSE. 
     1520      ENDIF 
     1521      IF (iom_use("OCN_PCO2")) THEN  
     1522         med_diag%OCN_PCO2%dgsave = .TRUE. 
     1523      ELSE  
     1524         med_diag%OCN_PCO2%dgsave = .FALSE. 
     1525      ENDIF 
     1526      IF (iom_use("OCNH2CO3")) THEN  
     1527         med_diag%OCNH2CO3%dgsave = .TRUE. 
     1528      ELSE  
     1529         med_diag%OCNH2CO3%dgsave = .FALSE. 
     1530      ENDIF 
     1531      IF (iom_use("OCN_HCO3")) THEN  
     1532         med_diag%OCN_HCO3%dgsave = .TRUE. 
     1533      ELSE  
     1534         med_diag%OCN_HCO3%dgsave = .FALSE. 
     1535      ENDIF 
     1536      IF (iom_use("OCN_CO3")) THEN  
     1537         med_diag%OCN_CO3%dgsave = .TRUE. 
     1538      ELSE  
     1539         med_diag%OCN_CO3%dgsave = .FALSE. 
     1540      ENDIF 
     1541      IF (iom_use("CO2FLUX")) THEN  
     1542         med_diag%CO2FLUX%dgsave = .TRUE. 
     1543      ELSE  
     1544         med_diag%CO2FLUX%dgsave = .FALSE. 
     1545      ENDIF 
     1546      IF (iom_use("OM_CAL")) THEN  
     1547         med_diag%OM_CAL%dgsave = .TRUE. 
     1548      ELSE  
     1549         med_diag%OM_CAL%dgsave = .FALSE. 
     1550      ENDIF 
     1551      IF (iom_use("OM_ARG")) THEN  
     1552         med_diag%OM_ARG%dgsave = .TRUE. 
     1553      ELSE  
     1554         med_diag%OM_ARG%dgsave = .FALSE. 
     1555      ENDIF 
     1556      IF (iom_use("TCO2")) THEN  
     1557         med_diag%TCO2%dgsave = .TRUE. 
     1558      ELSE  
     1559         med_diag%TCO2%dgsave = .FALSE. 
     1560      ENDIF 
     1561      IF (iom_use("TALK")) THEN  
     1562         med_diag%TALK%dgsave = .TRUE. 
     1563      ELSE  
     1564         med_diag%TALK%dgsave = .FALSE. 
     1565      ENDIF 
     1566      IF (iom_use("KW660")) THEN  
     1567         med_diag%KW660%dgsave = .TRUE. 
     1568      ELSE  
     1569         med_diag%KW660%dgsave = .FALSE. 
     1570      ENDIF 
     1571      IF (iom_use("ATM_PP0")) THEN  
     1572         med_diag%ATM_PP0%dgsave = .TRUE. 
     1573      ELSE  
     1574         med_diag%ATM_PP0%dgsave = .FALSE. 
     1575      ENDIF 
     1576      IF (iom_use("O2FLUX")) THEN  
     1577         med_diag%O2FLUX%dgsave = .TRUE. 
     1578      ELSE  
     1579         med_diag%O2FLUX%dgsave = .FALSE. 
     1580      ENDIF 
     1581      IF (iom_use("O2SAT")) THEN  
     1582         med_diag%O2SAT%dgsave = .TRUE. 
     1583      ELSE  
     1584         med_diag%O2SAT%dgsave = .FALSE. 
     1585      ENDIF 
     1586      IF (iom_use("CAL_CCD")) THEN  
     1587         med_diag%CAL_CCD%dgsave = .TRUE. 
     1588      ELSE  
     1589         med_diag%CAL_CCD%dgsave = .FALSE. 
     1590      ENDIF 
     1591      IF (iom_use("ARG_CCD")) THEN  
     1592         med_diag%ARG_CCD%dgsave = .TRUE. 
     1593      ELSE  
     1594         med_diag%ARG_CCD%dgsave = .FALSE. 
     1595      ENDIF 
     1596      IF (iom_use("SFR_OCAL")) THEN  
     1597         med_diag%SFR_OCAL%dgsave = .TRUE. 
     1598      ELSE  
     1599         med_diag%SFR_OCAL%dgsave = .FALSE. 
     1600      ENDIF 
     1601      IF (iom_use("SFR_OARG")) THEN  
     1602         med_diag%SFR_OARG%dgsave = .TRUE. 
     1603      ELSE  
     1604         med_diag%SFR_OARG%dgsave = .FALSE. 
     1605      ENDIF 
     1606      IF (iom_use("N_PROD")) THEN  
     1607         med_diag%N_PROD%dgsave = .TRUE. 
     1608      ELSE  
     1609         med_diag%N_PROD%dgsave = .FALSE. 
     1610      ENDIF 
     1611      IF (iom_use("N_CONS")) THEN  
     1612         med_diag%N_CONS%dgsave = .TRUE. 
     1613      ELSE  
     1614         med_diag%N_CONS%dgsave = .FALSE. 
     1615      ENDIF 
     1616      IF (iom_use("C_PROD")) THEN  
     1617         med_diag%C_PROD%dgsave = .TRUE. 
     1618      ELSE  
     1619         med_diag%C_PROD%dgsave = .FALSE. 
     1620      ENDIF 
     1621      IF (iom_use("C_CONS")) THEN  
     1622         med_diag%C_CONS%dgsave = .TRUE. 
     1623      ELSE  
     1624         med_diag%C_CONS%dgsave = .FALSE. 
     1625      ENDIF 
     1626      IF (iom_use("O2_PROD")) THEN  
     1627         med_diag%O2_PROD%dgsave = .TRUE. 
     1628      ELSE  
     1629         med_diag%O2_PROD%dgsave = .FALSE. 
     1630      ENDIF 
     1631      IF (iom_use("O2_CONS")) THEN  
     1632         med_diag%O2_CONS%dgsave = .TRUE. 
     1633      ELSE  
     1634         med_diag%O2_CONS%dgsave = .FALSE. 
     1635      ENDIF 
     1636      IF (iom_use("O2_ANOX")) THEN  
     1637         med_diag%O2_ANOX%dgsave = .TRUE. 
     1638      ELSE  
     1639         med_diag%O2_ANOX%dgsave = .FALSE. 
     1640      ENDIF 
     1641      IF (iom_use("RR_0100")) THEN  
     1642         med_diag%RR_0100%dgsave = .TRUE. 
     1643      ELSE  
     1644         med_diag%RR_0100%dgsave = .FALSE. 
     1645      ENDIF 
     1646      IF (iom_use("RR_0500")) THEN  
     1647         med_diag%RR_0500%dgsave = .TRUE. 
     1648      ELSE  
     1649         med_diag%RR_0500%dgsave = .FALSE. 
     1650      ENDIF 
     1651      IF (iom_use("RR_1000")) THEN  
     1652         med_diag%RR_1000%dgsave = .TRUE. 
     1653      ELSE  
     1654         med_diag%RR_1000%dgsave = .FALSE. 
     1655      ENDIF 
     1656      IF (iom_use("IBEN_N")) THEN  
     1657         med_diag%IBEN_N%dgsave = .TRUE. 
     1658      ELSE  
     1659         med_diag%IBEN_N%dgsave = .FALSE. 
     1660      ENDIF 
     1661      IF (iom_use("IBEN_FE")) THEN  
     1662         med_diag%IBEN_FE%dgsave = .TRUE. 
     1663      ELSE  
     1664         med_diag%IBEN_FE%dgsave = .FALSE. 
     1665      ENDIF 
     1666      IF (iom_use("IBEN_C")) THEN  
     1667         med_diag%IBEN_C%dgsave = .TRUE. 
     1668      ELSE  
     1669         med_diag%IBEN_C%dgsave = .FALSE. 
     1670      ENDIF 
     1671      IF (iom_use("IBEN_SI")) THEN  
     1672         med_diag%IBEN_SI%dgsave = .TRUE. 
     1673      ELSE  
     1674         med_diag%IBEN_SI%dgsave = .FALSE. 
     1675      ENDIF 
     1676      IF (iom_use("IBEN_CA")) THEN  
     1677         med_diag%IBEN_CA%dgsave = .TRUE. 
     1678      ELSE  
     1679         med_diag%IBEN_CA%dgsave = .FALSE. 
     1680      ENDIF 
     1681      IF (iom_use("OBEN_N")) THEN  
     1682         med_diag%OBEN_N%dgsave = .TRUE. 
     1683      ELSE  
     1684         med_diag%OBEN_N%dgsave = .FALSE. 
     1685      ENDIF 
     1686      IF (iom_use("OBEN_FE")) THEN  
     1687         med_diag%OBEN_FE%dgsave = .TRUE. 
     1688      ELSE  
     1689         med_diag%OBEN_FE%dgsave = .FALSE. 
     1690      ENDIF 
     1691      IF (iom_use("OBEN_C")) THEN  
     1692         med_diag%OBEN_C%dgsave = .TRUE. 
     1693      ELSE  
     1694         med_diag%OBEN_C%dgsave = .FALSE. 
     1695      ENDIF 
     1696      IF (iom_use("OBEN_SI")) THEN  
     1697         med_diag%OBEN_SI%dgsave = .TRUE. 
     1698      ELSE  
     1699         med_diag%OBEN_SI%dgsave = .FALSE. 
     1700      ENDIF 
     1701      IF (iom_use("OBEN_CA")) THEN  
     1702         med_diag%OBEN_CA%dgsave = .TRUE. 
     1703      ELSE  
     1704         med_diag%OBEN_CA%dgsave = .FALSE. 
     1705      ENDIF 
     1706      IF (iom_use("BEN_N")) THEN  
     1707         med_diag%BEN_N%dgsave = .TRUE. 
     1708      ELSE  
     1709         med_diag%BEN_N%dgsave = .FALSE. 
     1710      ENDIF 
     1711      IF (iom_use("BEN_FE")) THEN  
     1712         med_diag%BEN_FE%dgsave = .TRUE. 
     1713      ELSE  
     1714         med_diag%BEN_FE%dgsave = .FALSE. 
     1715      ENDIF 
     1716      IF (iom_use("BEN_C")) THEN  
     1717         med_diag%BEN_C%dgsave = .TRUE. 
     1718      ELSE  
     1719         med_diag%BEN_C%dgsave = .FALSE. 
     1720      ENDIF 
     1721      IF (iom_use("BEN_SI")) THEN  
     1722         med_diag%BEN_SI%dgsave = .TRUE. 
     1723      ELSE  
     1724         med_diag%BEN_SI%dgsave = .FALSE. 
     1725      ENDIF 
     1726      IF (iom_use("BEN_CA")) THEN  
     1727         med_diag%BEN_CA%dgsave = .TRUE. 
     1728      ELSE  
     1729         med_diag%BEN_CA%dgsave = .FALSE. 
     1730      ENDIF 
     1731      IF (iom_use("RUNOFF")) THEN  
     1732         med_diag%RUNOFF%dgsave = .TRUE. 
     1733      ELSE  
     1734         med_diag%RUNOFF%dgsave = .FALSE. 
     1735      ENDIF 
     1736      IF (iom_use("RIV_N")) THEN  
     1737         med_diag%RIV_N%dgsave = .TRUE. 
     1738      ELSE  
     1739         med_diag%RIV_N%dgsave = .FALSE. 
     1740      ENDIF 
     1741      IF (iom_use("RIV_SI")) THEN  
     1742         med_diag%RIV_SI%dgsave = .TRUE. 
     1743      ELSE  
     1744         med_diag%RIV_SI%dgsave = .FALSE. 
     1745      ENDIF 
     1746      IF (iom_use("RIV_C")) THEN  
     1747         med_diag%RIV_C%dgsave = .TRUE. 
     1748      ELSE  
     1749         med_diag%RIV_C%dgsave = .FALSE. 
     1750      ENDIF 
     1751      IF (iom_use("RIV_ALK")) THEN  
     1752         med_diag%RIV_ALK%dgsave = .TRUE. 
     1753      ELSE  
     1754         med_diag%RIV_ALK%dgsave = .FALSE. 
     1755      ENDIF 
     1756      IF (iom_use("DETC")) THEN  
     1757         med_diag%DETC%dgsave = .TRUE. 
     1758      ELSE  
     1759         med_diag%DETC%dgsave = .FALSE. 
     1760      ENDIF 
     1761      IF (iom_use("SDC__100")) THEN  
     1762         med_diag%SDC__100%dgsave = .TRUE. 
     1763      ELSE  
     1764         med_diag%SDC__100%dgsave = .FALSE. 
     1765      ENDIF 
     1766      IF (iom_use("SDC__200")) THEN  
     1767         med_diag%SDC__200%dgsave = .TRUE. 
     1768      ELSE  
     1769         med_diag%SDC__200%dgsave = .FALSE. 
     1770      ENDIF 
     1771      IF (iom_use("SDC__500")) THEN  
     1772         med_diag%SDC__500%dgsave = .TRUE. 
     1773      ELSE  
     1774         med_diag%SDC__500%dgsave = .FALSE. 
     1775      ENDIF 
     1776      IF (iom_use("SDC_1000")) THEN  
     1777         med_diag%SDC_1000%dgsave = .TRUE. 
     1778      ELSE  
     1779         med_diag%SDC_1000%dgsave = .FALSE. 
     1780      ENDIF 
     1781      IF (iom_use("INVTC")) THEN  
     1782         med_diag%INVTC%dgsave = .TRUE. 
     1783      ELSE  
     1784         med_diag%INVTC%dgsave = .FALSE. 
     1785      ENDIF 
     1786      IF (iom_use("INVTALK")) THEN  
     1787         med_diag%INVTALK%dgsave = .TRUE. 
     1788      ELSE  
     1789         med_diag%INVTALK%dgsave = .FALSE. 
     1790      ENDIF 
     1791      IF (iom_use("INVTO2")) THEN  
     1792         med_diag%INVTO2%dgsave = .TRUE. 
     1793      ELSE  
     1794         med_diag%INVTO2%dgsave = .FALSE. 
     1795      ENDIF 
     1796      IF (iom_use("LYSO_CA")) THEN  
     1797         med_diag%LYSO_CA%dgsave = .TRUE. 
     1798      ELSE  
     1799         med_diag%LYSO_CA%dgsave = .FALSE. 
     1800      ENDIF 
     1801      IF (iom_use("COM_RESP")) THEN  
     1802         med_diag%COM_RESP%dgsave = .TRUE. 
     1803      ELSE  
     1804         med_diag%COM_RESP%dgsave = .FALSE. 
     1805      ENDIF 
     1806      IF (iom_use("PN_LLOSS")) THEN  
     1807         med_diag%PN_LLOSS%dgsave = .TRUE. 
     1808      ELSE  
     1809         med_diag%PN_LLOSS%dgsave = .FALSE. 
     1810      ENDIF 
     1811      IF (iom_use("PD_LLOSS")) THEN  
     1812         med_diag%PD_LLOSS%dgsave = .TRUE. 
     1813      ELSE  
     1814         med_diag%PD_LLOSS%dgsave = .FALSE. 
     1815      ENDIF 
     1816      IF (iom_use("ZI_LLOSS")) THEN  
     1817         med_diag%ZI_LLOSS%dgsave = .TRUE. 
     1818      ELSE  
     1819         med_diag%ZI_LLOSS%dgsave = .FALSE. 
     1820      ENDIF 
     1821      IF (iom_use("ZE_LLOSS")) THEN  
     1822         med_diag%ZE_LLOSS%dgsave = .TRUE. 
     1823      ELSE  
     1824         med_diag%ZE_LLOSS%dgsave = .FALSE. 
     1825      ENDIF 
     1826      IF (iom_use("ZI_MES_N")) THEN  
     1827         med_diag%ZI_MES_N%dgsave = .TRUE. 
     1828      ELSE  
     1829         med_diag%ZI_MES_N%dgsave = .FALSE. 
     1830      ENDIF 
     1831      IF (iom_use("ZI_MES_D")) THEN  
     1832         med_diag%ZI_MES_D%dgsave = .TRUE. 
     1833      ELSE  
     1834         med_diag%ZI_MES_D%dgsave = .FALSE. 
     1835      ENDIF 
     1836      IF (iom_use("ZI_MES_C")) THEN  
     1837         med_diag%ZI_MES_C%dgsave = .TRUE. 
     1838      ELSE  
     1839         med_diag%ZI_MES_C%dgsave = .FALSE. 
     1840      ENDIF 
     1841      IF (iom_use("ZI_MESDC")) THEN  
     1842         med_diag%ZI_MESDC%dgsave = .TRUE. 
     1843      ELSE  
     1844         med_diag%ZI_MESDC%dgsave = .FALSE. 
     1845      ENDIF 
     1846      IF (iom_use("ZI_EXCR")) THEN  
     1847         med_diag%ZI_EXCR%dgsave = .TRUE. 
     1848      ELSE  
     1849         med_diag%ZI_EXCR%dgsave = .FALSE. 
     1850      ENDIF 
     1851      IF (iom_use("ZI_RESP")) THEN  
     1852         med_diag%ZI_RESP%dgsave = .TRUE. 
     1853      ELSE  
     1854         med_diag%ZI_RESP%dgsave = .FALSE. 
     1855      ENDIF 
     1856      IF (iom_use("ZI_GROW")) THEN  
     1857         med_diag%ZI_GROW%dgsave = .TRUE. 
     1858      ELSE  
     1859         med_diag%ZI_GROW%dgsave = .FALSE. 
     1860      ENDIF 
     1861      IF (iom_use("ZE_MES_N")) THEN  
     1862         med_diag%ZE_MES_N%dgsave = .TRUE. 
     1863      ELSE  
     1864         med_diag%ZE_MES_N%dgsave = .FALSE. 
     1865      ENDIF 
     1866      IF (iom_use("ZE_MES_D")) THEN  
     1867         med_diag%ZE_MES_D%dgsave = .TRUE. 
     1868      ELSE  
     1869         med_diag%ZE_MES_D%dgsave = .FALSE. 
     1870      ENDIF 
     1871      IF (iom_use("ZE_MES_C")) THEN  
     1872         med_diag%ZE_MES_C%dgsave = .TRUE. 
     1873      ELSE  
     1874         med_diag%ZE_MES_C%dgsave = .FALSE. 
     1875      ENDIF 
     1876      IF (iom_use("ZE_MESDC")) THEN  
     1877         med_diag%ZE_MESDC%dgsave = .TRUE. 
     1878      ELSE  
     1879         med_diag%ZE_MESDC%dgsave = .FALSE. 
     1880      ENDIF 
     1881      IF (iom_use("ZE_EXCR")) THEN  
     1882         med_diag%ZE_EXCR%dgsave = .TRUE. 
     1883      ELSE  
     1884         med_diag%ZE_EXCR%dgsave = .FALSE. 
     1885      ENDIF 
     1886      IF (iom_use("ZE_RESP")) THEN  
     1887         med_diag%ZE_RESP%dgsave = .TRUE. 
     1888      ELSE  
     1889         med_diag%ZE_RESP%dgsave = .FALSE. 
     1890      ENDIF 
     1891      IF (iom_use("ZE_GROW")) THEN  
     1892         med_diag%ZE_GROW%dgsave = .TRUE. 
     1893      ELSE  
     1894         med_diag%ZE_GROW%dgsave = .FALSE. 
     1895      ENDIF 
     1896      IF (iom_use("MDETC")) THEN  
     1897         med_diag%MDETC%dgsave = .TRUE. 
     1898      ELSE  
     1899         med_diag%MDETC%dgsave = .FALSE. 
     1900      ENDIF 
     1901      IF (iom_use("GMIDC")) THEN  
     1902         med_diag%GMIDC%dgsave = .TRUE. 
     1903      ELSE  
     1904         med_diag%GMIDC%dgsave = .FALSE. 
     1905      ENDIF 
     1906      IF (iom_use("GMEDC")) THEN  
     1907         med_diag%GMEDC%dgsave = .TRUE. 
     1908      ELSE  
     1909         med_diag%GMEDC%dgsave = .FALSE. 
     1910      ENDIF 
     1911      IF (iom_use("INT_ZMI")) THEN  
     1912         med_diag%INT_ZMI%dgsave = .TRUE. 
     1913      ELSE  
     1914         med_diag%INT_ZMI%dgsave = .FALSE. 
     1915      ENDIF 
     1916      IF (iom_use("INT_ZME")) THEN  
     1917         med_diag%INT_ZME%dgsave = .TRUE. 
     1918      ELSE  
     1919         med_diag%INT_ZME%dgsave = .FALSE. 
     1920      ENDIF 
     1921      IF (iom_use("INT_DET")) THEN  
     1922         med_diag%INT_DET%dgsave = .TRUE. 
     1923      ELSE  
     1924         med_diag%INT_DET%dgsave = .FALSE. 
     1925      ENDIF 
     1926      IF (iom_use("INT_DTC")) THEN  
     1927         med_diag%INT_DTC%dgsave = .TRUE. 
     1928      ELSE  
     1929         med_diag%INT_DTC%dgsave = .FALSE. 
     1930      ENDIF 
     1931      IF (iom_use("DMS_SURF")) THEN  
     1932         med_diag%DMS_SURF%dgsave = .TRUE. 
     1933      ELSE  
     1934         med_diag%DMS_SURF%dgsave = .FALSE. 
     1935      ENDIF 
     1936      IF (iom_use("DMS_ANDR")) THEN  
     1937         med_diag%DMS_ANDR%dgsave = .TRUE. 
     1938      ELSE  
     1939         med_diag%DMS_ANDR%dgsave = .FALSE. 
     1940      ENDIF 
     1941      IF (iom_use("DMS_SIMO")) THEN  
     1942         med_diag%DMS_SIMO%dgsave = .TRUE. 
     1943      ELSE  
     1944         med_diag%DMS_SIMO%dgsave = .FALSE. 
     1945      ENDIF 
     1946      IF (iom_use("DMS_ARAN")) THEN  
     1947         med_diag%DMS_ARAN%dgsave = .TRUE. 
     1948      ELSE  
     1949         med_diag%DMS_ARAN%dgsave = .FALSE. 
     1950      ENDIF 
     1951      IF (iom_use("DMS_HALL")) THEN  
     1952         med_diag%DMS_HALL%dgsave = .TRUE. 
     1953      ELSE  
     1954         med_diag%DMS_HALL%dgsave = .FALSE. 
     1955      ENDIF 
     1956      IF (iom_use("ATM_XCO2")) THEN  
     1957         med_diag%ATM_XCO2%dgsave = .TRUE. 
     1958      ELSE  
     1959         med_diag%ATM_XCO2%dgsave = .FALSE. 
     1960      ENDIF 
     1961      IF (iom_use("OCN_FCO2")) THEN  
     1962         med_diag%OCN_FCO2%dgsave = .TRUE. 
     1963      ELSE  
     1964         med_diag%OCN_FCO2%dgsave = .FALSE. 
     1965      ENDIF 
     1966      IF (iom_use("ATM_FCO2")) THEN  
     1967         med_diag%ATM_FCO2%dgsave = .TRUE. 
     1968      ELSE  
     1969         med_diag%ATM_FCO2%dgsave = .FALSE. 
     1970      ENDIF 
     1971      IF (iom_use("OCN_RHOSW")) THEN  
     1972         med_diag%OCN_RHOSW%dgsave = .TRUE. 
     1973      ELSE  
     1974         med_diag%OCN_RHOSW%dgsave = .FALSE. 
     1975      ENDIF 
     1976      IF (iom_use("OCN_SCHCO2")) THEN  
     1977         med_diag%OCN_SCHCO2%dgsave = .TRUE. 
     1978      ELSE  
     1979         med_diag%OCN_SCHCO2%dgsave = .FALSE. 
     1980      ENDIF 
     1981      IF (iom_use("OCN_KWCO2")) THEN  
     1982         med_diag%OCN_KWCO2%dgsave = .TRUE. 
     1983      ELSE  
     1984         med_diag%OCN_KWCO2%dgsave = .FALSE. 
     1985      ENDIF 
     1986      IF (iom_use("OCN_K0")) THEN  
     1987         med_diag%OCN_K0%dgsave = .TRUE. 
     1988      ELSE  
     1989         med_diag%OCN_K0%dgsave = .FALSE. 
     1990      ENDIF 
     1991      IF (iom_use("CO2STARAIR")) THEN  
     1992         med_diag%CO2STARAIR%dgsave = .TRUE. 
     1993      ELSE  
     1994         med_diag%CO2STARAIR%dgsave = .FALSE. 
     1995      ENDIF 
     1996      IF (iom_use("OCN_DPCO2")) THEN  
     1997         med_diag%OCN_DPCO2%dgsave = .TRUE. 
     1998      ELSE  
     1999         med_diag%OCN_DPCO2%dgsave = .FALSE. 
    21412000      ENDIF 
    21422001      !! 
    2143       IF  (iom_use("TPP3")) THEN  
    2144           med_diag%TPP3%dgsave = .TRUE. 
    2145       ELSE  
    2146           med_diag%TPP3%dgsave = .FALSE. 
    2147       ENDIF 
    2148       IF  (iom_use("DETFLUX3")) THEN  
    2149           med_diag%DETFLUX3%dgsave = .TRUE. 
    2150       ELSE  
    2151           med_diag%DETFLUX3%dgsave = .FALSE. 
    2152       ENDIF 
    2153       IF  (iom_use("REMIN3N")) THEN  
    2154           med_diag%REMIN3N%dgsave = .TRUE. 
    2155       ELSE  
    2156           med_diag%REMIN3N%dgsave = .FALSE. 
    2157       ENDIF 
    2158       IF  (iom_use("PH3")) THEN  
    2159           med_diag%PH3%dgsave = .TRUE. 
    2160       ELSE  
    2161           med_diag%PH3%dgsave = .FALSE. 
    2162       ENDIF 
    2163       IF  (iom_use("OM_CAL3")) THEN  
    2164           med_diag%OM_CAL3%dgsave = .TRUE. 
    2165       ELSE  
    2166           med_diag%OM_CAL3%dgsave = .FALSE. 
     2002      !! 3D fields 
     2003      IF (iom_use("TPP3")) THEN  
     2004         med_diag%TPP3%dgsave = .TRUE. 
     2005      ELSE  
     2006         med_diag%TPP3%dgsave = .FALSE. 
     2007      ENDIF 
     2008      IF (iom_use("DETFLUX3")) THEN  
     2009         med_diag%DETFLUX3%dgsave = .TRUE. 
     2010      ELSE  
     2011         med_diag%DETFLUX3%dgsave = .FALSE. 
     2012      ENDIF 
     2013      IF (iom_use("REMIN3N")) THEN  
     2014         med_diag%REMIN3N%dgsave = .TRUE. 
     2015      ELSE  
     2016         med_diag%REMIN3N%dgsave = .FALSE. 
     2017      ENDIF 
     2018      IF (iom_use("PH3")) THEN  
     2019         med_diag%PH3%dgsave = .TRUE. 
     2020      ELSE  
     2021         med_diag%PH3%dgsave = .FALSE. 
     2022      ENDIF 
     2023      IF (iom_use("OM_CAL3")) THEN  
     2024         med_diag%OM_CAL3%dgsave = .TRUE. 
     2025      ELSE  
     2026         med_diag%OM_CAL3%dgsave = .FALSE. 
    21672027      ENDIF 
    21682028      !! 
     
    21722032      !! 
    21732033      !! 2D fields 
    2174       IF  (iom_use("epC100")) THEN  
    2175           med_diag%epC100%dgsave = .TRUE. 
    2176       ELSE  
    2177           med_diag%epC100%dgsave = .FALSE. 
    2178       ENDIF 
    2179       IF  (iom_use("epCALC100")) THEN  
    2180           med_diag%epCALC100%dgsave = .TRUE. 
    2181       ELSE  
    2182           med_diag%epCALC100%dgsave = .FALSE. 
    2183       ENDIF 
    2184       IF  (iom_use("epN100")) THEN  
    2185           med_diag%epN100%dgsave = .TRUE. 
    2186       ELSE  
    2187           med_diag%epN100%dgsave = .FALSE. 
    2188       ENDIF 
    2189       IF  (iom_use("epSI100")) THEN  
    2190           med_diag%epSI100%dgsave = .TRUE. 
    2191       ELSE  
    2192           med_diag%epSI100%dgsave = .FALSE. 
    2193       ENDIF 
    2194       IF  (iom_use("FGCO2")) THEN  
    2195           med_diag%FGCO2%dgsave = .TRUE. 
    2196       ELSE  
    2197           med_diag%FGCO2%dgsave = .FALSE. 
    2198       ENDIF 
    2199       IF  (iom_use("INTDISSIC")) THEN  
    2200           med_diag%INTDISSIC%dgsave = .TRUE. 
    2201       ELSE  
    2202           med_diag%INTDISSIC%dgsave = .FALSE. 
    2203       ENDIF 
    2204       IF  (iom_use("INTDISSIN")) THEN  
    2205           med_diag%INTDISSIN%dgsave = .TRUE. 
    2206       ELSE  
    2207           med_diag%INTDISSIN%dgsave = .FALSE. 
    2208       ENDIF 
    2209       IF  (iom_use("INTDISSISI")) THEN  
    2210           med_diag%INTDISSISI%dgsave = .TRUE. 
    2211       ELSE  
    2212           med_diag%INTDISSISI%dgsave = .FALSE. 
    2213       ENDIF 
    2214       IF  (iom_use("INTTALK")) THEN  
    2215           med_diag%INTTALK%dgsave = .TRUE. 
    2216       ELSE  
    2217           med_diag%INTTALK%dgsave = .FALSE. 
    2218       ENDIF 
    2219       IF  (iom_use("O2min")) THEN  
    2220           med_diag%O2min%dgsave = .TRUE. 
    2221       ELSE  
    2222           med_diag%O2min%dgsave = .FALSE. 
    2223       ENDIF 
    2224       IF  (iom_use("ZO2min")) THEN  
    2225           med_diag%ZO2min%dgsave = .TRUE. 
    2226       ELSE  
    2227           med_diag%ZO2min%dgsave = .FALSE. 
    2228       ENDIF 
    2229       IF  (iom_use("FBDDTALK")) THEN  
    2230           med_diag%FBDDTALK%dgsave = .TRUE. 
    2231       ELSE  
    2232           med_diag%FBDDTALK%dgsave = .FALSE. 
    2233       ENDIF 
    2234       IF  (iom_use("FBDDTDIC")) THEN  
    2235           med_diag%FBDDTDIC%dgsave = .TRUE. 
    2236       ELSE  
    2237           med_diag%FBDDTDIC%dgsave = .FALSE. 
    2238       ENDIF 
    2239       IF  (iom_use("FBDDTDIFE")) THEN  
    2240           med_diag%FBDDTDIFE%dgsave = .TRUE. 
    2241       ELSE  
    2242           med_diag%FBDDTDIFE%dgsave = .FALSE. 
    2243       ENDIF 
    2244       IF  (iom_use("FBDDTDIN")) THEN  
    2245           med_diag%FBDDTDIN%dgsave = .TRUE. 
    2246       ELSE  
    2247           med_diag%FBDDTDIN%dgsave = .FALSE. 
    2248       ENDIF 
    2249       IF  (iom_use("FBDDTDISI")) THEN  
    2250           med_diag%FBDDTDISI%dgsave = .TRUE. 
    2251       ELSE  
    2252           med_diag%FBDDTDISI%dgsave = .FALSE. 
     2034      IF (iom_use("epC100")) THEN  
     2035         med_diag%epC100%dgsave = .TRUE. 
     2036      ELSE  
     2037         med_diag%epC100%dgsave = .FALSE. 
     2038      ENDIF 
     2039      IF (iom_use("epCALC100")) THEN  
     2040         med_diag%epCALC100%dgsave = .TRUE. 
     2041      ELSE  
     2042         med_diag%epCALC100%dgsave = .FALSE. 
     2043      ENDIF 
     2044      IF (iom_use("epN100")) THEN  
     2045         med_diag%epN100%dgsave = .TRUE. 
     2046      ELSE  
     2047         med_diag%epN100%dgsave = .FALSE. 
     2048      ENDIF 
     2049      IF (iom_use("epSI100")) THEN  
     2050         med_diag%epSI100%dgsave = .TRUE. 
     2051      ELSE  
     2052         med_diag%epSI100%dgsave = .FALSE. 
     2053      ENDIF 
     2054      IF (iom_use("FGCO2")) THEN  
     2055         med_diag%FGCO2%dgsave = .TRUE. 
     2056      ELSE  
     2057         med_diag%FGCO2%dgsave = .FALSE. 
     2058      ENDIF 
     2059      IF (iom_use("INTDISSIC")) THEN  
     2060         med_diag%INTDISSIC%dgsave = .TRUE. 
     2061      ELSE  
     2062         med_diag%INTDISSIC%dgsave = .FALSE. 
     2063      ENDIF 
     2064      IF (iom_use("INTDISSIN")) THEN  
     2065         med_diag%INTDISSIN%dgsave = .TRUE. 
     2066      ELSE  
     2067         med_diag%INTDISSIN%dgsave = .FALSE. 
     2068      ENDIF 
     2069      IF (iom_use("INTDISSISI")) THEN  
     2070         med_diag%INTDISSISI%dgsave = .TRUE. 
     2071      ELSE  
     2072         med_diag%INTDISSISI%dgsave = .FALSE. 
     2073      ENDIF 
     2074      IF (iom_use("INTTALK")) THEN  
     2075         med_diag%INTTALK%dgsave = .TRUE. 
     2076      ELSE  
     2077         med_diag%INTTALK%dgsave = .FALSE. 
     2078      ENDIF 
     2079      IF (iom_use("O2min")) THEN  
     2080         med_diag%O2min%dgsave = .TRUE. 
     2081      ELSE  
     2082         med_diag%O2min%dgsave = .FALSE. 
     2083      ENDIF 
     2084      IF (iom_use("ZO2min")) THEN  
     2085         med_diag%ZO2min%dgsave = .TRUE. 
     2086      ELSE  
     2087         med_diag%ZO2min%dgsave = .FALSE. 
     2088      ENDIF 
     2089      IF (iom_use("FBDDTALK")) THEN  
     2090         med_diag%FBDDTALK%dgsave = .TRUE. 
     2091      ELSE  
     2092         med_diag%FBDDTALK%dgsave = .FALSE. 
     2093      ENDIF 
     2094      IF (iom_use("FBDDTDIC")) THEN  
     2095         med_diag%FBDDTDIC%dgsave = .TRUE. 
     2096      ELSE  
     2097         med_diag%FBDDTDIC%dgsave = .FALSE. 
     2098      ENDIF 
     2099      IF (iom_use("FBDDTDIFE")) THEN  
     2100         med_diag%FBDDTDIFE%dgsave = .TRUE. 
     2101      ELSE  
     2102         med_diag%FBDDTDIFE%dgsave = .FALSE. 
     2103      ENDIF 
     2104      IF (iom_use("FBDDTDIN")) THEN  
     2105         med_diag%FBDDTDIN%dgsave = .TRUE. 
     2106      ELSE  
     2107         med_diag%FBDDTDIN%dgsave = .FALSE. 
     2108      ENDIF 
     2109      IF (iom_use("FBDDTDISI")) THEN  
     2110         med_diag%FBDDTDISI%dgsave = .TRUE. 
     2111      ELSE  
     2112         med_diag%FBDDTDISI%dgsave = .FALSE. 
    22532113      ENDIF 
    22542114      !! 
    2255       !! 3D 
    2256       IF  (iom_use("TPPD3")) THEN  
    2257           med_diag%TPPD3%dgsave = .TRUE. 
    2258       ELSE  
    2259           med_diag%TPPD3%dgsave = .FALSE. 
    2260       ENDIF 
    2261       IF  (iom_use("BDDTALK3")) THEN  
    2262           med_diag%BDDTALK3%dgsave = .TRUE. 
    2263       ELSE  
    2264           med_diag%BDDTALK3%dgsave = .FALSE. 
    2265       ENDIF 
    2266       IF  (iom_use("BDDTDIC3")) THEN  
    2267           med_diag%BDDTDIC3%dgsave = .TRUE. 
    2268       ELSE  
    2269           med_diag%BDDTDIC3%dgsave = .FALSE. 
    2270       ENDIF 
    2271       IF  (iom_use("BDDTDIFE3")) THEN  
    2272           med_diag%BDDTDIFE3%dgsave = .TRUE. 
    2273       ELSE  
    2274           med_diag%BDDTDIFE3%dgsave = .FALSE. 
    2275       ENDIF 
    2276       IF  (iom_use("BDDTDIN3")) THEN  
    2277           med_diag%BDDTDIN3%dgsave = .TRUE. 
    2278       ELSE  
    2279           med_diag%BDDTDIN3%dgsave = .FALSE. 
    2280       ENDIF 
    2281       IF  (iom_use("BDDTDISI3")) THEN  
    2282           med_diag%BDDTDISI3%dgsave = .TRUE. 
    2283       ELSE  
    2284           med_diag%BDDTDISI3%dgsave = .FALSE. 
    2285       ENDIF 
    2286       IF  (iom_use("FD_NIT3")) THEN  
    2287           med_diag%FD_NIT3%dgsave = .TRUE. 
    2288       ELSE  
    2289           med_diag%FD_NIT3%dgsave = .FALSE. 
    2290       ENDIF 
    2291       IF  (iom_use("FD_SIL3")) THEN  
    2292           med_diag%FD_SIL3%dgsave = .TRUE. 
    2293       ELSE  
    2294           med_diag%FD_SIL3%dgsave = .FALSE. 
    2295       ENDIF 
    2296       IF  (iom_use("FD_CAR3")) THEN  
    2297           med_diag%FD_CAR3%dgsave = .TRUE. 
    2298       ELSE  
    2299           med_diag%FD_CAR3%dgsave = .FALSE. 
    2300       ENDIF 
    2301       IF  (iom_use("FD_CAL3")) THEN  
    2302           med_diag%FD_CAL3%dgsave = .TRUE. 
    2303       ELSE  
    2304           med_diag%FD_CAL3%dgsave = .FALSE. 
    2305       ENDIF 
    2306       IF  (iom_use("CO33")) THEN  
    2307           med_diag%CO33%dgsave = .TRUE. 
    2308       ELSE  
    2309           med_diag%CO33%dgsave = .FALSE. 
    2310       ENDIF 
    2311       IF  (iom_use("CO3SATARAG3")) THEN  
    2312           med_diag%CO3SATARAG3%dgsave = .TRUE. 
    2313       ELSE  
    2314           med_diag%CO3SATARAG3%dgsave = .FALSE. 
    2315       ENDIF 
    2316       IF  (iom_use("CO3SATCALC3")) THEN  
    2317           med_diag%CO3SATCALC3%dgsave = .TRUE. 
    2318       ELSE  
    2319           med_diag%CO3SATCALC3%dgsave = .FALSE. 
    2320       ENDIF 
    2321       IF  (iom_use("DCALC3")) THEN  
    2322           med_diag%DCALC3%dgsave = .TRUE. 
    2323       ELSE  
    2324           med_diag%DCALC3%dgsave = .FALSE. 
    2325       ENDIF 
    2326       IF  (iom_use("EXPC3")) THEN  
    2327           med_diag%EXPC3%dgsave = .TRUE. 
    2328       ELSE  
    2329           med_diag%EXPC3%dgsave = .FALSE. 
    2330       ENDIF 
    2331       IF  (iom_use("EXPN3")) THEN  
    2332           med_diag%EXPN3%dgsave = .TRUE. 
    2333       ELSE  
    2334           med_diag%EXPN3%dgsave = .FALSE. 
    2335       ENDIF 
    2336       IF  (iom_use("FEDISS3")) THEN  
    2337           med_diag%FEDISS3%dgsave = .TRUE. 
    2338       ELSE  
    2339           med_diag%FEDISS3%dgsave = .FALSE. 
    2340       ENDIF 
    2341       IF  (iom_use("FESCAV3")) THEN  
    2342           med_diag%FESCAV3%dgsave = .TRUE. 
    2343       ELSE  
    2344           med_diag%FESCAV3%dgsave = .FALSE. 
    2345       ENDIF 
    2346       IF  (iom_use("MIGRAZP3")) THEN  
    2347           med_diag%MIGRAZP3%dgsave = .TRUE. 
    2348       ELSE  
    2349           med_diag%MIGRAZP3%dgsave = .FALSE. 
    2350       ENDIF 
    2351       IF  (iom_use("MIGRAZD3")) THEN  
    2352           med_diag%MIGRAZD3%dgsave = .TRUE. 
    2353       ELSE  
    2354           med_diag%MIGRAZD3%dgsave = .FALSE. 
    2355       ENDIF 
    2356       IF  (iom_use("MEGRAZP3")) THEN  
    2357           med_diag%MEGRAZP3%dgsave = .TRUE. 
    2358       ELSE  
    2359           med_diag%MEGRAZP3%dgsave = .FALSE. 
    2360       ENDIF 
    2361       IF  (iom_use("MEGRAZD3")) THEN  
    2362           med_diag%MEGRAZD3%dgsave = .TRUE. 
    2363       ELSE  
    2364           med_diag%MEGRAZD3%dgsave = .FALSE. 
    2365       ENDIF 
    2366       IF  (iom_use("MEGRAZZ3")) THEN  
    2367           med_diag%MEGRAZZ3%dgsave = .TRUE. 
    2368       ELSE  
    2369           med_diag%MEGRAZZ3%dgsave = .FALSE. 
    2370       ENDIF 
    2371       IF  (iom_use("O2SAT3")) THEN  
    2372           med_diag%O2SAT3%dgsave = .TRUE. 
    2373       ELSE  
    2374           med_diag%O2SAT3%dgsave = .FALSE. 
    2375       ENDIF 
    2376       IF  (iom_use("PBSI3")) THEN  
    2377           med_diag%PBSI3%dgsave = .TRUE. 
    2378       ELSE  
    2379           med_diag%PBSI3%dgsave = .FALSE. 
    2380       ENDIF 
    2381       IF  (iom_use("PCAL3")) THEN  
    2382           med_diag%PCAL3%dgsave = .TRUE. 
    2383       ELSE  
    2384           med_diag%PCAL3%dgsave = .FALSE. 
    2385       ENDIF 
    2386       IF  (iom_use("REMOC3")) THEN  
    2387           med_diag%REMOC3%dgsave = .TRUE. 
    2388       ELSE  
    2389           med_diag%REMOC3%dgsave = .FALSE. 
    2390       ENDIF 
    2391       IF  (iom_use("PNLIMJ3")) THEN  
    2392           med_diag%PNLIMJ3%dgsave = .TRUE. 
    2393       ELSE  
    2394           med_diag%PNLIMJ3%dgsave = .FALSE. 
    2395       ENDIF 
    2396       IF  (iom_use("PNLIMN3")) THEN  
    2397           med_diag%PNLIMN3%dgsave = .TRUE. 
    2398       ELSE  
    2399           med_diag%PNLIMN3%dgsave = .FALSE. 
    2400       ENDIF 
    2401       IF  (iom_use("PNLIMFE3")) THEN  
    2402           med_diag%PNLIMFE3%dgsave = .TRUE. 
    2403       ELSE  
    2404           med_diag%PNLIMFE3%dgsave = .FALSE. 
    2405       ENDIF 
    2406       IF  (iom_use("PDLIMJ3")) THEN  
    2407           med_diag%PDLIMJ3%dgsave = .TRUE. 
    2408       ELSE  
    2409           med_diag%PDLIMJ3%dgsave = .FALSE. 
    2410       ENDIF 
    2411       IF  (iom_use("PDLIMN3")) THEN  
    2412           med_diag%PDLIMN3%dgsave = .TRUE. 
    2413       ELSE  
    2414           med_diag%PDLIMN3%dgsave = .FALSE. 
    2415       ENDIF 
    2416       IF  (iom_use("PDLIMFE3")) THEN  
    2417           med_diag%PDLIMFE3%dgsave = .TRUE. 
    2418       ELSE  
    2419           med_diag%PDLIMFE3%dgsave = .FALSE. 
    2420       ENDIF 
    2421       IF  (iom_use("PDLIMSI3")) THEN  
    2422           med_diag%PDLIMSI3%dgsave = .TRUE. 
    2423       ELSE  
    2424           med_diag%PDLIMSI3%dgsave = .FALSE. 
    2425       ENDIF 
    2426  
     2115      !! 3D fields 
     2116      IF (iom_use("TPPD3")) THEN  
     2117         med_diag%TPPD3%dgsave = .TRUE. 
     2118      ELSE  
     2119         med_diag%TPPD3%dgsave = .FALSE. 
     2120      ENDIF 
     2121      IF (iom_use("BDDTALK3")) THEN  
     2122         med_diag%BDDTALK3%dgsave = .TRUE. 
     2123      ELSE  
     2124         med_diag%BDDTALK3%dgsave = .FALSE. 
     2125      ENDIF 
     2126      IF (iom_use("BDDTDIC3")) THEN  
     2127         med_diag%BDDTDIC3%dgsave = .TRUE. 
     2128      ELSE  
     2129         med_diag%BDDTDIC3%dgsave = .FALSE. 
     2130      ENDIF 
     2131      IF (iom_use("BDDTDIFE3")) THEN  
     2132         med_diag%BDDTDIFE3%dgsave = .TRUE. 
     2133      ELSE  
     2134         med_diag%BDDTDIFE3%dgsave = .FALSE. 
     2135      ENDIF 
     2136      IF (iom_use("BDDTDIN3")) THEN  
     2137         med_diag%BDDTDIN3%dgsave = .TRUE. 
     2138      ELSE  
     2139         med_diag%BDDTDIN3%dgsave = .FALSE. 
     2140      ENDIF 
     2141      IF (iom_use("BDDTDISI3")) THEN  
     2142         med_diag%BDDTDISI3%dgsave = .TRUE. 
     2143      ELSE  
     2144         med_diag%BDDTDISI3%dgsave = .FALSE. 
     2145      ENDIF 
     2146      IF (iom_use("FD_NIT3")) THEN  
     2147         med_diag%FD_NIT3%dgsave = .TRUE. 
     2148      ELSE  
     2149         med_diag%FD_NIT3%dgsave = .FALSE. 
     2150      ENDIF 
     2151      IF (iom_use("FD_SIL3")) THEN  
     2152         med_diag%FD_SIL3%dgsave = .TRUE. 
     2153      ELSE  
     2154         med_diag%FD_SIL3%dgsave = .FALSE. 
     2155      ENDIF 
     2156      IF (iom_use("FD_CAR3")) THEN  
     2157         med_diag%FD_CAR3%dgsave = .TRUE. 
     2158      ELSE  
     2159         med_diag%FD_CAR3%dgsave = .FALSE. 
     2160      ENDIF 
     2161      IF (iom_use("FD_CAL3")) THEN  
     2162         med_diag%FD_CAL3%dgsave = .TRUE. 
     2163      ELSE  
     2164         med_diag%FD_CAL3%dgsave = .FALSE. 
     2165      ENDIF 
     2166      IF (iom_use("CO33")) THEN  
     2167         med_diag%CO33%dgsave = .TRUE. 
     2168      ELSE  
     2169         med_diag%CO33%dgsave = .FALSE. 
     2170      ENDIF 
     2171      IF (iom_use("CO3SATARAG3")) THEN  
     2172         med_diag%CO3SATARAG3%dgsave = .TRUE. 
     2173      ELSE  
     2174         med_diag%CO3SATARAG3%dgsave = .FALSE. 
     2175      ENDIF 
     2176      IF (iom_use("CO3SATCALC3")) THEN  
     2177         med_diag%CO3SATCALC3%dgsave = .TRUE. 
     2178      ELSE  
     2179         med_diag%CO3SATCALC3%dgsave = .FALSE. 
     2180      ENDIF 
     2181      IF (iom_use("DCALC3")) THEN  
     2182         med_diag%DCALC3%dgsave = .TRUE. 
     2183      ELSE  
     2184         med_diag%DCALC3%dgsave = .FALSE. 
     2185      ENDIF 
     2186      IF (iom_use("EXPC3")) THEN  
     2187         med_diag%EXPC3%dgsave = .TRUE. 
     2188      ELSE  
     2189         med_diag%EXPC3%dgsave = .FALSE. 
     2190      ENDIF 
     2191      IF (iom_use("EXPN3")) THEN  
     2192         med_diag%EXPN3%dgsave = .TRUE. 
     2193      ELSE  
     2194         med_diag%EXPN3%dgsave = .FALSE. 
     2195      ENDIF 
     2196      IF (iom_use("FEDISS3")) THEN  
     2197         med_diag%FEDISS3%dgsave = .TRUE. 
     2198      ELSE  
     2199         med_diag%FEDISS3%dgsave = .FALSE. 
     2200      ENDIF 
     2201      IF (iom_use("FESCAV3")) THEN  
     2202         med_diag%FESCAV3%dgsave = .TRUE. 
     2203      ELSE  
     2204         med_diag%FESCAV3%dgsave = .FALSE. 
     2205      ENDIF 
     2206      IF (iom_use("MIGRAZP3")) THEN  
     2207         med_diag%MIGRAZP3%dgsave = .TRUE. 
     2208      ELSE  
     2209         med_diag%MIGRAZP3%dgsave = .FALSE. 
     2210      ENDIF 
     2211      IF (iom_use("MIGRAZD3")) THEN  
     2212         med_diag%MIGRAZD3%dgsave = .TRUE. 
     2213      ELSE  
     2214         med_diag%MIGRAZD3%dgsave = .FALSE. 
     2215      ENDIF 
     2216      IF (iom_use("MEGRAZP3")) THEN  
     2217         med_diag%MEGRAZP3%dgsave = .TRUE. 
     2218      ELSE  
     2219         med_diag%MEGRAZP3%dgsave = .FALSE. 
     2220      ENDIF 
     2221      IF (iom_use("MEGRAZD3")) THEN  
     2222         med_diag%MEGRAZD3%dgsave = .TRUE. 
     2223      ELSE  
     2224         med_diag%MEGRAZD3%dgsave = .FALSE. 
     2225      ENDIF 
     2226      IF (iom_use("MEGRAZZ3")) THEN  
     2227         med_diag%MEGRAZZ3%dgsave = .TRUE. 
     2228      ELSE  
     2229         med_diag%MEGRAZZ3%dgsave = .FALSE. 
     2230      ENDIF 
     2231      IF (iom_use("O2SAT3")) THEN  
     2232         med_diag%O2SAT3%dgsave = .TRUE. 
     2233      ELSE  
     2234         med_diag%O2SAT3%dgsave = .FALSE. 
     2235      ENDIF 
     2236      IF (iom_use("PBSI3")) THEN  
     2237         med_diag%PBSI3%dgsave = .TRUE. 
     2238      ELSE  
     2239         med_diag%PBSI3%dgsave = .FALSE. 
     2240      ENDIF 
     2241      IF (iom_use("PCAL3")) THEN  
     2242         med_diag%PCAL3%dgsave = .TRUE. 
     2243      ELSE  
     2244         med_diag%PCAL3%dgsave = .FALSE. 
     2245      ENDIF 
     2246      IF (iom_use("REMOC3")) THEN  
     2247         med_diag%REMOC3%dgsave = .TRUE. 
     2248      ELSE  
     2249         med_diag%REMOC3%dgsave = .FALSE. 
     2250      ENDIF 
     2251      IF (iom_use("PNLIMJ3")) THEN  
     2252         med_diag%PNLIMJ3%dgsave = .TRUE. 
     2253      ELSE  
     2254         med_diag%PNLIMJ3%dgsave = .FALSE. 
     2255      ENDIF 
     2256      IF (iom_use("PNLIMN3")) THEN  
     2257         med_diag%PNLIMN3%dgsave = .TRUE. 
     2258      ELSE  
     2259         med_diag%PNLIMN3%dgsave = .FALSE. 
     2260      ENDIF 
     2261      IF (iom_use("PNLIMFE3")) THEN  
     2262         med_diag%PNLIMFE3%dgsave = .TRUE. 
     2263      ELSE  
     2264         med_diag%PNLIMFE3%dgsave = .FALSE. 
     2265      ENDIF 
     2266      IF (iom_use("PDLIMJ3")) THEN  
     2267         med_diag%PDLIMJ3%dgsave = .TRUE. 
     2268      ELSE  
     2269         med_diag%PDLIMJ3%dgsave = .FALSE. 
     2270      ENDIF 
     2271      IF (iom_use("PDLIMN3")) THEN  
     2272         med_diag%PDLIMN3%dgsave = .TRUE. 
     2273      ELSE  
     2274         med_diag%PDLIMN3%dgsave = .FALSE. 
     2275      ENDIF 
     2276      IF (iom_use("PDLIMFE3")) THEN  
     2277         med_diag%PDLIMFE3%dgsave = .TRUE. 
     2278      ELSE  
     2279         med_diag%PDLIMFE3%dgsave = .FALSE. 
     2280      ENDIF 
     2281      IF (iom_use("PDLIMSI3")) THEN  
     2282         med_diag%PDLIMSI3%dgsave = .TRUE. 
     2283      ELSE  
     2284         med_diag%PDLIMSI3%dgsave = .FALSE. 
     2285      ENDIF 
     2286       
    24272287   END SUBROUTINE   trc_nam_iom_medusa 
    24282288    
Note: See TracChangeset for help on using the changeset viewer.