| 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 |
|---|
