find nens from ldas ICS when ladas coupling atmdet

   receive nens from fvsetup when ladas coupling atmens
   a minor fix to nml resolution insertion

Author: saraqzhang

GEOSldas_App/ldas_setup

@@ -196,15 +196,15 @@ class LDASsetup:
            assert self.bcs_version   is not None, "Error. Must have command line arg bcs_version for coupled land-atm DAS.\n"
            assert self.rstloc        is not None, "Error. Must have command line arg rstloc      for coupled land-atm DAS.\n"
            assert self.varwindow     is not None, "Error. Must have command line arg varwindow   for coupled land-atm DAS.\n"
-           assert self.nens          is not None, "Error. Must have command line arg nens        for coupled land-atm DAS.\n"
+           if self.ladas_cpl == 2:
+               assert self.nens          is not None, "Error. Must have command line arg nens        for coupled land-atmensDAS.\n"
 
            self.rqdExeInp['BEG_DATE']                 = f"{self.nymdb} {self.nhmsb}"
            rstloc_                                    = self.rstloc.rstrip('/')                  # remove trailing '/'
            assert os.path.isdir(rstloc_)                                                         # make sure rstloc_ is a valid directory
            self.rstloc                                = os.path.abspath(rstloc_)
            self.rqdExeInp['RESTART_PATH']             = os.path.dirname( self.rstloc) 
            self.rqdExeInp['RESTART_ID']               = os.path.basename(self.rstloc)
-
            self.adas_expdir                           = os.path.dirname( self.exphome)
            self.rqdExeInp['ADAS_EXPDIR']              = self.adas_expdir
            self.adas_expid                            = os.path.basename(self.adas_expdir)
@@ -214,16 +214,15 @@ class LDASsetup:
                # ldas coupled with determistic component of adas
                self.rqdExeInp['EXP_ID']               = self.adas_expid  + '_LDAS'
                self.rqdExeInp['MET_PATH']             = self.adas_expdir + '/recycle/holdpredout'
-               self.rqdExeInp['ENSEMBLE_FORCING']     = 'NO'
+               self.rqdExeInp['ENSEMBLE_FORCING']     = 'NO' 
            elif self.ladas_cpl == 2 :
                # ldas coupled with ensemble component of adas
                self.rqdExeInp['EXP_ID']               = self.adas_expid  + '_LDAS4ens'
                self.rqdExeInp['MET_PATH']             = self.adas_expdir + '/atmens/mem'
                self.rqdExeInp['ENSEMBLE_FORCING']     = 'YES' 
            else :
                exit("Error. Unknown value of self.ladas_cpl.\n")
-               
-           self.rqdExeInp['NUM_LDAS_ENSEMBLE']        = self.nens
+           
            self.first_ens_id                          = 1
            self.rqdExeInp['FIRST_ENS_ID']             = self.first_ens_id
 
@@ -261,6 +260,14 @@ class LDASsetup:
            _hours                                     = int(self.rqdExeInp['JOB_SGMT'][ 9:11])
            _end_date                                  = _beg_date + timedelta(hours=int(self.varwindow)/60) 
            self.rqdExeInp['END_DATE']                 = _end_date.strftime("%Y%m%d %H%M%S") 
+
+           # find nens  from  ldas ICS
+           rstdir=self.rqdExeInp['RESTART_PATH']+'/'+self.rqdExeInp['RESTART_ID']+'/output/'+ self.rqdExeInp['RESTART_DOMAIN'] 
+           enslist = glob.glob(os.path.join(rstdir, "rs", "ens????"))  
+           self.rqdExeInp['NUM_LDAS_ENSEMBLE']         = len(enslist) 
+           # if fvsetup send atm_nens via commandline, use that
+           if self.ladas_cpl == 2 :
+              self.rqdExeInp['NUM_LDAS_ENSEMBLE']      = self.nens 
 
            # end if self.ladas_cpl > 0 --------------------------------------------------------------------
 
@@ -1250,7 +1257,8 @@ class LDASsetup:
 
         if self.ladas_cpl > 0:
            # edit resolution info in ensupd nml file
-           sp.run(['sed', '-i', 's/<CFnnnn>/'+ self.agcm_res+'/g', self.rundir+'/LDASsa_SPECIAL_inputs_ensupd.nml'])
+           _agcm_res = 'CF'+self.agcm_res
+           sp.run(['sed', '-i', 's/<CFnnnn>/'+ _agcm_res+'/g', self.rundir+'/LDASsa_SPECIAL_inputs_ensupd.nml'])
 
         # get optimzed NX and IMS
         optimized_distribution_file = tempfile.NamedTemporaryFile(delete=False)