Changeset 10788 for NEMO/trunk/src/TOP/PISCES/trcice_pisces.F90

Timestamp:

2019-03-21T12:15:14+01:00 (5 years ago)

Author:

cetlod

Message:

trunk: Some bugfixes at ice/biogeochmeical interfaces, see ticket #2262

File:

• NEMO/trunk/src/TOP/PISCES/trcice_pisces.F90 (modified) (3 diffs)

NEMO/trunk/src/TOP/PISCES/trcice_pisces.F90

@@ -266 +266 @@
       zpisc(jpgon,4) =  1.05e-8_wp
       zpisc(jpgop,4) =  1.05e-8_wp
- 
+!
+!     ln_ironice and tracers in seaice are redundant. Thus, if tracers in ice
+!     is activated, ln_ironice should be set to false
+!     ------------------------------------------------------------------------
+      IF( nn_ice_tr >= 0 .AND. ln_ironice ) THEN
+         IF(lwp) THEN
+            WRITE(numout,*) '   ==>>>   ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr
+            WRITE(numout,*) '           Specify your sea ice iron concentration in nampisice instead '
+            WRITE(numout,*) '           ln_ironice is forced to .FALSE. '
+         ENDIF
+         ln_ironice = .FALSE.
+      ENDIF
+! 
       DO jn = jp_pcs0, jp_pcs1
          IF( cn_trc_o(jn) == 'GL ' ) trc_o(:,:,jn) = zpisc(jn,1)  ! Global case
@@ -310 +322 @@
             trc_i(:,:,jn) = trc_ice_prescr(jn)
          ENDIF
-       
          !-- Baltic
          IF( cn_cfg == "orca" .OR. cn_cfg == "ORCA" ) THEN     
@@ -318 +329 @@
                      trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn) 
                END WHERE
-            ELSE                                 ! prescribed tracer concentration in ice
-               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    &
-                   54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp )
-                     trc_i(:,:,jn) = trc_ice_prescr(jn)
-               END WHERE
-            ENDIF ! trc_ice_ratio
+            ENDIF
          ENDIF
       !