Changeset 7192 for branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.F90

Timestamp:

2016-11-04T12:23:02+01:00 (8 years ago)

Author:

cetlod

Message:

new top interface : minor bugfixes

File:

• branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.F90 (modified) (2 diffs)

branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.F90

@@ -529 +529 @@
       IF (jcpoc > 1) THEN
          !
-         remindelta = log(4. * 1000. ) / float(jcpoc-1)
-         reminup = 1./ 400. * exp(remindelta)
+         remindelta = LOG(4. * 1000. ) / REAL(jcpoc-1, wp)
+         reminup = 1./ 400. * EXP(remindelta)
          !
          ! Discretization based on incomplete gamma functions
@@ -540 +540 @@
          reminp(1) = gamain(reminup, rshape+1.0, ifault) * xremip / alphan(1)
          DO jn = 2, jcpoc-1
-            reminup = 1./ 400. * exp(float(jn) * remindelta)
-            remindown = 1. / 400. * exp(float(jn-1) * remindelta)
+            reminup = 1./ 400. * EXP( REAL(jn, wp) * remindelta)
+            remindown = 1. / 400. * EXP( REAL(jn-1, wp) * remindelta)
             alphan(jn) = gamain(reminup, rshape, ifault) - gamain(remindown, rshape, ifault)
             reminp(jn) = gamain(reminup, rshape+1.0, ifault) - gamain(remindown, rshape+1.0, ifault)
             reminp(jn) = reminp(jn) * xremip / alphan(jn)
          END DO
-         remindown = 1. / 400. * exp(float(jcpoc-1) * remindelta)
+         remindown = 1. / 400. * EXP( REAL(jcpoc-1, wp) * remindelta)
          alphan(jcpoc) = 1.0 - gamain(remindown, rshape, ifault)
          reminp(jcpoc) = 1.0 - gamain(remindown, rshape+1.0, ifault)