Changeset 72 for branches/ORCHIDEE_EXT/ORCHIDEE

Timestamp:

2011-01-18T16:27:06+01:00 (13 years ago)

Author:

didier.solyga

Message:

Add a default array pft_to_mtc for the case of a standard use of ORCHIDEE

Location:

branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters

Files:

• constantes.f90 (modified) (1 diff)
• pft_parameters.f90 (modified) (1 diff)

branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/constantes.f90

@@ -45 +45 @@
   ! One year in seconds
   REAL(r_std),SAVE :: one_year
+
+  ! undef integer for integer arrays
+  INTEGER(i_std), PARAMETER :: undef_integer = 999999999
 
   ! Specific value if no restart value

branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/pft_parameters.f90

@@ -405 +405 @@
 
       ! Initialisation of the correspondance table
-      pft_to_mtc (:) = zero_int
+      pft_to_mtc (:) = undef_integer
       
       ! Reading of the conrrespondance table in the .def file
       CALL getin('PFT_TO_MTC',pft_to_mtc)
+
+      ! Add the standard configuration
+      IF(nvm .EQ. 13 ) THEN
+         IF(pft_to_mtc(1) .EQ. undef_integer) THEN
+            WRITE(numout,*) 'Note to the user : we will use ORCHIDEE to its standard configuration'
+            pft_to_mtc(:) = (/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 /)
+         ENDIF
+      ELSE   
+         IF(pft_to_mtc(1) .EQ. undef_integer) THEN
+            WRITE(numout,*)' The array PFT_TO_MTC is empty : we stop'
+         ENDIF
+      ENDIF
 
      ! What happened if pft_to_mtc(j) > nvmc (if the mtc doesn't exist)?