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 7282 for branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90 – NEMO

Ignore:
Timestamp:
2016-11-21T12:13:57+01:00 (8 years ago)
Author:
flavoni
Message:

update CNRS 2016 to trunk revision 7161

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

    r7278 r7282  
    279279                  zprdia(ji,jj,jk) = zprdia(ji,jj,jk) * zmixdiat(ji,jj) 
    280280               ENDIF 
     281               zprbio(ji,jj,jk) = zprbio(ji,jj,jk) * ( 1. - fr_i(ji,jj) ) 
     282               zprdia(ji,jj,jk) = zprdia(ji,jj,jk) * ( 1. - fr_i(ji,jj) ) 
    281283            END DO 
    282284         END DO 
    … …  
    315317      END DO 
    316318 
    317       IF( ln_newprod ) THEN 
    318          DO jk = 1, jpkm1 
    319             DO jj = 1, jpj 
    320                DO ji = 1, jpi 
    321                   IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 
    322                      zprnch(ji,jj,jk) = zprnch(ji,jj,jk) * zmixnano(ji,jj) 
    323                      zprdch(ji,jj,jk) = zprdch(ji,jj,jk) * zmixdiat(ji,jj) 
    324                   ENDIF 
    325                   IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 
    326                      !  production terms for nanophyto. ( chlorophyll ) 
    327                      znanotot = enano(ji,jj,jk) * zstrn(ji,jj) 
    328                      zprod    = rday * zprorca(ji,jj,jk) * zprnch(ji,jj,jk) * xlimphy(ji,jj,jk) 
    329                      zprochln(ji,jj,jk) = chlcmin * 12. * zprorca (ji,jj,jk) 
    330                      zprochln(ji,jj,jk) = zprochln(ji,jj,jk) + (chlcnm-chlcmin) * 12. * zprod / & 
    331                                         & (  zpislopead(ji,jj,jk) * znanotot +rtrn) 
    332                      !  production terms for diatomees ( chlorophyll ) 
    333                      zdiattot = ediat(ji,jj,jk) * zstrn(ji,jj) 
    334                      zprod = rday * zprorcad(ji,jj,jk) * zprdch(ji,jj,jk) * xlimdia(ji,jj,jk) 
    335                      zprochld(ji,jj,jk) = chlcmin * 12. * zprorcad(ji,jj,jk) 
    336                      zprochld(ji,jj,jk) = zprochld(ji,jj,jk) + (chlcdm-chlcmin) * 12. * zprod / & 
    337                                         & ( zpislopead2(ji,jj,jk) * zdiattot +rtrn ) 
    338                   ENDIF 
    339                END DO 
    340             END DO 
    341          END DO 
    342       ELSE 
    343          DO jk = 1, jpkm1 
    344             DO jj = 1, jpj 
    345                DO ji = 1, jpi 
    346                   IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 
    347                      !  production terms for nanophyto. ( chlorophyll ) 
    348                      znanotot = enano(ji,jj,jk) 
    349                      zprod = rday * zprorca(ji,jj,jk) * zprnch(ji,jj,jk) * trb(ji,jj,jk,jpphy) * xlimphy(ji,jj,jk) 
    350                      zprochln(ji,jj,jk) = chlcmin * 12. * zprorca (ji,jj,jk) 
    351                      zprochln(ji,jj,jk) = zprochln(ji,jj,jk) + (chlcnm-chlcmin) * 144. * zprod            & 
    352                      &                    / ( zpislopead(ji,jj,jk) * trb(ji,jj,jk,jpnch) * znanotot +rtrn ) 
    353                      !  production terms for diatomees ( chlorophyll ) 
    354                      zdiattot = ediat(ji,jj,jk) 
    355                      zprod = rday * zprorcad(ji,jj,jk) * zprdch(ji,jj,jk) * trb(ji,jj,jk,jpdia) * xlimdia(ji,jj,jk) 
    356                      zprochld(ji,jj,jk) = chlcmin * 12. * zprorcad(ji,jj,jk) 
    357                      zprochld(ji,jj,jk) = zprochld(ji,jj,jk) + (chlcdm-chlcmin) * 144. * zprod             & 
    358                      &                    / ( zpislopead2(ji,jj,jk) * trb(ji,jj,jk,jpdch) * zdiattot +rtrn ) 
    359                   ENDIF 
    360                END DO 
    361             END DO 
    362          END DO 
    363       ENDIF 
     319      DO jk = 1, jpkm1 
     320         DO jj = 1, jpj 
     321            DO ji = 1, jpi 
     322               IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 
     323                  zprnch(ji,jj,jk) = zprnch(ji,jj,jk) * zmixnano(ji,jj) 
     324                  zprdch(ji,jj,jk) = zprdch(ji,jj,jk) * zmixdiat(ji,jj) 
     325               ENDIF 
     326               IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 
     327                  !  production terms for nanophyto. ( chlorophyll ) 
     328                  znanotot = enano(ji,jj,jk) * zstrn(ji,jj) 
     329                  zprod    = rday * zprorca(ji,jj,jk) * zprnch(ji,jj,jk) * xlimphy(ji,jj,jk) 
     330                  zprochln(ji,jj,jk) = chlcmin * 12. * zprorca (ji,jj,jk) 
     331                  zprochln(ji,jj,jk) = zprochln(ji,jj,jk) + (chlcnm-chlcmin) * 12. * zprod / & 
     332                                     & (  zpislopead(ji,jj,jk) * znanotot +rtrn) 
     333                  !  production terms for diatomees ( chlorophyll ) 
     334                  zdiattot = ediat(ji,jj,jk) * zstrn(ji,jj) 
     335                  zprod = rday * zprorcad(ji,jj,jk) * zprdch(ji,jj,jk) * xlimdia(ji,jj,jk) 
     336                  zprochld(ji,jj,jk) = chlcmin * 12. * zprorcad(ji,jj,jk) 
     337                  zprochld(ji,jj,jk) = zprochld(ji,jj,jk) + (chlcdm-chlcmin) * 12. * zprod / & 
     338                                     & ( zpislopead2(ji,jj,jk) * zdiattot +rtrn ) 
     339               ENDIF 
     340            END DO 
     341         END DO 
     342      END DO 
    364343 
    365344      !   Update the arrays TRA which contain the biological sources and sinks 
Note: See TracChangeset for help on using the changeset viewer.