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

gen be for ensemble-based perturbations

Do you mean attached to one of your messages? I cannot see it.
By the way I compiled gen_be_v3, now if I'm right the namelist for cv option 5 with perturbation from an ensemble should be:

&gen_be_nml
nc_list_file = 'flist.txt' ! flist.txt contains a list of wrfout files in full path
be_method = 'ENS'
varnames = 'PSI', 'CHI_U', 'T_U', 'RH', 'PSFC_U'
do_pert_calc = .true.
do_slen_calc = .true.
slen_opt = 2 ! 1: curve-fitting method (extremely slow, need to run with OPENMP on)
! 2: laplacian method (recommended)
do_eof_transform = .true. ! .true.: statistics on EOF space for WRFDA runs
! .false.: statistics on physical model levels for diagnostics
pert1_read_opt = 1 ! specify how the program accesses the pert1 data (output of do_pert_calc=.true.)
! 1: (default) read and store all cases in memory at once
! 2: read from pert1 file when need it
! This option writes out additional vertical-mode-projected fields when do_eof_transform=.true.
! Set pert1_read_opt=2 for large number of cases and large domain sizes that encounter memory insufficiency.
! Note: when unbalanced variables are requested in varnames, pert1_read_opt=2 is used regardless of the namelist setting.
/
&ens_nml
/
 
Do you mean attached to one of your messages? I cannot see it.
By the way I compiled gen_be_v3, now if I'm right the namelist for cv option 5 with perturbation from an ensemble should be:

&gen_be_nml
nc_list_file = 'flist.txt' ! flist.txt contains a list of wrfout files in full path
be_method = 'ENS'
varnames = 'PSI', 'CHI_U', 'T_U', 'RH', 'PSFC_U'
do_pert_calc = .true.
do_slen_calc = .true.
slen_opt = 2 ! 1: curve-fitting method (extremely slow, need to run with OPENMP on)
! 2: laplacian method (recommended)
do_eof_transform = .true. ! .true.: statistics on EOF space for WRFDA runs
! .false.: statistics on physical model levels for diagnostics
pert1_read_opt = 1 ! specify how the program accesses the pert1 data (output of do_pert_calc=.true.)
! 1: (default) read and store all cases in memory at once
! 2: read from pert1 file when need it
! This option writes out additional vertical-mode-projected fields when do_eof_transform=.true.
! Set pert1_read_opt=2 for large number of cases and large domain sizes that encounter memory insufficiency.
! Note: when unbalanced variables are requested in varnames, pert1_read_opt=2 is used regardless of the namelist setting.
/
&ens_nml
/
looks good, maybe exchange further by email: stephen.tjemkes@maatresearch.eu. Then i can also send the python module and exchange some experience if of interest
 
Top