Changeset 6022


Ignore:
Timestamp:
2015-12-08T13:33:34+01:00 (5 years ago)
Author:
jpalmier
Message:

JPALM — 8-12-2015 — diag with iom-use bugfixe - now they finally all work

File:
1 edited

Legend:

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

    r5982 r6022  
    6969# if defined key_iomput 
    7070      USE iom 
    71       USE trcnam_medusa                ! JPALM 13-11-2015 -- if iom_use for diag 
     71      USE trcnam_medusa         ! JPALM 13-11-2015 -- if iom_use for diag 
    7272      !!USE trc_nam_iom_medusa  ! JPALM 13-11-2015 -- if iom_use for diag 
    7373# endif 
     
    354354      REAL(wp), POINTER, DIMENSION(:,:  ) :: f_TDIC2d, f_TALK2d, f_kw6602d, f_pp02d, f_o2flux2d, f_o2sat2d 
    355355      REAL(wp), POINTER, DIMENSION(:,:  ) :: dms_surf2d, dms_andr2d, dms_simo2d, dms_aran2d, dms_hall2d 
     356      REAL(wp), POINTER, DIMENSION(:,:  ) :: iben_n2d, iben_fe2d, iben_c2d, iben_si2d, iben_ca2d, oben_n2d 
     357      REAL(wp), POINTER, DIMENSION(:,:  ) :: oben_fe2d, oben_c2d, oben_si2d, oben_ca2d, sfr_ocal2d 
     358      REAL(wp), POINTER, DIMENSION(:,:  ) :: sfr_oarg2d, lyso_ca2d  
    356359# endif 
    357360      !! 2D var for diagnostics. 
     
    434437      ffastc(:,:)  = 0.0        !! organic carbon 
    435438      ffastca(:,:) = 0.0        !! biogenic calcium carbonate 
     439      !! 
     440      fsedn(:,:)   = 0.0        !! Seafloor flux of N  
     441      fsedsi(:,:)  = 0.0        !! Seafloor flux of Si 
     442      fsedfe(:,:)  = 0.0        !! Seafloor flux of Fe 
     443      fsedc(:,:)   = 0.0        !! Seafloor flux of C 
     444      fsedca(:,:)  = 0.0        !! Seafloor flux of CaCO3 
    436445      !! 
    437446      fregenfast(:,:)   = 0.0   !! integrated  N regeneration (fast detritus) 
     
    522531      !! ----------------------------- 
    523532      !! Juju :: add kt condition !! 
    524       IF( lk_iomput .AND. .NOT.  ln_diatrc ) THEN    
     533      IF( lk_iomput .AND. .NOT.  ln_diatrc ) THEN  
    525534      !! 
    526535      if ( kt == nittrc000 )   CALL trc_nam_iom_medusa !! initialise iom_use test 
     
    831840         CALL wrk_alloc( jpi, jpj,    f_o2sat2d   ) 
    832841         f_o2sat2d(:,:)      = 0.0      !! 
    833       ENDIF    
     842      ENDIF  
     843      !! 
     844      IF( med_diag%IBEN_N%dgsave ) THEN 
     845         CALL wrk_alloc( jpi, jpj,    iben_n2d  ) 
     846         iben_n2d(:,:)      = 0.0      !! 
     847      ENDIF 
     848      IF( med_diag%IBEN_FE%dgsave ) THEN 
     849         CALL wrk_alloc( jpi, jpj,   iben_fe2d   ) 
     850         iben_fe2d(:,:)      = 0.0     !! 
     851      ENDIF 
     852      IF( med_diag%IBEN_C%dgsave ) THEN 
     853         CALL wrk_alloc( jpi, jpj,   iben_c2d   ) 
     854         iben_c2d(:,:)      = 0.0      !! 
     855      ENDIF 
     856      IF( med_diag%IBEN_SI%dgsave ) THEN 
     857         CALL wrk_alloc( jpi, jpj,   iben_si2d   ) 
     858         iben_si2d(:,:)      = 0.0      !! 
     859      ENDIF 
     860      IF( med_diag%IBEN_CA%dgsave ) THEN 
     861         CALL wrk_alloc( jpi, jpj,   iben_ca2d   ) 
     862         iben_ca2d(:,:)      = 0.0      !! 
     863      ENDIF 
     864      IF( med_diag%OBEN_N%dgsave ) THEN 
     865         CALL wrk_alloc( jpi, jpj,    oben_n2d  ) 
     866         oben_n2d(:,:)      = 0.0      !! 
     867      ENDIF 
     868      IF( med_diag%OBEN_FE%dgsave ) THEN 
     869         CALL wrk_alloc( jpi, jpj,    oben_fe2d  ) 
     870         oben_fe2d(:,:)      = 0.0      !! 
     871      ENDIF 
     872      IF( med_diag%OBEN_C%dgsave ) THEN 
     873         CALL wrk_alloc( jpi, jpj,    oben_c2d  ) 
     874         oben_c2d(:,:)      = 0.0      !! 
     875      ENDIF 
     876      IF( med_diag%OBEN_SI%dgsave ) THEN 
     877         CALL wrk_alloc( jpi, jpj,    oben_si2d  ) 
     878         oben_si2d(:,:)      = 0.0      !! 
     879      ENDIF 
     880      IF( med_diag%OBEN_CA%dgsave ) THEN 
     881         CALL wrk_alloc( jpi, jpj,    oben_ca2d  ) 
     882         oben_ca2d(:,:)      = 0.0      !! 
     883      ENDIF 
     884      IF( med_diag%SFR_OCAL%dgsave ) THEN 
     885         CALL wrk_alloc( jpi, jpj,    sfr_ocal2d  ) 
     886         sfr_ocal2d(:,:)      = 0.0      !! 
     887      ENDIF 
     888      IF( med_diag%SFR_OARG%dgsave ) THEN 
     889         CALL wrk_alloc( jpi, jpj,    sfr_oarg2d  ) 
     890         sfr_oarg2d(:,:)      = 0.0      !! 
     891      ENDIF 
     892      IF( med_diag%LYSO_CA%dgsave ) THEN 
     893         CALL wrk_alloc( jpi, jpj,    lyso_ca2d  ) 
     894         lyso_ca2d(:,:)      = 0.0      !! 
     895      ENDIF 
     896      !!   
    834897      IF (jdms .eq. 1) THEN 
    835898         IF( med_diag%DMS_SURF%dgsave ) THEN 
     
    10471110                     f3_omarg(ji,jj,jk) = f_omarg(ji,jj) 
    10481111                     !! 
    1049                      !! store 2D outputs 
    1050                      IF ( lk_iomput .and.(jk.eq.1) ) THEN 
    1051                         IF( med_diag%ATM_PCO2%dgsave ) THEN 
    1052                             f_pco2a2d(ji,jj) = f_pco2a 
    1053                         ENDIF 
    1054                         IF( med_diag%OCN_PCO2%dgsave ) THEN 
    1055                             f_pco2w2d(ji,jj) = f_pco2w 
    1056                         ENDIF 
    1057                         IF( med_diag%CO2FLUX%dgsave ) THEN 
    1058                             f_co2flux2d(ji,jj) = f_co2flux 
    1059                         ENDIF 
    1060                         IF( med_diag%TCO2%dgsave ) THEN 
    1061                              f_TDIC2d(ji,jj) = f_TDIC  
    1062                         ENDIF 
    1063                         IF( med_diag%TALK%dgsave ) THEN 
    1064                              f_TALK2d(ji,jj) = f_TALK  
    1065                         ENDIF 
    1066                         IF( med_diag%KW660%dgsave ) THEN 
    1067                              f_kw6602d(ji,jj) = f_kw660  
    1068                         ENDIF 
    1069                         IF( med_diag%ATM_PP0%dgsave ) THEN 
    1070                             f_pp02d(ji,jj) = f_pp0  
    1071                         ENDIF 
    1072                         IF( med_diag%O2FLUX%dgsave ) THEN 
    1073                              f_o2flux2d(ji,jj) = f_o2flux  
    1074                         ENDIF 
    1075                         IF( med_diag%O2SAT%dgsave ) THEN 
    1076                              f_o2sat2d(ji,jj) = f_o2sat  
    1077                         ENDIF                      
    1078                      ENDIF                      
    1079                      !! 
    10801112                     !! CCD calculation: calcite 
    10811113                     if (i2_omcal(ji,jj) .eq. 0 .and. f_omcal(ji,jj) .lt. 1.0) then 
     
    13141346                  ENDIF 
    13151347# endif 
    1316               ENDIF 
    1317                
     1348               ENDIF 
     1349 
    13181350               CALL flush(numout) 
    13191351 
     
    15061538                         dms_hall2d(ji,jj) = dms_hall 
    15071539                       ENDIF 
     1540#   if defined key_debug_medusa 
     1541                       IF (lwp) write (numout,*) 'trc_bio_medusa: finnish calculating dms' 
     1542                     CALL flush(numout) 
     1543#   endif  
    15081544                     ENDIF 
    15091545                     !! End iom 
    15101546                  ENDIF 
    15111547                  !! End DMS Loop 
     1548                  !! 
     1549                  !! store 2D outputs 
     1550                  IF ( lk_iomput ) THEN 
     1551                      IF( med_diag%ATM_PCO2%dgsave ) THEN 
     1552                         f_pco2a2d(ji,jj) = f_pco2a 
     1553                      ENDIF 
     1554                      IF( med_diag%OCN_PCO2%dgsave ) THEN 
     1555                         f_pco2w2d(ji,jj) = f_pco2w 
     1556                      ENDIF 
     1557                      IF( med_diag%CO2FLUX%dgsave ) THEN 
     1558                         f_co2flux2d(ji,jj) = f_co2flux 
     1559                      ENDIF 
     1560                      IF( med_diag%TCO2%dgsave ) THEN 
     1561                         f_TDIC2d(ji,jj) = f_TDIC 
     1562                      ENDIF 
     1563                      IF( med_diag%TALK%dgsave ) THEN 
     1564                         f_TALK2d(ji,jj) = f_TALK 
     1565                      ENDIF 
     1566                      IF( med_diag%KW660%dgsave ) THEN 
     1567                         f_kw6602d(ji,jj) = f_kw660 
     1568                      ENDIF 
     1569                      IF( med_diag%ATM_PP0%dgsave ) THEN 
     1570                         f_pp02d(ji,jj) = f_pp0 
     1571                      ENDIF 
     1572                      IF( med_diag%O2FLUX%dgsave ) THEN 
     1573                         f_o2flux2d(ji,jj) = f_o2flux 
     1574                      ENDIF 
     1575                      IF( med_diag%O2SAT%dgsave ) THEN 
     1576                         f_o2sat2d(ji,jj) = f_o2sat 
     1577                      ENDIF 
     1578                  ENDIF 
     1579                  !!  
    15121580               endif 
    15131581               !! End jk = 1 loop within ROAM key  
     
    34433511               endif                
    34443512 
     3513#   if defined key_debug_medusa 
     3514               IF (lwp) write (numout,*) '------' 
     3515               IF (lwp) write (numout,*) 'trc_bio_medusa: end all calculations' 
     3516               IF (lwp) write (numout,*) 'trc_bio_medusa: now outputs' 
     3517                     CALL flush(numout) 
     3518#   endif 
     3519 
    34453520# if defined key_axy_nancheck 
    34463521               !!---------------------------------------------------------------------- 
     
    35253600         !! 
    35263601         !! ** 2D diagnostics 
     3602#   if defined key_debug_medusa 
     3603                     IF (lwp) write (numout,*) 'trc_bio_medusa: diag in ij-jj-jk loop' 
     3604                     CALL flush(numout) 
     3605#   endif 
    35273606                  IF ( med_diag%PRN%dgsave ) THEN 
    35283607                      fprn2d(ji,jj) = fprn2d(ji,jj) + (fprn  * zphn * fthk)  
     
    36753754                         ffastca2d(ji,jj) =   & 
    36763755                            ffastca(ji,jj)/MAX(ffastc(ji,jj), rsmall) 
    3677                      ENDIF                      
     3756                     ENDIF 
     3757                 ELSE IF (jk.eq.(mbathy(ji,jj)-1)) THEN 
     3758                     IF( med_diag%IBEN_N%dgsave ) THEN 
     3759                         iben_n2d(ji,jj) = f_sbenin_n(ji,jj)  + f_fbenin_n(ji,jj) 
     3760                     ENDIF 
     3761                     IF( med_diag%IBEN_FE%dgsave ) THEN 
     3762                         iben_fe2d(ji,jj) = f_sbenin_fe(ji,jj) + f_fbenin_fe(ji,jj) 
     3763                     ENDIF 
     3764                     IF( med_diag%IBEN_C%dgsave ) THEN 
     3765                         iben_c2d(ji,jj) = f_sbenin_c(ji,jj)  + f_fbenin_c(ji,jj) 
     3766                     ENDIF 
     3767                     IF( med_diag%IBEN_SI%dgsave ) THEN 
     3768                         iben_si2d(ji,jj) = f_fbenin_si(ji,jj) 
     3769                     ENDIF 
     3770                     IF( med_diag%IBEN_CA%dgsave ) THEN 
     3771                         iben_ca2d(ji,jj) = f_fbenin_ca(ji,jj) 
     3772                     ENDIF 
     3773                     IF( med_diag%OBEN_N%dgsave ) THEN 
     3774                         oben_n2d(ji,jj) = f_benout_n(ji,jj) 
     3775                     ENDIF 
     3776                     IF( med_diag%OBEN_FE%dgsave ) THEN 
     3777                         oben_fe2d(ji,jj) = f_benout_fe(ji,jj) 
     3778                     ENDIF 
     3779                     IF( med_diag%OBEN_C%dgsave ) THEN 
     3780                         oben_c2d(ji,jj) = f_benout_c(ji,jj) 
     3781                     ENDIF 
     3782                     IF( med_diag%OBEN_SI%dgsave ) THEN 
     3783                         oben_si2d(ji,jj) = f_benout_si(ji,jj) 
     3784                     ENDIF 
     3785                     IF( med_diag%OBEN_CA%dgsave ) THEN 
     3786                         oben_ca2d(ji,jj) = f_benout_ca(ji,jj) 
     3787                     ENDIF 
     3788                     IF( med_diag%SFR_OCAL%dgsave ) THEN 
     3789                         sfr_ocal2d(ji,jj) = f3_omcal(ji,jj,jk) 
     3790                     ENDIF 
     3791                     IF( med_diag%SFR_OARG%dgsave ) THEN 
     3792                         sfr_oarg2d(ji,jj) =  f3_omarg(ji,jj,jk) 
     3793                     ENDIF 
     3794                     IF( med_diag%LYSO_CA%dgsave ) THEN 
     3795                         lyso_ca2d(ji,jj) = f_benout_lyso_ca(ji,jj) 
     3796                     ENDIF 
    36783797                 ENDIF 
    3679                 !! 
     3798                !! !! end bathy-1 diags 
    36803799                !! 
    36813800                  IF( med_diag%RIV_N%dgsave ) THEN 
     
    37803899                      !CALL iom_put( "TPP3"  , tpp3d ) 
    37813900                  ENDIF 
     3901 
    37823902                  IF( med_diag%REMIN3N%dgsave ) THEN 
    37833903                      remin3dn(ji,jj,jk) = fregen + (freminn * fthk)  !! remineralisation 
     
    37863906                  !! IF( med_diag%PH3%dgsave ) THEN 
    37873907                  !!     CALL iom_put( "PH3"  , f3_pH ) 
    3788                   !!  ENDIF 
     3908                  !! ENDIF 
    37893909                  !! IF( med_diag%OM_CAL3%dgsave ) THEN 
    37903910                  !!     CALL iom_put( "OM_CAL3"  , f3_omcal ) 
     
    37943914         !! ** Without using iom_use 
    37953915               ELSE IF( ln_diatrc ) THEN 
     3916#   if defined key_debug_medusa 
     3917                  IF (lwp) write (numout,*) 'trc_bio_medusa: diag in ij-jj-jk ln_diatrc' 
     3918                  CALL flush(numout) 
     3919#   endif 
    37963920                  !!---------------------------------------------------------------------- 
    37973921                  !! Prepare 2D diagnostics 
     
    41554279                 !!----------------------------------------- 
    41564280                 IF (jk.eq.1) THEN 
     4281#   if defined key_debug_medusa 
     4282                     IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 1' 
     4283                     CALL flush(numout) 
     4284#   endif 
    41574285                     IF( med_diag%MED_QSR%dgsave ) THEN 
    41584286                         CALL iom_put( "MED_QSR"  , qsr ) ! 
     
    42624390# endif                      
    42634391                 ELSE IF (jk.eq.i0100) THEN  
     4392#   if defined key_debug_medusa 
     4393                     IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 100' 
     4394                     CALL flush(numout) 
     4395#   endif 
    42644396                     IF( med_diag%SDT__100%dgsave ) THEN 
    42654397                         CALL iom_put( "SDT__100"  , fslownflux ) 
     
    42914423                     ENDIF                   
    42924424                 ELSE IF (jk.eq.i0150) THEN 
     4425#   if defined key_debug_medusa 
     4426                     IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 150' 
     4427                     CALL flush(numout) 
     4428#   endif 
    42934429                     IF( med_diag%BASIN_01%dgsave ) THEN 
    42944430                         CALL iom_put( "BASIN_01"  , fprn2d ) 
     
    44284564# endif                      
    44294565                 ELSE IF (jk.eq.i0200) THEN 
     4566#   if defined key_debug_medusa 
     4567                     IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 200' 
     4568                     CALL flush(numout) 
     4569#   endif 
    44304570                     IF( med_diag%SDT__200%dgsave ) THEN 
    44314571                         CALL iom_put( "SDT__200"  , fslownflux ) 
     
    44554595# endif                      
    44564596                 ELSE IF (jk.eq.i0500) THEN 
     4597#   if defined key_debug_medusa 
     4598                     IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 500' 
     4599                     CALL flush(numout) 
     4600#   endif 
    44574601                     IF( med_diag%SDT__500%dgsave ) THEN 
    44584602                         CALL iom_put( "SDT__500"  , fregen2d ) 
     
    44854629# endif                       
    44864630                 ELSE IF (jk.eq.i1000) THEN 
     4631#   if defined key_debug_medusa 
     4632                     IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 1000' 
     4633                     CALL flush(numout) 
     4634#   endif 
    44874635                     IF( med_diag%SDT_1000%dgsave ) THEN 
    44884636                         CALL iom_put( "SDT_1000"  , fslownflux ) 
     
    45154663                     ENDIF  
    45164664# endif                       
    4517                  ELSE IF (jk.eq.(mbathy(ji,jj)-1)) THEN 
    4518                      IF( med_diag%SEAFLRN%dgsave ) THEN 
    4519                          CALL iom_put( "SEAFLRN"  , fsedn ) 
    4520                      ENDIF 
    4521                      IF( med_diag%SEAFLRSI%dgsave ) THEN 
    4522                          CALL iom_put( "SEAFLRSI"  , fsedsi ) 
    4523                      ENDIF 
    4524                      IF( med_diag%SEAFLRFE%dgsave ) THEN 
    4525                          CALL iom_put( "SEAFLRFE"  , fsedfe ) 
    4526                      ENDIF 
    4527                      IF( med_diag%SEAFLRC%dgsave ) THEN 
    4528                          CALL iom_put( "SEAFLRC"  , fsedc ) 
    4529                      ENDIF 
    4530                      IF( med_diag%SEAFLRCA%dgsave ) THEN 
    4531                          CALL iom_put( "SEAFLRCA"  , fsedca ) 
    4532                      ENDIF 
    4533 # if defined key_roam                      
    4534                      IF( med_diag%IBEN_N%dgsave ) THEN 
    4535                          zw2d(:,:) = f_sbenin_n(:,:)  + f_fbenin_n(:,:)  
    4536                          CALL iom_put( "IBEN_N"  , zw2d ) 
    4537                      ENDIF 
    4538                      IF( med_diag%IBEN_FE%dgsave ) THEN 
    4539                          zw2d(:,:) = f_sbenin_fe(:,:) + f_fbenin_fe(:,:) 
    4540                          CALL iom_put( "IBEN_FE"  , zw2d ) 
    4541                      ENDIF 
    4542                      IF( med_diag%IBEN_C%dgsave ) THEN 
    4543                          zw2d(:,:) = f_sbenin_c(:,:)  + f_fbenin_c(:,:) 
    4544                          CALL iom_put( "IBEN_C"  , zw2d ) 
    4545                      ENDIF 
    4546                      IF( med_diag%IBEN_SI%dgsave ) THEN 
    4547                          CALL iom_put( "IBEN_SI"  , f_fbenin_si ) 
    4548                      ENDIF 
    4549                      IF( med_diag%IBEN_CA%dgsave ) THEN 
    4550                          CALL iom_put( "IBEN_CA"  , f_fbenin_ca ) 
    4551                      ENDIF 
    4552                      IF( med_diag%OBEN_N%dgsave ) THEN 
    4553                          CALL iom_put( "OBEN_N"  , f_benout_n ) 
    4554                      ENDIF 
    4555                      IF( med_diag%OBEN_FE%dgsave ) THEN 
    4556                          CALL iom_put( "OBEN_FE"  , f_benout_fe ) 
    4557                      ENDIF 
    4558                      IF( med_diag%OBEN_C%dgsave ) THEN 
    4559                          CALL iom_put( "OBEN_C"  , f_benout_c ) 
    4560                      ENDIF 
    4561                      IF( med_diag%OBEN_SI%dgsave ) THEN 
    4562                          CALL iom_put( "OBEN_SI"  , f_benout_si ) 
    4563                      ENDIF 
    4564                      IF( med_diag%OBEN_CA%dgsave ) THEN 
    4565                          CALL iom_put( "OBEN_CA"  , f_benout_ca ) 
    4566                      ENDIF 
    4567                      IF( med_diag%SFR_OCAL%dgsave ) THEN 
    4568                          zw2d(:,:) = f3_omcal(:,:,jk) 
    4569                          CALL iom_put( "SFR_OCAL"  , zw2d ) 
    4570                      ENDIF 
    4571                      IF( med_diag%SFR_OARG%dgsave ) THEN 
    4572                          zw2d(:,:) =  f3_omarg(:,:,jk) 
    4573                          CALL iom_put( "SFR_OARG"  , zw2d ) 
    4574                      ENDIF   
    4575                      IF( med_diag%LYSO_CA%dgsave ) THEN 
    4576                          CALL iom_put( "LYSO_CA"  , f_benout_lyso_ca ) 
    4577                      ENDIF  
    4578 # endif                      
    45794665                 ENDIF 
    45804666                 !! to do on every k loop : 
     
    49305016         !! 
    49315017         !! ** 2D diagnostics 
     5018#   if defined key_debug_medusa 
     5019         IF (lwp) write (numout,*) 'trc_bio_medusa: export all diag.' 
     5020         CALL flush(numout) 
     5021#   endif 
    49325022         IF ( med_diag%INVTN%dgsave ) THEN 
    49335023             CALL iom_put( "INVTN"  , ftot_n ) 
     
    51965286             CALL wrk_dealloc( jpi, jpj,   freminca2d  ) 
    51975287        ENDIF 
     5288        IF( med_diag%SEAFLRN%dgsave ) THEN 
     5289            CALL iom_put( "SEAFLRN"  , fsedn ) 
     5290        ENDIF 
     5291        IF( med_diag%SEAFLRSI%dgsave ) THEN 
     5292            CALL iom_put( "SEAFLRSI"  , fsedsi ) 
     5293        ENDIF 
     5294        IF( med_diag%SEAFLRFE%dgsave ) THEN 
     5295            CALL iom_put( "SEAFLRFE"  , fsedfe ) 
     5296        ENDIF 
     5297        IF( med_diag%SEAFLRC%dgsave ) THEN 
     5298            CALL iom_put( "SEAFLRC"  , fsedc ) 
     5299        ENDIF 
     5300        IF( med_diag%SEAFLRCA%dgsave ) THEN 
     5301            CALL iom_put( "SEAFLRCA"  , fsedca ) 
     5302        ENDIF 
     5303 
    51985304        !! 
    51995305        !! 
     
    53145420            CALL iom_put( "GMEDC"  , gmedc2d ) 
    53155421             CALL wrk_dealloc( jpi, jpj,    gmedc2d  ) 
     5422        ENDIF 
     5423        IF( med_diag%IBEN_N%dgsave ) THEN 
     5424           CALL iom_put( "IBEN_N"  , iben_n2d ) 
     5425             CALL wrk_dealloc( jpi, jpj,    iben_n2d  ) 
     5426        ENDIF 
     5427        IF( med_diag%IBEN_FE%dgsave ) THEN 
     5428           CALL iom_put( "IBEN_FE"  , iben_fe2d ) 
     5429             CALL wrk_dealloc( jpi, jpj,   iben_fe2d   ) 
     5430        ENDIF 
     5431        IF( med_diag%IBEN_C%dgsave ) THEN 
     5432           CALL iom_put( "IBEN_C"  , iben_c2d ) 
     5433             CALL wrk_dealloc( jpi, jpj,   iben_c2d   ) 
     5434        ENDIF 
     5435        IF( med_diag%IBEN_SI%dgsave ) THEN 
     5436           CALL iom_put( "IBEN_SI"  , iben_si2d ) 
     5437             CALL wrk_dealloc( jpi, jpj,   iben_si2d   ) 
     5438        ENDIF 
     5439        IF( med_diag%IBEN_CA%dgsave ) THEN 
     5440           CALL iom_put( "IBEN_CA"  , iben_ca2d ) 
     5441             CALL wrk_dealloc( jpi, jpj,   iben_ca2d   ) 
     5442        ENDIF 
     5443        IF( med_diag%OBEN_N%dgsave ) THEN 
     5444           CALL iom_put( "OBEN_N"  , oben_n2d ) 
     5445             CALL wrk_dealloc( jpi, jpj,    oben_n2d  ) 
     5446        ENDIF 
     5447        IF( med_diag%OBEN_FE%dgsave ) THEN 
     5448           CALL iom_put( "OBEN_FE"  , oben_fe2d ) 
     5449             CALL wrk_dealloc( jpi, jpj,    oben_fe2d  ) 
     5450        ENDIF 
     5451        IF( med_diag%OBEN_C%dgsave ) THEN 
     5452           CALL iom_put( "OBEN_C"  , oben_c2d ) 
     5453             CALL wrk_dealloc( jpi, jpj,    oben_c2d  ) 
     5454        ENDIF 
     5455        IF( med_diag%OBEN_SI%dgsave ) THEN 
     5456           CALL iom_put( "OBEN_SI"  , oben_si2d ) 
     5457             CALL wrk_dealloc( jpi, jpj,    oben_si2d  ) 
     5458        ENDIF 
     5459        IF( med_diag%OBEN_CA%dgsave ) THEN 
     5460           CALL iom_put( "OBEN_CA"  , oben_ca2d ) 
     5461             CALL wrk_dealloc( jpi, jpj, oben_ca2d     ) 
     5462        ENDIF 
     5463        IF( med_diag%SFR_OCAL%dgsave ) THEN 
     5464           CALL iom_put( "SFR_OCAL"  , sfr_ocal2d ) 
     5465             CALL wrk_dealloc( jpi, jpj,    sfr_ocal2d  ) 
     5466        ENDIF 
     5467        IF( med_diag%SFR_OARG%dgsave ) THEN 
     5468           CALL iom_put( "SFR_OARG"  , sfr_oarg2d ) 
     5469             CALL wrk_dealloc( jpi, jpj,    sfr_oarg2d  ) 
     5470        ENDIF 
     5471        IF( med_diag%LYSO_CA%dgsave ) THEN 
     5472           CALL iom_put( "LYSO_CA"  , lyso_ca2d ) 
     5473             CALL wrk_dealloc( jpi, jpj,    lyso_ca2d  ) 
    53165474        ENDIF 
    53175475# endif                    
Note: See TracChangeset for help on using the changeset viewer.