Get the update of JRA 55 - HiroyukiTsujino/JRA55-do GitHub Wiki

Getting the update of JRA-55 forecast and pre-processing

Getting the update of JRA-55

  • Get the following 3-hourly data obtained from the forecast phase of JRA-55 on a reduced TL319 grid.
    • fcst_surf (Hist/Daily/fcst_surf/YYYYMM/fcst_surf.YYYYMMDDhh)
    • fcst_phy2m (Hist/Daily/fcst_phy2m/YYYYMM/fcst_phy2m.YYYYMMDDhh)
    • ice (Hist/Daily/ice/YYYYMM/ice.YYYYMMDDhh)
    • fcst_phyland (Hist/Daily/fcst_phyland/YYYYMM/fcst_phyland.YYYYMMDDhh).
      • This is used as input for the river routing model. If you do not intend to create river discharge data, fcst_phyland is not necessary.
  • After getting updates of JRA-55, edit the following files in the extract_jra55 directory to specify the path (orgdir) to the obtained data.
    • grib2flatbin_surf.sh
    • grib2flatbin_phy2m.sh
    • grib2flatbin_ice.sh
  • Also get the land-sea mask data from JRA-55 (Const/TL319.grib) if you haven't and put it in the ${DATADIR}/work/const. The land-sea mask can be decoded by the following procedures.
   $ cd ${DATADIR}/work
   $ mkdir const
   $ cd ${HOMEDIR}/extract_jra55
   $ ./grib2flatbin_mask.sh (Make sure that the target directory (orgdir) is correct and give necessary edits)
  • ToDo: Except for assessing the original JRA-55 and COBE-SST data, all scripts should be written with relative paths, using symbolic links from linkdir.

Decode GRIB and make flat binary

  • In extract_jra55 directory, execute
  • $ ./process_year_last.sh 2016 5 29 (When data is available to 29 May 2016. Do not use process_year.sh in the same directory, because it processes unnecessary data)
    • Output ⇒ ../linkdir/work/jra55fcst_3hr_TL319r (little_endian)
  • Note that versions of wgrib may be using distinct grib tables. That is, a particular variable may have been given a different name depending on the version. For example,
    • "total precipitation": "APCP" for wgrib v1.8.1.2a, "TPRAT" for wgrib v1.8.2
    • "runoff": "WATR" for wgrib v1.8.1.2a, "ROF" for wgrib v1.8.2.

Apply zonal 1-2-1 filter to winds

  • In jra55_org_grid_calib directory,
    • compile src/filter_wind_interannual.F90
    • execute $ ./filter_wind_interannual_latest_exec.sh 2020 1 1 2020 7 16 (If data is available to 15 June 2020, specify the day next of this date.)
    • Output ⇒ ../linkdir/work/jra55fcst_filt_3hr_TL319r

Regridding from reduced TL319 to normal TL319

  • In jra55_org_grid_product directory, execute

    • '$ ./red2reg_exec_raw_v1_0.sh 2018 2018 1 32` (except for wind vectors)
      • Note: Process 1 year by one operation. Specify the first day (1) and the last day (32 for 01Feb) of the year to which data is available.
    • Output ⇒ ../linkdir/work/jra55fcst_filt_3hr_TL319
    • `$ ./red2reg_exec_raw_v1_0_wind.sh 2018 2018 1 32' (wind vectors)
      • Note: Process 1 year by one operation. Specify the first day (1) and the last day (32 for 01Feb) of the year to which data is available.
    • Output ⇒ ../linkdir/work/jra55fcst_filt_3hr_TL319

  • In ${DATADIR}/work, make symbolic link from jra55fcst_3hr_TL319 to jra55fcst_filt_3hr_TL319

⚠️ **GitHub.com Fallback** ⚠️