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

Ignore:
Timestamp:
2018-12-03T12:45:01+01:00 (5 years ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: merge with trunk@10365, see #2133

File:
1 edited

Legend:

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

    r10069 r10368  
    109109         DO jk = 1, nksrp       
    110110            etot_ndcy(:,:,jk) =        ze1(:,:,jk) +        ze2(:,:,jk) +       ze3(:,:,jk) 
    111             enano    (:,:,jk) =  2.1 * ze1(:,:,jk) + 0.42 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 
    112             ediat    (:,:,jk) =  1.6 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.7 * ze3(:,:,jk) 
     111            enano    (:,:,jk) =  1.85 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.46 * ze3(:,:,jk) 
     112            ediat    (:,:,jk) =  1.62 * ze1(:,:,jk) + 0.74 * ze2(:,:,jk) + 0.63 * ze3(:,:,jk) 
    113113         END DO 
    114114         IF( ln_p5z ) THEN 
    115115            DO jk = 1, nksrp       
    116               epico  (:,:,jk) =  2.1 * ze1(:,:,jk) + 0.42 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 
     116              epico  (:,:,jk) =  1.94 * ze1(:,:,jk) + 0.66 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 
    117117            END DO 
    118118         ENDIF 
     
    134134         DO jk = 1, nksrp       
    135135            etot (:,:,jk) =        ze1(:,:,jk) +        ze2(:,:,jk) +       ze3(:,:,jk) 
    136             enano(:,:,jk) =  2.1 * ze1(:,:,jk) + 0.42 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 
    137             ediat(:,:,jk) =  1.6 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.7 * ze3(:,:,jk) 
     136            enano(:,:,jk) =  1.85 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.46 * ze3(:,:,jk) 
     137            ediat(:,:,jk) =  1.62 * ze1(:,:,jk) + 0.74 * ze2(:,:,jk) + 0.63 * ze3(:,:,jk) 
    138138         END DO 
    139139         IF( ln_p5z ) THEN 
    140140            DO jk = 1, nksrp       
    141               epico(:,:,jk) =  2.1 * ze1(:,:,jk) + 0.42 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 
     141              epico(:,:,jk) =  1.94 * ze1(:,:,jk) + 0.66 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 
    142142            END DO 
    143143         ENDIF 
     
    182182      zetmp1 (:,:)   = 0.e0 
    183183      zetmp2 (:,:)   = 0.e0 
    184       zetmp3 (:,:)   = 0.e0 
    185       zetmp4 (:,:)   = 0.e0 
    186184 
    187185      DO jk = 1, nksrp 
     
    191189                  zetmp1 (ji,jj) = zetmp1 (ji,jj) + etot     (ji,jj,jk) * e3t_n(ji,jj,jk) ! remineralisation 
    192190                  zetmp2 (ji,jj) = zetmp2 (ji,jj) + etot_ndcy(ji,jj,jk) * e3t_n(ji,jj,jk) ! production 
    193                   zetmp3 (ji,jj) = zetmp3 (ji,jj) + enano    (ji,jj,jk) * e3t_n(ji,jj,jk) ! production 
    194                   zetmp4 (ji,jj) = zetmp4 (ji,jj) + ediat    (ji,jj,jk) * e3t_n(ji,jj,jk) ! production 
    195191                  zdepmoy(ji,jj) = zdepmoy(ji,jj) +                       e3t_n(ji,jj,jk) 
    196192               ENDIF 
     
    209205                  emoy (ji,jj,jk) = zetmp1(ji,jj) * z1_dep 
    210206                  zpar (ji,jj,jk) = zetmp2(ji,jj) * z1_dep 
    211                   enano(ji,jj,jk) = zetmp3(ji,jj) * z1_dep 
    212                   ediat(ji,jj,jk) = zetmp4(ji,jj) * z1_dep 
     207               ENDIF 
     208            END DO 
     209         END DO 
     210      END DO 
     211      ! 
     212      zdepmoy(:,:)   = 0.e0 
     213      zetmp3 (:,:)   = 0.e0 
     214      zetmp4 (:,:)   = 0.e0 
     215      ! 
     216      DO jk = 1, nksrp 
     217         DO jj = 1, jpj 
     218            DO ji = 1, jpi 
     219               IF( gdepw_n(ji,jj,jk+1) <= MIN(hmld(ji,jj), heup_01(ji,jj)) ) THEN 
     220                  zetmp3 (ji,jj) = zetmp3 (ji,jj) + enano    (ji,jj,jk) * e3t_n(ji,jj,jk) ! production 
     221                  zetmp4 (ji,jj) = zetmp4 (ji,jj) + ediat    (ji,jj,jk) * e3t_n(ji,jj,jk) ! production 
     222                  zdepmoy(ji,jj) = zdepmoy(ji,jj) +                       e3t_n(ji,jj,jk) 
     223               ENDIF 
     224            END DO 
     225         END DO 
     226      END DO 
     227      enanom(:,:,:) = enano(:,:,:) 
     228      ediatm(:,:,:) = ediat(:,:,:) 
     229      ! 
     230      DO jk = 1, nksrp 
     231         DO jj = 1, jpj 
     232            DO ji = 1, jpi 
     233               IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 
     234                  z1_dep = 1. / ( zdepmoy(ji,jj) + rtrn ) 
     235                  enanom(ji,jj,jk) = zetmp3(ji,jj) * z1_dep 
     236                  ediatm(ji,jj,jk) = zetmp4(ji,jj) * z1_dep 
    213237               ENDIF 
    214238            END DO 
     
    221245            DO jj = 1, jpj 
    222246               DO ji = 1, jpi 
    223                   IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN  
     247                  IF( gdepw_n(ji,jj,jk+1) <= MIN(hmld(ji,jj), heup_01(ji,jj)) ) THEN 
     248                     zetmp5(ji,jj)  = zetmp5 (ji,jj) + epico(ji,jj,jk) * e3t_n(ji,jj,jk) ! production 
     249                  ENDIF 
     250               END DO 
     251            END DO 
     252         END DO 
     253         ! 
     254         epicom(:,:,:) = epico(:,:,:) 
     255         ! 
     256         DO jk = 1, nksrp 
     257            DO jj = 1, jpj 
     258               DO ji = 1, jpi 
     259                  IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 
    224260                     z1_dep = 1. / ( zdepmoy(ji,jj) + rtrn ) 
    225                      zetmp5(ji,jj)  = zetmp5 (ji,jj) + epico(ji,jj,jk) * e3t_n(ji,jj,jk) ! production 
    226                      epico(ji,jj,jk) = zetmp5(ji,jj) * z1_dep 
     261                     epicom(ji,jj,jk) = zetmp5(ji,jj) * z1_dep 
    227262                  ENDIF 
    228263               END DO 
Note: See TracChangeset for help on using the changeset viewer.