source: trunk/SOURCES/limit_file.f90 @ 23

Last change on this file since 23 was 4, checked in by dumas, 10 years ago

initial import GRISLI trunk

File size: 2.0 KB
Line 
1!> \file limit_file.f90
2!! Lecture dans un petit fichier au cas ou les parametres du run seraient changes
3!<
4
5!> SUBROUTINE: limit_file()
6!!
7!!Ecrit le fichier limit
8!! \author ...
9!! \date ...
10!! @param  nt               
11!! @param  time               
12!! @param  dt           
13!! @param  tend                       
14!! @param  dtsortie
15!! @param  dtcpt
16!! @param  testdiag
17!! @param  dtt
18!! @param  runname
19!! @return
20
21
22subroutine limit_file(nt,time,dt,tend,dtsortie,dtcpt,testdiag,dtt,runname)
23
24! ecrit le fichier limit
25
26implicit none
27INTEGER :: nt
28REAL :: dt
29REAL :: tend
30REAL :: dtsortie
31REAL :: dtcpt
32REAL :: testdiag
33REAL :: dtt
34REAL :: time
35CHARACTER(len=8) :: RUNNAME 
36      integer :: num_limit = 99
37
38
39
40!--------------------------------------------------------------------!
41!     lecture dans un petit fichier au cas ou les parametres du run  !
42!     seraient changes                                               !
43!--------------------------------------------------------------------!
44      if (time.ge.5000.)  DTSORTIE=10.
45      if (time.ge.6000.)  DTSORTIE=100.
46      if (time.ge.20000.) DTSORTIE=500.
47
48     if (nt.eq.1) then
49!      if ((nt.eq.1).or.(abs(time-5000.).lt.5.).or.(abs(time-6000.).lt.5.)   &
50!           .or.(abs(time-20000.).lt.5.)) then
51        open(num_limit,file=runname//'.limit')
52        write(num_limit,*) TEND,DTSORTIE,DTCPT,testdiag,dtt
53        write(num_limit,*) 'TEND,DTSORTIE,DTCPT,testdiag,dtt'
54        write(num_limit,*)'time=', time,'   dt=',dt
55        close(num_limit)
56     endif
57
58
59! lecture au cas ou cela ait change
60     open(num_limit,file=runname//'.limit')
61     read(num_limit,*,end=07,err=07) TEND,DTSORTIE,dtcpt, &
62          testdiag,dtt
6307     close(num_limit)
64
65! ecriture
66     open(num_limit,file=runname//'.limit')
67     write(num_limit,*) TEND,DTSORTIE,DTCPT,testdiag,dtt
68     write(num_limit,*) 'TEND,DTSORTIE,DTCPT,testdiag,dtt'
69     write(num_limit,*)'time=', time,'   dt=',dt
70     close(num_limit)
71
72
73     return
74     end subroutine 
75!--------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.