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 8074 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 – NEMO

Ignore:
Timestamp:
2017-05-25T18:43:39+02:00 (7 years ago)
Author:
jpalmier
Message:

JPALM -- reverse MEDUSA cleaning and update MOCSY

File:
1 edited

Legend:

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

    r7894 r8074  
    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 
    1514   !!---------------------------------------------------------------------- 
    1615#if defined key_medusa 
     
    6867      ! 
    6968      NAMELIST/natbio/ xxi,xaln,xald,jphy,xvpn,xvpd,          & 
    70       &    xsin0,xnsi0,xuif,jliebig, xq10,                    & 
     69      &    xsin0,xnsi0,xuif,jliebig, jq10,                    & 
    7170      &    xthetam,xthetamd,xnln,xnld,xsld,xfln,xfld,         & 
    7271      &  xgmi,xgme,xkmi,xkme,xphi,xbetan,xbetac,xkc,          & 
     
    137136# endif 
    138137 
    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 
     138     IF( ( .NOT.lk_iomput .AND. ln_diatrc ) .OR. ( ln_diatrc .AND. lk_medusa ) ) THEN 
    141139         ! 
    142140         ! Namelist nammeddia 
     
    223221      ! ------------------------------------------- 
    224222       
    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 
     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 
    363343             
    364344      !REWIND(numnatm) 
    365345      !READ(numnatm,natbio) 
    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  
     346         ! Namelist natbio 
     347         ! ------------------- 
     348         REWIND( numnatp_ref )              ! Namelist natbio in reference namelist : MEDUSA diagnostics 
     349         READ  ( numnatp_ref, natbio, IOSTAT = ios, ERR = 903) 
     350903      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 ) 
     354904      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!! 
    377497      IF(lwp) THEN 
    378  
     498!! 
    379499!! AXY (08/11/13): compilation key notification 
    380500         WRITE(numout,*) '=== Compilation keys' 
     
    392512         WRITE(numout,*)     & 
    393513         &   ' 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' 
    394521#endif         
    395522#if defined key_avgqsr_medusa 
     
    474601         &   ' maximum Chl to C ratio for diatoms,                        xthetamd    = ', xthetamd 
    475602         WRITE(numout,*)     & 
    476          &   ' specific Q10 value (jphy==2),                                  xq10    = ', xq10 
     603         &   ' specific Q10 value (jphy==2),                                  jq10    = ', jq10 
    477604!! 
    478605!! Diatom silicon parameters 
     
    664791            WRITE(numout,*)     & 
    665792            &   ' 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 
    666802         else 
    667803            WRITE(numout,*)     & 
     
    8681004!! 
    8691005!! UKESM1 - new diagnostics  !! Jpalm; AXY (08/07/15) 
    870          WRITE(numout,*) '=== UKESM1 DMS parameters' 
     1006         WRITE(numout,*) '=== UKESM1-related parameters' 
    8711007         WRITE(numout,*)     & 
    8721008         &   ' include DMS diagnostic?,                                   jdms        = ', jdms 
     
    9481084 
    9491085      !READ(numnatm,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 ) 
     1086         ! Namelist natroam 
     1087         ! ------------------- 
     1088         REWIND( numnatp_ref )              ! Namelist natroam in reference namelist : MEDUSA diagnostics 
     1089         READ  ( numnatp_ref, natroam, IOSTAT = ios, ERR = 905) 
     1090905      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 ) 
     1094906      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in configuration namelist', lwp ) 
     1095         IF(lwm) WRITE ( numonp, natroam ) 
    9601096 
    9611097!! ROAM carbon, alkalinity and oxygen cycle parameters 
     
    9651101!!       xthetarem :  oxygen consumption by carbon remineralisation 
    9661102!!       xo2min    :  oxygen minimum concentration 
     1103 
    9671104      IF(lwp) THEN 
    9681105          WRITE(numout,*) 'natroam' 
    9691106          WRITE(numout,*) ' ' 
     1107!! 
     1108!! ROAM carbon, alkalinity and oxygen cycle parameters 
    9701109          WRITE(numout,*) '=== ROAM carbon, alkalinity and oxygen cycle parameters' 
    9711110          WRITE(numout,*)     & 
     
    9971136 
    9981137      !READ(numnatm,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 ) 
     1138         ! Namelist natopt 
     1139         ! ------------------- 
     1140         REWIND( numnatp_ref )              ! Namelist natopt in reference namelist : MEDUSA diagnostics 
     1141         READ  ( numnatp_ref, natopt, IOSTAT = ios, ERR = 907) 
     1142907      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 ) 
     1146908      IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in configuration namelist', lwp ) 
     1147         IF(lwm) WRITE ( numonp, natopt ) 
    10091148 
    10101149      IF(lwp) THEN 
     
    10191158         WRITE(numout,*) ' chla/chla+phea ratio            rpig  = ',rpig 
    10201159         WRITE(numout,*) ' ' 
     1160 
    10211161      ENDIF 
    10221162 
     
    10341174         endif 
    10351175         WRITE(numout,*) ' ' 
    1036       ENDIF 
     1176       ENDIF 
    10371177 
    10381178   END SUBROUTINE trc_nam_medusa 
     
    10501190      !!--------------------------------------------------------------------- 
    10511191      !! 
     1192      !! 
    10521193      !!----------------------------------------------------------------------             
    10531194      !! Variable conventions 
     
    10551196      !! 
    10561197      IF (iom_use("INVTN")) THEN  
    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. 
     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. 
    20002141      ENDIF 
    20012142      !! 
    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. 
     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. 
    20272167      ENDIF 
    20282168      !! 
     
    20322172      !! 
    20332173      !! 2D fields 
    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. 
     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. 
    21132253      ENDIF 
    21142254      !! 
    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        
     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 
    22872427   END SUBROUTINE   trc_nam_iom_medusa 
    22882428    
Note: See TracChangeset for help on using the changeset viewer.