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	jscott	1.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