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 15127 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zlim.F90 – NEMO

Ignore:
Timestamp:
2021-07-16T20:00:12+02:00 (3 years ago)
Author:
cetlod
Message:

dev_PISCO : merge with trunk@15119

Location:
NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO

    • Property svn:externals

      •  

        old new  
        99 
        1010# SETTE 
        11 ^/utils/CI/sette@14244        sette 
         11^/utils/CI/sette@HEAD        sette 
         12 

  • NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zlim.F90

    r14786 r15127  
    105105      sizena(:,:,:) = 1.0  ;  sizeda(:,:,:) = 1.0 
    106106      ! 
    107       DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
     107      DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1) 
    108108          
    109109         ! Computation of a variable Ks for iron on diatoms taking into account 
    … …  
    232232      ! This fraction depends on nutrient limitation, light, temperature 
    233233      ! -------------------------------------------------------------------- 
    234       DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 
     234      DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1) 
    235235         zlim1  = xnanonh4(ji,jj,jk) + xnanono3(ji,jj,jk)  
    236236         zlim2  = tr(ji,jj,jk,jppo4,Kbb) / ( tr(ji,jj,jk,jppo4,Kbb) + concnnh4 ) 
    … …  
    249249         xfracal(ji,jj,jk) = MIN( 0.8 , xfracal(ji,jj,jk) ) 
    250250         xfracal(ji,jj,jk) = MAX( 0.02, xfracal(ji,jj,jk) ) 
     251      END_3D 
     252      ! 
     253      DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1) 
     254         ! denitrification factor computed from O2 levels 
     255         nitrfac(ji,jj,jk) = MAX(  0.e0, 0.4 * ( 6.e-6  - tr(ji,jj,jk,jpoxy,Kbb) )    & 
     256            &                                / ( oxymin + tr(ji,jj,jk,jpoxy,Kbb) )  ) 
     257         nitrfac(ji,jj,jk) = MIN( 1., nitrfac(ji,jj,jk) ) 
     258         ! 
     259         ! denitrification factor computed from NO3 levels 
     260         nitrfac2(ji,jj,jk) = MAX( 0.e0,       ( 1.E-6 - tr(ji,jj,jk,jpno3,Kbb) )  & 
     261            &                                / ( 1.E-6 + tr(ji,jj,jk,jpno3,Kbb) ) ) 
     262         nitrfac2(ji,jj,jk) = MIN( 1., nitrfac2(ji,jj,jk) ) 
    251263      END_3D 
    252264      ! 
Note: See TracChangeset for help on using the changeset viewer.