Ignore:
Timestamp:
2019-12-16T14:36:58+01:00 (11 months ago)
Author:
cetlod
Message:

dev_r12072_MERGE_OPTION2_2019: Cleaning PISCES diagnostics output ; fully sette tested

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/TOP/PISCES/P4Z/p4zflx.F90

    r12210 r12257  
    8080      CHARACTER (len=25) ::   charout 
    8181      REAL(wp), DIMENSION(jpi,jpj) ::   zkgco2, zkgo2, zh2co3, zoflx,  zpco2atm   
    82       REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zw2d 
    8382      !!--------------------------------------------------------------------- 
    8483      ! 
     
    186185 
    187186      IF( lk_iomput .AND. knt == nrdttrc ) THEN 
    188          CALL iom_put( "AtmCo2" , satmco2(:,:) * tmask(:,:,1) )   ! Atmospheric CO2 concentration 
    189          !  
    190          ALLOCATE( zw2d(jpi,jpj) )   
    191          IF( iom_use( "Cflx"  ) )  THEN 
    192             zw2d(:,:) = oce_co2(:,:) * 1000.  ! conversion in molC/m2/s 
    193             CALL iom_put( "Cflx"     , zw2d )  
    194          ENDIF 
    195          IF( iom_use( "Oflx"  ) )  THEN 
    196             zw2d(:,:) =  zoflx(:,:) * 1000. 
    197             CALL iom_put( "Oflx" , zw2d ) 
    198          ENDIF 
    199          IF( iom_use( "Kg"    ) )  THEN 
    200             zw2d(:,:) =  zkgco2(:,:) * tmask(:,:,1) 
    201             CALL iom_put( "Kg"   , zw2d ) 
    202          ENDIF 
    203          IF( iom_use( "Dpco2" ) ) THEN 
    204            zw2d(:,:) = ( zpco2atm(:,:) - zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1) 
    205            CALL iom_put( "Dpco2" ,  zw2d ) 
    206          ENDIF 
    207          IF( iom_use( "pCO2sea" ) ) THEN 
    208            zw2d(:,:) = ( zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1) 
    209            CALL iom_put( "pCO2sea" ,  zw2d ) 
    210          ENDIF 
    211  
    212          IF( iom_use( "Dpo2" ) )  THEN 
    213            zw2d(:,:) = ( atcox * patm(:,:) - atcox * trb(:,:,1,jpoxy) / ( chemo2(:,:,1) + rtrn ) ) * tmask(:,:,1) 
    214            CALL iom_put( "Dpo2"  , zw2d ) 
    215          ENDIF 
    216          CALL iom_put( "tcflx"    , t_oce_co2_flx     )   ! molC/s 
    217          CALL iom_put( "tcflxcum" , t_oce_co2_flx_cum )   ! molC 
    218          ! 
    219          DEALLOCATE( zw2d ) 
     187         CALL iom_put( "AtmCo2"  , satmco2(:,:) * tmask(:,:,1) )   ! Atmospheric CO2 concentration 
     188         CALL iom_put( "Cflx"    , oce_co2(:,:) * 1000. )  
     189         CALL iom_put( "Oflx"    , zoflx(:,:) * 1000.  ) 
     190         CALL iom_put( "Kg"      , zkgco2(:,:) * tmask(:,:,1)  ) 
     191         CALL iom_put( "Dpco2"   , ( zpco2atm(:,:) - zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1) ) 
     192         CALL iom_put( "pCO2sea" , ( zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1) ) 
     193         CALL iom_put( "Dpo2"    , ( atcox * patm(:,:) - atcox * trb(:,:,1,jpoxy) / ( chemo2(:,:,1) + rtrn ) ) * tmask(:,:,1) ) 
     194         CALL iom_put( "tcflx"   , t_oce_co2_flx     )   ! molC/s 
     195         CALL iom_put( "tcflxcum", t_oce_co2_flx_cum )   ! molC 
    220196      ENDIF 
    221197      ! 
Note: See TracChangeset for help on using the changeset viewer.