地波雷达质控与同化模块
################################

.. toctree::
   :titlesonly:

沈阳地波雷达资料转化程序
----------------------

将地波雷达资料处理为little_r格式，按10m风观测同化进入系统即可：


.. code-block:: python

   #!/usr/bin/env python
   
   #----------------------------------------------------------------------------
   
   if __name__ == '__main__':    
       import sys
       import datetime
       import argparse
       import fortranformat as ff
       parser = argparse.ArgumentParser(description='Example of usage')
       parser.add_argument('-d', '--date',
                           dest='start_date',
                           default=None,
                           required=True,
                           help='date in format "YYMMDDHHmm"')
       args = parser.parse_args()
       start_date_object = args.start_date
       filename = start_date_object + '.000'
       
       datetime = '20'+start_date_object
       # Constants
       rpt_format =  ff.FortranRecordWriter ('( 2f20.5, 2a40, 2a40, 1f20.5, 5i10, 3L10, 2i10, a20, 13( f13.5 , i7 ) ) ')
       mesa_format = ff.FortranRecordWriter ('( 10( f13.5 , i7 ) )')
       end_format = ff.FortranRecordWriter (' ( 3 ( i7 ) ) ')
   
       lr_missing  = -888888
       lr_end_data = -777777.0
       lr_missing_r = lr_missing
   
       platform  = "FM-12 SYNOP                             "
   
       # Counter
       counter = 1
   
       print("Reading " + filename)
       with open(filename) as lndradar:
           for row in lndradar:
               [id, lon, lat, wdir, wspd] = map(float, row.split())
   
               print(rpt_format.write([lat           , lon           , ""               , ""          ,
                           platform         , 'ShenYang Land Radar'       , 0.          , 6           ,
                           0                , 0                , counter          , 0           ,
                           False            , False            , False            , lr_missing  ,
                           lr_missing       , datetime         , lr_missing_r     , 0           ,
                           lr_missing_r     , 0                , lr_missing_r     , 0           ,
                           lr_missing_r     , 0                , lr_missing_r     , 0           ,
                           lr_missing_r     , 0                , lr_missing_r     , 0           ,
                           lr_missing_r     , 0                , lr_missing_r     , 0           ,
                           lr_missing_r     , 0                , lr_missing_r     , 0           ,
                           lr_missing_r     , 0                , lr_missing_r     , 0]))
               print(mesa_format.write([lr_missing_r        , 0 , 0.                  , 0 , lr_missing_r        , 0 ,
                                        lr_missing_r        , 0 , wspd                , 0 , wdir                , 0 ,
                                        lr_missing_r        , 0 , lr_missing_r        , 0 , lr_missing_r        , 0 , lr_missing_r , 0
                                       ]))
               print(mesa_format.write([lr_end_data , 0 , lr_end_data , 0 , 1., 0 , 
                                       lr_missing_r , 0 , lr_missing_r , 0 , lr_missing_r , 0 ,
                                       lr_missing_r , 0 , lr_missing_r , 0 , lr_missing_r , 0 ,
                                       lr_missing_r , 0]))
               print(end_format.write([1, 0, 0]))
               counter = counter + 1