slide1
Download
Skip this Video
Download Presentation
Observation pre-processing in the ARPEGE/ALADIN model

Loading in 2 Seconds...

play fullscreen
1 / 21

Observation pre-processing in the ARPEGE/ALADIN model - PowerPoint PPT Presentation


  • 104 Views
  • Uploaded on

Observation pre-processing in the ARPEGE/ALADIN model. Roger Randriamampianina Hungarian Meteorological Service ALADIN 3D-Var/ODB training course, Budapest, 6-10 June, 2006. Outline. Local pre-processing:. OULAN . BATOR. 1C RAD. BIAS CORRECTION. BATOR .

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about ' Observation pre-processing in the ARPEGE/ALADIN model' - zev


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
slide1

Observation pre-processing in theARPEGE/ALADIN model

Roger Randriamampianina

Hungarian Meteorological Service

ALADIN 3D-Var/ODB training course, Budapest, 6-10 June, 2006

slide2

Outline

Local pre-processing:

OULAN

BATOR

1C RAD. BIAS CORRECTION

slide4

The bator procedure

1 – Initialisation of the observation errors It contains the observation error for several observation types

2 –Blacklisting of the observations

1 – Creation of (an) ODB database from the OBSOUL file(s)

slide5

Bator and the observation errors

ECTERO(NSYNOP,1,1:5,1)=(/Z_VAL,1.4_JPRB,2.0_JPRB,0.1_JPRB,1.5_JPRB /) ! synop

ECTERO(NSYNOP,2,1:5,1)=(/Z_VAL,1.4_JPRB,3.0_JPRB,0.1_JPRB,1.5_JPRB /) ! ship

ECTERO(NDRIBU,1,1:5,1)=(/Z_VAL,1.4_JPRB,3.0_JPRB,Z_V99,1.5_JPRB /) ! buoy

ECTERO(NSATOB,1,3,1:19)=REF_STAT(1:19,4)*1.1_JPRB ! Meteosat

ECTERO(NSATOB,2,3,1:19)=REF_STAT(1:19,4)*1.1_JPRB ! Goes

ECTERO(NSATOB,3,3,1:19)=REF_STAT(1:19,4)*1.3_JPRB ! sat. japonais

ECTERO(NSATOB,4,3,1:19)=REF_STAT(1:19,4)*1.4_JPRB ! sat. indien

ECTERO(NSATOB,5,3,1:19)=REF_STAT(1:19,4)*1.5_JPRB ! autre

ECTERO(NAIREP,1,3, 1:19)=(/ 2.5_JPRB, 2.5_JPRB, 2.5_JPRB, 2.6_JPRB, &! airep vent

& 2.7_JPRB, 2.8_JPRB, 3.2_JPRB, 3.4_JPRB, 3.8_JPRB, 4.0_JPRB, &

& 4.1_JPRB, 4.2_JPRB, 4.2_JPRB, 4.2_JPRB, 4.2_JPRB, 4.3_JPRB, &

& 4.6_JPRB, 4.8_JPRB, 5.0_JPRB /)

ECTERO(NAIREP,1,2, 1:19)=(/ 1.9_JPRB, 1.8_JPRB, 1.7_JPRB, 1.6_JPRB, &! airep tempe

& 1.5_JPRB, 1.4_JPRB, 1.4_JPRB, 1.4_JPRB, 1.4_JPRB, 1.5_JPRB, &

& 1.6_JPRB, 1.7_JPRB, 1.8_JPRB, 1.9_JPRB, 2.0_JPRB, 2.1_JPRB, &

& 2.2_JPRB, 2.3_JPRB, 2.4_JPRB /)

ECTERO(NAIREP,1,2,1:19) = ECTERO(NAIREP,1,2,1:19) / 1.1_JPRB

ECTERO(NTEMP,1,3, 1:19) =(/ 2.3_JPRB, 2.3_JPRB, 2.3_JPRB, 2.4_JPRB, &! temp vent

& 2.5_JPRB, 2.5_JPRB, 2.8_JPRB, 3.0_JPRB, 3.3_JPRB, 3.6_JPRB, &

& 3.7_JPRB, 3.8_JPRB, 3.8_JPRB, 3.8_JPRB, 3.8_JPRB, 3.9_JPRB, &

& 4.1_JPRB, 4.3_JPRB, 4.5_JPRB /)

ECTERO(NTEMP,1,2, 1:19 )=(/ 1.7_JPRB, 1.6_JPRB, 1.5_JPRB, 1.4_JPRB, &! temp tempe

& 1.4_JPRB, 1.3_JPRB, 1.3_JPRB, 1.3_JPRB, 1.3_JPRB, 1.4_JPRB, &

& 1.5_JPRB, 1.5_JPRB, 1.6_JPRB, 1.6_JPRB, 1.6_JPRB, 1.7_JPRB, &

& 1.8_JPRB, 1.9_JPRB, 2.0_JPRB /)

ECTERO(NTEMP,1,2,1:19) = ECTERO(NTEMP,1,2,1:19) / 1.2_JPRB

ECTERO(NTEMP,1,1, 1:19)=(/ 8.0_JPRB, 8.2_JPRB, 8.6_JPRB, &! temp geop

& 9.0_JPRB, 9.4_JPRB, 9.9_JPRB, 11.4_JPRB, 12.7_JPRB, &

& 14.0_JPRB, 16.0_JPRB, 15.7_JPRB, 17.2_JPRB, 20.1_JPRB, &

& 22.0_JPRB, 24.4_JPRB, 27.0_JPRB, 30.0_JPRB, 31.5_JPRB, &

bator_util.F90

slide6

Bator and blacklisting

ival=-1

if (LN_DIAP(jst)%code == \'PROF\') ival=jst

if ( (CLSID(kobs) == LN_DIAP(jst)%cident).and. &

& ((kcodmes == LN_DIAP(jst)%icodmes).or.(icmd == LN_DIAP(jst)%icodmes)) ) then

klndiap( LN_DIAP(jst)%ivarno ) = ival

if ( kotp == NSYNOP ) then !------- SYNOP -------------

if ( LN_DIAP(jst)%ivarno == NVNUMB(10) ) then ! temp 2m

klndiap( NVNUMB( 7) ) = ival ! hum 2m

klndiap( NVNUMB(63) ) = ival ! hum spec.

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(7) ) then ! hum 2m

klndiap( NVNUMB(63) ) = ival

endif

ELSEIF ( kotp == NTEMP ) then !------- NTEMP -------------

if ( LN_DIAP(jst)%ivarno == NVNUMB(3) ) then ! geop.

klndiap( NVNUMB( 5) ) = ival ! hum

klndiap( NVNUMB( 8) ) = ival ! temp

klndiap( NVNUMB( 9) ) = ival ! temp pt de rosee

klndiap( NVNUMB(58) ) = ival ! hum. de couche

klndiap( NVNUMB(63) ) = ival ! hum spec.

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(8) ) then ! temp

klndiap( NVNUMB( 5) ) = ival

klndiap( NVNUMB(63) ) = ival

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(5) ) then ! hum

klndiap( NVNUMB(63) ) = ival

endif

endif

bator_saisies.F90

1 SHIP 21 1 62301 01022003

1 SYNOP 14 1 71094 01032004

1 SYNOP 11 39 03590 20050718

2 ACAR 145 2 1RYFVQBA 01092003

2 AMDAR 144 2 EU3781 01072001

4 BATHY 63 39 ZSAF 01032004

4 BUOY 165 1 17546 01122000

5 TEMP 35 2 42314 10062003

5 TEMP 35 3 42339 01121997

5 TEMP 35 1 01001 18052005

6 EUROPROFIL 134 3 ABWWP 01112003

6 PILOT 32 3 07162 01112001

6 PROFILER 34 4 74630 22032004

slide7

bator_saisies.F90

1 SHIP 21 1 62301 01022003

1 SYNOP 14 1 71094 01032004

1 SYNOP 11 39 03590 20050718

2 ACAR 145 2 1RYFVQBA 01092003

2 AMDAR 144 2 EU3781 01072001

4 BATHY 63 39 ZSAF 01032004

4 BUOY 165 1 17546 01122000

5 TEMP 35 2 42314 10062003

5 TEMP 35 3 42339 01121997

5 TEMP 35 1 01001 18052005

6 EUROPROFIL 134 3 ABWWP 01112003

6 PILOT 32 3 07162 01112001

6 PROFILER 34 4 74630 22032004

if (LN_DIAP(jst)%code == \'PROF\') ival=jst

if ( (CLSID(kobs) == LN_DIAP(jst)%cident).and. &

& ((kcodmes == LN_DIAP(jst)%icodmes).or.(icmd == LN_DIAP(jst)%icodmes)) ) then

klndiap( LN_DIAP(jst)%ivarno ) = ival

if ( kotp == NSYNOP ) then !------- SYNOP -------------

if ( LN_DIAP(jst)%ivarno == NVNUMB(10) ) then ! temp 2m

klndiap( NVNUMB( 7) ) = ival ! hum 2m

klndiap( NVNUMB(63) ) = ival ! hum spec.

klndiap( NVNUMB(56) ) = ival ! 10M U COMPONENT

klndiap( NVNUMB(3) ) = ival ! geop.

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(7) ) then ! hum 2m

klndiap( NVNUMB(63) ) = ival

endif

ELSEIF ( kotp == NTEMP ) then !------- NTEMP -------------

if ( LN_DIAP(jst)%ivarno == NVNUMB(3) ) then ! geop.

!klndiap( NVNUMB( 5) ) = ival ! hum

!klndiap( NVNUMB( 8) ) = ival ! temp

!klndiap( NVNUMB( 9) ) = ival ! temp pt de rosee

!klndiap( NVNUMB(58) ) = ival ! hum. de couche

!klndiap( NVNUMB(63) ) = ival ! hum spec.

!klndiap( NVNUMB(10) ) = ival ! temp 2m

!klndiap( NVNUMB( 7) ) = ival ! hum 2m

!klndiap( NVNUMB(63) ) = ival ! hum spec.

!klndiap( NVNUMB(56) ) = ival ! 10M U COMPONENT

!klndiap( NVNUMB(1) ) = ival ! UPPER AIR U COMPONENT

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(8) ) then ! temp

klndiap( NVNUMB( 5) ) = ival

klndiap( NVNUMB(63) ) = ival

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(5) ) then ! hum

klndiap( NVNUMB(63) ) = ival

endif

if ( (CLSID(kobs) == LN_DIAP(jst)%cident).and. &

& ((kcodmes == LN_DIAP(jst)%icodmes).or.(icmd == LN_DIAP(jst)%icodmes)) ) then

klndiap( LN_DIAP(jst)%ivarno ) = ival

if ( kotp == NSYNOP ) then !------- SYNOP -------------

if ( LN_DIAP(jst)%ivarno == NVNUMB(10) ) then ! temp 2m

klndiap( NVNUMB( 7) ) = ival ! hum 2m

klndiap( NVNUMB(63) ) = ival ! hum spec.

klndiap( NVNUMB(56) ) = ival ! 10M U COMPONENT

klndiap( NVNUMB(3) ) = ival ! geop.

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(7) ) then ! hum 2m

klndiap( NVNUMB(63) ) = ival

endif

ELSEIF ( kotp == NTEMP ) then !------- NTEMP -------------

if ( LN_DIAP(jst)%ivarno == NVNUMB(3) ) then ! geop.

klndiap( NVNUMB( 5) ) = ival ! hum

klndiap( NVNUMB( 8) ) = ival ! temp

klndiap( NVNUMB( 9) ) = ival ! temp pt de rosee

klndiap( NVNUMB(58) ) = ival ! hum. de couche

klndiap( NVNUMB(63) ) = ival ! hum spec.

klndiap( NVNUMB(10) ) = ival ! temp 2m

klndiap( NVNUMB( 7) ) = ival ! hum 2m

klndiap( NVNUMB(63) ) = ival ! hum spec.

klndiap( NVNUMB(56) ) = ival ! 10M U COMPONENT

klndiap( NVNUMB(1) ) = ival ! UPPER AIR U COMPONENT

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(8) ) then ! temp

klndiap( NVNUMB( 5) ) = ival

klndiap( NVNUMB(63) ) = ival

ELSEIF ( LN_DIAP(jst)%ivarno == NVNUMB(5) ) then ! hum

klndiap( NVNUMB(63) ) = ival

endif

slide8

The bator procedure

- 2 files are needed for blacklisting: LISTE_NOIRE_DIAP and LISTE_LOC

 Blacklisting of the observations: The LISTE_NOIRE_DIAP file

Obstype

1 SHIP 21 1 62301 01022003

1 SYNOP 14 1 71094 01032004

1 SYNOP 11 39 03590 20050718

2 ACAR 145 2 1RYFVQBA 01092003

2 AMDAR 144 2 EU3781 01072001

4 BATHY 63 39 ZSAF 01032004

4 BUOY 165 1 17546 01122000

5 TEMP 35 2 42314 10062003

5 TEMP 35 3 42339 01121997

6 EUROPROFIL 134 3 ABWWP 01112003

6 PILOT 32 3 07162 01112001

6 PROFILER 34 4 74630 22032004

Character describing the obs subtype

A number describing the obs subtype

A number describing the parameter (varno)

Station identifier

Date: start of theblacklisting

slide9

Action:

N: blacklist

O: force to use

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

slide10

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

obstype:

(Arpege/Aladin)

slide11

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

Obs. subtype

slide12

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

Satellite indicator

slide13

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

Center of production for

Satellite data

slide14

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

Parameter orsensor for Sat data

slide15

!---------------------------------------------------------------------------------------------------!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

max

ZONB4

min

Key word: Blacklisting Satob in different zones

ZONx4min max min max

x=B: < min or > max

or  < min or > max

x=C: < min or > max

or (  > min and < max)

max

min

max

ZONC4

min

max

min

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

slide16

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

Key word: blacklisting of

Satellite channels

TOVSn C1, C2, … Cn

Where 1 n  9

Hirs: 1,20

Amsua:1,15

Airs: 1,2377

Seviri: 1,8

slide17

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PPPP2700 400

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

Key word: blacklisting of pressure levels

PPPPn P1, P2, …Pn

Where 1 n  9

slide18

The bator procedure

Blacklisting of the observations: The LISTE_LOC file

!---------------------------------------------------------------------------------------------------

! action / obstype / [ codetype / ident / code prod. / varno / code / liste ] / ! commentaire

! i1 i2 i3 a8 i3 i3 (a3i2)

x xx xxx xxxxxxxx xxx xxx XXXXx xxxxxxx xxxxxxx (1x,x7)

!---------------------------------------------------------------------------------------------------

N 1 16

N 2 141 29

N 2 144 29

N 2 145 29

N 3 88 052

N 3 88 054

N 3 90 052 ZONB4 -50 50 13 113

N 3 90 054 ZONB4 -50 50 -50 50

N 3 88 253 ZONC4 -50 50 -155 105

N 3 88 254 ZONC4 -50 50 -85 175

N 3 88 256 ZONB4 -50 50 -125 -25

N 6 34 4 PROF2 700 400 0 0 1

N 6 134 3 PROF2 700 400 1 0 1

N 7 210 206 3 TOVS2 6 11

N 9 122 ! ERS

N 9 210 ! NSCAT

N 9 300 ! Quikscat before 02/2000

Key word: blacklisting of pressure thicknesses

PROFn P1, P2, …Pn I1,I2, …In+1

Where 1 n  9 and I1,I2,…,In arethe index of the blacklisting associated to thicknesses

P[1000,P1[,… [Pj-1,Pj[ …

Blacklisting Ij = 1

slide19

To be continued …

Any questions?

slide20

The bator procedure- compilation

  • You can find the bator directory under ~Preproc

 for simplicity we will work in this directory

 bator executable will not use MPI facilities

we had problem with multiprocessor executable

 few empty routines are put in bator.F90

slide21

55 º

55º

30 º

30°

20º

10 º

20°

10°

Exercises

1 (a,b)-Please blacklist Meteosat-8 geowind (AMV) data

a)

b)

2-please blacklist the US profiler data in 2 thicknesses 850 - 600 and 400 – 200 hPa

3- blacklist AMSU-A channels 7, 9,10 and 11 from NOAA-16

4- blacklist only humidity from Budapest (12843) radiosonde(this needs re-compilation of the code)

3- Reduce the obs. errors for the radiosonde temperature data by half (this needs re-compilation of the code)

ad