RANDOMCV IN WRFDA 3.9

[liaofanlin]

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

 

[hclin2018]

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.

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

 

[liaofanlin]

Thank you hclin2018 for the quick and correct answer! The problem is resolved.