diff --git a/components.yaml b/components.yaml index 2c02c3ec..14f76997 100644 --- a/components.yaml +++ b/components.yaml @@ -28,7 +28,7 @@ NCEP_Shared: GMAO_Shared: local: ./src/Shared/@GMAO_Shared remote: ../GMAO_Shared.git - tag: rt1_4_10_stoch_grid_fix + tag: rt1_4_10_geosit develop: main MAPL: @@ -46,13 +46,13 @@ FMS: GEOSana_GridComp: local: ./src/Components/@GEOSana_GridComp remote: ../GEOSana_GridComp.git - tag: v1.5.4.1 + tag: v1.5.4.2 develop: develop GEOSgcm_GridComp: local: ./src/Components/@GEOSgcm_GridComp remote: ../GEOSgcm_GridComp.git - tag: rt1_12_4_tskinice + tag: rt1_12_4_geosit sparse: ./config/GEOSgcm_GridComp.sparse develop: develop @@ -83,7 +83,7 @@ fvdycore: GEOSchem_GridComp: local: ./src/Components/@GEOSgcm_GridComp/GEOSagcm_GridComp/GEOSphysics_GridComp/@GEOSchem_GridComp remote: ../GEOSchem_GridComp.git - tag: rt1.6.2 + tag: rt1.6.2_geosit develop: develop HEMCO: @@ -101,7 +101,7 @@ geos-chem: GOCART: local: ./src/Components/@GEOSgcm_GridComp/GEOSagcm_GridComp/GEOSphysics_GridComp/@GEOSchem_GridComp/@GOCART remote: ../GOCART.git - tag: rt1.0.1_CEDS + tag: rt1.0.1_geosit sparse: ./config/GOCART.sparse develop: develop @@ -121,7 +121,7 @@ mom6: GEOSgcm_App: local: ./src/Applications/@GEOSgcm_App remote: ../GEOSgcm_App.git - tag: rt1.5.6_geosit_hist1 + tag: rt1.5.6_geosit_hist2 develop: develop UMD_Etc: diff --git a/src/Applications/GEOSdas_App/fvsetup b/src/Applications/GEOSdas_App/fvsetup index 9ad7eb59..94ccefce 100755 --- a/src/Applications/GEOSdas_App/fvsetup +++ b/src/Applications/GEOSdas_App/fvsetup @@ -1216,11 +1216,12 @@ sub get_specific_info { $r21c = 0; $ans = query(" Is this a MERRA2 (1), GEOSIT (2), R21C (3) experiment?", "0"); - return 0 unless ( $ans ); +# return 0 unless ( $ans ); if ($ans == 1) {$merra2 = 1}; if ($ans == 2) {$geosit = 1}; if ($ans == 3) {$r21c = 1}; + return 0; } @@ -2607,6 +2608,7 @@ EOF #======================================================================= sub get_lsmodel { + my $lndbcs_r21c = "Icarus-NLv3"; my $lndbcs_def = "Icarus_Updated"; my $lsdflt = 1; my $rrdflt = 0; @@ -2629,8 +2631,11 @@ LAND BCS Version EOF ; - - $landbcs = query(" Land Boundary Conditions?", $lndbcs_def); + if ( $r21c) { + $landbcs = query(" Land Boundary Conditions?", $lndbcs_r21c); + } else { + $landbcs = query(" Land Boundary Conditions?", $lndbcs_def); + } if ( "$landbcs" eq "Icarus-NLv3" ) { @@ -3138,6 +3143,7 @@ sub get_times { $asynbkg_min = 60; # time frequency of background fields (min) } else { $asynbkg_min = 180; # time frequency of background fields (min) + if ( $r21c) { $asynbkg_min = 60}; } } else { $splite = 0; $splitexe = 0; @@ -3256,8 +3262,7 @@ EOF $asynbkg_hms = sprintf("%6.6d","$asynbkg_hms"); #print "asynbkg_hms = $asynbkg_hms \n"; } - $pcp_forced = $merra2; # for now, we attach these two choices (it can be that PCP-forcing will be applied - # for other non-MERRA2-like runs, but not for a while ... + $pcp_forced = ($merra2 or $r21c ); $beg_date = "$nymdb $nhmsb" ; # GEOS-5 GCM starting date parameter $end_date = "$nymde $nhmse" ; # GEOS-5 GCM ending date parameter @@ -3321,6 +3326,7 @@ sub set_radcor { if ($obClass =~ "merra_upa_pre-qc_bufr" || $obClass =~ "merra_cdas_pre-qc_bufr" || $obClass =~ "merra2_cdas0_pre-qc_bufr" || $obClass =~ "merra2_upa_pre-qc_bufr" || $obClass =~ "merra2_cdas_pre-qc_bufr" ) { $radcor = "HAIMB_HRAD" }; + if($geosit || $r21c) {$radcor = "HAIMB_HRAD"}; # query user for radcor choice #----------------------------- @@ -3438,6 +3444,7 @@ sub get_setgsi { $anaexec = "GSIsa.x"; $hybrid = ".false."; $hyb_ens = -1; + if ($r21c) { $hyb_ens = 4} if ($merra2) { $nosfcana = 0 } # will apply similarity to produce xana.sfc else { $nosfcana = 1 } # will not generation xana.sfc @@ -3557,7 +3564,11 @@ sub set_atmens { # defaults (here for now) # ----------------------- - $aens_replay = "yes"; + if ($r21c) { + $aens_replay = "no"; + } else { + $aens_replay = "yes"; + } $aens_sppt = "yes"; $aens_res = "C90"; $aens_lev = "72"; @@ -3630,7 +3641,10 @@ EOF # If so, setup EnADAS # ------------------- - $flags = "-nlevs $aens_lev $setacftbc $rcorrarg $setradbc $aens_lsmodel -expdir $fvhome/.. $sppt_flag"; + if ($r21c) { + $bcopt = "-r21c"; + } + $flags = "-nlevs $aens_lev $setacftbc $rcorrarg $setradbc $aens_lsmodel -expdir $fvhome/.. $sppt_flag $bcopt"; $flags .= " -fvhome $fvhome" if $checkFLG; $params = " $aens_ana $expid $aens_im $aens_jm $aens_ocn $landbcs"; @@ -3805,11 +3819,13 @@ sub get_obsys { } elsif ($ans == 4) { # GEOS-IT $obsysrc = "$fvetc/obsys-${loc}-geosit.rc"; -# $rflags .= " -stem geosit"; + $obsysGrc = "$fvetc/obsys-${loc}-gaas-geosit.rc"; + $rflags .= " -stem geosit"; } elsif ($ans == 5) { # R21C $obsysrc = "$fvetc/obsys-${loc}-r21c.rc"; -# $rflags .= " -stem r21c"; + $obsysGrc = "$fvetc/obsys-${loc}-gaas-r21c.rc"; + $rflags .= " -stem r21c"; } $reqobs = $fvbin . "/require_obsys.pl"; @@ -4263,11 +4279,13 @@ EOF $emission{"4"} = "OPS"; $emission{"5"} = "PIESA"; $emission{"6"} = "GEOSIT"; + $emission{"7"} = "R21C"; %remission = reverse %emission; if ($nrt) { $dflt = "OPS" } elsif ($merra2) { $dflt = "MERRA2" } elsif ($geosit) { $dflt = "GEOSIT" } + elsif ($r21c) { $dflt = "R21C" } else { $dflt = "PIESA" } if ($gocart_tracers) { @@ -4670,6 +4688,9 @@ sub get_dimsg5gcm { $res = "C48"; $vres = 72; if ( $merra2 ) { $ores = "e34" } + elsif ($r21c) { $ores = "C"; + $res = "C360"; + } else { $ores = "f" } print <<"EOF"; @@ -5448,6 +5469,7 @@ sub archiving_rules { \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_msu_%c%c%c_%c%c%c.%y4%m2%d2_%h2z.bin \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_msu_tirosn_%c%c%c.%y4%m2%d2_%h2z.bin \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_omi_aura_%c%c%c.%y4%m2%d2_%h2z.bin +\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslp_g_npp_%c%c%c.%y4%m2%d2_%h2z.bin \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpuv_npp_%c%c%c.%y4%m2%d2_%h2z.bin \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpvis_npp_%c%c%c.%y4%m2%d2_%h2z.bin \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompsnm_npp_%c%c%c.%y4%m2%d2_%h2z.bin @@ -5505,6 +5527,7 @@ sub archiving_rules { \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_msu_%c%c%c.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_msu_tirosn.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_omi_aura.%y4%m2%d2_%h2z.ods +\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslp_g_npp.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpuv_npp.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpvis_npp.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompsnm_npp.%y4%m2%d2_%h2z.ods @@ -5561,6 +5584,7 @@ sub archiving_rules { \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp0hr_diag_msu_%c%c%c.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp0hr_diag_msu_tirosn.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp0hr_diag_omi_aura.%y4%m2%d2_%h2z.ods +\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp0hr_diag_ompslp_g_npp.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp0hr_diag_ompslpuv_npp.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp0hr_diag_ompslpvis_npp.%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp0hr_diag_ompsnm_npp.%y4%m2%d2_%h2z.ods @@ -5601,10 +5625,11 @@ sub archiving_rules { #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_seviri_m%c%c_%c%c%c_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_sbuv%c_%c%c%c_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_omi_aura_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslp_g_npp_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpuv_npp_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpvis_npp_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin -#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpnm_npp_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin -#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompslpnp_npp_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompsnm_npp_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_ompsnp_npp_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_hirs%c_%c%c%c_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_hirs%c_metop-%c_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.diag_hirs%c_tirosn_%c%c%c_%c%c%c%c.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.bin @@ -5637,10 +5662,11 @@ sub archiving_rules { #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_pcp_tmi_%c%c%c%c_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_sbuv%c_%c%c%c_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_omi_aura_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslp_g_npp_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpuv_npp_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpvis_npp_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods -#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnm_npp_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods -#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnp_npp_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompsnm_npp_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompsnp_npp_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_hirs%c_%c%c%c_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_hirs%c_metop-%c_%c%c%c_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_hirs%c_tirosn_%c%c%c_%c%c%c_%c%c%c%c.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods @@ -5674,10 +5700,11 @@ sub archiving_rules { #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_pcp_tmi_%c%c%c%c_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_sbuv%c_%c%c%c_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_omi_aura_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslp_g_npp_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpuv_npp_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpvis_npp_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods -#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnm_npp_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods -#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnp_npp_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompsnm_npp_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompsnp_npp_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_hirs%c_%c%c%c_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_hirs%c_metop-%c_%c%c%c_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods #\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_hirs%c_tirosn_%c%c%c_%c%c%c_%c%c%c%c.sigo.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods @@ -6868,10 +6895,11 @@ sub arch_asens { \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_msu_tirosn.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_omi_aura.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_omieff_aura.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslp_g_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpuv_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpvis_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods -\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnm_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods -\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnp_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnm_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods +#\${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompslpnp_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompsnm_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompsnmeff_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods \${PESTOROOT}%s/obs/Y%y4/M%m2/D%d2/H%h2/%s.imp%c_%c%c%c_ompsnp_npp.obs.%y4%m2%d2_%h2z+%y4%m2%d2_%h2z-%y4%m2%d2_%h2z.ods @@ -7451,12 +7479,16 @@ close(SCRIPT); #======================================================================= sub build_lnbcs { - $mymerra2 = ""; - if ( $merra2 ) { $mymerra2 = "-merra2" }; + $bcopt = ""; + if ( $merra2 ) { + $bcopt = "-merra2" ; + }elsif ( $r21c ) { + $bcopt = "-r21c"; + } if ( $cubed ) { - $cmd = "$fvbin/gen_lnbcs.pl -fvhome $fvhome -o lnbcs $mymerra2 -cubed $agcm_im $agcm_jm $ores $landbcs"; + $cmd = "$fvbin/gen_lnbcs.pl -fvhome $fvhome -o lnbcs $bcopt -cubed $agcm_im $agcm_jm $ogcm $landbcs"; } else { - $cmd = "$fvbin/gen_lnbcs.pl -fvhome $fvhome -o lnbcs $mymerra2 $agcm_im $agcm_jm $ores $landbcs"; + $cmd = "$fvbin/gen_lnbcs.pl -fvhome $fvhome -o lnbcs $bcopt $agcm_im $agcm_jm $ogcm $landbcs"; } print "$cmd\n"; system($cmd); @@ -8685,8 +8717,9 @@ sub write_FVDAS_Run_Config { $ENV{"VTXLEVS"} = $vtxlevs; $ENV{"VTXRELOC"} = $vtxreloc; - if ($geosit | $r21c) { + if ($geosit || $r21c) { $ENV{"MKSI_SIDB"} = "\$FVHOME/run/gmao_satinfo.db"; + $ENV{"MKSI_OZDB"} = "\$FVHOME/run/gmao_ozinfo.db"; } unless ( $ENV{"ARCHIVE"} ) { @@ -9734,6 +9767,7 @@ sub init_agcm_rc { # default is based on NCA settings (but line is commented since we don't want this done) if ( $pcp_forced ) { $pcp_fntmpl = "d5_merra/Y%y4/M%m2/d5_merra.tavg1_2d_lfo_Nx_corr.%y4%m2%d2_%h230z.nc"; + #if (r21c) { $pcp_fntmpl = " "}; # ADD when pcp files are ready } else { if ( $nymd < 19890101 ) { $pcp_fntmpl = "d5_merra_jan79/diag/Y%y4/d5_merra_jan79.tavg1_2d_lfo_Nx_corr.%y4%m2%d2_%h230z.nc"; @@ -10062,6 +10096,26 @@ sub copy_resources { # Create daotovs_fv.rc in FVHOME/run # ed_daotovs_fv(); + # When applicable, overwrite w/ templated RC files for GSI + # ------------------------------------------------------ + if ( $r21c ){ + if ( -d "$fvetc/gsi/R21C" ) { + my @files = glob("$fvetc/gsi/R21C" . "/*"); + foreach my $fn ( @files ) { + cp("$fn", "$fvhome/run"); + } + } + if ( -d "$fvetc/atmens/R21C" ){ + my @files = glob("$fvetc/atmens/R21C" . "/*"); + foreach my $fn ( @files ) { + cp("$fn","$fvhome/run/atmens"); + } + } + } + # Update gocart files for pre-2000 case + # To Be Done + # + # Edit psas.rc ed_psas_rc(); @@ -10166,34 +10220,29 @@ sub copy_resources { cp("$fn", "$fvhome/run"); } } + } + if ( $r21c ) { $casedir = "R21C" }; + if ( $geosit ) { $casedir = "GEOSIT"}; + if ( $r21c || $geosit ){ mkdir ("$fvhome/run/gmao_satinfo.db"); - if ( -d "$fvetc/gmao_satinfo.db/GEOSIT" ) { - my @files = glob("$fvetc/gmao_satinfo.db/GEOSIT" . "/*"); - foreach my $fn ( @files ) { - cp("$fn", "$fvhome/run/gmao_satinfo.db"); - } + if ( -d "$fvetc/gmao_satinfo.db/$casedir" ) { + my @files = glob("$fvetc/gmao_satinfo.db/$casedir" . "/*"); + foreach my $fn ( @files ) { + cp("$fn", "$fvhome/run/gmao_satinfo.db"); + } } else { - die "Cannot find GEOSIT sat.db under $fvetc, aborting ..."; - } - } - if ( $r21c ) { - if ( -d "$fvetc/gsi/R21C" ) { - my @files = glob("$fvetc/gsi/R21C" . "/*"); - foreach my $fn ( @files ) { - cp("$fn", "$fvhome/run"); - } + die "Cannot find $casedir sat.db under $fvetc, aborting ..."; } - mkdir ("$fvhome/run/gmao_satinfo.db"); - if ( -d "$fvetc/gmao_satinfo.db/R21C" ) { - my @files = glob("$fvetc/gmao_satinfo.db/R21C" . "/*"); + mkdir ("$fvhome/run/gmao_ozinfo.db"); + if ( -d "$fvetc/gmao_ozinfo.db/$casedir" ) { + my @files = glob("$fvetc/gmao_ozinfo.db/$casedir" . "/*"); foreach my $fn ( @files ) { - cp("$fn", "$fvhome/run/gmao_satinfo.db"); + cp("$fn", "$fvhome/run/gmao_ozinfo.db"); } } else { - die "Cannot find R21C sat.db under $fvetc, aborting ..."; + die "Cannot find $casedir oz.db under $fvetc, aborting ..."; } } - } #======================================================================= diff --git a/src/Applications/GEOSdas_App/gen_lnbcs.pl b/src/Applications/GEOSdas_App/gen_lnbcs.pl index dbefbc6e..6dcca493 100755 --- a/src/Applications/GEOSdas_App/gen_lnbcs.pl +++ b/src/Applications/GEOSdas_App/gen_lnbcs.pl @@ -32,6 +32,7 @@ "sstdir=s", "cubed", "merra2", + "r21c", "h" ); usage() if $opt_h; @@ -80,6 +81,8 @@ sub init { if ( $opt_merra2 ) { $pcp_loc = "/discover/nobackup/projects/gmao/share/dao_ops/fvInput/merra_land/precip_CPCUexcludeAfrica-CMAP_corrected_MERRA/GEOSdas-2_1_4"; +} elsif ( $opt_r21c ) { # This will updated + $pcp_loc = "/discover/nobackup/projects/gmao/share/dao_ops/fvInput/merra_land/precip_CPCUexcludeAfrica-CMAP_corrected_MERRA/GEOSdas-2_1_4"; } else { $pcp_loc = "/gpfsm/dnb51/projects/p15/iau/merra_land/precip_CPCU-CMAP_corrected_MERRA/GEOSdas-2_1_4"; } @@ -425,6 +428,7 @@ sub usage { -fvhome location of FVHOME (default: write script locally) -cubed needed for cubed GCM -merra2 specify to set related BCs + -r21c specify to set related BCs -h prints this usage notice EXAMPLE COMMAND LINE diff --git a/src/Applications/GEOSdas_App/testsuites/geos_it.input b/src/Applications/GEOSdas_App/testsuites/geos_it.input index 699ccad9..ae7301ea 100644 --- a/src/Applications/GEOSdas_App/testsuites/geos_it.input +++ b/src/Applications/GEOSdas_App/testsuites/geos_it.input @@ -45,7 +45,7 @@ The directory /discover/nobackup/projects/gmao/dadev/rtodling/geos_it does not e > Processing nodes (1:Westmere, 2:SandyBridge, 3:Ivy Bridge, 4:Haswell, 5:Skylake, 6:Cascase)? [4] -> 5 +> 6 Which case of variational analysis? [1] > @@ -138,7 +138,7 @@ Which main class of ObsSys (1: NRT; 2: MERRA; 3: MERRA-2; 4: GEOS-IT; 5: R21C)? > 4 OBSERVING SYSTEM CLASSES? -> merra2_cdas0_pre-qc_bufr,merra2_cdas_pre-qc_bufr,merra2_avhrrwnd_pre-qc_bufr,merra2_ascat_pre-qc_bufr,merra2_ers1_pre-qc_bufr,merra2_repro_ers2_pre-qc_bufr,merra2_qscat_jpl_pre-qc_bufr,merra2_wspd_pre-qc_bufr,merra2_nmodis_pre-qc_bufr,merra2_prof_pre-qc_bufr,merra2_cdas0_pre-qc_bufr,merra2_cdas_pre-qc_bufr,merra2_avhrrwnd_pre-qc_bufr,merra2_ascat_pre-qc_bufr,merra2_ers1_pre-qc_bufr,merra2_repro_ers2_pre-qc_bufr,merra2_qscat_jpl_pre-qc_bufr,merra2_wspd_pre-qc_bufr,merra2_nmodis_pre-qc_bufr,merra2_prof_pre-qc_bufr,merra2_ncep_tcvitals,ncep_ssmis_bufr,ncep_1bamua_bufr,ncep_mhs_bufr,ncep_1bhrs4_bufr,ncep_goesfv_bufr,ncep_gpsro_bufr,ncep_mtiasi_bufr,ncep_atms_bufr,ncep_cris_bufr,ncep_crisfsr_bufr,ncep_satwnd_bufr,ncep_sevcsr_bufr,ncep_avcsam_bufr,ncep_avcspm_bufr,disc_airs_bufr,disc_amsua_bufr,ncep_acftpfl_bufr,aura_omieff_nc,npp_ompsnp_nc,npp_ompsnmeff_nc,gmao_gmi_bufr,gmao_amsr2_bufr +> geosit_cdas_raob_pre-qc_bufr,geosit_cdas_conv_pre-qc_bufr,geosit_prof_pre-qc_bufr,geosit_avhrrwnd_pre-qc_bufr,geosit_nmodis_pre-qc_bufr,geosit_goes_pre-qc_prep_bufr,geosit_metsat_pre-qc_prep_bufr,geosit_jma_pre-qc_prep_bufr,geosit_ascat_pre-qc_bufr,geosit_repro_ers2_pre-qc_bufr,geosit_qscat_jpl_pre-qc_bufr,geosit_wspd_pre-qc_bufr,geosit_satwnd_bufr,geosit_avhrr_satwnd_bufr,geosit_ncep_tcvitals,geosit_tmi_bufr,geosit_gpsro_bufr,geosit_sevcsr_bufr,geosit_1bamua_bufr,geosit_1bamub_bufr,geosit_1bhrs2_bufr,geosit_1bhrs3_bufr,geosit_1bhrs4_bufr,geosit_1bmsu_bufr,geosit_1bmhs_bufr,geosit_1bssu_bufr,geosit_eosairs_bufr,geosit_eosamsua_bufr,geosit_mtiasi_bufr,geosit_atms_bufr,geosit_ssmit11_bufr,geosit_ssmit13_bufr,geosit_ssmit14_bufr,geosit_ssmit15_bufr,geosit_amsre_bufr,geosit_osbuv8_bufr,geosit_npp_ompsnp_nc,geosit_aura_omieff_nc,geosit_npp_ompsnmeff_nc,geosit_avcsam_bufr,geosit_avcspm_bufr,geosit_acftpfl_bufr,geosit_amsr2_bufr,geosit_crisfsr_bufr,geosit_gmi_bufr CHECKING OBSYSTEM? [2] > 1 @@ -195,7 +195,7 @@ Do Aerosol Analysis (y/n)? [y] > AOD OBSERVING CLASSES [or type 'none']? -> +> mod04_061_his,myd04_061_his Enable GAAS feedback to model (y/n)? [y] > diff --git a/src/Applications/GEOSdas_App/write_FVDAS_Run_Config.pl b/src/Applications/GEOSdas_App/write_FVDAS_Run_Config.pl index 12598afd..5e9a5b8a 100755 --- a/src/Applications/GEOSdas_App/write_FVDAS_Run_Config.pl +++ b/src/Applications/GEOSdas_App/write_FVDAS_Run_Config.pl @@ -14,7 +14,7 @@ my ($ARCH, $HOST); my ($FVHOME, $FVROOT, $RUNDIR); my ($AOD_OBSCLASS, $BERROR, $DO_ECS_OUT, $DO_REM_SYNC, $EXPID, $FVARCH, - $FVBCS, $GID, $MONTHLY_MEANS, $MKSI_SIDB, $MP_SET_NUMTHREADS, $NCEPINPUT, $NOBACKUP, + $FVBCS, $GID, $MONTHLY_MEANS, $MKSI_SIDB, $MKSI_OZDB, $MP_SET_NUMTHREADS, $NCEPINPUT, $NOBACKUP, $OBSCLASS, $OBSCLASS_NOAIRS, $OMP_NUM_THREADS, $RUN_QUADS, $PYRADMON, $VTRACK, $VTXLEVS, $VTXRELOC); my ($BASEDIR, $FCSTID, $FVDMGET, $G5MODULES, $PLOTS_LOC, $GEOSUTIL, $GTAG); @@ -93,6 +93,7 @@ sub init { $GID = $ENV{"GID"}; $MONTHLY_MEANS = $ENV{"MONTHLY_MEANS"}; $MKSI_SIDB = $ENV{"MKSI_SIDB"}; + $MKSI_OZDB = $ENV{"MKSI_OZDB"}; $MP_SET_NUMTHREADS = $ENV{"MP_SET_NUMTHREADS"}; $NCEPINPUT = $ENV{"NCEPINPUT"}; $OBSCLASS = $ENV{"OBSCLASS"}; @@ -323,6 +324,7 @@ sub writefile { print RUNCONF "setenv FVDOLMS $FVDOLMS\n" if $FVDOLMS; print RUNCONF "setenv CASE $CASE\n" if $CASE; print RUNCONF "setenv MKSI_SIDB $MKSI_SIDB\n" if $MKSI_SIDB; + print RUNCONF "setenv MKSI_OZDB $MKSI_OZDB\n" if $MKSI_OZDB; print RUNCONF "setenv MP_SET_NUMTHREADS $MP_SET_NUMTHREADS\n" if $MP_SET_NUMTHREADS; print RUNCONF "setenv OMP_NUM_THREADS $OMP_NUM_THREADS\n" if $OMP_NUM_THREADS; print RUNCONF "setenv ARCH_QUEUE \"$ARCH_QUEUE\"\n" if $ARCH_QUEUE; diff --git a/src/Applications/GEOSdas_App/write_monthly_rc_arc.pl b/src/Applications/GEOSdas_App/write_monthly_rc_arc.pl index 1d494f23..3e1422ca 100755 --- a/src/Applications/GEOSdas_App/write_monthly_rc_arc.pl +++ b/src/Applications/GEOSdas_App/write_monthly_rc_arc.pl @@ -32,13 +32,17 @@ # # => %list: keys = analysis products, e.g. "ana.eta", and names of model # collections found in $histfile, e.g. "inst1_2d_asm_Nx" -# value = a single code containing one or more of the following +# values = a single code containing one or more of the following # letters, {C, M, T, P}, with the following meanings: # C => comment entry in monthly.rc # M => create monthly products # T => tar this product # P => create monthly plots for this product -# the code, except for the 'C', is written to the monthly.rc file +# the code, except for the 'C', is written to the monthly.rc file +# => %freq: keys = names of collections found in $histfile +# values = hour frequency (e.g. 1, 3, or 6) as specified in $histfile +# => %mode: keys = names of collections found in $histfile +# values = mode ("inst" or "tavg") as specified in $histfile # # key local variables in main program: # => @outRcArr: lines to write to the $outRc file (monthly.rc) @@ -56,9 +60,11 @@ # global variables #----------------- -my ($histfile, $siloarc, $ncana, %list); +my ($histfile, $siloarc, $ncana); my ($outRc, $outArc, $outArk); my ($rcFLG, $arcFLG, $arkFLG); +my (%list, %freq, %mode); + my $script = basename $0; # analysis and aod output names with default hour types @@ -100,6 +106,7 @@ "prog.sfc" => 1, "ptrj.prs" => 1, "traj.lcv" => 1, + "traj_lcv_rst" => 1, "vtx.mix" => 1, "vtx.prs" => 1 ); # main program @@ -183,15 +190,8 @@ sub init { #======================================================================= # name - get_list_from_HIST # purpose - add collection names from $histfile to list -# -# key local variables -# => %freq: keys = names of collections found in $histfile -# values = hour frequency (e.g. 1, 3, or 6) as specified in $histfile -# => %mode: keys = names of collections found in $histfile -# values = mode ("inst" or "tavg") as specified in $histfile #======================================================================= sub get_list_from_HIST { - my (%freq, %mode); my ($colFLG, $extFLG); my ($colon, $doublecolon); @@ -227,8 +227,8 @@ sub get_list_from_HIST { # get frequency and mode information #----------------------------------- - extract_freq($_, \%freq) if /\.frequency\s*$colon/; - extract_mode($_, \%mode) if /\.mode\s*$colon/; + extract_freq($_) if /\.frequency\s*$colon/; + extract_mode($_) if /\.mode\s*$colon/; } close HIST; } @@ -302,18 +302,12 @@ sub plots_type { # # input parameter # => $line: line from $histfile containing frequency information for a collection -# -# input/output parameter -# => $freqAddr: address of %freq hash containing collection frequency info #======================================================================= sub extract_freq { - my ($line, $freqAddr); - + my ($line); $line = shift; - $freqAddr = shift; - if ( $line =~ /^\s*(\S*)\.frequency\s*:\s*(\d{6})\s*,/ ) { - $$freqAddr{$1} = $2 / 10000; + $freq{$1} = $2 / 10000; } } @@ -323,22 +317,14 @@ sub extract_freq { # # input parameter # => $line: line from $histfile containing mode information for a collection -# -# input/output parameter -# => $modeAddr: address of %mode hash containing collection mode info #======================================================================= sub extract_mode { - my ($line, $modeAddr); - my ($name); - + my ($line, $name); $line = shift; - $modeAddr = shift; - - $name = undef; - $name = $1 if $line =~ /^\s*(\S*)\.mode\s*:/; - if ($name) { - $$modeAddr{$name} = "inst" if $line =~ /instantaneous/; - $$modeAddr{$name} = "tavg" if $line =~ /time-averaged/; + if ( $line =~ /^\s*(\S*)\.mode\s*:/ ) { + $name = $1; + $mode{$name} = "inst" if $line =~ /instantaneous/; + $mode{$name} = "tavg" if $line =~ /time-averaged/; } } @@ -401,6 +387,12 @@ sub get_info_from_SILO { } close SILO; + # add additional hourtype requirements + #------------------------------------- + foreach (qw(traj_lcv_rst)) { + $hourtype{$_} = "$mode{$_}$freq{$_}"; + } + # entries for monthly.rc file #---------------------------- @$outRcArrAddr = (); diff --git a/src/Applications/NCEP_Etc/NCEP_bias/readme_init_coeffs.txt b/src/Applications/NCEP_Etc/NCEP_bias/readme_init_coeffs.txt index 87f5655a..ccf5ce84 100644 --- a/src/Applications/NCEP_Etc/NCEP_bias/readme_init_coeffs.txt +++ b/src/Applications/NCEP_Etc/NCEP_bias/readme_init_coeffs.txt @@ -45,6 +45,13 @@ directory with names in the format %s.diag_(dtype)_(dplat)_ges.%y4%m2%d2_%h2z.bin where the %s is substituted with the expid specified on the command line +Note: Since this program uses the "gsi.rc.tmpl" file to get values for 'dtype', +'dplat', and 'dsis' for each satellite instrument being fitted, if you are +using the 'generic' value for these variables in your "gsi.rc.tmpl" to +configure your instrument you will need to supply another file with an +'OBS_INPUT::' table like in the "gsi.rc.tmpl" which specifies the actual +'dtype', 'dplat' and 'dsis' values for your satellite(s). + --------------------- EXAMPLES -------------------------------------------- 1) Using archived diag_*_ges.*.bin files as input diff --git a/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/CMakeLists.txt b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/CMakeLists.txt index 949528ab..6a0f7409 100644 --- a/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/CMakeLists.txt +++ b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/CMakeLists.txt @@ -39,3 +39,14 @@ install ( FILES ${ALLETC} DESTINATION etc/atmens ) + +set (R21C_files + R21C/atmos_enkf.nml.tmpl + R21C/obs1gsi_mean.rc + R21C/obs1gsi_member.rc + ) + +install ( + FILES ${R21C_files} + DESTINATION etc/atmens/R21C + ) diff --git a/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/atmos_enkf.nml.tmpl b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/atmos_enkf.nml.tmpl new file mode 100644 index 00000000..04172e78 --- /dev/null +++ b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/atmos_enkf.nml.tmpl @@ -0,0 +1,142 @@ +&nam_enkf + datestring=">>>YYYYMMDDHH<<<",datapath="./", + expid=">>>EXPID<<<", + analpertwtnh=0.85,analpertwtsh=0.85,analpertwttr=0.85, + covinflatemax=1.e2,covinflatemin=1,pseudo_rh=.true., + corrlengthnh=2000,corrlengthsh=2000,corrlengthtr=2000, + lnsigcutoffnh=2.0,lnsigcutoffsh=2.0,lnsigcutofftr=2.0, + lnsigcutoffsatnh=2.0,lnsigcutoffsatsh=2.0,lnsigcutoffsattr=2.0, + lnsigcutoffpsnh=2.0,lnsigcutoffpssh=2.0,lnsigcutoffpstr=2.0, + obtimelnh=1.e30,obtimelsh=1.e30,obtimeltr=1.e30, + saterrfact=1.0,numiter=1, + sprd_tol=1.e30,paoverpb_thresh=0.990,iassim_order=0, + nlevs=>>>ENS_NLEVS<<<,nanals=>>>NMEM<<<, + nsfcvars=2, +! ntrac=3, + nvars=9, + stats_usedob_only=.true., + deterministic=.true.,sortinc=.true.,lupd_satbiasc=.false., + readin_localization=.true., + nlats=>>>ENS_NLATS<<<,nlons=>>>ENS_NLONS<<<, +! nhr_anal(1)=3, +! nhr_anal(2)=4, +! nhr_anal(3)=5, +! nhr_anal(4)=6, +! nhr_anal(5)=7, +! nhr_anal(6)=8, +! nhr_anal(7)=9, +@RADBC newpc4pred=.true.,adp_anglebc=.true.,angord=4, +@RADBC use_edges=.false.,emiss_bc=.true., +! +! forward EFSO + fso_flag = .false., + fso_cycling = .true., + / +&satobs_enkf + sattypes_rad(1) = 'amsua_n15', dsis(1) = 'amsua_n15', + sattypes_rad(2) = 'amsub_n15', dsis(2) = 'amsub_n15', + sattypes_rad(3) = 'hirs3_n16', dsis(3) = 'hirs3_n16', + sattypes_rad(4) = 'amsua_n16', dsis(4) = 'amsua_n16', + sattypes_rad(5) = 'amsub_n16', dsis(5) = 'amsub_n16', + sattypes_rad(6) = 'hirs3_n17', dsis(6) = 'hirs3_n17', + sattypes_rad(7) = 'amsua_n17', dsis(7) = 'amsua_n17', + sattypes_rad(8) = 'amsub_n17', dsis(8) = 'amsub_n17', + sattypes_rad(9) = 'hirs4_n18', dsis(9) = 'hirs4_n18', + sattypes_rad(10)= 'amsua_n18', dsis(10)= 'amsua_n18', +! sattypes_rad(11)= 'amsub_n18', dsis(11)= 'amsub_n18', + sattypes_rad(12)= 'mhs_n18', dsis(12)= 'mhs_n18', + sattypes_rad(13)= 'hirs4_n19', dsis(13)= 'hirs4_n19', + sattypes_rad(14)= 'amsua_n19', dsis(14)= 'amsua_n19', +! sattypes_rad(15)= 'amsub_n19', dsis(15)= 'amsub_n19', + sattypes_rad(16)= 'mhs_n19', dsis(16)= 'mhs_n19', + sattypes_rad(17)= 'hirs4_metop-a', dsis(17)= 'hirs4_metop-a', + sattypes_rad(18)= 'amsua_metop-a', dsis(18)= 'amsua_metop-a', + sattypes_rad(19)= 'amsua_metop-b', dsis(19)= 'amsua_metop-b', + sattypes_rad(20)= 'iasi_metop-a', dsis(20)= 'iasi_metop-a', + sattypes_rad(21)= 'airs_aqua', dsis(21)= 'airs_aqua', + sattypes_rad(22)= 'amsua_aqua', dsis(22)= 'amsua_aqua', +! sattypes_rad(15)= 'goes_img_g11', dsis(15)= 'imgr_g11', +! sattypes_rad(16)= 'goes_img_g12', dsis(16)= 'imgr_g12', +! sattypes_rad(17)= 'goes_img_g13', dsis(17)= 'imgr_g13', +! sattypes_rad(18)= 'goes_img_g14', dsis(18)= 'imgr_g14', +! sattypes_rad(19)= 'goes_img_g15', dsis(19)= 'imgr_g15', + sattypes_rad(23)= 'amsre_aqua', dsis(23)= 'amsre_aqua', +! sattypes_rad(24)= 'ssmis_f16', dsis(24)= 'ssmis_f16', + sattypes_rad(25)= 'ssmis_f17', dsis(25)= 'ssmis_f17', + sattypes_rad(26)= 'ssmis_f18', dsis(26)= 'ssmis_f18', +! sattypes_rad(27)= 'ssmis_f19', dsis(27)= 'ssmis_f19', +! sattypes_rad(28)= 'ssmis_f20', dsis(28)= 'ssmis_f20', +! sattypes_rad(29)= 'sndrd1_g11', dsis(29)= 'sndrD1_g11', +! sattypes_rad(30)= 'sndrd2_g11', dsis(30)= 'sndrD2_g11', +! sattypes_rad(31)= 'sndrd3_g11', dsis(31)= 'sndrD3_g11', +! sattypes_rad(32)= 'sndrd4_g11', dsis(32)= 'sndrD4_g11', +! sattypes_rad(33)= 'sndrd1_g12', dsis(33)= 'sndrD1_g12', +! sattypes_rad(34)= 'sndrd2_g12', dsis(34)= 'sndrD2_g12', +! sattypes_rad(35)= 'sndrd3_g12', dsis(35)= 'sndrD3_g12', +! sattypes_rad(36)= 'sndrd4_g12', dsis(36)= 'sndrD4_g12', +! sattypes_rad(37)= 'sndrd1_g13', dsis(37)= 'sndrD1_g13', +! sattypes_rad(38)= 'sndrd2_g13', dsis(38)= 'sndrD2_g13', +! sattypes_rad(39)= 'sndrd3_g13', dsis(39)= 'sndrD3_g13', +! sattypes_rad(40)= 'sndrd4_g13', dsis(40)= 'sndrD4_g13', +! sattypes_rad(41)= 'sndrd1_g14', dsis(41)= 'sndrD1_g14', +! sattypes_rad(42)= 'sndrd2_g14', dsis(42)= 'sndrD2_g14', +! sattypes_rad(43)= 'sndrd3_g14', dsis(43)= 'sndrD3_g14', +! sattypes_rad(44)= 'sndrd4_g14', dsis(44)= 'sndrD4_g14', +! sattypes_rad(45)= 'sndrd1_g15', dsis(45)= 'sndrD1_g15', +! sattypes_rad(46)= 'sndrd2_g15', dsis(46)= 'sndrD2_g15', +! sattypes_rad(47)= 'sndrd3_g15', dsis(47)= 'sndrD3_g15', +! sattypes_rad(48)= 'sndrd4_g15', dsis(48)= 'sndrD4_g15', + sattypes_rad(49)= 'mhs_metop-a', dsis(49)= 'mhs_metop-a', + sattypes_rad(50)= 'atms_npp', dsis(50)= 'atms_npp', + sattypes_rad(51)= 'cris_npp', dsis(51)= 'cris_npp', +! sattypes_rad(52)= 'seviri_m08', dsis(52)= 'seviri_m08', +! sattypes_rad(53)= 'seviri_m09', dsis(53)= 'seviri_m09', +! sattypes_rad(54)= 'seviri_m10', dsis(54)= 'seviri_m10', + sattypes_rad(55)= 'mhs_metop-b', dsis(55)= 'mhs_metop-b', + sattypes_rad(56)= 'iasi_metop-b', dsis(56)= 'iasi_metop-b', + sattypes_rad(57)= 'avhrr_n18', dsis(57)= 'avhrr3_n18', + sattypes_rad(58)= 'avhrr_n19', dsis(58)= 'avhrr3_n19', + sattypes_rad(59)= 'avhrr_metop-a', dsis(59)= 'avhrr3_metop-a', + sattypes_rad(60)= 'avhrr_metop-b', dsis(60)= 'avhrr3_metop-b', + sattypes_rad(61)= 'gmi_gpm', dsis(61)= 'gmi_gpm', + sattypes_rad(62)= 'amsua_metop-c', dsis(62)= 'amsua_metop-c', + sattypes_rad(63)= 'mhs_metop-c', dsis(63)= 'mhs_metop-c', + sattypes_rad(64)= 'cris-fsr_n20', dsis(64)= 'cris-fsr_n20', + sattypes_rad(65)= 'cris-fsr_npp', dsis(65)= 'cris-fsr_npp', + sattypes_rad(66)= 'iasi_metop-c', dsis(66)= 'iasi_metop-c', + sattypes_rad(67)= 'msu_n11', dsis(67)= 'msu_n11', + sattypes_rad(68)= 'msu_n12', dsis(68)= 'msu_n12', + sattypes_rad(69)= 'msu_n14', dsis(69)= 'msu_n14', + sattypes_rad(70)= 'ssmi_f13', dsis(70)= 'ssmi_f13', + sattypes_rad(71)= 'ssmi_f14', dsis(71)= 'ssmi_f14', + sattypes_rad(72)= 'ssmi_f15', dsis(72)= 'ssmi_f15', + sattypes_rad(73)= 'tmi_trmm', dsis(73)= 'tmi_trmm', + sattypes_rad(74)= 'avhrr_n15', dsis(74)= 'avhrr3_n15', + sattypes_rad(75)= 'avhrr_n16', dsis(75)= 'avhrr3_n16', + sattypes_rad(76)= 'avhrr_n17', dsis(76)= 'avhrr3_n17', + sattypes_rad(77)= 'atms_n20', dsis(77)= 'atms_n20', + sattypes_rad(78)= 'amsr2_gcom-w1', dsis(78)= 'amsr2_gcom-w1', +/ +&ozobs_enkf + sattypes_oz(1) = 'sbuv8_n11', + sattypes_oz(2) = 'sbuv8_n14', + sattypes_oz(3) = 'sbuv8_n16', + sattypes_oz(4) = 'sbuv8_n17', + sattypes_oz(5) = 'sbuv8_n18', + sattypes_oz(6) = 'sbuv8_n19', + sattypes_oz(7) = 'mls55_aura', +!sattypes_oz(8) = 'omi_aura', +!sattypes_oz(9) = 'gome_metop-a', +!sattypes_oz(10) = 'gome_metop-b', +!sattypes_oz(11) = 'o3lev_aura', +/ + +spectral_truncation: 96 + +eps_for_log_transform_aod: -1.0 + +fso_ana_increment_test: no + +analysis_type4fso: null + +verification_type4fso: null diff --git a/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/obs1gsi_mean.rc b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/obs1gsi_mean.rc new file mode 100644 index 00000000..2e5adc68 --- /dev/null +++ b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/obs1gsi_mean.rc @@ -0,0 +1,195 @@ + &SETUP + miter=0,niter(1)=1,niter(2)=1, + niter_no_qc(1)=999,niter_no_qc(2)=999, + write_diag(1)=.true.,write_diag(2)=.false.,write_diag(3)=.false., + jiterstart=1, + gencode=82,qoption=2, + factqmin=0.005,factqmax=0.005,deltim=300, +!--gencode=82,qoption=1, +!--factqmin=0.8,factqmax=1.2,deltim=300, + ifact10=0, + pseudo_q2=.true., + use_prepb_satwnd=>>>USE_PREPB_SATWND<<<, + ec_amv_qc=.false., + iguess=-1, + diag_version=30303, + id_drifter=.true., + id_ship=.false., + tzr_qc=1, + oneobtest=.false.,retrieval=.false., + biascor=-0.10,bcoption=0,diurnalbc=1.0, + crtm_coeffs_path="CRTM_Coeffs/", + print_diag_pcg=.false., + use_compress=.true.,nsig_ext=@NSIG_EXT,gpstop=60.,lgpsbnd_revint=.true., + use_sp_eqspace=.true., + reduce_diag=.true., + luse_obsdiag=.true., + lread_obs_save=.true.,lread_obs_skip=.false.,lwrite_predterms=.true.,lwrite_peakwt=.true., + ens_nstarthr=3, +@RADBC newpc4pred=.true.,adp_anglebc=.true.,angord=4, +@RADBC passive_bc=.true.,use_edges=.false., +@RADBC diag_precon=.true.,step_start=1.e-3,emiss_bc=.true., + lrun_subdirs=.false., + / + &GRIDOPTS + JCAP=@GSI_JCAP,NLAT=@GSI_JM,NLON=@GSI_IM,nsig=@GSI_LM, + regional=.false., + nvege_type=13, + / + &BKGERR + vs=0.6, + hzscl=0.588,1.25,2.0, + hswgt=0.45,0.3,0.25, + bw=0.0,norsp=4, + bkgv_flowdep=.true.,bkgv_rewgtfct=1.5, + fpsproj=.false., + adjustozvar=.true., + / + &ANBKGERR + anisotropic=.false., + / + &JCOPTS + ljcpdry=.true.,bamp_jcpdry=2.5e7, + / + &STRONGOPTS +! tlnmc_option=1,nstrong=1,nvmodes_keep=24,period_max=6.,period_width=1.5, +! baldiag_full=.true.,baldiag_inc=.true., + / + &OBSQC + dfact=0.75,dfact1=3.0,noiqc=.false.,oberrflg=.true.,c_varqc=0.02,blacklst=.true., + use_poq7=.true.,qc_noirjaco3=.true.,qc_satwnds=.false.,cld_det_dec2bin=.true., +! tcp_ermin=0.75,tcp_ermax=0.75, + >>>AIRCFT_BIAS<<< + / + &OBS_INPUT + dmesh(1)=145.0,dmesh(2)=150.0,dmesh(3)=180.0,time_window_max=3.0, + / +OBS_INPUT:: +! dfile dtype dplat dsis dval dthin dsfcalc obclass + prepbufr ps null ps 0.0 0 0 gmao_prep_bufr + prepbufr t null t 0.0 0 0 gmao_prep_bufr + prepbufr q null q 0.0 0 0 gmao_prep_bufr + prepbufr uv null uv 0.0 0 0 gmao_prep_bufr + prepbufr_profl t prof t 0.0 0 0 gmao_acftpfl_bufr + prepbufr_profl uv prof uv 0.0 0 0 gmao_acftpfl_bufr + mlstbufr t aura t 0.0 0 0 r21c_gmao_mlst_bufr + gpsrobufr gps_bnd null gps 0.0 0 0 r21c_gpsro_bufr + tcvitl tcp null tcp 0.0 0 0 r21c_ncep_tcvitals + sbuvbufr sbuv2 n11 sbuv8_n11 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n14 sbuv8_n14 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n16 sbuv8_n16 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n17 sbuv8_n17 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n18 sbuv8_n18 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n19 sbuv8_n19 0.0 0 0 r21c_osbuv8_bufr + ompsnmeffnc ompsnmeff npp ompsnmeff_npp 0.0 2 0 r21c_npp_ompsnmeff_nc +! ompsnpnc ompsnpnc npp ompsnpnc_npp 0.0 0 0 r21c_npp_ompsnp_nc + ompslpgnc ompslp_g npp ompslp_g_npp 1.0 0 0 r21c_ompslp_g_nc + ompslpvisnc ompslpvis npp ompslpvis_npp 0.0 0 0 r21c_ompslpvis_nc + mlsnc mls55 aura mls55_aura 0.0 0 0 r21c_mls_nc + omieffnc omieff aura omieff_aura 0.0 2 0 r21c_aura_omieff_nc + hirs2bufr hirs2 n11 hirs2_n11 0.0 1 0 r21c_1bhrs2_bufr + hirs2bufr hirs2 n12 hirs2_n12 0.0 1 0 r21c_1bhrs2_bufr + hirs2bufr hirs2 n14 hirs2_n14 0.0 1 0 r21c_1bhrs2_bufr + hirs3bufr hirs3 n15 hirs3_n15 0.0 1 0 r21c_1bhrs3_bufr + hirs3bufr hirs3 n16 hirs3_n16 0.0 1 0 r21c_1bhrs3_bufr + hirs3bufr hirs3 n17 hirs3_n17 0.0 1 0 r21c_1bhrs3_bufr + hirs4bufr hirs4 n18 hirs4_n18 0.0 1 0 r21c_1bhrs4_bufr + hirs4bufr hirs4 n19 hirs4_n19 0.0 1 0 r21c_1bhrs4_bufr + hirs4bufr hirs4 metop-a hirs4_metop-a 0.0 1 0 r21c_1bhrs4_bufr + hirs4bufr hirs4 metop-b hirs4_metop-b 0.0 1 0 r21c_1bhrs4_bufr + airsbufr airs aqua airs_aqua 0.0 3 0 r21c_eosairs_bufr + eosamsua amsua aqua amsua_aqua 0.0 1 0 r21c_eosamsua_bufr + msubufr msu n11 msu_n11 0.0 1 0 r21c_1bmsu_bufr + msubufr msu n12 msu_n12 0.0 1 0 r21c_1bmsu_bufr + msubufr msu n14 msu_n14 0.0 1 0 r21c_1bmsu_bufr + ssubufr ssu n14 ssu_n14 0.0 1 0 r21c_1bssu_bufr + amsuabufr amsua n15 amsua_n15 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n16 amsua_n16 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n17 amsua_n17 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n18 amsua_n18 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n19 amsua_n19 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua metop-a amsua_metop-a 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua metop-b amsua_metop-b 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua metop-c amsua_metop-c 0.0 1 0 r21c_1bamua_bufr + amsubbufr amsub n15 amsub_n15 0.0 1 0 r21c_1bamub_bufr + amsubbufr amsub n16 amsub_n16 0.0 1 0 r21c_1bamub_bufr + amsubbufr amsub n17 amsub_n17 0.0 1 0 r21c_1bamub_bufr + mhsbufr mhs n18 mhs_n18 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs n19 mhs_n19 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs metop-a mhs_metop-a 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs metop-b mhs_metop-b 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs metop-c mhs_metop-c 0.0 1 0 r21c_mhs_bufr + seviribufr seviri m08 seviri_m08 0.0 1 0 r21c_sevcsr_bufr + seviribufr seviri m09 seviri_m09 0.0 1 0 r21c_sevcsr_bufr + seviribufr seviri m10 seviri_m10 0.0 1 0 r21c_sevcsr_bufr + smit11bufr ssmi f11 ssmi_f11 0.0 1 0 r21c_ssmit11_bufr + smit13bufr ssmi f13 ssmi_f13 0.0 1 0 r21c_ssmit13_bufr + smit14bufr ssmi f14 ssmi_f14 0.0 1 0 r21c_ssmit14_bufr + smit15bufr ssmi f15 ssmi_f15 0.0 1 0 r21c_ssmit15_bufr + iasibufr iasi metop-a iasi_metop-a 0.0 3 0 r21c_mtiasi_bufr + iasibufr iasi metop-b iasi_metop-b 0.0 3 0 r21c_mtiasi_bufr + iasibufr iasi metop-c iasi_metop-c 0.0 3 0 r21c_mtiasi_bufr + atmsbufr atms npp atms_npp 0.0 1 0 r21c_atms_bufr + atmsbufr atms n20 atms_n20 0.0 1 0 r21c_atms_bufr + crisfsrbufr cris-fsr npp cris-fsr_npp 0.0 3 0 r21c_crisfsr_bufr + crisfsrbufr cris-fsr n20 cris-fsr_n20 0.0 3 0 r21c_crisfsr_bufr + tmibufr tmi trmm tmi_trmm 0.0 1 0 r21c_tmi_bufr + gmibufr gmi gpm gmi_gpm 0.0 1 0 r21c_gmi_bufr + amsr2bufr amsr2 gcom-w1 amsr2_gcom-w1 0.0 1 0 r21c_amsr2_bufr + amsregmao amsre aqua amsre_aqua 0.0 1 0 r21c_amsre_bufr + satwndbufr uv null uv 0.0 0 0 r21c_satwnd_bufr + satwndavhr uv null uv 0.0 0 0 r21c_avhrr_satwnd_bufr + avcsambufr avhrr n15 avhrr3_n15 0.0 1 0 r21c_avcsam_bufr + avcsambufr avhrr n17 avhrr3_n17 0.0 1 0 r21c_avcsam_bufr + avcsambufr avhrr metop-a avhrr3_metop-a 0.0 1 0 r21c_avcsam_bufr + avcsambufr avhrr metop-b avhrr3_metop-b 0.0 1 0 r21c_avcsam_bufr + avcspmbufr avhrr n16 avhrr3_n16 0.0 1 0 r21c_avcspm_bufr + avcspmbufr avhrr n18 avhrr3_n18 0.0 1 0 r21c_avcspm_bufr + avcspmbufr avhrr n19 avhrr3_n19 0.0 1 0 r21c_avcspm_bufr +!!!!! +! prepbufr spd null spd 0.0 0 0 gmao_prep_bufr +! prepbufr pw null pw 0.0 0 0 gmao_prep_bufr +! ascatbufr uv null uv 0.0 0 0 merra2_ascat_bufr +! ssmisbufr ssmis f17 ssmis_f17 0.0 1 0 hist_ssmis_bufr +! ssmisbufr ssmis f18 ssmis_f18 0.0 1 0 hist_ssmis_bufr +! crisbufr cris npp cris_npp 0.0 3 0 hist_cris_bufr +! ompsnmbufr ompsnm npp ompsnm_npp 0.0 2 0 npp_ompsnm_bufr +! ompsnpbufr ompsnp npp ompsnp_npp 0.0 0 0 npp_ompsnp_bufr +! ompslpuvnc ompslpuv npp ompslpuv_npp 1.0 0 0 ompslpuv_nc +! sbuvbufr sbuv2 nim07 sbuv8_nim07 0.0 0 0 osbuv8_bufr +:: + &SUPEROB_RADAR + / + &LAG_DATA +! lag_accur=1e-6, +! infile_lag='inistate_lag.dat', +! lag_stepduration=900., +! lag_nmax_bal=100, +! lag_vorcore_stderr_a=2e3, +! lag_vorcore_stderr_b=0, + / + &HYBRID_ENSEMBLE + l_hyb_ens=.false.,n_ens=32,beta_s0=0.50,generate_ens=.false.,uv_hyb_ens=.true., + s_ens_h=800.,s_ens_v=-0.5, + jcap_ens=126,nlat_ens=181,nlon_ens=288,aniso_a_en=.false., + jcap_ens_test=126, + oz_univ_static=.true., +! use_localization_grid=.true., + readin_localization=.true., + readin_beta=.true., + use_gfs_ens=.false., + eqspace_ensgrid=.true., + / + &RAPIDREFRESH_CLDSURF +! dfi_radar_latent_heat_time_period=30.0, + / + &SINGLEOB_TEST +! maginnov=0.1,magoberr=0.1,oneob_type='t', +! oblat=45.,oblon=180.,obpres=1000.,obdattim=2004041512, +! obhourset=0., + / + &CHEM + / + &NST + nst_gsi=3, + / diff --git a/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/obs1gsi_member.rc b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/obs1gsi_member.rc new file mode 100644 index 00000000..5f2c8075 --- /dev/null +++ b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/etc/R21C/obs1gsi_member.rc @@ -0,0 +1,185 @@ + &SETUP + miter=0,niter(1)=1,niter(2)=1, + niter_no_qc(1)=999,niter_no_qc(2)=999, + write_diag(1)=.true.,write_diag(2)=.false.,write_diag(3)=.false., + jiterstart=1, + gencode=82,qoption=2, + factqmin=0.005,factqmax=0.005,deltim=300, +!--gencode=82,qoption=1, +!--factqmin=0.8,factqmax=1.2,deltim=300, + ifact10=0, + pseudo_q2=.true., + use_prepb_satwnd=>>>USE_PREPB_SATWND<<<, + ec_amv_qc=.false., + iguess=-1, + diag_version=30303, + id_drifter=.true., + id_ship=.false., + tzr_qc=1, + oneobtest=.false.,retrieval=.false., + biascor=-0.10,bcoption=0,diurnalbc=1.0, + crtm_coeffs_path="CRTM_Coeffs/", + print_diag_pcg=.false., + use_compress=.true.,nsig_ext=@NSIG_EXT,gpstop=60.,lgpsbnd_revint=.true., + use_sp_eqspace=.true., + reduce_diag=.true., + luse_obsdiag=.true., + lread_obs_save=.false.,lread_obs_skip=.true.,lwrite_predterms=.true.,lwrite_peakwt=.true., + ens_nstarthr=3, +@RADBC newpc4pred=.true.,adp_anglebc=.true.,angord=4, +@RADBC passive_bc=.true.,use_edges=.false., +@RADBC diag_precon=.true.,step_start=1.e-3,emiss_bc=.true., + lrun_subdirs=.false., + / + &GRIDOPTS + JCAP=@GSI_JCAP,NLAT=@GSI_JM,NLON=@GSI_IM,nsig=@GSI_LM, + regional=.false., + nvege_type=13, + / + &BKGERR + vs=0.6, + hzscl=0.588,1.25,2.0, + hswgt=0.45,0.3,0.25, + bw=0.0,norsp=4, + bkgv_flowdep=.true.,bkgv_rewgtfct=1.5, + fpsproj=.false., + adjustozvar=.true., + / + &ANBKGERR + anisotropic=.false., + / + &JCOPTS + ljcpdry=.true.,bamp_jcpdry=2.5e7, + / + &STRONGOPTS +! tlnmc_option=1,nstrong=1,nvmodes_keep=24,period_max=6.,period_width=1.5, +! baldiag_full=.true.,baldiag_inc=.true., + / + &OBSQC + dfact=0.75,dfact1=3.0,noiqc=.false.,oberrflg=.true.,c_varqc=0.02,blacklst=.true., + use_poq7=.true.,qc_noirjaco3=.true.,qc_satwnds=.false.,cld_det_dec2bin=.true., +! tcp_ermin=0.75,tcp_ermax=0.75, + >>>AIRCFT_BIAS<<< + / + &OBS_INPUT + dmesh(1)=145.0,dmesh(2)=150.0,dmesh(3)=180.0,time_window_max=3.0, + / +OBS_INPUT:: +! dfile dtype dplat dsis dval dthin dsfcalc obclass + prepbufr ps null ps 0.0 0 0 gmao_prep_bufr + prepbufr t null t 0.0 0 0 gmao_prep_bufr + prepbufr q null q 0.0 0 0 gmao_prep_bufr + prepbufr uv null uv 0.0 0 0 gmao_prep_bufr + prepbufr_profl t prof t 0.0 0 0 gmao_acftpfl_bufr + prepbufr_profl uv prof uv 0.0 0 0 gmao_acftpfl_bufr + mlstbufr t aura t 0.0 0 0 r21c_gmao_mlst_bufr + gpsrobufr gps_bnd null gps 0.0 0 0 r21c_gpsro_bufr + tcvitl tcp null tcp 0.0 0 0 r21c_ncep_tcvitals + sbuvbufr sbuv2 n11 sbuv8_n11 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n14 sbuv8_n14 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n16 sbuv8_n16 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n17 sbuv8_n17 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n18 sbuv8_n18 0.0 0 0 r21c_osbuv8_bufr + sbuvbufr sbuv2 n19 sbuv8_n19 0.0 0 0 r21c_osbuv8_bufr + ompsnmeffnc ompsnmeff npp ompsnmeff_npp 0.0 2 0 r21c_npp_ompsnmeff_nc +! ompsnpnc ompsnpnc npp ompsnpnc_npp 0.0 0 0 r21c_npp_ompsnp_nc + ompslpgnc ompslp_g npp ompslp_g_npp 1.0 0 0 r21c_ompslp_g_nc + ompslpvisnc ompslpvis npp ompslpvis_npp 0.0 0 0 r21c_ompslpvis_nc + mlsnc mls55 aura mls55_aura 0.0 0 0 r21c_mls_nc + omieffnc omieff aura omieff_aura 0.0 2 0 r21c_aura_omieff_nc + hirs2bufr hirs2 n11 hirs2_n11 0.0 1 0 r21c_1bhrs2_bufr + hirs2bufr hirs2 n12 hirs2_n12 0.0 1 0 r21c_1bhrs2_bufr + hirs2bufr hirs2 n14 hirs2_n14 0.0 1 0 r21c_1bhrs2_bufr + hirs3bufr hirs3 n15 hirs3_n15 0.0 1 0 r21c_1bhrs3_bufr + hirs3bufr hirs3 n16 hirs3_n16 0.0 1 0 r21c_1bhrs3_bufr + hirs3bufr hirs3 n17 hirs3_n17 0.0 1 0 r21c_1bhrs3_bufr + hirs4bufr hirs4 n18 hirs4_n18 0.0 1 0 r21c_1bhrs4_bufr + hirs4bufr hirs4 n19 hirs4_n19 0.0 1 0 r21c_1bhrs4_bufr + hirs4bufr hirs4 metop-a hirs4_metop-a 0.0 1 0 r21c_1bhrs4_bufr + hirs4bufr hirs4 metop-b hirs4_metop-b 0.0 1 0 r21c_1bhrs4_bufr + airsbufr airs aqua airs_aqua 0.0 3 0 r21c_eosairs_bufr + eosamsua amsua aqua amsua_aqua 0.0 1 0 r21c_eosamsua_bufr + msubufr msu n11 msu_n11 0.0 1 0 r21c_1bmsu_bufr + msubufr msu n12 msu_n12 0.0 1 0 r21c_1bmsu_bufr + msubufr msu n14 msu_n14 0.0 1 0 r21c_1bmsu_bufr + ssubufr ssu n14 ssu_n14 0.0 1 0 r21c_1bssu_bufr + amsuabufr amsua n15 amsua_n15 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n16 amsua_n16 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n17 amsua_n17 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n18 amsua_n18 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua n19 amsua_n19 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua metop-a amsua_metop-a 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua metop-b amsua_metop-b 0.0 1 0 r21c_1bamua_bufr + amsuabufr amsua metop-c amsua_metop-c 0.0 1 0 r21c_1bamua_bufr + amsubbufr amsub n15 amsub_n15 0.0 1 0 r21c_1bamub_bufr + amsubbufr amsub n16 amsub_n16 0.0 1 0 r21c_1bamub_bufr + amsubbufr amsub n17 amsub_n17 0.0 1 0 r21c_1bamub_bufr + mhsbufr mhs n18 mhs_n18 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs n19 mhs_n19 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs metop-a mhs_metop-a 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs metop-b mhs_metop-b 0.0 1 0 r21c_1bmhs_bufr + mhsbufr mhs metop-c mhs_metop-c 0.0 1 0 r21c_mhs_bufr + seviribufr seviri m08 seviri_m08 0.0 1 0 r21c_sevcsr_bufr + seviribufr seviri m09 seviri_m09 0.0 1 0 r21c_sevcsr_bufr + seviribufr seviri m10 seviri_m10 0.0 1 0 r21c_sevcsr_bufr + smit11bufr ssmi f11 ssmi_f11 0.0 1 0 r21c_ssmit11_bufr + smit13bufr ssmi f13 ssmi_f13 0.0 1 0 r21c_ssmit13_bufr + smit14bufr ssmi f14 ssmi_f14 0.0 1 0 r21c_ssmit14_bufr + smit15bufr ssmi f15 ssmi_f15 0.0 1 0 r21c_ssmit15_bufr + iasibufr iasi metop-a iasi_metop-a 0.0 3 0 r21c_mtiasi_bufr + iasibufr iasi metop-b iasi_metop-b 0.0 3 0 r21c_mtiasi_bufr + iasibufr iasi metop-c iasi_metop-c 0.0 3 0 r21c_mtiasi_bufr + atmsbufr atms npp atms_npp 0.0 1 0 r21c_atms_bufr + atmsbufr atms n20 atms_n20 0.0 1 0 r21c_atms_bufr + crisfsrbufr cris-fsr npp cris-fsr_npp 0.0 3 0 r21c_crisfsr_bufr + crisfsrbufr cris-fsr n20 cris-fsr_n20 0.0 3 0 r21c_crisfsr_bufr + tmibufr tmi trmm tmi_trmm 0.0 1 0 r21c_tmi_bufr + gmibufr gmi gpm gmi_gpm 0.0 1 0 r21c_gmi_bufr + amsr2bufr amsr2 gcom-w1 amsr2_gcom-w1 0.0 1 0 r21c_amsr2_bufr + amsregmao amsre aqua amsre_aqua 0.0 1 0 r21c_amsre_bufr + satwndbufr uv null uv 0.0 0 0 r21c_satwnd_bufr + satwndavhr uv null uv 0.0 0 0 r21c_avhrr_satwnd_bufr + avcsambufr avhrr n15 avhrr3_n15 0.0 1 0 r21c_avcsam_bufr + avcsambufr avhrr n17 avhrr3_n17 0.0 1 0 r21c_avcsam_bufr + avcsambufr avhrr metop-a avhrr3_metop-a 0.0 1 0 r21c_avcsam_bufr + avcsambufr avhrr metop-b avhrr3_metop-b 0.0 1 0 r21c_avcsam_bufr + avcspmbufr avhrr n16 avhrr3_n16 0.0 1 0 r21c_avcspm_bufr + avcspmbufr avhrr n18 avhrr3_n18 0.0 1 0 r21c_avcspm_bufr + avcspmbufr avhrr n19 avhrr3_n19 0.0 1 0 r21c_avcspm_bufr +!!!!! +! prepbufr spd null spd 0.0 0 0 gmao_prep_bufr +! prepbufr pw null pw 0.0 0 0 gmao_prep_bufr +! ascatbufr uv null uv 0.0 0 0 merra2_ascat_bufr +! ssmisbufr ssmis f17 ssmis_f17 0.0 1 0 hist_ssmis_bufr +! ssmisbufr ssmis f18 ssmis_f18 0.0 1 0 hist_ssmis_bufr +! crisbufr cris npp cris_npp 0.0 3 0 hist_cris_bufr +! ompsnmbufr ompsnm npp ompsnm_npp 0.0 2 0 npp_ompsnm_bufr +! ompsnpbufr ompsnp npp ompsnp_npp 0.0 0 0 npp_ompsnp_bufr +! ompslpuvnc ompslpuv npp ompslpuv_npp 1.0 0 0 ompslpuv_nc +! sbuvbufr sbuv2 nim07 sbuv8_nim07 0.0 0 0 osbuv8_bufr +:: + &SUPEROB_RADAR + / + &LAG_DATA +! lag_accur=1e-6, +! infile_lag='inistate_lag.dat', +! lag_stepduration=900., +! lag_nmax_bal=100, +! lag_vorcore_stderr_a=2e3, +! lag_vorcore_stderr_b=0, + / + &HYBRID_ENSEMBLE + / + &RAPIDREFRESH_CLDSURF +! dfi_radar_latent_heat_time_period=30.0, + / + &SINGLEOB_TEST +! maginnov=0.1,magoberr=0.1,oneob_type='t', +! oblat=45.,oblon=180.,obpres=1000.,obdattim=2004041512, +! obhourset=0., + / + &CHEM + / + &NST + nst_gsi=3, + / diff --git a/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/setup_atmens.pl b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/setup_atmens.pl index e5741fd0..c8b42bff 100755 --- a/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/setup_atmens.pl +++ b/src/Applications/NCEP_Etc/NCEP_enkf/scripts/gmao/setup_atmens.pl @@ -49,6 +49,7 @@ "nosppt", "ose", "rcorr", + "r21c", "h" ); usage() if $opt_h; @@ -121,6 +122,11 @@ sub init { } } + $bcopt = ""; + if ( $opt_r21c ) { + $bcopt = "-r21c"; + } + $setradbc = 0; if ( $opt_radbc ) { $setradbc = 1; @@ -420,7 +426,8 @@ sub install { cp("$FVROOT/bin/atm_ens.j","$FVHOME/run"); # generate boundary condition script -$cmd = "$FVROOT/bin/gen_lnbcs.pl $cubed -o $FVHOME/run/lnbcs_ens $aim $ajm $ogrid $lndbcs"; +# $cmd = "$FVROOT/bin/gen_lnbcs.pl $cubed -o $FVHOME/run/lnbcs_ens $aim $ajm $ogrid $lndbcs"; +$cmd = "$FVROOT/bin/gen_lnbcs.pl $cubed $bcopt -o $FVHOME/run/lnbcs_ens $aim $ajm $ogrid $lndbcs"; $rc = system($cmd); # make sure .no_archiving exists in ATMENS