Ignore:
Timestamp:
2018-12-19T20:46:30+01:00 (22 months 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/p4zsed.F90

    r10368 r10419  
    6060      REAL(wp) ::  ztrfer, ztrpo4s, ztrdp, zwdust, zmudia, ztemp 
    6161      REAL(wp) ::  xdiano3, xdianh4 
    62       REAL(wp) ::  zwssfep 
    6362      ! 
    6463      CHARACTER (len=25) :: charout 
     
    6867      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zsoufer, zlight 
    6968      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrpo4, ztrdop, zirondep, zpdep 
    70       REAL(wp), ALLOCATABLE, DIMENSION(:,:  ) :: zsidep, zwsfep, zironice 
     69      REAL(wp), ALLOCATABLE, DIMENSION(:,:  ) :: zsidep, zironice 
    7170      !!--------------------------------------------------------------------- 
    7271      ! 
     
    8786      ALLOCATE( ztrpo4(jpi,jpj,jpk) ) 
    8887      IF( ln_p5z )    ALLOCATE( ztrdop(jpi,jpj,jpk) ) 
    89       IF( ln_ligand ) ALLOCATE( zwsfep(jpi,jpj) ) 
    9088 
    9189      zdenit2d(:,:) = 0.e0 
     
    131129         ELSE 
    132130            zirondep(:,:,1) = dustsolub  * dust(:,:) * mfrac * rfact2 / e3t_n(:,:,1) / 55.85 + 3.e-10 * r1_ryyss  
    133          ENDIF 
    134          IF ( ln_ligand ) THEN 
    135             IF( ln_solub ) THEN 
    136                tra(:,:,1,jpfep) = tra(:,:,1,jpfep) + rdustfep * (1.0 - solub(:,:)) * dust(:,:) * mfrac * rfact2 / e3t_n(:,:,1) / 55.85 
    137             ELSE 
    138                tra(:,:,1,jpfep) = tra(:,:,1,jpfep) + rdustfep * (1.0 - dustsolub) * dust(:,:) * mfrac * rfact2 / e3t_n(:,:,1) / 55.85 
    139             ENDIF 
    140131         ENDIF 
    141132         zsidep(:,:)   = 8.8 * 0.075 * dust(:,:) * mfrac * rfact2 / e3t_n(:,:,1) / 28.1  
     
    216207            tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + hydrofe(:,:,:) * rfact2 
    217208         IF( ln_ligand ) THEN 
    218             tra(:,:,:,jpfep) = tra(:,:,:,jpfep) + ( hydrofe(:,:,:) * fep_rath ) * rfact2 
    219209            tra(:,:,:,jplgw) = tra(:,:,:,jplgw) + ( hydrofe(:,:,:) * lgw_rath ) * rfact2 
    220210         ENDIF 
     
    235225      END DO 
    236226      ! 
    237       IF( ln_ligand ) THEN 
    238          DO jj = 1, jpj 
    239             DO ji = 1, jpi 
    240                ikt  = mbkt(ji,jj) 
    241                zdep = e3t_n(ji,jj,ikt) / xstep 
    242                zwsfep(ji,jj)  = MIN( 0.99 * zdep, wsfep(ji,jj,ikt)  ) 
    243             END DO 
    244          ENDDO 
    245       ENDIF 
    246  
    247227      IF( .NOT.lk_sed ) THEN 
    248228! 
     
    251231         IF( ln_ironsed ) THEN 
    252232                            tra(:,:,:,jpfer) = tra(:,:,:,jpfer) + ironsed(:,:,:) * rfact2 
    253             IF( ln_ligand ) tra(:,:,:,jpfep) = tra(:,:,:,jpfep) + ( ironsed(:,:,:) * fep_rats ) * rfact2 
    254233            ! 
    255234            IF( lk_iomput .AND. knt == nrdttrc .AND. iom_use( "Ironsed" ) )   & 
     
    334313         END DO 
    335314      END DO 
    336       ! 
    337       IF( ln_ligand ) THEN 
    338          DO jj = 1, jpj 
    339             DO ji = 1, jpi 
    340                ikt     = mbkt(ji,jj) 
    341                zdep    = xstep / e3t_n(ji,jj,ikt)  
    342                zwssfep = zwsfep(ji,jj) * zdep 
    343                tra(ji,jj,ikt,jpfep) = tra(ji,jj,ikt,jpfep) - trb(ji,jj,ikt,jpfep) * zwssfep 
    344             END DO 
    345          END DO 
    346       ENDIF 
    347315      ! 
    348316      IF( ln_p5z ) THEN 
     
    524492      ! 
    525493      IF( ln_p5z )    DEALLOCATE( ztrpo4, ztrdop ) 
    526       IF( ln_ligand ) DEALLOCATE( zwsfep ) 
    527494      ! 
    528495      IF( ln_timing )  CALL timing_stop('p4z_sed') 
Note: See TracChangeset for help on using the changeset viewer.