Changeset 72


Ignore:
Timestamp:
01/18/11 16:27:06 (11 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:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/constantes.f90

    r64 r72  
    4545  ! One year in seconds 
    4646  REAL(r_std),SAVE :: one_year 
     47 
     48  ! undef integer for integer arrays 
     49  INTEGER(i_std), PARAMETER :: undef_integer = 999999999 
    4750 
    4851  ! Specific value if no restart value 
  • branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/pft_parameters.f90

    r64 r72  
    405405 
    406406      ! Initialisation of the correspondance table 
    407       pft_to_mtc (:) = zero_int 
     407      pft_to_mtc (:) = undef_integer 
    408408       
    409409      ! Reading of the conrrespondance table in the .def file 
    410410      CALL getin('PFT_TO_MTC',pft_to_mtc) 
     411 
     412      ! Add the standard configuration 
     413      IF(nvm .EQ. 13 ) THEN 
     414         IF(pft_to_mtc(1) .EQ. undef_integer) THEN 
     415            WRITE(numout,*) 'Note to the user : we will use ORCHIDEE to its standard configuration' 
     416            pft_to_mtc(:) = (/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 /) 
     417         ENDIF 
     418      ELSE    
     419         IF(pft_to_mtc(1) .EQ. undef_integer) THEN 
     420            WRITE(numout,*)' The array PFT_TO_MTC is empty : we stop' 
     421         ENDIF 
     422      ENDIF 
    411423 
    412424     ! What happened if pft_to_mtc(j) > nvmc (if the mtc doesn't exist)? 
Note: See TracChangeset for help on using the changeset viewer.