1 | program elc |
---|
2 | c Header(1)=Rlat |
---|
3 | c Header(2)=Rlon |
---|
4 | c Header(3)=IASI_Satellite_Zenith(J,I) |
---|
5 | c Header(4)=Land_Sea_Alti(J,I) |
---|
6 | c Header(5)=cloud_cover(J,I) |
---|
7 | c Header(6)=hh*1.e4+mn*1.e2+ss |
---|
8 | c Header(7)=float(indtm*10+indjn) |
---|
9 | c Header(8)=yyyymmdd |
---|
10 | real head(8),mots(339) |
---|
11 | real ang_sol,norb |
---|
12 | integer indi(339) |
---|
13 | character*200 filein |
---|
14 | c header : rlat,rlon,angle_zen, |
---|
15 | call getarg(1,filein) |
---|
16 | open(1,file=filein,form='unformatted') |
---|
17 | read(1) nbmots,(indi(i),i=1,nbmots) |
---|
18 | print *,'nbmots',nbmots,(indi(i),i=1,nbmots) |
---|
19 | |
---|
20 | ! nbre de temperatures de brillance extraites = nbmots - nbre mots header - angle solaire - numero orbite |
---|
21 | nb_TBs=nbmots-8-2 |
---|
22 | write(*,*) nb_TBs |
---|
23 | icpt=0 |
---|
24 | ic=0 |
---|
25 | imn=0 |
---|
26 | imj=0 |
---|
27 | itn=0 |
---|
28 | itj=0 |
---|
29 | ! ixx calcule le nombre de points dans 4 classes combinant les indices terre/mer/jour/nuit |
---|
30 | * xx=mn pour mer/nuit |
---|
31 | 10 read(1,end=100) head,(mots(i),i=1,nb_TBs),ang_sol,norb |
---|
32 | print * ,"head=",nb_TBs |
---|
33 | sectta=1./cosd(head(3)) |
---|
34 | !write(*,1) head(1),head(6),head(8),mots(1),ang_sol,int(norb) |
---|
35 | 1 format(f8.2,2f10.0,2f8.2,i7) |
---|
36 | if(head(7).eq.0) imn=imn+1 |
---|
37 | if(head(7).eq.1) imj=imj+1 |
---|
38 | if(head(7).eq.10) itn=itn+1 |
---|
39 | if(head(7).eq.11) itj=itj+1 |
---|
40 | icpt=icpt+1 |
---|
41 | ic=ic+1 |
---|
42 | goto 10 |
---|
43 | 100 continue |
---|
44 | 14 format(8f11.3,2i4) |
---|
45 | 15 format(3f11.3,2i4) |
---|
46 | print *,icpt,imn,imj,itn,itj |
---|
47 | end |
---|