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 10419 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zsink.F90 – NEMO

Ignore:
Timestamp:
2018-12-19T20:46:30+01:00 (5 years ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: merge with trunk@10418, see #2133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zsink.F90

    r10377 r10419  
    3535   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sinkfer            !: Small BFe sinking fluxes 
    3636   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sinkfer2           !: Big iron sinking fluxes 
    37    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sinkfep      !: Fep sinking fluxes 
    3837 
    3938   INTEGER  :: ik100 
     
    126125      ENDIF 
    127126 
    128       IF( ln_ligand ) THEN 
    129          wsfep (:,:,:) = wfep 
    130          ! 
    131          sinkfep(:,:,:) = 0.e0 
    132          CALL trc_sink( kt, wsfep, sinkfep , jpfep, rfact2 ) 
    133       ENDIF 
    134  
    135127     ! Total carbon export per year 
    136128     IF( iom_use( "tcexp" ) .OR. ( ln_check_mass .AND. kt == nitend .AND. knt == nrdttrc )  )  & 
     
    214206      !!                     ***  ROUTINE p4z_sink_alloc  *** 
    215207      !!---------------------------------------------------------------------- 
    216       INTEGER :: ierr(3) 
     208      INTEGER :: ierr(2) 
    217209      !!---------------------------------------------------------------------- 
    218210      ! 
     
    224216         &      sinkfer(jpi,jpj,jpk)                                            , STAT=ierr(1) )                 
    225217         ! 
    226       IF( ln_ligand ) ALLOCATE( sinkfep(jpi,jpj,jpk)                            , STAT=ierr(2) )   
    227           
    228218      IF( ln_p5z    ) ALLOCATE( sinkingn(jpi,jpj,jpk), sinking2n(jpi,jpj,jpk)   ,     & 
    229          &                      sinkingp(jpi,jpj,jpk), sinking2p(jpi,jpj,jpk)   , STAT=ierr(3) ) 
     219         &                      sinkingp(jpi,jpj,jpk), sinking2p(jpi,jpj,jpk)   , STAT=ierr(2) ) 
    230220      ! 
    231221      p4z_sink_alloc = MAXVAL( ierr ) 
Note: See TracChangeset for help on using the changeset viewer.