CDXWRFuse

Use

These are the steps to use the module

1. One need to add to the 'namelist.input' the auxiliar output number 19 [NOTE: this has changed from previous version !!] (e.g. for output every 3 hours and 1-day files) at the `&history' section:

auxhist19_outname = "wrfcdx_d<domain>_<date>"
auxhist19_interval = 180, 180,
frames_per_auxhist19 = 8, 8,
io_form_auxhist19 = 2

1. One need to add to the 'namelist.input' the auxiliar output number 18 [NOTE: this has changed from previous version !!] for high-frequency output [by now, only COLMAX and PRFLUX] (e.g. for output every 10 minutes and 1-day files) at the `&history' section:

auxhist18_outname = "wrfhfcdx_d<domain>_<date>"
auxhist18_interval = 10, 10,
frames_per_auxhist18 = 144, 144,
io_form_auxhist18 = 2

1. Also a new section should be added (assuming it will get complex and different implementations of the diagnostics might be necessary...) NOTE: do not copy directly this text. Fortran do not understand ':' as namelist input.

&cordex
 output_cordex      = 1
 psl_diag           = 1: sea-level pressure diagnostic following hydrostatic Shuell correction 
                      (Stackpole and Cooley, 1970)
                    = 2: psl diagnostic following a target pressure (Benjamin and Miller, 1990)
                    = 3: psl diagnostic following ECMWF method (Yesad, 2015) [default]
 psmooth            = 5: passes of neighborgh filtering (3x3-grid point mean) of psfc for psl_diag=2 
                      [default 5]
 ptarget            = 70000.: pressure [Pa] target to be used by psl_diag=2 [default 70000.]
 wsgs_diag          = 1: wind-gust diagnostic following (Brasseur, 2001, MWR) [default]
                    = 2: wsgs folllowing heavy precipitation method
 output_wb          = 1: whether water-budget variables (Fita and Flaounas, 2018, QJRMS) have to 
                      computed (1) or not [0, default] (only available with CDXWRF=4)
 wsz100_diag        = 1: wind extraoplation at z100m_wind using power-law method [default]
                    = 2: wind extraoplation at z100m_wind using logarithmic method
                    = 3: wind extraoplation at z100m_wind using Monin-Obukhov theory
 z100m_wind         = 100.: height to extraplate winds [100. default]
 zmlagen_dqv        = 0.1: percentage of variation of mixing ratio to determine mixed layer depth used in 
                        zmla computation (Nielsen-Gammon et al., 2008) [0.1 default]
 zmlagen_dtheta     = 1.5: increment in K of potantial temperature from its minimum within the MLD used in 
                        zmla computation (1.5 default)
 potevap_diag       = 1 potential evapotranspiration using bulk computation (Manabe, 1969)
 potevap_diag       = 2 potential evapotranspiration using Milly92 correction (Milly, 1992) after ORCHIDEE 
                        implementation [default]
 convxtrm_diag      = 0: diagnostic of extremes from convection indices: 0: No [default]; 1: yes
 fogvisibility_diag = 1: diagnostic of visibility inside fog following Kunkel (1984)
                    = 2: RUC method (Smirnova et al., 2000)
                    = 3: FRAML 50% prob Gultepe and Milbrandt, (2010) [default]
 fogvars            = 1: variables to use to diagnose fog using 3D [hur] [default]
                    = 2: sfc [hus] (not available for fogvisibility_diag = 1)

 zmlagen_diag       = 1: diagnostic of generic zmla after Nielsen-Gammon et al. (2008)
                    = 2: bulk-Richardson (Vogelezang and Holtslag, 1996; Seidel et al., 2004) (default)
 timedbg            = frequency in minutes of time-debugging (720, default)
 londbg             = longitude in degrees of grid-point debugging (-999., not used) (-999., default)
 latdbg             = latitude in degrees of grid-point debugging (-999., not used) (-999., default)
 wgsts_wbdswsstres  = whether wind gusts should be used for wbds-wss residence-time (0: no, 1: yes (1, default)

#### ONLY for CDXWRF>=3
 nhtasrng           = amount + 1 of high-resolution bins for temporal residence of tas (38, default)
 nhhursrng          = amount + 1 of high-resolution bins for temporal residence of hurs (25, default)
 nltasrng           = amount + 1 of low-resolution bins for temporal residence of tas (21, default)
 nlhursrng          = amount + 1 of low-resolution bins for temporal residence of hurs (25, default)
 nwbdsrng           = amount + 1 of  bins for temporal residence of wbds (33, default)
 nwssrng            = amount + 1 of  bins for temporal residence of wss (22, default)
 nhtashursrng       = amount + 1 of high-resolution bins for temporal residence of tas and hurs (nhtas+1)*(nhhurs+1) (950, default)
 nltashursrng       = amount + 1 of low-resolution bins for temporal residence of wss and wbds (nwbds+1)*(nwss+1) (525, default)
 nwbdswssrng        = amount + 1 of bins for temporal residence of wbds and wss (nwbds+1)*(nwss+1) (726, default)
 htasrng            = high-resolution bins for temporal residence of tas (/173.15,268.15,269.15,270.15,271.15,272.15,273.15,
    274.15,275.15,276.15,277.15,278.15,279.15,280.15,281.15,282.15,283.15,284.15,285.15,286.15,287.15,288.15,289.15,
    290.15,291.15,292.15,293.15,294.15,295.15,296.15,297.15,298.15,299.15,300.15,301.15,302.15,303.15,0.00,0.00,0.00,
    0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00/)
 hhursrng           = high-resolution bins for temporal residence of hurs (/0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,
    0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0.85,0.90,0.95,1.00,1.05,1.10,1.15,0.00,0.00,0.00,0.00,0.00,0.00,0.00,
    0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00/)
 ltasrng            = low-resolution bins for temporal residence of tas (/173.15,233.15,238.15,243.15,248.15,253.15,
    258.15,263.15,268.15,273.15,278.15,283.15,288.15,293.15,298.15,303.15,308.15,313.15,318.15,323.15,0.00,0.00,0.00,
    0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,
    0.00,0.00,0.00,0.00,0.00/)
 lhursrng           = low-resolution bins for temporal residence of hurs (/0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,
    0.50,0.55,0.60,0.65,0.70,0.75,0.80,0.85,0.90,0.95,1.00,1.05,1.10,1.15,0.00,0.00,0.00,0.00,0.00,0.00,
    0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00/)
 wbdsrng            = bins for temporal residence of wbds (/-0.0982,0.0982,0.2945,0.4909,0.6872,0.8836,   
    1.0799,1.2763,1.4726,1.669,1.8653,2.0617,2.258,2.4544,2.6507,2.8471,3.0434,3.2398,3.4361,3.6325,3.8288,
    4.0252,4.2215,4.4179,4.6142,4.8106,5.0069,5.2033,5.3996,5.596,5.7923,5.9887,6.185,0.0,0.0,0.0,0.0,0.0,
    0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/)
 wssrng             = bins for temporal residence of wss (/0.0,0.5,1.0,2.5,5.0,10.0,15.0,20.0,25.0,30.0,35.0,40.0,45.0,
    50.0,60.0,70.0,80.0,90.0,100.0,125.0,150.0,200.0,95.2375,107.6506,120.5606,133.9489,147.7992,0.0,0.0,0.0,
    0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/)
/

Optional values for wbds & wss residence-time

In the code, user is provided with 2 different complementary set-ups for the wbds & wss residence-time

kt-based The references for wind-speed are based on knots (kt)

• wss: [0., 2.5, 5., 10., 15., 20., 25., 30., 35., 40.. 45., 50., 60., 70., 80., 100., 150., 200., 250., 300. ]
• Transforming to ms-1 as wss = kt * 1852. / 3600.

 nwssrng            = amount + 1 of  bins for temporal residence of wss kt-based (20, default)
 nwbdswssrng        = amount + 1 of bins for temporal residence of wbds and wss (nwbds+1)*(nwss+1) kt-based (660, default)
 wssrng             = bins for temporal residence of wss kt based (/0.0,1.2861,2.5722,5.1444,7.7167,10.2889,12.8611,
   15.4333,18.0056,20.5778,23.15,25.7222,30.8667,36.0111,41.1556,51.4444,77.1667,102.8889,128.6111,154.3333,0.0,0.0,0.0,0.0,
   0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0/)        

Beafourt scale-based The references of wind speed are based on the Beafourt scale

• wss: [0., 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5, 12.5, 13.5, 14.5, 15.5, 16.5, 17.5, 18.5, 19.5, 21.5, 23.5, 25.5, 27.5, 29.5, 31.5]
• Transforming to ms-1 as wss = 0.836*Bnum**(3./2.)

 nwssrng            = amount + 1 of  bins for temporal residence of wss Beafourt-based (27, default)
 nwbdswssrng        = amount + 1 of bins for temporal residence of wbds and wss (nwbds+1)*(nwss+1) Beafourt-based (891, default)
 wssrng             = bins for temporal residence of wss Beafourt-based (/0.0,0.2956,1.5358,3.3046,5.474,7.9804,10.7833,
   13.854,17.1711,20.7174,24.4789,28.444,32.6027,36.9463,41.4674,46.1592,51.0157,56.0315,61.2017,66.5218,71.9877,83.342,
   95.2375,107.6506,120.5606,133.9489,147.7992,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
   0.0,0.0,0.0,0.0,0.0/)

Pressure interpolation

Remember to activate section &diags in order to get pressure-level vertical interpolation of state variables (g.e.: assuming 6 levels only and output every 3 hours)

&time_control
(...)
auxhist23_outname="wrfpress_d<domain>_<date>"
io_form_auxhist23 = 2,
auxhist23_interval = 180, 180,
frames_per_auxhist23 = 100, 100,
(...)
/

&diags
p_lev_diags = 1,
num_press_levels = 6,
press_levels = 100000, 92500, 85000, 70000, 50000, 20000
use_tot_or_hyd_p = 1
p_lev_missing = -999.
/

• Optional for 23 p-level values:

 num_press_levels = 23,
 press_levels = 100000,97500,95000,92500,90000,87500,85000,82500,80000,75000,70000,65000,60000,55000,50000,45000,40000,35000,
   30000,25000,20000,15000,10000,

The module will provide additional mass-conservative vertically integrated values for the water species

Back to the main page CDXWRF