Changeset 12759 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/sms_pisces.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/sms_pisces.F90 (modified) (3 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/sms_pisces.F90

@@ -132 +132 @@
       !!----------------------------------------------------------------------
       USE lib_mpp , ONLY: ctl_stop
-      INTEGER ::   ierr(10)        ! Local variables
+      INTEGER ::   ierr(11)        ! Local variables
       !!----------------------------------------------------------------------
       ierr(:) = 0
@@ -165 +165 @@
          !* Temperature dependency of SMS terms
          ALLOCATE( tgfunc (jpi,jpj,jpk) , tgfunc2(jpi,jpj,jpk),     &
-            &      tgfunc3(jpi,jpj,jpk) , STAT=ierr(6) )
+            &      tgfunc3(jpi,jpj,jpk) ,                         STAT=ierr(6) )
          !
          !* Sinking speed
-         ALLOCATE( wsbio3 (jpi,jpj,jpk) , wsbio4 (jpi,jpj,jpk),     &
-            &                             STAT=ierr(7) )   
+         ALLOCATE( wsbio3 (jpi,jpj,jpk) , wsbio4 (jpi,jpj,jpk),   STAT=ierr(7) )
+
+         !*  Size of phytoplankton cells
+         ALLOCATE( sizen (jpi,jpj,jpk), sized (jpi,jpj,jpk),        &
+           &       sizena(jpi,jpj,jpk), sizeda(jpi,jpj,jpk),      STAT=ierr(8) )
          ! 
          !* Prognostic ligand
          IF( ln_ligand ) THEN
-           ALLOCATE( plig(jpi,jpj,jpk)  ,                         STAT=ierr(8) )
+           ALLOCATE( plig(jpi,jpj,jpk)  ,                         STAT=ierr(9) )
          ENDIF
       ENDIF
@@ -179 +182 @@
       IF( ln_p5z ) THEN
          ! PISCES-QUOTA specific part      
-         ALLOCATE( epico(jpi,jpj,jpk)   , epicom(jpi,jpj,jpk) ,   STAT=ierr(9) ) 
+         ALLOCATE( epico(jpi,jpj,jpk)   , epicom(jpi,jpj,jpk) ,   STAT=ierr(10) ) 
 
          !*  Size of phytoplankton cells
-         ALLOCATE( sizen(jpi,jpj,jpk), sizep(jpi,jpj,jpk),         &
-           &       sized(jpi,jpj,jpk), sizena(jpi,jpj,jpk),        &
-           &       sizepa(jpi,jpj,jpk), sizeda(jpi,jpj,jpk),       STAT=ierr(10) )
+         ALLOCATE( sizep(jpi,jpj,jpk), sizepa(jpi,jpj,jpk),       STAT=ierr(11) )
       ENDIF
       !