Changeset 12759 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsms.F90

Timestamp:

2020-04-17T00:17:29+02:00 (4 years ago)

Author:

aumont

Message:

make parameterizations in PISCES-operationnal more similar to thos of PISCES-QUOTA (prey switching, optimal allocation, size, ...)

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsms.F90 (modified) (2 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsms.F90

@@ -320 +320 @@
          ENDIF
          !
+         ! PISCES size proxy
+         IF( iom_varid( numrtr, 'sized', ldstop = .FALSE. ) > 0 ) THEN
+            CALL iom_get( numrtr, jpdom_autoglo, 'sized' , sized(:,:,:)  )
+         ELSE
+            sized(:,:,:) = 1.
+         ENDIF
+         IF( iom_varid( numrtr, 'sizen', ldstop = .FALSE. ) > 0 ) THEN
+            CALL iom_get( numrtr, jpdom_autoglo, 'sizen' , sizen(:,:,:)  )
+         ELSE
+            sizen(:,:,:) = 1.
+         ENDIF
+
          ! PISCES-QUOTA specific part
          IF( ln_p5z ) THEN
             ! Read the size of the different phytoplankton groups
             ! If not in the restart file, they are set to 1
-            IF( iom_varid( numrtr, 'sized', ldstop = .FALSE. ) > 0 ) THEN
+            IF( iom_varid( numrtr, 'sizep', ldstop = .FALSE. ) > 0 ) THEN
                CALL iom_get( numrtr, jpdom_autoglo, 'sizep' , sizep(:,:,:)  )
-               CALL iom_get( numrtr, jpdom_autoglo, 'sizen' , sizen(:,:,:)  )
-               CALL iom_get( numrtr, jpdom_autoglo, 'sized' , sized(:,:,:)  )
             ELSE
                sizep(:,:,:) = 1.
-               sizen(:,:,:) = 1.
-               sized(:,:,:) = 1.
             ENDIF
         ENDIF
@@ -346 +354 @@
          CALL iom_rstput( kt, nitrst, numrtw, 'Silicamax', xksimax(:,:) ) ! Si max concentration
          CALL iom_rstput( kt, nitrst, numrtw, 'tcflxcum', t_oce_co2_flx_cum ) ! Cumulative CO2 flux
+         CALL iom_rstput( kt, nitrst, numrtw, 'sizen', sizen(:,:,:) )  ! Size of nanophytoplankton
+         CALL iom_rstput( kt, nitrst, numrtw, 'sized', sized(:,:,:) )  ! Size of diatoms
          IF( ln_p5z ) THEN
             CALL iom_rstput( kt, nitrst, numrtw, 'sizep', sizep(:,:,:) )  ! Size of picophytoplankton
-            CALL iom_rstput( kt, nitrst, numrtw, 'sizen', sizen(:,:,:) )  ! Size of nanophytoplankton
-            CALL iom_rstput( kt, nitrst, numrtw, 'sized', sized(:,:,:) )  ! Size of diatoms
          ENDIF
       ENDIF