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

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

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


Revision 1.1 - (hide annotations) (download)
Fri Aug 11 19:35:32 2006 UTC (18 years, 11 months ago) by jscott
Branch: MAIN
CVS Tags: HEAD
atm2d package

1 jscott 1.1 #include "ctrparam.h"
2    
3     subroutine read_staer (nydata,amn)
4     #include "BD2G04.COM"
5     parameter( naery=2000)
6     character * 40 title
7     character * 80 title1
8     character * 120 fl_volaer
9     character *5 xxx
10     character *41 fordir
11     dimension amn(jm0,12,naery),amnglav(12,naery)
12     common/vaerosol/fl_volaer
13     SS=0.0
14     do j=1,jm
15     SS=SS+DXYP(j)
16     enddo
17     do ii =1,25
18     print *,fl_volaer
19     enddo
20     open (560,file=fl_volaer,form='formatted',status='old')
21     c read(560,'(a40')title
22     c print *,title
23     c read(560,'(a80')title1
24     c print *,title1
25     do ny=1,naery
26     do j=1,jm
27     read (560,101,end=500),(amn(j,m,ny),m=1,12)
28     enddo
29     enddo
30     500 continue
31     nydata=ny-1
32     print *,'Volcanic aerosol data for year =',1850+ny-2
33     do n=1,nydata
34     do nm=1,12
35     amnglav(nm,n)=0.0
36     do j=1,jm
37     amnglav(nm,n)=amnglav(nm,n)+amn(j,nm,n)*DXYP(j)
38     enddo
39     enddo
40     print 102,1849+n,(amnglav(nm,n)/SS,nm=1,12)
41     enddo
42     close (560)
43     101 format(12f7.4)
44     102 format(i5,12f7.4)
45     return
46     end
47    

  ViewVC Help
Powered by ViewVC 1.1.22