!!! $Header$
!!!
!!!                       trclsm.hamocc3.h
!!!                       ****************
!!!
!!!  PURPOSE :
!!!  ---------
!!!     READs and PRINT options for HAMOCC3 or P3ZD NAMELIST
!!!
!!   METHOD :                   : no
!!   -------
!!
!!   MODIFICATIONS:
!!   --------------
!!      original  : 99-10 (M.A. Foujols, M. Levy) passive tracer
!!      addition  : 00-01 (L. Bopp) hamocc3,p3zd
!!     
!!----------------------------------------------------------------------
!! local declarations
!! ==================

#if defined key_passivetrc && defined key_trc_hamocc3
      INTEGER :: ji
      CHARACTER (len=32) :: clname

!!!---------------------------------------------------------------------
!!!  OPA8, LODYC (15/11/96)
!!!---------------------------------------------------------------------

! 0. initializations
! ------------------

       NAMELIST/natgas/ gasfac, igaswind, icice
       NAMELIST/natext/ ilecvit,atcco2
#    if defined key_trc_p3zd
       NAMELIST/natbio/silic0,pendec,pendes,fluexp,caco3r,dispo0,conc0,   &
       &               oxymin,grosip,jpkb,sedlam,sedlostpoc,sedlostsil,   &
       &               sedlostcal,pislope,betslope,ekc,ekw1,ekw2,excret,  &
       &               epsbio,wsbio,wchl,resrat,mprat,mzrat,grazrat,      &
       &               xprefc,xprefp,unass,xkgraz,xkmort,xksi1,xksi2,sicmax, &
       &               xremip,xremik,smax,xkdoc1,xkdoc2
#else
       NAMELIST/natbio/silic0,pendec,pendes,fluexp,caco3r,dispo0,conc0,   &
       &               oxymin,grosip,jpkb,sedlam,sedlostpoc,sedlostsil,   &
       &               sedlostcal
#endif


#    if defined key_trc_p3zd
      IF(lwp) THEN
          WRITE(numout,*) ' '
          WRITE(numout,*) ' ROUTINE trclec'
          WRITE(numout,*) ' **************'
          WRITE(numout,*) ' '
          WRITE(numout,*) ' NAMELIST for P3ZD model'
          WRITE(numout,*) ' ***********************'
          WRITE(numout,*) ' '
      ENDIF
#    else
      IF(lwp) THEN
          WRITE(numout,*) ' '
          WRITE(numout,*) ' ROUTINE trclec'
          WRITE(numout,*) ' **************'
          WRITE(numout,*) ' '
          WRITE(numout,*) ' NAMELIST for HAMOCC3 model'
          WRITE(numout,*) ' ***********************'
          WRITE(numout,*) ' '
      ENDIF
#    endif

      numnat=80
      clname ='namelist.trc.sms'
      OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old')

! 1 Namelist natgas :
! -------------------

      READ(numnat,natgas)

      IF(lwp) THEN
          WRITE(numout,*) ' '
          WRITE(numout,*) 'natgas'
          write(numout,*) 'gasfac = ',gasfac
          WRITE(numout,*) ' '
          write(numout,*) 'igaswind = ',igaswind
          WRITE(numout,*) ' '
          write(numout,*) 'icice = ',icice
          WRITE(numout,*) ' '

      ENDIF

! 2 Namelist natext :
! -------------------

      READ(numnat,natext)

      IF(lwp) THEN
          WRITE(numout,*) ' '
          WRITE(numout,*) 'natext'
          write(numout,*)'ilecvit = ',ilecvit
          WRITE(numout,*) ' '
          WRITE(numout,*) 'atmospheric pCO2= ',atcco2
          WRITE(numout,*) ' '
      ENDIF

      READ(numnat,natbio)
      IF(lwp) THEN
          WRITE(numout,*) 'natbio'
          WRITE(numout,*) ' '
          WRITE(numout,*)    &
          &   ' Silicon half saturation constant           =', silic0
          WRITE(numout,*)    &
          &   ' export length for calcite                  =', pendec
          WRITE(numout,*)    &
          &   ' export length for biogenic silicon         =', pendes
          WRITE(numout,*)    &
          &   ' power constant for POC export              =', fluexp
          WRITE(numout,*)    &
          &   ' mean rainratio                             =', caco3r
          WRITE(numout,*)    &
          &   ' Calcite dissolution half saturation        =', dispo0
          WRITE(numout,*)    &
          &   ' Phosphate half saturation                  =', conc0
          WRITE(numout,*)    &
          &   ' Productive zone vertical indice            =', jpkb
          WRITE(numout,*)    &
          &   ' Sediment bioturbation constant             =', sedlam
          WRITE(numout,*)    &
          &   ' Sediment geol loss for POC            =', sedlostpoc
          WRITE(numout,*)    &
          &   ' Sediment geol loss for SIL             =', sedlostsil
          WRITE(numout,*)    &
          &   ' Sediment geol loss for CAL             =', sedlostcal
#    if defined key_trc_p3zd
          WRITE(numout,*)    &
          &   ' P-I slope                                  =', pislope
          WRITE(numout,*)    &
          &   ' Light inhibition factor                    =', betslope
          WRITE(numout,*)    &
          &   ' pigment light absorption coeff.            =', ekc
          WRITE(numout,*)    &
          &   ' green light absorption coeff. of water     =', ekw1
          WRITE(numout,*)    &
          &   ' red light absorption coeff. of water       =', ekw2
          WRITE(numout,*)    &
          &   ' excretion ratio of phytoplankton           =', excret
          WRITE(numout,*)    &
          &   ' fraction of dead zooplankton --> POC       =', epsbio
          WRITE(numout,*)    &
          &   ' POC sinking speed                          =', wsbio
          WRITE(numout,*)    &
          &   ' exsudation rate of zooplankton             =', resrat
          WRITE(numout,*)    &
          &   ' phytoplankton mortality rate               =', mprat
          WRITE(numout,*)    &
          &   ' zooplankton mortality rate                 =', mzrat
          WRITE(numout,*)    &
          &   ' zoo preference for phyto                   =', xprefc
          WRITE(numout,*)    &
          &   ' zoo preference for POC                     =', xprefp
          WRITE(numout,*)    &
          &   ' maximal zoo grazing rate                   =', grazrat
          WRITE(numout,*)    &
          &   ' non assimilated fraction of phyto by zoo   =', unass
          WRITE(numout,*)    &
          &   ' half saturation constant for grazing        =', xkgraz
          WRITE(numout,*)    &
          &   ' half saturation constant for mortality     =', xkmort
          WRITE(numout,*)    &
          &   ' half saturation constant for Si uptake     =', xksi1
          WRITE(numout,*)    &
          &   ' half saturation constant for Si/C          =', xksi2
          WRITE(numout,*)    &
          &   ' maximum Si/C                               =', sicmax
          WRITE(numout,*)    &
          &   ' remineralisation rate of POC               =', xremip
          WRITE(numout,*)    &
          &   ' remineralization rate of DOC               =', xremik
          WRITE(numout,*)    &
          &   ' maximum sinking rate of phyto              =', smax
          WRITE(numout,*)    &
          &   ' 1st half-sat. of DOC remineralization      =', xkdoc1
          WRITE(numout,*)    &
          &   ' 2nd half-sat. of DOC remineralization      =', xkdoc2
          WRITE(numout,*) ' '
#    endif
      ENDIF

#else

! no passive tracers

#endif