1 / 21

Observation pre-processing in the ARPEGE/ALADIN model

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 .

zev
Download Presentation

Observation pre-processing in the ARPEGE/ALADIN model

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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Observation pre-processing in theARPEGE/ALADIN model Roger Randriamampianina Hungarian Meteorological Service ALADIN 3D-Var/ODB training course, Budapest, 6-10 June, 2006

  2. Outline Local pre-processing: OULAN BATOR 1C RAD. BIAS CORRECTION

  3. BATOR

  4. 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)

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

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

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

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

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

  10. 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)

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

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

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

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

  15. !---------------------------------------------------------------------------------------------------!--------------------------------------------------------------------------------------------------- ! 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

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

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

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

  19. To be continued … Any questions?

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

  21. 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)

More Related