Changeset 14957 for NEMO/branches/2021

Timestamp:

2021-06-07T14:26:34+02:00 (3 years ago)

Author:

cetlod

Message:

NEWDEV_PISCO : minor changes to avoid compilation error when using AGRIF

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zfechem.F90 (modified) (3 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zfechem.F90

@@ -59 +59 @@
       REAL(wp) ::   zscave, zaggdfea, zaggdfeb, ztrc, zdust, zklight
       REAL(wp) ::   ztfe, zhplus, zxlam, zaggliga, zaggligb
-      REAL(wp) ::   zprecip, zprecipno3, zvar
+      REAL(wp) ::   zprecip, zprecipno3,  zconsfe, za1
       REAL(wp) ::   zrfact2
       CHARACTER (len=25) :: charout
@@ -98 +98 @@
           zTL1(ji,jj,jk)  = ztotlig(ji,jj,jk)
           zkeq            = fekeq(ji,jj,jk)
-          zklight         = 4.77E-7 * etot(ji,jj,jk) * 0.5 / 10**-6.3
+          zklight         = 4.77E-7 * etot(ji,jj,jk) * 0.5 / ( 10**(-6.3) )
+          zconsfe         = consfe3(ji,jj,jk) / ( 10**(-6.3) )
           zfesatur        = zTL1(ji,jj,jk) * 1E-9
           ztfe            = (1.0 + zklight) * tr(ji,jj,jk,jpfer,Kbb) 
           ! Fe' is the root of a 2nd order polynom
-          zvar =  1. + zfesatur * zkeq + zklight + consfe3(ji,jj,jk)/10**-6.3 - zkeq * tr(ji,jj,jk,jpfer,Kbb)
-          zFe3 (ji,jj,jk) = ( -1 * zvar + SQRT( zvar**2 + 4. * ztfe * zkeq) ) / ( 2. * zkeq + rtrn )
+          za1 =  1. + zfesatur * zkeq + zklight +  zconsfe - zkeq * tr(ji,jj,jk,jpfer,Kbb)
+          zFe3 (ji,jj,jk) = ( -1 * za1 + SQRT( za1**2 + 4. * ztfe * zkeq) ) / ( 2. * zkeq + rtrn )
           zFeL1(ji,jj,jk) = MAX( 0., tr(ji,jj,jk,jpfer,Kbb) - zFe3(ji,jj,jk) )
       END_3D
@@ -264 +265 @@
          WRITE(numout,*) '      Scavenged iron that is added to POFe      scaveff      =', scaveff
       ENDIF
-      ! 
+      !
    END SUBROUTINE p4z_fechem_init