Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

updates for running C90C_ens and C90C_replay test cases #93

Merged
merged 1 commit into from
May 26, 2021
Merged

updates for running C90C_ens and C90C_replay test cases #93

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
17 changes: 8 additions & 9 deletions src/Applications/GEOSdas_App/fvsetup
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -341,7 +341,7 @@
# 25Jul2017 Holdaway Added resource files associated with NGGPS FV3 tlm/adjoint, sens and svec
# 19Jul2018 Wargan Additions for OMPS
#-------------------------------------------------------------------------
use Cwd qw(cwd);
use Cwd qw(abs_path cwd);
use English;
use Env; # make env vars readily available
use File::Basename qw(basename dirname);
Expand Down Expand Up @@ -2278,29 +2278,28 @@ EOF
#=======================================================================
sub ed_aens_das_replay_acq {
my($mydir) = @_;
my($replay_arcdir, $aens_replay_expid, $acq);
my($replay_arcdir, $replay_expid, $acq);

$aens_replay_expid = "x0044";
$replay_arcdir = query("Replay exp name?", $aens_replay_expid);
$replay_arcdir = "/discover/nobackup/projects/gmao/advda/rtodling/archive/x0044";
$replay_arcdir = "/discover/nobackup/projects/gmao/dadev/rtodling/archive/x0044";
$replay_arcdir = query("Replay archive directory?", $replay_arcdir);
$replay_expid = basename(abs_path($replay_arcdir));

$acq = "$fvhome/$mydir/atmens_replay.acq";
open(SCRIPT,">$acq") or die ">>> ERROR <<< cannot write $acq";
open(SCRIPT,">$acq") or die ">>> ERROR <<< cannot write $acq";
print SCRIPT <<"EOF";
$replay_arcdir/atmens/Y%y4/M%m2/$aens_replay_expid.atmens_ebkg.%y4%m2%d2_%h2z.tar => $expid.atmens_ebkg.%y4%m2%d2_%h2z.tar
$replay_arcdir/atmens/Y%y4/M%m2/$replay_expid.atmens_ebkg.%y4%m2%d2_%h2z.tar => $expid.atmens_ebkg.%y4%m2%d2_%h2z.tar
EOF

if ( $mydir eq "run") {
$acq = "$fvhome/anasa/atmens_replay.acq";
open(SCRIPT,">$acq") or die ">>> ERROR <<< cannot write $acq";
print SCRIPT <<"EOF";
$replay_arcdir/atmens/Y%y4/M%m2/$aens_replay_expid.atmens_ebkg.%y4%m2%d2_%h2z.tar => $expid.atmens_ebkg.%y4%m2%d2_%h2z.tar
$replay_arcdir/atmens/Y%y4/M%m2/$replay_expid.atmens_ebkg.%y4%m2%d2_%h2z.tar => $expid.atmens_ebkg.%y4%m2%d2_%h2z.tar
EOF
$acq = "$fvhome/asens/atmens_asens.acq";
open(SCRIPT,">$acq") or die ">>> ERROR <<< cannot write $acq";
print SCRIPT <<"EOF";
$replay_arcdir/atmens/Y%y4/M%m2/$aens_replay_expid.atmens_ebkg.%y4%m2%d2_%h2z.tar => $expid.atmens_ebkg.%y4%m2%d2_%h2z.tar
$replay_arcdir/atmens/Y%y4/M%m2/$replay_expid.atmens_ebkg.%y4%m2%d2_%h2z.tar => $expid.atmens_ebkg.%y4%m2%d2_%h2z.tar
EOF
}

Expand Down
8 changes: 4 additions & 4 deletions src/Applications/GEOSdas_App/testsuites/C90C_ens.input
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -74,7 +74,7 @@ AeroCom? [/discover/nobackup/projects/gmao/share/gmao_ops/fvInput_4dvar/AeroCom]
>

FVICS? [/archive/u/jstassi/restarts/GEOSadas-5_24_0]
> /archive/u/jstassi/restarts/GEOSadas-5_25_0/C90CS_x0039_p6.rst.20190729_21z.tar
> /archive/u/jstassi/restarts/GEOSadas-5_27_0/C90CS_x0044.rst.20201215_21z.tar

Run model-adjoint-related applications (0=no,1=yes)? [0]
> 1
Expand All @@ -95,7 +95,7 @@ Verifying experiment id: [C90C_ens]
>

Ending year-month-day? [20190731]
> 20190801
> 20201220

Length of FORECAST run segments (in hours)? [123]
>
Expand Down Expand Up @@ -143,7 +143,7 @@ Which main class of ObsSys (1: NRT; 2: MERRA; 3: MERRA-2)? [1]
>

OBSERVING SYSTEM CLASSES?
> ncep_prep_bufr,ncep_1bamua_bufr,disc_airs_bufr,disc_amsua_bufr,ncep_mhs_bufr,ncep_1bhrs4_bufr,ncep_goesfv_bufr,ncep_mtiasi_bufr,ncep_gpsro_bufr,ncep_aura_omi_bufr,ncep_satwnd_bufr,ncep_atms_bufr,ncep_sevcsr_bufr,ncep_ssmis_bufr,mls_nrt_nc,ncep_avcsam_bufr,ncep_avcspm_bufr,ncep_tcvitals,gmao_gmi_bufr,ncep_crisfsr_bufr,npp_ompsnm_bufr,ncep_acftpfl_bufr
> ncep_ssmis_bufr,ncep_1bamua_bufr,ncep_prep_bufr,ncep_tcvitals,ncep_mhs_bufr,ncep_1bhrs4_bufr,ncep_gpsro_bufr,ncep_mtiasi_bufr,ncep_atms_bufr,ncep_crisfsr_bufr,ncep_satwnd_bufr,ncep_sevcsr_bufr,ncep_avcsam_bufr,ncep_avcspm_bufr,ncep_aura_omi_bufr,disc_airs_bufr,disc_amsua_bufr,gmao_gmi_bufr,mls_nrt_nc,gmao_amsr2_bufr,npp_ompsnm_bufr,ncep_acftpfl_bufr

CHECKING OBSYSTEM? [2]
>
Expand Down Expand Up @@ -227,7 +227,7 @@ Ensemble Vertical Levels? [72]
>

Experiment archive directory for ensemble restarts or 'later': [/archive/u/jstassi/C90C_ens]
> /archive/u/dao_it/x0039_p6
> /discover/nobackup/projects/gmao/dadev/rtodling/archive/x0044

Edit COLLECTIONS list in run/HISTORY.rc.tmpl (y/n)? [n]
>
Expand Down
5 changes: 1 addition & 4 deletions src/Applications/GEOSdas_App/testsuites/C90C_replay.input
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -217,11 +217,8 @@ Select group: [g0613]
Replayed Ensemble (from OPS)? [yes]
>

Replay exp name? [x0044]
>

Replay archive directory? [/discover/nobackup/projects/gmao/advda/rtodling/archive/x0044]
>
> /discover/nobackup/projects/gmao/dadev/rtodling/archive/x0044

Edit COLLECTIONS list in run/HISTORY.rc.tmpl (y/n)? [n]
>
Expand Down
26 changes: 21 additions & 5 deletions src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/get_atmens_rst.pl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,10 +1,12 @@
#!/usr/bin/env perl
use strict;
use warnings;
use Cwd qw(cwd);
use Cwd qw(abs_path cwd);
use File::Basename qw(basename dirname);
use File::Copy qw(cp mv);
use File::Path qw(mkpath rmtree);
use Getopt::Long qw(GetOptions);

use FindBin qw($Bin);
use lib "$Bin";
use Manipulate_time qw(tick);
Expand All @@ -30,7 +32,7 @@
my ($atmens_stat_dir, $atmens_ebkg_dir, $atmens_erst_dir, $atmens_ecbkg_dir);
my ($tarfile, $tarpath, $label, $pid);
my ($ens, $mem, $mfile, $mfile_new);
my (@tarList);
my (@archList, @tarList);

init();
chdir($atmens_dir);
Expand All @@ -52,12 +54,12 @@
foreach $label ("stat", "ebkg", "ecbkg", "erst") {
$tarfile = "$expid.atmens_$label.${yyyymmdd}_${hh}z.tar";
$tarpath = "$atmens_date_dir/$tarfile";
push @archList, $tarpath if archFile($tarpath);
push @tarList, $tarpath;
}

defined($pid = fork) or die "Error while attempting to fork;";
unless ($pid) {
system "dmget @tarList";
system "dmget @archList";
exit;
}
foreach $tarpath (@tarList) { system_("tar xvf $tarpath") }
Expand Down Expand Up @@ -88,7 +90,6 @@
# purpose - get runtime parameters and flags
#=======================================================================
sub init {
use Getopt::Long qw(GetOptions);
my ($fvhome, $help, $exparcdir);

GetOptions("fvhome=s" => \$fvhome,
Expand All @@ -102,6 +103,7 @@ sub init {
($exparcdir, $newid, $yyyymmdd, $hh) = @ARGV;

$exparcdir =~ s/[\s\/]*$//;
$exparcdir = abs_path($exparcdir);
$arcdir = dirname($exparcdir);
$expid = basename($exparcdir);

Expand All @@ -123,6 +125,20 @@ sub init {

}

#=======================================================================
# name - archFile
# purpose - Return true (1) if $file is an archive file;
# Return false (0) if not
#=======================================================================
sub archFile {
my ($file, $dmLine, $dmFLG);
$file = shift @_;
$dmLine = `dmls -l $file`;
$dmFLG = 1;
$dmFLG = 0 if $dmLine =~ m|(N/A)|;
return $dmFLG;
}

#=======================================================================
# name - rename_new
#=======================================================================
Expand Down