source: LMDZ.3.3/branches/rel-LF/libf/dyn3d/defrun.F @ 517

      SUBROUTINE defrun(tapedef,creinit)
      IMPLICIT NONE
c-----------------------------------------------------------------------
c   declarations:
c   -------------

#include "dimensions.h"
#include "paramet.h"
#include "control.h"
#include "logic.h"
#include "serre.h"
#include "comdissip.h"
      LOGICAL creinit

c   local:
c   ------

      CHARACTER ch1*72,ch2*72,ch3*72,ch4*8
      INTEGER tapedef,tapeout
      REAL clonn,clatt,alphaxx,alphayy
c
c-----------------------------------------------------------------------
c   initialisations:
c   ----------------

      tapeout=6

c-----------------------------------------------------------------------
c  Parametres de controle du run:
c-----------------------------------------------------------------------

      READ (tapedef,9000) ch1,ch2,ch3
      WRITE(tapeout,9000) ch1,ch2,ch3

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    nday
      WRITE(tapeout,9001) ch1,'nday'
      WRITE(tapeout,*)    nday

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    day_step
      WRITE(tapeout,9001) ch1,'day_step'
      WRITE(tapeout,*)    day_step

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    iperiod
      WRITE(tapeout,9001) ch1,'iperiod'
      WRITE(tapeout,*)    iperiod

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    iconser
      WRITE(tapeout,9001) ch1,'iconser'
      WRITE(tapeout,*)    iconser

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    iecri
      WRITE(tapeout,9001) ch1,'iecri'
      WRITE(tapeout,*)    iecri

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    periodav
      WRITE(tapeout,9001) ch1,'periodav'
      WRITE(tapeout,*)    periodav

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    idissip
      WRITE(tapeout,9001) ch1,'idissip'
      WRITE(tapeout,*)    idissip

ccc  ....   P. Le Van , ajout  le 5/09/95 .pour la dissipation  ...
ccc
      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    lstardis
      WRITE(tapeout,9001) ch1,'lstardis'
      WRITE(tapeout,*)    lstardis

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    niterdis
      WRITE(tapeout,9001) ch1,'niterdis'
      WRITE(tapeout,*)    niterdis

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    coefdis
      WRITE(tapeout,9001) ch1,'coefdis'
      WRITE(tapeout,*)    coefdis

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    tetavel
      WRITE(tapeout,9001) ch1,'tetavel'
      WRITE(tapeout,*)    tetavel

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    tetatemp
      WRITE(tapeout,9001) ch1,'tetatemp'
      WRITE(tapeout,*)    tetatemp

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    iphysiq
      WRITE(tapeout,9001) ch1,'iphysiq'
      WRITE(tapeout,*)    iphysiq

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    purmats
      WRITE(tapeout,9001) ch1,'purmats'
      WRITE(tapeout,*)    purmats

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    physic
      WRITE(tapeout,9001) ch1,'physic'
      WRITE(tapeout,*)    physic

ccc   .... P.Le Van, ajout le 03/01/96 pour l'ecriture phys ...
c
      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    ecritphy
      WRITE(tapeout,9001) ch1,'ecritphy'
      WRITE(tapeout,*)    ecritphy

ccc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...

      IF( creinit ) go to 100

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    clonn
      WRITE(tapeout,9001) ch1,'clon'
      WRITE(tapeout,*)    clonn
      IF( clon - clonn .gt.0.001)  THEN
        PRINT *,' La valeur de clon passee par run.def est differente de 
     *  celle lue sur le fichier  start '
        STOP
      ENDIF

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    clatt
      WRITE(tapeout,9001) ch1,'clat'
      WRITE(tapeout,*)    clatt

      IF( clat -clatt.gt.0.001 )  THEN
        PRINT *,' La valeur de clat passee par run.def est differente de 
     *  celle lue sur le fichier  start '
        STOP
      ENDIF

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    alphaxx
      WRITE(tapeout,9001) ch1,'alphax'
      WRITE(tapeout,*)    alphaxx

      IF( alphax-alphaxx.gt.0.001 )  THEN
        PRINT *,' La valeur de alphax passee par run.def est differente  
     *  de celle lue sur le fichier  start '
        STOP
      ENDIF

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    alphayy
      WRITE(tapeout,9001) ch1,'alphay'
      WRITE(tapeout,*)    alphayy

      IF( alphay-alphayy.gt.0.001 )  THEN
        PRINT *,' La valeur de alphay passee par run.def est differente  
     * de celle lue sur le fichier  start '
        STOP
      ENDIF

      RETURN
c
c    
100   CONTINUE
      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    clon
      WRITE(tapeout,9001) ch1,'clon'
      WRITE(tapeout,*)    clon

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    clat
      WRITE(tapeout,9001) ch1,'clat'
      WRITE(tapeout,*)    clat

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    alphax
      WRITE(tapeout,9001) ch1,'alphax'
      WRITE(tapeout,*)    alphax

      READ (tapedef,9001) ch1,ch4
      READ (tapedef,*)    alphay
      WRITE(tapeout,9001) ch1,'alphay'
      WRITE(tapeout,*)    alphay

c-----------------------------------------------------------------------

      RETURN
9000  FORMAT(3(/,a72))
9001  FORMAT(/,a72,/,a8)
      END