var/convertor/kma_netcdf/netcdf2kma.f90

   1 !NETCDF2KMA:DRIVER_LAYER:MAIN
   2 !
   3
   4 PROGRAM netcdf2kma
   5
   6    use da_control, only : debug_level
   7
   8    USE module_machine
   9    USE module_domain
  10    USE module_driver_constants
  11    USE module_configure
  12
  13    USE module_timing
  14    USE module_wrf_error
  15 #ifdef DM_PARALLEL
  16    USE module_dm
  17 #endif
  18
  19    USE da_wrfvar_io
  20    USE module_netcdf2kma_interface
  21    USE module_wave2grid_kma
  22
  23    IMPLICIT NONE
  24
  25    REAL    :: time
  26
  27    TYPE (domain) , POINTER :: keep_grid, grid_ptr, null_domain
  28    TYPE (grid_config_rec_type)              :: config_flags
  29
  30    INTEGER :: domain_id , fid , oid , idum1 , idum2
  31
  32 !--------------
  33  INTERFACE
  34 SUBROUTINE Setup_Timekeeping ( grid )
  35    USE module_domain
  36    TYPE(domain), POINTER :: grid
  37 END SUBROUTINE Setup_Timekeeping
  38 END INTERFACE
  39 !--------------
  40
  41 #ifdef DM_PARALLEL
  42    INTEGER                 :: nbytes
  43    INTEGER, PARAMETER      :: configbuflen = 4*8192
  44    INTEGER                 :: configbuf( configbuflen )
  45    LOGICAL , EXTERNAL      :: wrf_dm_on_monitor
  46 #endif
  47
  48 !--Define the name of this program (program_name defined in module_domain)
  49
  50    program_name = "NETCDF2KMA"
  51
  52 !--Get the NAMELIST data for input.
  53
  54    CALL init_modules(1)
  55 !rizvi's add  start
  56      ! Initialize utilities (time manager, etc.)
  57       call wrfu_initialize(defaultCalendar=WRFU_CAL_GREGORIAN)
  58 !rizvi's add over
  59
  60 #ifdef DM_PARALLEL
  61    IF ( wrf_dm_on_monitor() ) THEN
  62      CALL start_timing
  63      CALL initial_config
  64    ENDIF
  65    CALL get_config_as_buffer( configbuf, configbuflen, nbytes )
  66    CALL wrf_dm_bcast_bytes( configbuf, nbytes )
  67    CALL set_config_as_buffer( configbuf, configbuflen )
  68    CALL wrf_dm_initialize
  69 #else
  70    CALL start_timing
  71    CALL initial_config
  72 #endif
  73
  74 !shc-wei start
  75 !  CALL get_debug_level ( debug_level )
  76    CALL nl_get_debug_level ( 1,debug_level )
  77 !shc-wei end
  78    CALL set_wrf_debug_level ( debug_level )
  79
  80 !--allocated and configure the mother domain
  81
  82    NULLIFY( null_domain )
  83
  84    CALL       wrf_debug ( 100 , 'wrf: calling alloc_and_configure_domain ' )
  85    CALL alloc_and_configure_domain ( domain_id  = 1 ,                  &
  86 !                                    local_time = 0 ,                  &     !shc-wei
  87                                      grid       = head_grid ,          &
  88                                      parent     = null_domain ,        &
  89                                      kid        = -1                   )
  90
  91    CALL       wrf_debug ( 100 , 'wrf: calling model_to_grid_config_rec ' )
  92    CALL model_to_grid_config_rec ( head_grid%id , model_config_rec , config_flags )
  93    CALL       wrf_debug ( 100 , 'wrf: calling set_scalar_indices_from_config ' )
  94    CALL set_scalar_indices_from_config ( head_grid%id , idum1, idum2 )
  95    CALL       wrf_debug ( 100 , 'wrf: calling init_wrfio' )
  96    CALL Setup_Timekeeping ( head_grid )
  97    if(config_flags%real_data_init_type == 1 .or. &
  98       config_flags%real_data_init_type == 3) then
  99      CALL init_wrfio
 100
 101 #ifdef DM_PARALLEL
 102      CALL get_config_as_buffer( configbuf, configbuflen, nbytes )
 103      CALL wrf_dm_bcast_bytes( configbuf, nbytes )
 104      CALL set_config_as_buffer( configbuf, configbuflen )
 105 #endif
 106
 107    endif
 108    call med_add_config_info_to_grid ( head_grid )
 109
 110
 111 !-- read input in netcdf format
 112      CALL da_med_initialdata_input( head_grid , config_flags, 'rizvi_to_choose' )
 113
 114 !--write in grid-array in kma-ascii format
 115
 116    CALL wrf_debug ( 100 , 'da_3dvar: calling netcdf2kma_interface ' )
 117    CALL netcdf2kma_interface ( head_grid, config_flags )
 118    CALL wrf_debug ( 100 , 'da_3dvar: back from netcdf2kma_interface' )
 119
 120    CALL med_shutdown_io ( head_grid , config_flags )
 121
 122 #ifdef DM_PARALLEL
 123    IF ( wrf_dm_on_monitor() ) THEN
 124       message(1) = 'Time elapsed in netcdf2kma'
 125       CALL end_timing ( TRIM(message(1)) )
 126    END IF
 127 #else
 128    message(1) = 'Time elapsed in netcdf2kma'
 129    CALL end_timing ( TRIM(message(1)) )
 130 #endif
 131
 132    CALL wrf_shutdown
 133
 134    STOP "Stopped normally. in netcdf2kma"
 135
 136 END PROGRAM netcdf2kma
 137