Skip to content

Commit

Permalink
Merge remote-tracking branch 'EMC/develop' into feature/toff_fix
Browse files Browse the repository at this point in the history
  • Loading branch information
@XuLu-NOAA
XuLu-NOAA committed Mar 24, 2024
2 parents 100c488 + 2167bc9 commit d6a4b09
Show file tree
Hide file tree
Showing 29 changed files with 975 additions and 452 deletions.
2 changes: 1 addition & 1 deletion fix
Submodule fix updated 2 files
+1 −0 .gitignore
+1 −0 gsi_binary_files.cmake
2 changes: 1 addition & 1 deletion modulefiles/gsi_gaea.intel.lua
View file
Original file line number Diff line number Diff line change
Expand Up @@ -31,7 +31,7 @@ local MKLROOT="/opt/intel/oneapi/mkl/2022.0.2/"
prepend_path("LD_LIBRARY_PATH",pathJoin(MKLROOT,"lib/intel64"))
pushenv("MKLROOT", MKLROOT)

pushenv("GSI_BINARY_SOURCE_DIR", "/lustre/f2/dev/role.epic/contrib/GSI_data/fix/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/lustre/f2/dev/role.epic/contrib/GSI_data/fix/20240208")

setenv("CC","cc")
setenv("FC","ftn")
Expand Down
2 changes: 1 addition & 1 deletion modulefiles/gsi_hera.gnu.lua
View file
Original file line number Diff line number Diff line change
Expand Up @@ -22,6 +22,6 @@ load("gsi_common")
load(pathJoin("prod_util", prod_util_ver))
load(pathJoin("openblas", openblas_ver))

pushenv("GSI_BINARY_SOURCE_DIR", "/scratch1/NCEPDEV/global/glopara/fix/gsi/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/scratch1/NCEPDEV/global/glopara/fix/gsi/20240208")

whatis("Description: GSI environment on Hera with GNU Compilers")
4 changes: 2 additions & 2 deletions modulefiles/gsi_hera.intel.lua
View file
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
help([[
]])

prepend_path("MODULEPATH", "/scratch1/NCEPDEV/nems/role.epic/spack-stack/spack-stack-1.6.0/envs/gsi-addon-dev/install/modulefiles/Core")
prepend_path("MODULEPATH", "/scratch1/NCEPDEV/nems/role.epic/spack-stack/spack-stack-1.6.0/envs/gsi-addon-dev-rocky8/install/modulefiles/Core")

local python_ver=os.getenv("python_ver") or "3.11.6"
local stack_intel_ver=os.getenv("stack_intel_ver") or "2021.5.0"
Expand All @@ -20,6 +20,6 @@ load(pathJoin("prod_util", prod_util_ver))
pushenv("CFLAGS", "-xHOST")
pushenv("FFLAGS", "-xHOST")

pushenv("GSI_BINARY_SOURCE_DIR", "/scratch1/NCEPDEV/global/glopara/fix/gsi/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/scratch1/NCEPDEV/global/glopara/fix/gsi/20240208")

whatis("Description: GSI environment on Hera with Intel Compilers")
2 changes: 1 addition & 1 deletion modulefiles/gsi_hercules.intel.lua
View file
Original file line number Diff line number Diff line change
Expand Up @@ -21,6 +21,6 @@ load("intel-oneapi-mkl/2022.2.1")
pushenv("CFLAGS", "-xHOST")
pushenv("FFLAGS", "-xHOST")

pushenv("GSI_BINARY_SOURCE_DIR", "/work/noaa/global/glopara/fix/gsi/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/work/noaa/global/glopara/fix/gsi/20240208")

whatis("Description: GSI environment on Hercules with Intel Compilers")
2 changes: 1 addition & 1 deletion modulefiles/gsi_jet.intel.lua
View file
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,6 @@ load(pathJoin("prod_util", prod_util_ver))
pushenv("CFLAGS", "-axSSE4.2,AVX,CORE-AVX2")
pushenv("FFLAGS", "-axSSE4.2,AVX,CORE-AVX2")

pushenv("GSI_BINARY_SOURCE_DIR", "/mnt/lfs4/HFIP/hfv3gfs/glopara/git/fv3gfs/fix/gsi/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/mnt/lfs4/HFIP/hfv3gfs/glopara/git/fv3gfs/fix/gsi/20240208")

whatis("Description: GSI environment on Jet with Intel Compilers")
2 changes: 1 addition & 1 deletion modulefiles/gsi_orion.intel.lua
View file
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,6 @@ load(pathJoin("prod_util", prod_util_ver))
pushenv("CFLAGS", "-xHOST")
pushenv("FFLAGS", "-xHOST")

pushenv("GSI_BINARY_SOURCE_DIR", "/work/noaa/global/glopara/fix/gsi/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/work/noaa/global/glopara/fix/gsi/20240208")

whatis("Description: GSI environment on Orion with Intel Compilers")
2 changes: 1 addition & 1 deletion modulefiles/gsi_s4.intel.lua
View file
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,6 @@ load(pathJoin("prod_util", prod_util_ver))
pushenv("CFLAGS", "-march=ivybridge")
pushenv("FFLAGS", "-march=ivybridge")

pushenv("GSI_BINARY_SOURCE_DIR", "/data/prod/glopara/fix/gsi/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/data/prod/glopara/fix/gsi/20240208")

whatis("Description: GSI environment on S4 with Intel Compilers")
2 changes: 1 addition & 1 deletion modulefiles/gsi_wcoss2.intel.lua
View file
Original file line number Diff line number Diff line change
Expand Up @@ -46,6 +46,6 @@ load(pathJoin("ncio", ncio_ver))
load(pathJoin("crtm", crtm_ver))
load(pathJoin("ncdiag",ncdiag_ver))

pushenv("GSI_BINARY_SOURCE_DIR", "/lfs/h2/emc/global/noscrub/emc.global/FIX/fix/gsi/20230911")
pushenv("GSI_BINARY_SOURCE_DIR", "/lfs/h2/emc/global/noscrub/emc.global/FIX/fix/gsi/20240208")

whatis("Description: GSI environment on WCOSS2")
55 changes: 28 additions & 27 deletions regression/global_4denvar.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -55,14 +55,15 @@ cycg=`echo $gdate | cut -c9-10`
dumpobs=gdas
prefix_obs=${dumpobs}.t${cyca}z
prefix_ges=gdas.t${cycg}z
prefix_ens=gdas.t${cycg}z
prefix_ens=enkfgdas.t${cycg}z
suffix=tm00.bufr_d

dumpges=gdas
COMROOTgfs=$casesdir/gfs/prod
datobs=$COMROOTgfs/$dumpobs.$PDYa/${cyca}/atmos
datges=$COMROOTgfs/$dumpges.$PDYg/${cycg}/atmos
datens=$COMROOTgfs/enkfgdas.$PDYg/${cycg}/atmos
datobs=$COMROOTgfs/$dumpobs.$PDYa/${cyca}/obs
dathis=$COMROOTgfs/$dumpges.$PDYg/${cycg}/model_data/atmos/history
datanl=$COMROOTgfs/gdas.$PDYg/${cycg}/analysis/atmos
datens=$COMROOTgfs/enkfgdas.$PDYg/${cycg}


# Set up $tmpdir
Expand Down Expand Up @@ -265,28 +266,28 @@ $nln $datobs/${prefix_obs}.esatms.${suffix} ./atmsbufrears
## $nln $datobs/${prefix_obs}.amsr2.tm00.bufr_d ./amsr2bufr

# Copy bias correction, atmospheric and surface files
$nln $datges/${prefix_ges}.abias ./satbias_in
$nln $datges/${prefix_ges}.abias_pc ./satbias_pc
$nln $datges/${prefix_ges}.abias_air ./aircftbias_in
$nln $datges/${prefix_ges}.radstat ./radstat.gdas

$nln $datges/${prefix_ges}.sfcf003.nc ./sfcf03
$nln $datges/${prefix_ges}.sfcf004.nc ./sfcf04
$nln $datges/${prefix_ges}.sfcf005.nc ./sfcf05
$nln $datges/${prefix_ges}.sfcf006.nc ./sfcf06
$nln $datges/${prefix_ges}.sfcf007.nc ./sfcf07
$nln $datges/${prefix_ges}.sfcf008.nc ./sfcf08
$nln $datges/${prefix_ges}.sfcf009.nc ./sfcf09

$nln $datges/${prefix_ges}.atmf003.nc ./sigf03
$nln $datges/${prefix_ges}.atmf004.nc ./sigf04
$nln $datges/${prefix_ges}.atmf005.nc ./sigf05
$nln $datges/${prefix_ges}.atmf006.nc ./sigf06
$nln $datges/${prefix_ges}.atmf007.nc ./sigf07
$nln $datges/${prefix_ges}.atmf008.nc ./sigf08
$nln $datges/${prefix_ges}.atmf009.nc ./sigf09

$nln $datens/${prefix_ens}.sfcf006.ensmean.nc ./sfcf06_anlgrid
$nln $datanl/${prefix_ges}.abias ./satbias_in
$nln $datanl/${prefix_ges}.abias_pc ./satbias_pc
$nln $datanl/${prefix_ges}.abias_air ./aircftbias_in
$nln $datanl/${prefix_ges}.radstat ./radstat.gdas

$nln $dathis/${prefix_ges}.sfcf003.nc ./sfcf03
$nln $dathis/${prefix_ges}.sfcf004.nc ./sfcf04
$nln $dathis/${prefix_ges}.sfcf005.nc ./sfcf05
$nln $dathis/${prefix_ges}.sfcf006.nc ./sfcf06
$nln $dathis/${prefix_ges}.sfcf007.nc ./sfcf07
$nln $dathis/${prefix_ges}.sfcf008.nc ./sfcf08
$nln $dathis/${prefix_ges}.sfcf009.nc ./sfcf09

$nln $dathis/${prefix_ges}.atmf003.nc ./sigf03
$nln $dathis/${prefix_ges}.atmf004.nc ./sigf04
$nln $dathis/${prefix_ges}.atmf005.nc ./sigf05
$nln $dathis/${prefix_ges}.atmf006.nc ./sigf06
$nln $dathis/${prefix_ges}.atmf007.nc ./sigf07
$nln $dathis/${prefix_ges}.atmf008.nc ./sigf08
$nln $dathis/${prefix_ges}.atmf009.nc ./sigf09

$nln $datens/ensstat/model_data/atmos/history/${prefix_ens}.sfcf006.ensmean.nc ./sfcf06_anlgrid

export ENS_PATH='./ensemble_data/'
mkdir -p ${ENS_PATH}
Expand All @@ -296,7 +297,7 @@ for fh in $flist; do
imem=1
while [[ $imem -le $NMEM_ENKF ]]; do
member="mem"`printf %03i $imem`
$nln $datens/$member/$sigens ${ENS_PATH}sigf${fh}_ens_${member}
$nln $datens/$member/model_data/atmos/history/$sigens ${ENS_PATH}sigf${fh}_ens_${member}
(( imem = $imem + 1 ))
done
done
Expand Down
19 changes: 8 additions & 11 deletions regression/global_enkf.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -51,17 +51,14 @@ cyca=`echo $global_adate | cut -c9-10`
PDYg=`echo $gdate | cut -c1-8`
cycg=`echo $gdate | cut -c9-10`

dumpobs=gdas
prefix_obs=${dumpobs}.t${cyca}z
prefix_ges=gdas.t${cycg}z
prefix_ens=gdas.t${cycg}z
prefix_obs=enkfgdas.t${cyca}z
prefix_ens=enkfgdas.t${cycg}z
suffix=tm00.bufr_d

dumpges=gdas
COMROOTgfs=$casesdir/gfs/prod
datobs=$COMROOTgfs/enkfgdas.$PDYa/${cyca}/atmos
datges=$COMROOTgfs/$dumpges.$PDYg/${cycg}/atmos
datens=$COMROOTgfs/enkfgdas.$PDYg/${cycg}/atmos
datobs=$COMROOTgfs/enkfgdas.$PDYa/${cyca}/ensstat/analysis/atmos
datens=$COMROOTgfs/enkfgdas.$PDYg/${cycg}


# Set up $tmpdir
Expand Down Expand Up @@ -166,19 +163,19 @@ nfhrs=`echo $IAUFHRS_ENKF | sed 's/,/ /g'`
for fhr in $nfhrs; do
for imem in $(seq 1 $NMEM_ENKF); do
memchar="mem"$(printf %03i $imem)
$nln $datens/$memchar/${prefix_ens}.atmf00${fhr}.nc sfg_${global_adate}_fhr0${fhr}_${memchar}
$nln $datens/$memchar/model_data/atmos/history/${prefix_ens}.atmf00${fhr}.nc sfg_${global_adate}_fhr0${fhr}_${memchar}
if [ $cnvw_option = ".true." ]; then
$nln $datens/$memchar/${prefix_ens}sfcf00${fhr}.nc sfgsfc_${global_adate}_fhr0${fhr}_${memchar}
$nln $datens/$memchar/model_data/atmos/history/${prefix_ens}sfcf00${fhr}.nc sfgsfc_${global_adate}_fhr0${fhr}_${memchar}
fi
(( imem = $imem + 1 ))
done
$nln $datens/${prefix_ens}.atmf00${fhr}.ensmean.nc sfg_${global_adate}_fhr0${fhr}_ensmean
$nln $datens/ensstat/model_data/atmos/history/${prefix_ens}.atmf00${fhr}.ensmean.nc sfg_${global_adate}_fhr0${fhr}_ensmean
if [ $cnvw_option = ".true." ]; then
$nln $datens/${prefix_ens}.sfcf00${fhr}.ensmean.nc sfgsfc_${global_adate}_fhr0${fhr}_ensmean
fi
done

$nln $datobs/${prefix_obs}.abias_int ./satbias_in
$nln $datobs/${prefix_obs}.abias_int.ensmean ./satbias_in


cd $tmpdir
Expand Down
2 changes: 2 additions & 0 deletions regression/regression_driver.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -36,10 +36,12 @@ for jn in `seq ${RSTART} ${REND}`; do
export scripts=${scripts_updat:-$scripts}
export fixgsi=${fixgsi_updat:-$fixgsi}
export modulefiles=${modulefiles_updat:-$modulefiles}
export ush=${ush_update:-$ush}
else
export scripts=${scripts_contrl:-$scripts}
export fixgsi=${fixgsi_contrl:-$fixgsi}
export modulefiles=${modulefiles_contrl:-$modulefiles}
export ush=${ush_cntrl:-$ush}
fi
rm -f ${job[$jn]}.out

Expand Down
2 changes: 1 addition & 1 deletion regression/regression_namelists.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -181,7 +181,7 @@ OBS_INPUT::
sstviirs viirs-m j1 viirs-m_j1 0.0 4 0
abibufr abi g18 abi_g18 0.0 1 0
ahibufr ahi himawari9 ahi_himawari9 0.0 1 0
atmsbufr atms n21 atms_n21 0.0 1 1
atmsbufr atms n21 atms_n21 0.0 1 0
crisfsbufr cris-fsr n21 cris-fsr_n21 0.0 1 0
sstviirs viirs-m j2 viirs-m_j2 0.0 4 0
ompsnpbufr ompsnp n21 ompsnp_n21 0.0 0 0
Expand Down
48 changes: 24 additions & 24 deletions regression/regression_param.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -87,17 +87,17 @@ case $regtest in
rrfs_3denvar_glbens)

if [[ "$machine" = "Hera" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/1"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/1"
elif [[ "$machine" = "Orion" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/2"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/2"
elif [[ "$machine" = "Hercules" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/2"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/2"
elif [[ "$machine" = "Jet" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/1"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/1"
elif [[ "$machine" = "Gaea" ]]; then
topts[1]="0:15:00" ; popts[1]="18/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="18/2/" ; ropts[2]="/1"
Expand All @@ -117,17 +117,17 @@ case $regtest in
hafs_3denvar_hybens)

if [[ "$machine" = "Hera" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/1"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/1"
elif [[ "$machine" = "Orion" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/2"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/2"
elif [[ "$machine" = "Hercules" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/2"
topts[2]="0:15:00" ; popts[2]="5/8/" ; ropts[2]="/2"
elif [[ "$machine" = "Jet" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/1"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/1"
elif [[ "$machine" = "Gaea" ]]; then
topts[1]="0:15:00" ; popts[1]="18/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="18/2/" ; ropts[2]="/1"
Expand All @@ -146,17 +146,17 @@ case $regtest in

hafs_4denvar_glbens)
if [[ "$machine" = "Hera" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/1"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/1"
elif [[ "$machine" = "Orion" ]]; then
topts[1]="0:20:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[1]="0:20:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:20:00" ; popts[2]="20/2/" ; ropts[2]="/1"
elif [[ "$machine" = "Hercules" ]]; then
topts[1]="0:20:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:20:00" ; popts[2]="20/2/" ; ropts[2]="/1"
topts[1]="0:20:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:20:00" ; popts[2]="10/4/" ; ropts[2]="/1"
elif [[ "$machine" = "Jet" ]]; then
topts[1]="0:15:00" ; popts[1]="20/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="20/2/" ; ropts[2]="/1"
topts[1]="0:15:00" ; popts[1]="5/4/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="10/4/" ; ropts[2]="/1"
elif [[ "$machine" = "Gaea" ]]; then
topts[1]="0:15:00" ; popts[1]="18/1/" ; ropts[1]="/1"
topts[2]="0:15:00" ; popts[2]="18/2/" ; ropts[2]="/1"
Expand Down Expand Up @@ -300,10 +300,10 @@ if [[ "$machine" = "Hera" ]]; then
export APRUN="srun"
elif [[ "$machine" = "Orion" ]]; then
export OMP_STACKSIZE=2048M
export APRUN="srun -n \$ntasks --cpus-per-task=\$threads"
export APRUN="srun -n \$ntasks --mem=0 --cpus-per-task=\$threads"
elif [[ "$machine" = "Hercules" ]]; then
export OMP_STACKSIZE=2048M
export APRUN="srun -n \$ntasks --cpus-per-task=\$threads"
export APRUN="srun -n \$ntasks --mem=0 --cpus-per-task=\$threads"
elif [[ "$machine" = "Jet" ]]; then
export OMP_STACKSIZE=1024M
export MPI_BUFS_PER_PROC=256
Expand Down
2 changes: 1 addition & 1 deletion regression/regression_var.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -190,7 +190,7 @@ export savdir="$ptmp"
export JCAP="62"

# Case Study analysis dates
export global_adate="2022110900"
export global_adate="2024022300"
export rtma_adate="2020022420"
export fv3_netcdf_adate="2017030100"
export rrfs_3denvar_glbens_adate="2021072518"
Expand Down
16 changes: 10 additions & 6 deletions src/enkf/controlvec.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -191,7 +191,7 @@ subroutine read_control()
! read ensemble members on IO tasks
implicit none
real(r_double) :: t1,t2
integer(i_kind) :: nb,ne
integer(i_kind) :: nb,nlev,ne
integer(i_kind) :: q_ind
integer(i_kind) :: ierr

Expand All @@ -218,19 +218,23 @@ subroutine read_control()
if (nproc == 0) t1 = mpi_wtime()
call readgriddata_pnc(cvars3d,cvars2d,nc3d,nc2d,clevels,ncdim,nbackgrounds, &
fgfileprefixes,fgsfcfileprefixes,reducedgrid,grdin,qsat)
if (nproc == 0) then
t2 = mpi_wtime()
print *,'time in readgrid_pnc on root',t2-t1,'secs'
end if
end if
if (nproc <= ntasks_io-1) then
if (.not. paranc) then
if (nproc == 0) t1 = mpi_wtime()
call readgriddata(nanal1(nproc),nanal2(nproc),cvars3d,cvars2d,nc3d,nc2d,clevels,ncdim,nbackgrounds, &
fgfileprefixes,fgsfcfileprefixes,reducedgrid,grdin,qsat)
if (nproc == 0) then
t2 = mpi_wtime()
print *,'time in readgrid on root',t2-t1,'secs'
end if
end if
!print *,'min/max qsat',nanal,'=',minval(qsat),maxval(qsat)
q_ind = getindex(cvars3d, 'q')
if (nproc == 0) then
t2 = mpi_wtime()
print *,'time in readgridata on root',t2-t1,'secs'
end if
if (pseudo_rh .and. q_ind > 0) then
do ne=1,nanals_per_iotask
do nb=1,nbackgrounds
Expand Down Expand Up @@ -357,7 +361,7 @@ subroutine write_control(no_inflate_flag)
endif
deallocate(grdin_mean)
t2 = mpi_wtime()
print *,'time in write_control on root',t2-t1,'secs'
print *,'time in write_control paranc on root',t2-t1,'secs'
endif
end if

Expand Down
Loading

0 comments on commit d6a4b09

Please sign in to comment.