igsm/src/robso3.F


#include "ctrparam.h"

      subroutine  robso3(iyearr1)
#include "BD2G04.COM" 
#include "chem_para"
#include "TEM.COM"
      character * 120 flo3,o3datadir
      character * 4 year
      character * 6 name
      data ifirst /1/
      iyearr=iyearr1
      if(ifirst.eq.1) then
       print *,'FROM ROBSO3'
       print *,'O3DATADIR=',o3datadir
       if(CLIMO3) print *,'With fixed O3 for 1860'
c      stop
       id2=index(o3datadir," ")
       print *,'O3DATADIR=',o3datadir(1:id2-1)
       ifirst=0
      endif
        if(CLIMO3) iyearr=1860
        write (year,'(i4)'),iyearr
        flo3=o3datadir(1:id2-1)//'o3bau'//year//'.dat46'
        print *,flo3
        open(561,file=flo3,status='old')
        do mn=1,12
        read (561,*),iiyear,nm1,name
        print *,'O3 for ',iiyear,nm1
        do j=1,jm0
         read (561,*),(obso3(n,j,mn),n=1,8)
        enddo
        enddo
        close (561)
      return
      end
1
2	#include "ctrparam.h"
3
4	subroutine robso3(iyearr1)
5	#include "BD2G04.COM"
6	#include "chem_para"
7	#include "TEM.COM"
8	character * 120 flo3,o3datadir
9	character * 4 year
10	character * 6 name
11	data ifirst /1/
12	iyearr=iyearr1
13	if(ifirst.eq.1) then
14	print *,'FROM ROBSO3'
15	print *,'O3DATADIR=',o3datadir
16	if(CLIMO3) print *,'With fixed O3 for 1860'
17	c stop
18	id2=index(o3datadir," ")
19	print *,'O3DATADIR=',o3datadir(1:id2-1)
20	ifirst=0
21	endif
22	if(CLIMO3) iyearr=1860
23	write (year,'(i4)'),iyearr
24	flo3=o3datadir(1:id2-1)//'o3bau'//year//'.dat46'
25	print *,flo3
26	open(561,file=flo3,status='old')
27	do mn=1,12
28	read (561,*),iiyear,nm1,name
29	print *,'O3 for ',iiyear,nm1
30	do j=1,jm0
31	read (561,*),(obso3(n,j,mn),n=1,8)
32	enddo
33	enddo
34	close (561)
35	return
36	end