Index: unk/NEMO/TOP_SRC/trclec.F
===================================================================
--- /trunk/NEMO/TOP_SRC/trclec.F	(revision 265)
+++ 	(revision )
@@ -1,283 +1,0 @@
-CCC $Header$ 
-CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-C ---------------------------------------------------------------------------
-CDIR$ LIST
-      SUBROUTINE trclec
-CCC---------------------------------------------------------------------
-CCC
-CCC                       ROUTINE trclec
-CCC                     ******************
-CCC
-CCC  PURPOSE :
-CCC  ---------
-CCC     READ and PRINT options for the passive tracer run (namelist)
-CCC
-CC   METHOD :                   : no
-CC   -------
-CC
-CC   INPUT :
-CC   -----
-CC      the namelist FILE ( UNIT numnat ) :
-CC            &nattrc           : general
-CC            &natnum           : numerical schemes
-CC            &natrtd           : dynamical trends 
-CC            &natadd           : additional arrays for diagnostics
-CC
-CC   OUTPUT :
-CC   ------
-CC      COMMON                  : (#ifdef "key_passivetrc")
-CC            /cittrc/          : tracer names
-CC            /cot3ad/          : options for smolarkiewicz scheme
-CC            /cit3ad/          : other options for smolarkiewicz scheme
-CC            /cottrd/          : passive tracer trends (#ifef
-CC                                 "key_trc_diatrd")
-CC
-CC   WORKSPACE :                : no
-CC   ---------
-CC
-CC   MODIFICATIONS:
-CC   --------------
-CC      original  : 96-11 (M.A. Foujols, M. Levy) passive tracer
-CC      modification : 98-04 (M.A Foujols, L. Bopp) ahtrb0 for isopycnal
-CC                                                  diffusion
-CC      modification : 99-10(M.A. Foujols, M. Levy) separation of sms
-CC      additions : 00-05(A. Estublier) TVD Limiter Scheme : Tests 
-CC                                      on ndttrc
-CC      additions : 00-06(A. Estublier) MUSCL Scheme : Tests 
-CC                                      on ndttrc
-CC      additions : 00-07(A. Estublier) PPM Scheme : Tests on ndttrc
-CC      modification : 00-11 (M.A Foujols, E Kestenare) trcrat, ahtrc0 and aeivtr0
-CC      modification : 01-01 (E Kestenare) suppress ndttrc=1 
-CC                                         for Arakawa and TVD schemes
-CC----------------------------------------------------------------------
-CC parameters and commons
-CC ======================
-CDIR$ NOLIST
-      USE oce_trc
-      USE trc
-      USE trctrp_lec
-      IMPLICIT NONE
-CDIR$ LIST
-CC----------------------------------------------------------------------
-CC local declarations
-CC ==================
-
-#if defined key_passivetrc
-      INTEGER ji
-      CHARACTER*32 clname
-CC----------------------------------------------------------------------
-CC statement functions
-CC ===================
-CDIR$ NOLIST
-#include "passivetrc_substitute.h90"
-CDIR$ LIST
-CC
-C
-C 0. initializations
-C ------------------
-C
-      namelist/nattrc/nwritetrc,lrsttr,nrsttr, ctrcnm,ctrcnl,ctrcun !general
-     $    ,lutini,ijulian
-
-      namelist/natnum/rsc,rtrn,ncortrc,ndttrc,crosster
-
-#    if defined key_trc_diatrd
-      namelist/natrtd/luttrd,nwritetrd                      ! dynamical trends
-#    endif
-
-#    if defined key_trc_diaadd
-      namelist/natadd/ctrc3d,ctrc3l,ctrc2d,ctrc2l, ctrc3u, ctrc2u,
-     $    nwriteadd             !additional diagnostics
-#    endif
-
-C
-C
-
-      IF(lwp) THEN
-          WRITE(numout,*) ' '
-          WRITE(numout,*) ' ROUTINE trclec'
-          WRITE(numout,*) ' **************'
-          WRITE(numout,*) ' '
-          WRITE(numout,*) ' namelist for passive tracer'
-          WRITE(numout,*) ' ***************************'
-          WRITE(numout,*) ' '
-      ENDIF
-C
-      numnat=80
-      REWIND (numnat)
-C
-      clname='namelist.passivetrc'
-      OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old')
-
-C
-C
-C 1., 2. & 3. initialization with namelist files
-C ----------------------------------------------
-C
-C 1.0 namelist nattrc :
-C
-      nwritetrc = 10
-      lrsttr=.FALSE.
-      nrsttr = 0
-
-      DO ji=1,jptra
-        WRITE (ctrcnm(ji),'("TR_",I1)') ji
-        WRITE (ctrcnl(ji),'("TRACER NUMBER ",I1)') ji
-        ctrcun(ji)='mmole/m3'
-        lutini(ji)=.FALSE. 
-      END DO
-C
-
-      REWIND(numnat)
-      READ(numnat,nattrc)
-C
-      IF(lwp) THEN
-          WRITE(numout,*) ' '
-          WRITE(numout,*) 'nattrc'
-          WRITE(numout,*) ' '
-          WRITE(numout,*)
-     $        ' frequency of outputs for passive tracers nwritetrc = '
-     $        ,nwritetrc  
-          WRITE(numout,*) ' restart LOGICAL for passive tr. lrsttr = ',
-     $        lrsttr
-          WRITE(numout,*) ' control of time step for p. tr. nrsttr = ',
-     $        nrsttr
-          DO ji=1,jptra
-            WRITE(numout,*) ' tracer nb: ',ji,' name = ',ctrcnm(ji)
-     $          ,ctrcnl(ji) 
-            WRITE(numout,*) ' in unit = ',ctrcun(ji)
-            WRITE(numout,*) ' initial value in FILE : ',lutini(ji) 
-            WRITE(numout,*) ' '
-          END DO
-          WRITE(numout,*) ' '
-      ENDIF
-
-# if defined key_trc_diatrd
-C
-C 1.2 namelist nattrd : passive tracers dynamical trends
-C
-      nwritetrd=10
-C
-C default : no dynamical trend recording
-C
-      DO ji=1,jptra
-        luttrd(ji) = .FALSE.
-      END DO 
-
-      REWIND(numnat)
-      READ(numnat,natrtd)
-C
-      IF(lwp) THEN
-          WRITE(numout,*) 'natrtd'
-          WRITE(numout,*) ' '
-          WRITE(numout,*)
-     $        ' frequency of outputs for dynamical trends nwritetrd = '
-     $        ,nwritetrd
-          DO ji=1,jptra
-            WRITE(numout,*)
-     $          ' keep dynamical trends for tracer number :',ji
-     $          ,luttrd(ji)   
-          END DO 
-      ENDIF
-# endif
-C
-C 1.3 namelist natadd : passive tracers diagnostics
-C
-# if defined key_trc_diaadd
-C
-      nwriteadd = 10
-C
-C default value for 3D output arrays : short and long name, units
-C
-      DO ji=1,jpdia3d
-        WRITE (ctrc3d(ji),'("3D_",I1)') ji
-        WRITE (ctrc3l(ji),'("3D DIAGNOSTIC NUMBER ",I2)') ji
-        ctrc3u(ji)=' '
-      END DO 
-C
-C
-C default value for 2D output arrays : short and long name, units
-C
-      DO ji=1,jpdia2d
-        WRITE (ctrc2d(ji),'("2D_",I1)') ji
-        WRITE (ctrc2l(ji),'("2D DIAGNOSTIC NUMBER ",I2)') ji
-        ctrc2u(ji)=' '
-      END DO 
-C
-      REWIND(numnat)
-      READ(numnat,natadd)
-C
-      IF(lwp) THEN
-          WRITE(numout,*) ' natadd'
-          WRITE(numout,*) ' '
-          WRITE(numout,*)
-     $        ' frequency of outputs for additional arrays nwriteadd = '
-     $        ,nwriteadd
-          DO ji=1,jpdia3d
-            WRITE(numout,*)
-     $          'name of 3d output field number :',ji,' : ',ctrc3d(ji)  
-            WRITE(numout,*) ctrc3l(ji)  
-            WRITE(numout,*) ' in unit = ',ctrc3u(ji)
-          END DO 
-          WRITE(numout,*) ' '
-          DO ji=1,jpdia2d
-            WRITE(numout,*)
-     $          'name of 2d output field number :',ji,' : ',ctrc2d(ji)  
-            WRITE(numout,*) ctrc2l(ji)  
-            WRITE(numout,*) ' in unit = ',ctrc2u(ji)
-          END DO 
-          WRITE(numout,*) ' '
-      ENDIF
-#endif
-
-C 1.1 namelist natnum :
-C
-      rsc=1.
-      rtrn=1.e-15
-      ncortrc=1
-      ndttrc=4
-      crosster=.FALSE.
-C
-      REWIND(numnat)
-      READ(numnat,natnum)
-
-Chris  computes the first time step of tracer model
-      nittrc000 = nit000 + ndttrc - 1
-Chris
-
-C
-      IF(lwp) THEN
-          WRITE(numout,*) ' '
-          WRITE(numout,*) 'natnum'
-          WRITE(numout,*) ' '
-          WRITE(numout,*) ' tuning coefficient              rsc     = ',
-     $        rsc
-          WRITE(numout,*) ' truncation value                rtrn    = ',
-     $        rtrn
-          WRITE(numout,*) ' number of corrective phase      ncortrc = ',
-     $        ncortrc
-          WRITE(numout,*) ' time step freq. for pass. trac. ndttrc  = ',
-     $        ndttrc
-          WRITE(numout,*) ' 1st time step for pass. trac. nittrc000 = ',
-     $        nittrc000
-          WRITE(numout,*) ' computes or not crossterms    crosster  = ',
-     $        crosster
-      ENDIF
-
-
-C
-C namelist of transport
-C
-      CALL trc_trp_lec
-C
-C namelist of SMS
-      CALL trclsm
-
-#else
-C
-C no passive tracers
-C
-#endif
-      RETURN
-      END
Index: unk/NEMO/TOP_SRC/trclsm.F
===================================================================
--- /trunk/NEMO/TOP_SRC/trclsm.F	(revision 265)
+++ 	(revision )
@@ -1,64 +1,0 @@
-CCC $Header$ 
-CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-C ---------------------------------------------------------------------------
-CDIR$ LIST
-      SUBROUTINE trclsm
-CCC---------------------------------------------------------------------
-CCC
-CCC                       ROUTINE trclsm
-CCC                     ******************
-CCC
-CCC  PURPOSE :
-CCC  ---------
-CCC     READS specific NAMELIST for sms terms 
-CCC
-CC   METHOD :                   : no
-CC   -------
-CC
-CC   INPUT :
-CC   -----
-CC
-CC   OUTPUT :
-CC   ------
-CC
-CC   WORKSPACE :                : no
-CC   ---------
-CC
-CC   MODIFICATIONS:
-CC   --------------
-CC      original  : 99-10 (M.A. Foujols, M. Levy) passive tracer
-CC                  00-04 (O. Aumont, M.A. Foujols) HAMOCC3 and P3ZD
-CC                  02-05 (O. Aumont) AGE TRACER
-CC----------------------------------------------------------------------
-CC parameters and commons
-CC ======================
-CDIR$ NOLIST
-      USE oce_trc
-      USE trc
-      USE sms
-      IMPLICIT NONE
-CDIR$ LIST
-
-#if defined key_passivetrc
-
-#    if defined key_trc_npzd
-#    include "trclsm.npzd.h"
-#    elif defined key_trc_lobster1
-#    include "trclsm.lobster1.h"
-#    elif defined key_trc_hamocc3
-#    include "trclsm.hamocc3.h"
-#    elif defined key_trc_age
-#    include "trclsm.age.h"
-#    elif defined key_trc_pisces
-#    include "trclsm.pisces.h"
-#    else
-C      user can define here his own sms namelist
-#    endif
-#else
-C
-C no passive tracers
-C
-#endif
-      RETURN
-      END
Index: unk/NEMO/TOP_SRC/trcrst.F
===================================================================
--- /trunk/NEMO/TOP_SRC/trcrst.F	(revision 265)
+++ 	(revision )
@@ -1,241 +1,0 @@
-CCC $Header$ 
-CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-C ---------------------------------------------------------------------------
-      SUBROUTINE trcrst
-CCC---------------------------------------------------------------------
-C
-CCC
-CCC                       ROUTINE trcrst
-CCC                     ******************
-CCC
-CCC  PURPOSE :
-CCC  ---------
-CCC     READ files for restart for passive tracer
-CCC
-CC   METHOD :
-CC   -------
-CC      READ the previous fields on the FILE nutrst
-CC      the first record indicates previous characterics
-CC      after control with the present run, we READ :
-CC      - prognostic variables on the second and more record
-CC
-CC   INPUT :			  no
-CC   -----
-CC      COMMON
-CC            /citrst/          : restart PARAMETER
-CC            /comctl/          : parameters for the control
-CC
-CC   OUTPUT :
-CC   ------
-CC      argument                : no
-CC      COMMON
-CC            /cottrc/          : passive tracers fields now
-CC                                and before 
-CC
-CC   WORKSPACE :		  no
-CC   ---------
-CC      ji,jj,jk,jl
-CC      ino0,it0,ipcg0,isor0,itke0,ino1,it1,isor1,ipcg1,itke1,idast1
-CC      iarak0,iarak1
-CC
-CC   MODIFICATIONS:
-CC   --------------
-CC      original  : 96-11
-CC      additions : 00-05 (A. Estublier) TVD Limiter Scheme
-CC                                       key_trc_tvd
-CC                  00-12 (O. Aumont, E. Kestenare) read restart file 
-CC                                                  for sediments                
-CC                  01-05 (O. Aumont, E. Kestenare) read restart file 
-CC                                   for calcite and silicate sediments
-CC----------------------------------------------------------------------
-CC parameters and commons
-CC ======================
-CDIR$ NOLIST
-#if defined key_mpp
-c no ioipsl
-#include "trcrst.mpp.h"
-#    else
-CDIR$ NOLIST
-      USE ioipsl
-      USE oce_trc
-      USE trc
-      USE sms
-      USE trctrp_lec
-      IMPLICIT NONE
-CDIR$ LIST
-#if defined key_passivetrc
-CC----------------------------------------------------------------------
-CC local declarations
-CC ==================
-      INTEGER ino0,ino1,it0,it1,jn,iarak0,iarak1
-      INTEGER ji, jj, jk
-      INTEGER itime, ibvar
-      REAL zdate0, zinfo(3)
-      CHARACTER*32 clname,clname1,clname2
-      LOGICAL clog
-      CHARACTER*8 clvnames(30)
-      REAL caralk,bicarb
-
-C
-C
-C 0. initialisations
-C ------------------
-C
-      ino0 =no
-      it0  =nit000
-C
-      IF( ln_trcadv_cen2 .OR. ln_trcadv_tvd ) THEN
-          iarak0 = 1
-      ELSE
-          iarak0=0
-      ENDIF
-
-C
-      IF(lwp) THEN
-          WRITE(numout,*) ' '
-          WRITE(numout,*) ' *** trcrst beginning of restart for'
-          WRITE(numout,*) ' passive tracer'
-          WRITE(numout,*) ' the present run :'
-          WRITE(numout,*) '   number job is  : ',no
-          WRITE(numout,*) '   with the time nit000 : ',nit000
-          IF(iarak0.eq.1) then
-              WRITE(numout,*) '   and before fields for Arakawa sheme '
-          ENDIF 
-          WRITE(numout,*) ' '
-      ENDIF
-C
-C 1. READ nutrst
-C --------------
-C
-C ... first information
-C
-      clname='restart.trc'
-      itime=0
-      clog=.false.
-      CALL restini(clname,jpi,jpj,glamt,gphit,jpk,gdept,clname
-     $    ,itime,zdate0,rdt,nutrst)
-C
-      CALL ioget_vname(nutrst, ibvar, clvnames)
-      CALL restget(nutrst,'info',1,1,3,0,clog,zinfo)
-      ino1  = nint(zinfo(1))
-      it1   = nint(zinfo(2))
-      iarak1 = nint(zinfo(3))
-C
-      IF(lwp) THEN
-          WRITE(numout,*) ' '
-          WRITE(numout,*) ' READ nutrst with '
-          WRITE(numout,*) '   number job is  : ',ino1
-          WRITE(numout,*) '   with the time it : ',it1
-          IF(iarak1.eq.1) then
-              WRITE(numout,*) '   and before fields for Arakawa sheme '
-          ENDIF 
-          WRITE(numout,*) '   number of variables   : ', ibvar
-          WRITE(numout,*) '   NetCDF variables      : '
-          WRITE(numout,*) ' ',clvnames (:ibvar)
-          WRITE(numout,*) ' '
-      ENDIF
-C
-C 1.2 control of date
-C
-      IF((it0-it1).ne.1.and.abs(nrsttr).eq.1) THEN
-          IF(lwp) THEN
-              WRITE(numout,*) ' ===>>>> : problem with nit000 for the',
-     $            ' passive tracer restart'
-              WRITE(numout,*) ' =======                              ',
-     $            ' ======================'
-              WRITE(numout,*) ' we stop. verify the FILE'
-              WRITE(numout,*) ' or rerun with the value  0 for the'
-              WRITE(numout,*) ' control of time PARAMETER   nrstdt'
-              WRITE(numout,*) ' '
-          ENDIF
-          STOP 'trcrst'
-      ENDIF
-C
-C 1.3 Control of the sheme
-C
-      IF(iarak0.ne.iarak1) THEN
-          IF(lwp) THEN
-              WRITE(numout,*) ' ===>>>> : problem with the',
-     $            ' passive tracer restart file'
-              WRITE(numout,*) ' =======                              ',
-     $            ' ==========================='
-              WRITE(numout,*) ' we stop. verify the FILE'
-              WRITE(numout,*) ' before field required IF 1=',iarak0
-              WRITE(numout,*) ' before field present in file IF 1=',
-     $            iarak1
-              WRITE(numout,*) ' '
-          ENDIF
-          STOP 'trcrst'
-      ENDIF
-
-C
-C ... READ prognostic variables and computes diagnostic variable
-C
-       DO jn=1,jptra
-         clname='TRN'//ctrcnm(jn)
-         CALL restget(nutrst,clname,jpi,jpj,jpk,0,clog,trn(:,:,:,jn))
-       END DO
-
-       DO jn=1,jptra
-         clname='TRB'//ctrcnm(jn)
-         CALL restget(nutrst,clname,jpi,jpj,jpk,0,clog,trb(:,:,:,jn))
-       END DO
-
-#     if defined key_trc_hamocc3
-         clname='SED'//ctrcnm(jppoc)
-         clname1='SED'//ctrcnm(jpcal)
-         clname2='SED'//ctrcnm(jpsil)
-         CALL restget(nutrst,clname1,jpi,jpj,1,0,clog,sedcal(:,:))
-         CALL restget(nutrst,clname2,jpi,jpj,1,0,clog,sedsil(:,:))
-         CALL restget(nutrst,clname,jpi,jpj,1,0,clog,sedpoc(:,:))
-#     elif defined key_trc_npzd || defined key_trc_lobster1
-         clname='SED'//ctrcnm(jpdet)
-         CALL restget(nutrst,clname,jpi,jpj,1,0,clog,sedpoc(:,:))
-#     elif defined key_trc_pisces
-C
-         clname='Silicalim'
-         CALL restget(nutrst,clname,jpi,jpj,1,0,clog,xksi)
-         xksimax=xksi
-C
-         clname='SED'//ctrcnm(jppoc)
-         clname1='SED'//ctrcnm(jpcal)
-         clname2='SED'//ctrcnm(jpsil)
-         CALL restget(nutrst,clname1,jpi,jpj,1,0,clog,sedcal(:,:))
-         CALL restget(nutrst,clname2,jpi,jpj,1,0,clog,sedsil(:,:))
-         CALL restget(nutrst,clname,jpi,jpj,1,0,clog,sedpoc(:,:))
-#     endif
-
-#if defined key_trc_hamocc3 || defined key_trc_pisces 
-
-C
-C  Initialization of chemical variables of the carbon cycle
-C  --------------------------------------------------------
-C
-        DO jk = 1,jpk
-          DO jj = 1,jpj
-            DO ji = 1,jpi
-              caralk = trn(ji,jj,jk,jptal)-
-     &                 borat(ji,jj,jk)/(1.+1.E-8/(rtrn+akb3(ji,jj,jk)))
-              co3(ji,jj,jk)=(caralk-trn(ji,jj,jk,jpdic))*tmask(ji,jj,jk)
-     &           +(1.-tmask(ji,jj,jk))*.5e-3
-              bicarb = (2.*trn(ji,jj,jk,jpdic)-caralk)
-              hi(ji,jj,jk) = (ak23(ji,jj,jk)*bicarb/co3(ji,jj,jk))
-     &           *tmask(ji,jj,jk)+(1.-tmask(ji,jj,jk))*1.e-9
-              h2co3(ji,jj) = 1.e-5
-            ENDDO
-          ENDDO
-        ENDDO
-#endif
-
-#else
-C
-C no passive tracers
-C
-#endif
-C
-#endif
-C
-      RETURN
-      END
-
Index: unk/NEMO/TOP_SRC/trcrst.mpp.h
===================================================================
--- /trunk/NEMO/TOP_SRC/trcrst.mpp.h	(revision 265)
+++ 	(revision )
@@ -1,110 +1,0 @@
-CCC $Header$ 
-CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-C ---------------------------------------------------------------------------
-CCC
-CCC                       trcrst.mpp.h
-CCC                     ***************
-CCC
-CDIR$ NOLIST
-      USE oce_trc
-      USE trc
-      USE sms
-      IMPLICIT NONE
-CDIR$ LIST
-#if defined key_passivetrc
-CC----------------------------------------------------------------------
-CC local declarations
-CC ==================
-      INTEGER ino0,ino1,it0,it1,idast1,jn,iarak0,iarak1
-
-C
-C
-C 0. initialisations
-C ------------------
-C
-      ino0 =no
-      it0  =nit000
-
-#if defined key_trc_cen2
-      iarak0=1
-#else
-      iarak0=0
-#endif
-C
-      IF(lwp) THEN
-          WRITE(numout,*) ' '
-          WRITE(numout,*) ' *** trcrst beginning of restart for'
-          WRITE(numout,*) ' passive tracer'
-          WRITE(numout,*) ' the present run :'
-          WRITE(numout,*) '   number job is  : ',no
-          WRITE(numout,*) '   with the time nit000 : ',nit000
-          IF(iarak0.eq.1) then
-              WRITE(numout,*) '   and before fields for Arakawa sheme '
-          ENDIF 
-          WRITE(numout,*) ' '
-      ENDIF
-C
-C 1. READ
-C -------
-C
-C 1.1 first record
-C
-      READ(nutrst,REC=1) ino1,it1,iarak1
-C
-      IF(lwp) THEN
-          WRITE(numout,*) ' '
-          WRITE(numout,*) ' READ nutrst with '
-          WRITE(numout,*) '   number job is  : ',ino1
-          WRITE(numout,*) '   with the time it : ',it1
-          IF(iarak1.eq.1) then
-              WRITE(numout,*) '   and before fields for Arakawa sheme '
-          ENDIF 
-          WRITE(numout,*) ' '
-      ENDIF
-C
-C 1.2 control of date
-C
-      IF((it0-it1).ne.1.and.abs(nrsttr).eq.1) THEN
-          IF(lwp) THEN
-              WRITE(numout,*) ' ===>>>> : problem with nit000 for the',
-     $            ' passive tracer restart'
-              WRITE(numout,*) ' =======                              ',
-     $            ' ======================'
-              WRITE(numout,*) ' we stop. verify the FILE'
-              WRITE(numout,*) ' or rerun with the value  0 for the'
-              WRITE(numout,*) ' control of time PARAMETER   nrstdt'
-              WRITE(numout,*) ' '
-          ENDIF
-          STOP 'trcrst'
-      ENDIF
-C
-C 1.3 Control of the sheme
-C
-      IF(iarak0.ne.iarak1) THEN
-          IF(lwp) THEN
-              WRITE(numout,*) ' ===>>>> : problem with the',
-     $            ' passive tracer restart file'
-              WRITE(numout,*) ' =======                              ',
-     $            ' ==========================='
-              WRITE(numout,*) ' we stop. verify the FILE'
-              WRITE(numout,*) ' before field required IF 1=',iarak0
-              WRITE(numout,*) ' before field present in file IF 1=',
-     $            iarak1
-              WRITE(numout,*) ' '
-          ENDIF
-          STOP 'trcrst'
-      ENDIF
-
-C
-C 1.4 READ prognostic variables
-C
-      DO jn=1,jptra
-        CALL read3(nutrst,trn(1,1,1,jn), jn + 1 )
-      END DO
-
-      DO jn=1,jptra
-        CALL read3(nutrst,trb(1,1,1,jn), jptra + jn + 1 )
-      END DO
-
-
Index: unk/NEMO/TOP_SRC/trcsms.F
===================================================================
--- /trunk/NEMO/TOP_SRC/trcsms.F	(revision 265)
+++ 	(revision )
@@ -1,152 +1,0 @@
-CCC $Header$ 
-CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-C ---------------------------------------------------------------------------
-CDIR$ LIST
-      SUBROUTINE trcsms( kt )
-CCC---------------------------------------------------------------------
-CCC
-CCC                       ROUTINE trcsms
-CCC                     *****************
-CCC
-CCC  PURPOSE :
-CCC  ---------
-CCC	time loop of opa for passive tracer
-CCC
-CC   METHOD :
-CC   -------
-CC      compute the well/spring evolution
-CC
-CC   INPUT :
-CC   -----
-CC      argument
-CC              ktask           : task identificator
-CC              kt              : time step
-CC      COMMON
-CC		all the COMMON defined in opa
-CC
-CC
-CC   OUTPUT :			: no
-CC   ------
-CC
-CC   WORKSPACE :
-CC   ---------
-CC
-CC   EXTERNAL :
-CC   --------
-CC      trcbio, trcsed, trcopt for NPZD or LOBSTER1 models
-CC
-CC      h3cprg for HAMOC3 and P3ZD models
-CC
-CC
-CC   MODIFICATIONS:
-CC   --------------
-CC      original  : 96-11
-CC      additions : 99-07 (M. Levy)
-CC                  04-00 (O. Aumont, M.A. Foujols) HAMOCC3 and P3ZD
-CC                  12-00 (O. Aumont, E. Kestenare) add trcexp for 
-CC                                                 instantaneous export
-CC----------------------------------------------------------------------
-CC parameters and commons
-CC ======================
-CDIR$ NOLIST
-      USE oce_trc
-      USE trc
-      IMPLICIT NONE
-CDIR$ LIST
-CC----------------------------------------------------------------------
-CC local declarations
-CC ==================
-      INTEGER kt
-
-#if defined key_passivetrc
-C
-C
-C
-C this ROUTINE is called only every ndttrc time step
-C
-
-      IF (mod(kt,ndttrc) .NE. 0) RETURN
-C
-C this first routines are parallelized on vertical slab
-C
-C
-#    if defined key_trc_npzd || defined key_trc_lobster1
-C
-C tracers: optical model
-C ----------------------
-C
-       CALL trcopt( kt)
-C
-C tracers: biological model
-C -------------------------
-C
-       CALL trcbio( kt)
-C
-C tracers: sedimentation model
-C ----------------------------
-C
-       CALL trcsed(kt)
-       CALL trcexp
-C
-#    elif defined key_trc_hamocc3
-C
-C    h3cprg: main PROGRAM for HAMOCC3 or P3ZD models
-C    -----------------------------------------------
-       CALL h3cprg(kt)
-#    elif defined key_trc_pisces
-C
-C    p4zprg: main PROGRAM for PISCES 
-C    -------------------------------
-       CALL p4zprg(kt)
-#    else
-C
-C SMS to DO
-C
-#if defined key_cfc
-C
-C   CFC's code taken from K. Rodgers
-C
-C   This part is still experimental
-C   -------------------------------
-C
-       CALL trcfreons(kt)
-C
-#elif defined key_gosac
-C
-C    C14 tracer modeled according to the OCMIP protocol
-C    This part has not been really tested in this configuration
-C
-C Compute the alkalinity
-C
-          CALL alkal(kt)
-C
-C  Compute the radioactive decay
-C
-          CALL decay
-C
-C Tracers: surface flux
-C
-          CALL dicgasflux(kt)
-
-#    endif
-
-#    endif
-
-# if defined key_trc_age
-C
-C Age Tracer 
-C
-          CALL trcage(kt)
-
-# endif
-
-#else
-C
-C no passive tracers
-C
-#endif
-C
-
-      RETURN
-      END
Index: unk/NEMO/TOP_SRC/trcwri.F
===================================================================
--- /trunk/NEMO/TOP_SRC/trcwri.F	(revision 265)
+++ 	(revision )
@@ -1,223 +1,0 @@
-CCC $Header$ 
-CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-C ---------------------------------------------------------------------------
-CDIR$ LIST
-      SUBROUTINE trcwri(kt)
-CCC---------------------------------------------------------------------
-CCC
-CCC                       ROUTINE trcwri
-CCC                     ******************
-CCC
-CCC  PURPOSE :
-CCC  ---------
-CCC     WRITE restart fields in nutwrs
-CCC
-CC   METHOD :
-CC   -------
-CC
-CC	nutwrs FILE:
-CC	each nstock time step , SAVE fields which are necessary for
-CC	passive tracer restart
-CC
-CC
-CC   INPUT :
-CC   -----
-CC      argument
-CC              kt              : time step
-CC      COMMON
-CC            /cottrc/          : passive tracers fields (before,now
-CC                                  ,after)
-CC
-CC   OUTPUT :
-CC   ------
-CC      FILE
-CC		nutwrs		: standard restart fields OUTPUT
-CC
-CC
-CC   WORKSPACE :
-CC   ---------
-CC      ji,jj,jk,jl,ino0,it0,iarak0
-CC
-CC   MODIFICATIONS:
-CC   --------------
-CC      original : 96-12
-CC      addition : 99-12 (M.-A. Foujols) NetCDF FORMAT with ioipsl
-CC      additions : 00-05 (A. Estublier) 
-CC                  TVD Limiter Scheme : key_trc_tvd
-CC      additions : 01-01 (M.A Foujols, E. Kestenare) bug fix: restclo
-CC      additions : 01-01 (O. Aumont, E. Kestenare) 
-CC                  write restart file for sediments
-CC      additions : 01-05 (O. Aumont, E. Kestenare) 
-CC                  write restart file for calcite and silicate sediments
-CC----------------------------------------------------------------------
-CC parameters and commons
-CC ======================
-#if defined key_mpp
-c no ioipsl
-#include "trcwri.mpp.h"
-#    else
-CDIR$ NOLIST
-      USE ioipsl 
-      USE oce_trc
-      USE trc
-      USE sms
-      IMPLICIT NONE
-CDIR$ LIST
-CC----------------------------------------------------------------------
-CC local declarations
-CC ==================
-      INTEGER kt
-
-#if defined key_passivetrc
-
-      INTEGER jn
-      INTEGER ino0,it0,iarak0
-      INTEGER ic,jc,ji,jj,jk
-      INTEGER iyear, imonth, iday, iy, itime
-      REAL zsec, zdate0, zdt, zinfo(3),zdiag_var,zdiag_varmin,
-     $     zdiag_varmax
-      CHARACTER*50 clname,clname1,clname2,cln,clnet
-      LOGICAL clbon
-
-C
-C
-C 1. OUTPUT of restart fields (nutwrs)
-C ---------------------------
-C
-      IF( (mod(kt,nstock).eq.0) .OR. (kt.eq.nitend) ) THEN
-C
-C 0. initialisations
-C ------------------
-C
-          IF(lwp) THEN
-              WRITE(numout,*) ' '
-              WRITE(numout,*)
-     $            ' trcwri: restart OUTPUT done in nutwrs = ',nutwrs
-     $            ,' at it= ',kt,' date= ',ndastp
-              WRITE(numout,*) ' -------'
-          ENDIF
-C
-          ino0 =no
-          it0  =kt
-
-#if defined key_trc_cen2 || defined key_trc_tvd
-          iarak0=1
-#else
-          iarak0=0
-#endif
-C
-C 1. WRITE in nutwrs
-C
-C ... first information
-C
-          INQUIRE (FILE=trestart,EXIST=clbon)
-          IF(clbon) THEN
-              OPEN(UNIT=numwrs,FILE=trestart,STATUS='old')
-              CLOSE(numwrs,STATUS='delete')
-          ENDIF
-
-          iyear = ndastp/10000
-          imonth = ndastp/100 - iyear*100
-          iday = ndastp - imonth*100 - iyear*10000
-          iyear = ndastp/10000
-          zsec=0.
-          iy=iyear-(iyear/100)*100
-          ic=1
-          DO jc=1,16
-            IF(cexper(jc:jc).ne.' ') ic=jc
-          END DO
-          WRITE(cln,'("_",i2.2,i2.2,i2.2,"_restart.trc")') iy,imonth
-     $        ,iday 
-          clname=cexper(1:ic)//cln
-          ic=1
-          DO jc=1,48
-            IF(clname(jc:jc).ne.' ') ic=jc
-          END DO
-          trestart=clname(1:ic)//".nc"
-          itime=0
-          CALL ymds2ju(iyear,imonth,iday,zsec,zdate0)
-          CALL restini('NONE',jpi,jpj,glamt,gphit,jpk,gdept,clname
-     $        ,itime,zdate0,rdt,nutwrs)
-          zinfo(1)=FLOAT(ino0)
-          zinfo(2)=FLOAT(it0)
-          zinfo(3)=FLOAT(iarak0)
-          CALL restput(nutwrs,'info',1,1,3,0,zinfo)
-C
-C prognostic variables
-C
-
-         DO jn=1,jptra
-            clname='TRN'//ctrcnm(jn)
-            CALL restput(nutwrs,clname,jpi,jpj,jpk,0,trn(:,:,:,jn))
-
-         zdiag_var=0.
-         zdiag_varmin=0.
-         zdiag_varmax=0.
-         WRITE(numout,*) '----TRACER STAT----'
-
-         DO ji=1,jpi
-            DO jj=1,jpj
-              DO jk=1,jpk
-
-              zdiag_var=zdiag_var+tmask(ji,jj,jk)*trn(ji,jj,jk,jn)
-              
-              IF (tmask(ji,jj,jk).EQ.1.) THEN
-
-                  IF (zdiag_varmin.GT.trn(ji,jj,jk,jn)) 
-     $                zdiag_varmin =  trn(ji,jj,jk,jn)
-                  IF (zdiag_varmax.LT.trn(ji,jj,jk,jn))
-     $                zdiag_varmax =  trn(ji,jj,jk,jn)
-            
-              ENDIF
-
-              END DO
-            END DO
-         END DO
-      
-         zdiag_var=zdiag_var/(jpi*jpj*jpk)
-         
-         WRITE(numout,*) 'MEAN NO ',jn,' =',zdiag_var,'MIN= '
-     $       ,zdiag_varmin,'MAX= ',zdiag_varmax
-
-         END DO
-
-
-
-          DO jn=1,jptra
-            clname='TRB'//ctrcnm(jn)
-            CALL restput(nutwrs,clname,jpi,jpj,jpk,0,trb(:,:,:,jn))
-          END DO
-
-#         if defined key_trc_hamocc3
-             clname='SED'//ctrcnm(jppoc)
-             clname1='SED'//ctrcnm(jpcal)
-             clname2='SED'//ctrcnm(jpsil)
-             CALL restput(nutwrs,clname1,jpi,jpj,1,0,sedcal(:,:))
-             CALL restput(nutwrs,clname2,jpi,jpj,1,0,sedsil(:,:))
-             CALL restput(nutwrs,clname,jpi,jpj,1,0,sedpoc(:,:))
-#         elif defined key_trc_npzd || defined key_trc_lobster1
-             clname='SED'//ctrcnm(jpdet)
-             CALL restput(nutwrs,clname,jpi,jpj,1,0,sedpoc(:,:))
-#         elif defined key_trc_pisces
-             clname='SED'//ctrcnm(jppoc)
-             clname1='SED'//ctrcnm(jpcal)
-             clname2='SED'//ctrcnm(jpsil)
-             CALL restput(nutwrs,clname1,jpi,jpj,1,0,sedcal(:,:))
-             CALL restput(nutwrs,clname2,jpi,jpj,1,0,sedsil(:,:))
-             CALL restput(nutwrs,clname,jpi,jpj,1,0,sedpoc(:,:))
-#         endif
-
-          CALL restclo(nutwrs)
-
-      ENDIF
-C
-#else
-C
-C no passive tracers
-C
-#endif
-C
-#endif
-      RETURN
-      END
Index: unk/NEMO/TOP_SRC/trcwri.mpp.h
===================================================================
--- /trunk/NEMO/TOP_SRC/trcwri.mpp.h	(revision 265)
+++ 	(revision )
@@ -1,68 +1,0 @@
-CCC $Header$ 
-CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-C ---------------------------------------------------------------------------
-CCC
-CCC                       trcwri.mpp
-CCC                     **************
-CCC
-CDIR$ NOLIST
-      USE oce_trc
-      USE trc
-      IMPLICIT NONE
-
-CDIR$ LIST
-CC----------------------------------------------------------------------
-CC local declarations
-CC ==================
-      INTEGER kt
-
-#if defined key_passivetrc
-
-      INTEGER jn
-      INTEGER ino0,it0,iarak0
-
-C
-C
-C 1. OUTPUT of restart fields (nutwrs)
-C ---------------------------
-C
-      IF( (mod(kt,nstock).eq.0) .OR. (kt.eq.nitend) ) THEN
-C
-C 1.0 initializations
-C
-          IF(lwp) THEN
-              WRITE(numout,*) ' '
-              WRITE(numout,*)
-     $            ' trcwri: restart OUTPUT done in nutwrs = ',nutwrs
-     $            ,' at it= ',kt,' date= ',ndastp
-              WRITE(numout,*) ' -------'
-          ENDIF
-C
-          ino0 =no
-          it0  =kt
-
-#if defined key_trc_cen2
-          iarak0=1
-#else
-          iarak0=0
-#endif
-C
-C 1.1 WRITE
-C
-C first record
-C
-          IF(lwp) WRITE (nutwrs,REC=1) ino0,it0,iarak0
-C
-C prognostic variables
-C
-
-          DO jn=1,jptra
-            CALL write3(nutwrs,trn(1,1,1,jn), jn + 1 )
-          END DO
-
-          DO jn=1,jptra
-            CALL write3(nutwrs,trb(1,1,1,jn), jptra + jn + 1 )
-          END DO
-      ENDIF
-C
