Changeset 6453 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90

Timestamp:

2016-04-08T10:10:11+02:00 (8 years ago)

Author:

aumont

Message:

New developments of PISCES (QUOTA, ligands, lability, ...)

File:

• branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90 (modified) (5 diffs)

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90

@@ -37 +37 @@
       !!----------------------------------------------------------------------
 
-      IF( lk_p4z ) THEN  ;   CALL p4z_ice_ini   !  PISCES
-      ELSE               ;   CALL p2z_ice_ini   !  LOBSTER
+      IF( nn_p4z == 1 ) THEN  ;   CALL p2z_ice_ini   !  PISCES
+      ELSE                    ;   CALL p4z_ice_ini   !  LOBSTER
       ENDIF
 
@@ -128 +128 @@
       zpisc(jpno3,1) =  5.79e-6_wp / rno3 
       zpisc(jpnh4,1) =  3.22e-7_wp / rno3
+#if defined key_pisces_quota
+      zpisc(jppic,1) =  9.57e-8_wp
+      zpisc(jpnpi,1) =  9.57e-8_wp
+      zpisc(jpppi,1) =  9.57e-8_wp
+      zpisc(jppfe,1) =  1.76e-11_wp
+      zpisc(jppch,1) =  1.67e-7_wp
+      zpisc(jpnph,1) =  9.57e-8_wp
+      zpisc(jppph,1) =  9.57e-8_wp
+      zpisc(jpndi,1) =  4.24e-7_wp
+      zpisc(jppdi,1) =  4.24e-7_wp
+      zpisc(jppon,1) =  9.57e-8_wp
+      zpisc(jppop,1) =  9.57e-8_wp
+      zpisc(jpdon,1) =  2.04e-5_wp
+      zpisc(jpdop,1) =  2.04e-5_wp
+      zpisc(jpgon,1) =  5.23e-8_wp
+      zpisc(jpgop,1) =  5.23e-8_wp
+#endif
 
       !--- Arctic specificities (dissolved inorganic & DOM)
@@ -158 +175 @@
       zpisc(jpno3,2) =  3.51e-06_wp / rno3 
       zpisc(jpnh4,2) =  6.15e-08_wp / rno3 
+#if defined key_pisces_quota
+      zpisc(jppic,1) =  5.25e-7_wp
+      zpisc(jpnpi,1) =  5.25e-7_wp
+      zpisc(jpppi,1) =  5.25e-7_wp
+      zpisc(jppfe,1) =  1.75e-11_wp
+      zpisc(jppch,1) =  1.46e-07_wp
+      zpisc(jpnph,1) =  5.25e-7_wp
+      zpisc(jppph,1) =  5.25e-7_wp
+      zpisc(jpndi,1) =  7.75e-7_wp
+      zpisc(jppdi,1) =  7.75e-7_wp
+      zpisc(jppon,1) =  4.05e-7_wp
+      zpisc(jppop,1) =  4.05e-7_wp
+      zpisc(jpdon,1) =  6.00e-6_wp
+      zpisc(jpdop,1) =  6.00e-6_wp
+      zpisc(jpgon,1) =  2.84e-8_wp
+      zpisc(jpgop,1) =  2.84e-8_wp
+#endif
+
 
       !--- Antarctic specificities (dissolved inorganic & DOM)
@@ -188 +223 @@
       zpisc(jpno3,3) =  2.64e-5_wp / rno3  
       zpisc(jpnh4,3) =  3.39e-7_wp / rno3  
+#if defined key_pisces_quota
+      zpisc(jppic,1) =  8.10e-7_wp
+      zpisc(jpnpi,1) =  8.10e-7_wp
+      zpisc(jpppi,1) =  8.10e-7_wp 
+      zpisc(jppfe,1) =  1.48e-11_wp
+      zpisc(jppch,1) =  2.02e-7_wp
+      zpisc(jpnph,1) =  9.57e-8_wp
+      zpisc(jppph,1) =  9.57e-8_wp
+      zpisc(jpndi,1) =  5.77e-7_wp
+      zpisc(jppdi,1) =  5.77e-7_wp
+      zpisc(jppon,1) =  1.13e-6_wp
+      zpisc(jppop,1) =  1.13e-6_wp
+      zpisc(jpdon,1) =  7.02e-6_wp
+      zpisc(jpdop,1) =  7.02e-6_wp
+      zpisc(jpgon,1) =  2.89e-8_wp
+      zpisc(jpgop,1) =  2.89e-8_wp
+#endif
+
 
       !--- Baltic Sea particular case for ORCA configurations
@@ -218 +271 @@
       zpisc(jpno3,4) = 5.36e-5_wp / rno3
       zpisc(jpnh4,4) = 7.18e-7_wp / rno3
+#if defined key_pisces_quota
+      zpisc(jppic,1) =  6.64e-7_wp
+      zpisc(jpnpi,1) =  6.64e-7_wp
+      zpisc(jpppi,1) =  6.64e-7_wp
+      zpisc(jppfe,1) =  3.89e-11_wp
+      zpisc(jppch,1) =  1.17e-7_wp
+      zpisc(jpnph,1) =  6.64e-7_wp
+      zpisc(jppph,1) =  6.64e-7_wp
+      zpisc(jpndi,1) =  3.41e-7_wp
+      zpisc(jppdi,1) =  3.41e-7_wp
+      zpisc(jppon,1) =  4.84e-7_wp
+      zpisc(jppop,1) =  4.84e-7_wp
+      zpisc(jpdon,1) =  1.06e-5_wp
+      zpisc(jpdop,1) =  1.06e-5_wp
+      zpisc(jpgon,1) =  1.05e-8_wp
+      zpisc(jpgop,1) =  1.05e-8_wp
+#endif
+
  
       DO jn = jp_pcs0, jp_pcs1