Changeset 10223


Ignore:
Timestamp:
2018-10-25T12:05:18+02:00 (19 months ago)
Author:
aumont
Message:

changes to account for the new sediment module

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/TOP/PISCES/P4Z/p4zsed.F90

    r10127 r10223  
    1717   USE p4zsbc          !  External source of nutrients  
    1818   USE p4zint          !  interpolation and computation of various fields 
     19   USE sed             !  Sediment module 
    1920   USE iom             !  I/O manager 
    2021   USE prtctl_trc      !  print control for debugging 
     
    2930   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:  ) :: sdenit     !: Nitrate reduction in the sediments 
    3031   REAL(wp) :: r1_rday                  !: inverse of rday 
     32   LOGICAL, SAVE :: lk_sed 
    3133 
    3234   !!---------------------------------------------------------------------- 
     
    7072      ! 
    7173      IF( ln_timing )  CALL timing_start('p4z_sed') 
     74      ! 
     75      IF( kt == nittrc000 .AND. knt == 1 )   THEN 
     76          r1_rday  = 1. / rday 
     77          IF (ln_sediment .AND. ln_sed_2way) THEN 
     78             lk_sed = .TRUE. 
     79          ELSE 
     80             lk_sed = .FALSE. 
     81          ENDIF 
     82      ENDIF 
    7283      ! 
    7384      IF( kt == nittrc000 .AND. knt == 1 )   r1_rday  = 1. / rday 
     
    185196      ENDIF 
    186197 
    187       ! Add the external input of iron from sediment mobilization 
    188       ! ------------------------------------------------------ 
    189       IF( ln_ironsed ) THEN 
    190                          tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + ironsed(:,:,:) * rfact2 
    191          IF( ln_ligand ) tra(:,:,:,jpfep) = tra(:,:,:,jpfep) + ( ironsed(:,:,:) * fep_rats ) * rfact2 
    192          ! 
    193          IF( lk_iomput .AND. knt == nrdttrc .AND. iom_use( "Ironsed" ) )   & 
    194             &   CALL iom_put( "Ironsed", ironsed(:,:,:) * 1.e+3 * tmask(:,:,:) ) ! iron inputs from sediments 
    195       ENDIF 
    196  
    197198      ! Add the external input of iron from hydrothermal vents 
    198199      ! ------------------------------------------------------ 
     
    231232 
    232233      IF( .NOT.lk_sed ) THEN 
     234! 
     235         ! Add the external input of iron from sediment mobilization 
     236         ! ------------------------------------------------------ 
     237         IF( ln_ironsed ) THEN 
     238                            tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + ironsed(:,:,:) * rfact2 
     239            IF( ln_ligand ) tra(:,:,:,jpfep) = tra(:,:,:,jpfep) + ( ironsed(:,:,:) * fep_rats ) * rfact2 
     240            ! 
     241            IF( lk_iomput .AND. knt == nrdttrc .AND. iom_use( "Ironsed" ) )   & 
     242               &   CALL iom_put( "Ironsed", ironsed(:,:,:) * 1.e+3 * tmask(:,:,:) ) ! iron inputs from sediments 
     243         ENDIF 
     244 
    233245         ! Computation of the sediment denitrification proportion: The metamodel from midlleburg (2006) is being used 
    234246         ! Computation of the fraction of organic matter that is permanently buried from Dunne's model 
Note: See TracChangeset for help on using the changeset viewer.