/[MITgcm]/MITgcm_contrib/jscott/igsm/src/robso3.F
ViewVC logotype

Annotation of /MITgcm_contrib/jscott/igsm/src/robso3.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.2 - (hide annotations) (download)
Mon Apr 23 21:20:18 2007 UTC (18 years, 3 months ago) by jscott
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +1 -1 lines
bring igsm atmos code up to date

1 jscott 1.1
2     #include "ctrparam.h"
3    
4     subroutine robso3(iyearr1)
5     #include "BD2G04.COM"
6     #include "chem_para"
7 jscott 1.2 #include "TEM.h"
8 jscott 1.1 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

  ViewVC Help
Powered by ViewVC 1.1.22