Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

RANDOMCV IN WRFDA 3.9

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

liaofanlin

New member
Hello!

I have a namelist script that used to work fine with RANDOMCV. However, it could not work anymore. Could you help me to take a quick look at this issue? I am including the name list and the last few lines of the rsl outputs for your information. Thank you very much!

Best,
Liaofan


============
NAMELIST
============
&wrfvar1
print_detail_outerloop=false,
print_detail_grad=true,
/
&wrfvar2
/
&wrfvar3
ob_format=2,
/
&wrfvar4
/
&wrfvar5
/
&wrfvar6
/
&wrfvar7
cv_options = 3,
as1 = 0.25, 1.0, 1.5,
as2 = 0.25, 1.0, 1.5,
as3 = 0.25, 1.0, 1.5,
as4 = 0.25, 1.0, 1.5,
as5 = 0.25, 1.0, 1.5,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
/
&wrfvar11
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
analysis_type = "RANDOMCV",
/
&wrfvar18
analysis_date="2016-07-01_00:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
/
&wrfvar22
/
&wrfvar23
/
&time_control
start_year = 2016,
start_month = 07,
start_day = 01,
start_hour = 00,
start_minute = 00,
start_second = 00,
end_year = 2016,
end_month = 07,
end_day = 01,
end_hour = 06,
end_minute = 00,
end_second = 00,
/
&domains
time_step = 200,
e_we = 121,
e_sn = 151,
e_vert = 41,
p_top_requested = 5000,
dx = 9000,
dy = 9000,
smooth_option = 0,
nproc_x = 0,
/
&fdda
/
&dfi_control
/
&tc
/
&physics
mp_physics = 3,
ra_lw_physics = 1,
ra_sw_physics = 1,
radt = 10,
sf_sfclay_physics = 1,
sf_surface_physics = 2,
bl_pbl_physics = 1,
cu_physics = 1,
cudt = 5,
num_soil_layers = 4,
mp_zero_out = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&scm
/
&dynamics
/
&bdy_control
/
&grib2
/
&namelist_quilt
/
&perturbation
/





===============
rsl.out.0000
===============
----------------------------------------------------------
[3.0] Set up background errors (be) for cv_option: 3
----------------------------------------------------------

DA_Setup_Background_Errors: 3DVAR dry control variables are:
psi, chi_u, t_u and psfc
DA_Setup_Background_Errors: 3DVAR humidity control variable is q/qsg

Number of vertical level for stats = 42
Number of latitude nlath = 96

Number of vertical level for WRFVar= 40

===============
rsl.err.0000
===============
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1 , 1559397164 bytes allocated
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 120 JS 1 JE 150
WRF NUMBER OF TILES = 1
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 61
Error 2 opening gts obs file ob.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 100
Error 2 opening gts obs file ob.ascii
---------------------------------------------------------------
At line 2938 of file da_define_structures.f
Fortran runtime error: End of file
 
Please try modifying line 80 of var/da/da_define_structures/da_random_seed.inc as shown below.
With gnu/7.2.0 compiler, the size of the seed_array is 33, larger than the specified format 10i16.

Code:
diff --git a/var/da/da_define_structures/da_random_seed.inc b/var/da/da_define_structures/da_random_seed.inc
index 3ab61ccf..4d7f7e37 100644
--- a/var/da/da_define_structures/da_random_seed.inc
+++ b/var/da/da_define_structures/da_random_seed.inc
@@ -28,6 +28,8 @@ subroutine da_random_seed
 #endif
 
    call random_seed(size=seed_size)              ! Get size of seed array.
+   write(unit=message(1),fmt='(a,i6)') 'Size of the random_seed array is ', seed_size
+   call da_message(message(1:1))
    allocate(seed_array(1:seed_size))
    seed_array(1:seed_size) = 1
 
@@ -77,7 +79,7 @@ subroutine da_random_seed
    else                                 ! Random seed set "randomly"
       call random_seed
       call random_seed(get=seed_array(1:seed_size))
-      write(unit=message(1),fmt='(a,10i16)') 'Random number seed array = ', seed_array
+      write(unit=message(1),fmt='(a,100i16)') 'Random number seed array = ', seed_array
       call da_message(message(1:1))
    end if
 
Top