MetClim

WRFV3.9.1 with ERA-Interim data on ORR2

Before trying to run WRF with ERA-Interim data, make sure you have completed the following:

    1. WPS: Use this as a template: namelist.wps, then run geogrid.exe to generate geo_em.d01.nc
    2. WPS: Copy the Vtable for ERA-Interim data into the WPS directory 
      cp ungrib/Variable_Tables/Vtable.ERA-interim.ml Vtable

      Now link to your data files

      ./link_grib.csh ../../GRIB/ERAI/ERA

      and run ungrib.exe:

      [conor@orr2 WPS]$ ls -hl FILE*
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-01_12
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-01_18
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-02_00
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-02_06
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-02_12
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-02_18
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-03_00
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-03_06
-rw-rw-r-- 1 conor conor 115M Mar 2 18:01 FILE:2017-12-03_12
    3. We need to convert the ERA-Interim data from model level to pressure levels. To do this, we need to create a file with information about the ERA-Interim model levels: (still in the WPS directory) create a file called ecmwf_coeffs file and put in the following data: 
          0         0.000000    0.00000000
    1        20.000000    0.00000000
    2        38.425343    0.00000000
    3        63.647804    0.00000000
    4        95.636963    0.00000000
    5       134.483307    0.00000000
    6       180.584351    0.00000000
    7       234.779053    0.00000000
    8       298.495789    0.00000000
    9       373.971924    0.00000000
   10       464.618134    0.00000000
   11       575.651001    0.00000000
   12       713.218079    0.00000000
   13       883.660522    0.00000000
   14      1094.834717    0.00000000
   15      1356.474609    0.00000000
   16      1680.640259    0.00000000
   17      2082.273926    0.00000000
   18      2579.888672    0.00000000
   19      3196.421631    0.00000000
   20      3960.291504    0.00000000
   21      4906.708496    0.00000000
   22      6018.019531    0.00000000
   23      7306.631348    0.00000000
   24      8765.053711    0.00007582
   25     10376.126953    0.00046139
   26     12077.446289    0.00181516
   27     13775.325195    0.00508112
   28     15379.805664    0.01114291
   29     16819.474609    0.02067788
   30     18045.183594    0.03412116
   31     19027.695313    0.05169041
   32     19755.109375    0.07353383
   33     20222.205078    0.09967469
   34     20429.863281    0.13002251
   35     20384.480469    0.16438432
   36     20097.402344    0.20247594
   37     19584.330078    0.24393314
   38     18864.750000    0.28832296
   39     17961.357422    0.33515489
   40     16899.468750    0.38389215
   41     15706.447266    0.43396294
   42     14411.124023    0.48477158
   43     13043.218750    0.53570992
   44     11632.758789    0.58616841
   45     10209.500977    0.63554746
   46      8802.356445    0.68326861
   47      7438.803223    0.72878581
   48      6144.314941    0.77159661
   49      4941.778320    0.81125343
   50      3850.913330    0.84737492
   51      2887.696533    0.87965691
   52      2063.779785    0.90788388
   53      1385.912598    0.93194032
   54       855.361755    0.95182151
   55       467.333588    0.96764523
   56       210.393890    0.97966272
   57        65.889244    0.98827010
   58         7.367743    0.99401945
   59         0.000000    0.99763012
   60         0.000000    1.00000000

      Then copy the WRF conversion program into the WPS directory:

      cp ./util/calc_ecmwf_p.exe .

      and run it:

      ./calc_ecmwf_p.exe

      This should produce the files:

      [conor@orr2 WPS]$ ls -hl PRES*
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-01_12
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-01_18
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-02_00
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-02_06
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-02_12
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-02_18
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-03_00
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-03_06
-rw-rw-r-- 1 conor conor 54M Mar 2 18:01 PRES:2017-12-03_12
    4. Run metgrid.exe, and you should see the following files: 
      [conor@orr2 WPS]$ ls -hl met_em*
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-01_12:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-01_18:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-02_00:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-02_06:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-02_12:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-02_18:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-03_00:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-03_06:00:00.nc
-rw-rw-r-- 1 conor conor 13M Mar 2 18:01 met_em.d01.2017-12-03_12:00:00.nc
    5.  Change into your WRFV3/run directory and use this as a template: namelist.input. Copy or link the met_em files to your run directory, then run real.exe. You should get files like this: 
      [conor@orr2 run]$ ls -hl wrfinput_d01 wrfbdy_d01 
-rw-rw-r-- 1 conor conor  20M Mar  2 18:02 wrfbdy_d01
-rw-rw-r-- 1 conor conor 8.6M Mar  2 18:02 wrfinput_d01
    6. Now run wrf: qsub RUN_WRF.qsub and you forecast should produce output files like this: 
      [conor@orr2 run]$ ls -hl wrfout*
-rw-rw-r-- 1 conor conor 395M Mar  2 18:05 wrfout_d01_2017-12-01_12:00:00

    Done!

Posted

in

by

Conor Sweeney

Tags: