Observation pre-processing Observation pre-processing in in thethe ARPEGE/ALADIN modelARPEGE/ALADIN model
Roger RandriamampianinaRoger Randriamampianina Hungarian Meteorological ServiceHungarian Meteorological Service
ALADIN 3D-Var/ODB training course, Budapest, 6-10 June, 2006ALADIN 3D-Var/ODB training course, Budapest, 6-10 June, 2006
Local pre-processing:
OULAN
BATOR
1C RAD. BIAS CORRECTION
OutlineOutline
BATOR
The bator procedureThe bator procedure
2 – Blacklisting of the observations2 – Blacklisting of the observations
1 – 1 – CreationCreation of (an) ODB database from the OBSOUL file(s) of (an) ODB database from the OBSOUL file(s)
1 – Initialisation of the observation errors1 – Initialisation of the observation errors It contains the observation error for several observation types It contains the observation error for several observation types
Bator and the observation errorsBator and the observation errorsECTERO(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
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
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
bator_saisies.F90
Bator and blacklistingBator and blacklisting
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
bator_saisies.F901 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
- 2 files are needed for blacklisting: LISTE_NOIRE_DIAP and LISTE_LOC LISTE_LOC Blacklisting of the observations: The LISTE_NOIRE_DIAP file
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
Obstype
Character describing the obs subtype
A number describing the obs subtype
A number describing the parameter (varno)
Station identifier
Date: start of theblacklisting
The bator procedureThe bator procedure
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
The bator procedureThe bator procedureBlacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Action:
N: blacklistO: force to use
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
The bator procedureThe bator procedureBlacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
obstype: (Arpege/Aladin)
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Blacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Obs. subtype
The bator procedureThe bator procedure
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Blacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Satellite indicator
The bator procedureThe bator procedure
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Blacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Center of production forSatellite data
The bator procedure
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Blacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Parameter orsensor for Sat data
The bator procedureThe bator procedure
Blacklisting of the observations: The LISTE_LOC fileBlacklisting 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Key word: Blacklisting Satob in different zonesZONx4 min max min maxx=B: < min or > max or < min or > maxx=C: < min or > maxor ( > min and < max)
The bator procedureThe bator procedure
ZONB4
min max
max
min
minmax
max
min
ZONC4
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Blacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Key word: blacklisting of Satellite channelsTOVSn C1, C2, … CnWhere 1 n 9Hirs: 1,20Amsua:1,15Airs: 1,2377Seviri: 1,8
The bator procedureThe bator procedure
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PPPP2 700 400 N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Blacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Key word: blacklisting of pressure levelsPPPPn P1, P2, …PnWhere 1 n 9
The bator procedureThe bator procedure
!---------------------------------------------------------------------------------------------------! 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 16N 2 141 29N 2 144 29N 2 145 29N 3 88 052N 3 88 054N 3 90 052 ZONB4 -50 50 13 113N 3 90 054 ZONB4 -50 50 -50 50N 3 88 253 ZONC4 -50 50 -155 105N 3 88 254 ZONC4 -50 50 -85 175N 3 88 256 ZONB4 -50 50 -125 -25N 6 34 4 PROF2 700 400 0 0 1N 6 134 3 PROF2 700 400 1 0 1N 7 210 206 3 TOVS2 6 11N 9 122 ! ERSN 9 210 ! NSCATN 9 300 ! Quikscat before 02/2000
Blacklisting of the observations: The LISTE_LOC fileBlacklisting of the observations: The LISTE_LOC file
Key word: blacklisting of pressure thicknessesPROFn P1, P2, …Pn I1,I2, …In+1Where 1 n 9 and I1,I2,…,In arethe index of the blacklisting associated to thicknesses P[1000,P1[,… [Pj-1,Pj[ …Blacklisting Ij = 1
The bator procedureThe bator procedure
To be continued …
Any questions?
The bator procedure- compilationThe bator procedure- compilation
You can find the bator directory under ~PreprocYou can find the bator directory under ~Preproc for simplicity we will work in this directoryfor simplicity we will work in this directory
bator executable will not use MPI facilitiesbator executable will not use MPI facilitieswe had problem with multiprocessor executablewe had problem with multiprocessor executable
few empty routines are put in bator.F90few empty routines are put in bator.F90
Exercises Exercises 1 (a,b)-1 (a,b)- Please blacklist Meteosat-8 geowind (AMV) data Please blacklist Meteosat-8 geowind (AMV) data
2-2- please blacklist the US profiler data in 2 thicknesses please blacklist the US profiler data in 2 thicknesses 850 - 600 and 400 – 200 hPa850 - 600 and 400 – 200 hPa
3- 3- blacklist AMSU-A channels 7, 9,10 and 11 from NOAA-16blacklist AMSU-A channels 7, 9,10 and 11 from NOAA-16
10° 20°
55º
30°
44- - blacklist only humidity from Budapest (12843) radiosondeblacklist only humidity from Budapest (12843) radiosonde((this this needneedss re-compilation of the code) re-compilation of the code)
3- 3- Reduce the obs. errors for the radiosonde temperatureReduce the obs. errors for the radiosonde temperature data by half ( data by half (this this needneedss re-compilation of the code) re-compilation of the code)
10 º 20º
55 º
30 º
a) b)