source: IOIPSL/trunk/example/testrest.f90 @ 122

Last change on this file since 122 was 122, checked in by bellier, 15 years ago

JB: some cleaning (-> fortran 90)

  • Property svn:keywords set to Id
File size: 3.6 KB
Line 
1PROGRAM testrest
2!-
3!$Id$
4!---------------------------------------------------------------------
5!- This program provide a an example of the basic usage of REST.
6!- Here the test the time sampling and averaging. Thus a long
7!- time-series is produced and sampled in different ways.
8!---------------------------------------------------------------------
9  USE ioipsl
10!
11  IMPLICIT NONE
12!
13  INTEGER :: iim,jjm,llm
14  PARAMETER (iim=12,jjm=10,llm=2)
15!
16  REAL :: champ1(iim,jjm,llm),champ2(iim,jjm,llm+1),champ3(iim,jjm,llm)
17  REAL :: champ4(iim,jjm)
18  REAL :: champ_read(iim,jjm,llm)
19  REAL :: lon(iim,jjm),lat(iim,jjm),lev(llm)
20  REAL :: x
21!
22  INTEGER :: i,j,l,fid,t,ij,sig_id,hori_id,it
23  INTEGER :: day=1,month=1,year=1997
24  INTEGER :: itau=1,start,INDEX(1)
25!
26  REAL :: julday,un_mois,un_an
27  REAL :: deltat=86400,dt_wrt,dt_op,dt_wrt2,dt_op2
28  CHARACTER*20 :: fnamein,fnameout,keyword,value
29!
30  REAL :: pi=3.1415
31!---------------------------------------------------------------------
32!-
33! 0.0 Choose a  gregorian calendar
34!-
35  CALL ioconf_calendar ('gregorian')
36!-
37! 1.0 Define a few variables we will need.
38!     These are the coordinates the file name and the date.
39!-
40  DO i=1,iim
41    DO j=1,jjm
42      lon(i,j) = &
43 &     ((float(iim/2)+0.5)-float(i))*pi/float(iim/2)*(-1.)*180./pi
44      lat(i,j) = &
45 &     (180./pi)*ASIN(((float(jjm/2)+0.5)-float(j))/float(jjm/2))
46    ENDDO
47  ENDDO
48!-
49  DO l=1,llm
50    lev(l) = float(l)/llm
51  ENDDO
52!-
53! 1.1 The chosen date is 15 Feb. 1997 as stated above. It has to be
54!     transformed into julian days using the calendar provided by
55!     IOIPSL.
56!-
57  CALL ymds2ju (year,month,day,0.,julday)
58  CALL ioget_calendar (un_an)
59  un_mois = un_an/12.
60  dt_wrt = un_mois*deltat
61  dt_op = deltat
62  dt_wrt2 = -1.
63  dt_op2 = deltat
64!-
65  fnamein = 'NONE'
66  fnameout = 'restfile'
67!-
68! 2.0 Create a restart file from nothing !
69!-
70  CALL restini (fnamein,iim,jjm,lon,lat,llm,lev,fnameout, &
71 &              itau,julday,deltat,fid)
72!-
73  champ1(:,:,:) = ASIN(1.0)
74  champ2(:,:,:) = EXP(ASIN(1.0))
75!-
76  CALL ioconf_setatt ('units','?')
77  CALL ioconf_setatt ('long_name','Tests 1 for a real variable')
78  CALL restput (fid,'test1',iim,jjm,llm,itau,champ1)
79!-
80  CALL ioconf_setatt ('units','?')
81  CALL ioconf_setatt ('long_name','Tests 2 for a real variable')
82  CALL restput (fid,'test2',iim,jjm,llm+1,itau,champ2)
83!-
84  CALL restclo ()
85!-
86  WRITE(*,*) '============== FIRST FILE CLOSED =============='
87!-
88!  3.0 Reopen the restart file and check that the values read are correct
89!-
90  fnamein = 'restfile'
91  fnameout = 'restfilebis'
92!-
93  CALL restini (fnamein,iim,jjm,lon,lat,llm,lev,fnameout, &
94 &              itau,julday,deltat,fid)
95!-
96  CALL restget (fid,'test1',iim,jjm,llm,itau,.FALSE.,champ_read)
97!-
98  itau = itau+10
99  CALL restput (fid,'test1',iim,jjm,llm,itau,champ_read)
100  CALL restput (fid,'test2',iim,jjm,llm+1,itau,champ2)
101!-
102  itau = itau+10
103  champ3(:,:,:) = champ_read(:,:,:)+champ2(:,:,1:llm) 
104  CALL restput (fid,'test1',iim,jjm,llm,itau,champ3)
105!-
106  CALL restclo ()
107!-
108  WRITE(*,'(a25,e36.30)') 'The input data : ',champ1(1,1,1)
109  WRITE(*,'(a25,e36.30)') 'The restart data : ',champ_read(1,1,1)
110!-
111!  4.0 Reopen the restart file and add another time step
112!-
113  fnamein = 'restfilebis'
114  fnameout = 'restfilebis'
115!-
116  CALL restini (fnamein,iim,jjm,lon,lat,llm,lev,fnameout, &
117 &              itau,julday,deltat,fid)
118!-
119  itau = itau+10
120  CALL restput (fid,'test1',iim,jjm,llm,itau,champ1)
121  CALL ioconf_setatt ('units','?')
122  CALL ioconf_setatt ('long_name', &
123 &                    'Test a variable with another dimension')
124  CALL restput (fid,'test4',iim,jjm,0,itau,champ4)
125!-
126  CALL restclo ()
127!-------------------
128END PROGRAM testrest
Note: See TracBrowser for help on using the repository browser.