CDXWRFuse - Historial de revisiones

Lluis.fita: /* Use */

2025-07-30T16:00:43Z

Use

← Revisión anterior		Revisión del 16:00 30 jul 2025
Línea 32:		Línea 32:
	= 2: wsgs folllowing heavy precipitation method		= 2: wsgs folllowing heavy precipitation method
	output_wb = 1: whether water-budget variables (Fita and Flaounas, 2018, QJRMS) have to		output_wb = 1: whether water-budget variables (Fita and Flaounas, 2018, QJRMS) have to
	computed (1) or not [0, default]		computed (1) or not [0, default] (only available with CDXWRF=4)
	wsz100_diag = 1: wind extraoplation at z100m_wind using power-law method [default]		wsz100_diag = 1: wind extraoplation at z100m_wind using power-law method [default]
	= 2: wind extraoplation at z100m_wind using logarithmic method		= 2: wind extraoplation at z100m_wind using logarithmic method

Lluis.fita: /* Optional values for wbds & wss residence-time */

2025-07-30T15:59:21Z

Optional values for wbds & wss residence-time

← Revisión anterior		Revisión del 15:59 30 jul 2025
Línea 95:		Línea 95:
	In the code, user is provided with 2 different complementary set-ups for the wbds & wss residence-time		In the code, user is provided with 2 different complementary set-ups for the wbds & wss residence-time

	[[kt-based]]		<ins>kt-based</ins>
	The references for wind-speed are based on knots (kt)		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. ]		* wss: [0., 2.5, 5., 10., 15., 20., 25., 30., 35., 40.. 45., 50., 60., 70., 80., 100., 150., 200., 250., 300. ]

Lluis.fita: /* Use */

2025-07-30T15:57:21Z

Use

@@ Línea 2: / Línea 2: @@
 These are the steps to use the module
-# One need to add to the <CODE>'namelist.input'</CODE> the auxiliar output number 9 (e.g. for output every 3 hours and 1-day files) at the <CODE>`&history'</CODE> section:
+# One need to add to the <CODE>'namelist.input'</CODE> 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 <CODE>`&history'</CODE> section:
 <PRE>
-auxhist9_outname = "wrfcdx_d<domain>_<date>"
+auxhist19_outname = "wrfcdx_d<domain>_<date>"
-auxhist9_interval = 180, 180,
+auxhist19_interval = 180, 180,
-frames_per_auxhist9 = 8, 8,
+frames_per_auxhist19 = 8, 8,
-io_form_auxhist9 = 2
+io_form_auxhist19 = 2
 </PRE>
 # 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.
 <PRE>
@@ Línea 41: / Línea 50: @@
   fogvars            = 1: variables to use to diagnose fog using 3D [hur] [default]
                      = 2: sfc [hus] (not available for fogvisibility_diag = 1)
+/
 </PRE>
@@ Línea 55: / Línea 131: @@
 (...)
+/
 &diags
@@ Línea 66: / Línea 140: @@
+/
 </PRE>
 Back to the main page [[CDXWRF]]

Lluis.fita en 16:48 27 feb 2019

2019-02-27T16:48:16Z

← Revisión anterior		Revisión del 16:48 27 feb 2019
Línea 66:		Línea 66:
	/		/
	</PRE>		</PRE>

			Back to the main page [[CDXWRF]]

Lluis.fita: Creando la página de uso

2019-02-27T15:18:56Z

Creando la página de uso

Página nueva

= Use =

These are the steps to use the module
# One need to add to the <CODE>'namelist.input'</CODE> the auxiliar output number 9 (e.g. for output every 3 hours and 1-day files) at the <CODE>`&history'</CODE> section:
<PRE>
auxhist9_outname = "wrfcdx_d<domain>_<date>"
auxhist9_interval = 180, 180,
frames_per_auxhist9 = 8, 8,
io_form_auxhist9 = 2
</PRE>
# 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.
<PRE>
&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]
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)
/
</PRE>

== Pressure interpolation ==
Remember to activate section <CODE>&diags</CODE> in order to get pressure-level vertical interpolation of state variables (g.e.: assuming 6 levels only and output every 3 hours)
<PRE>
&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.
/
</PRE>