La versión para imprimir ya no se admite y puede contener errores de representación. Actualiza los marcadores del navegador y utiliza en su lugar la función de impresión predeterminada del navegador.

Use

These are the steps to use the module

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

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

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