add option not to write anal members

UMD-AOSC · Feb 15, 2024 · 580ec30 · 580ec30
1 parent fa66205
commit 580ec30
Show file tree

Hide file tree

Showing 3 changed files with 16 additions and 12 deletions.
diff --git a/src/letkf/letkf.f90 b/src/letkf/letkf.f90
@@ -228,20 +228,22 @@ PROGRAM letkf
   !----------------------------------------------------------------------------
   ! Write analysis ensemble members
   !----------------------------------------------------------------------------
-  CALL MPI_BARRIER(MPI_COMM_WORLD,ierr)
-  CALL write_ens_mpi('anal',nbv,anal3d,anal2d)
+  if (DO_WRITE_ANAL_MEMBERS) then
+    CALL MPI_BARRIER(MPI_COMM_WORLD,ierr)
+    CALL write_ens_mpi('anal',nbv,anal3d,anal2d)
 
-  !-----------------------------------------------------------------------------
-  ! Check timer for writing analysis ensemble members
-  !-----------------------------------------------------------------------------
-  CALL MPI_BARRIER(MPI_COMM_WORLD,ierr)  !STEVE: to make sure timers are aligned (mostly for debugging)
-  CALL CPU_TIME(rtimer)
-  WRITE(6,'(A,2F10.2)') '### TIMER(WRITE_ENS_MPI):',rtimer,rtimer-rtimer00
-  if (dortout) then !STEVE: force output to file
-    CLOSE(6)
-    OPEN(6,FILE=stdoutf,POSITION='APPEND',STATUS = 'OLD')
+    !-----------------------------------------------------------------------------
+    ! Check timer for writing analysis ensemble members
+    !-----------------------------------------------------------------------------
+    CALL MPI_BARRIER(MPI_COMM_WORLD,ierr)  !STEVE: to make sure timers are aligned (mostly for debugging)
+    CALL CPU_TIME(rtimer)
+    WRITE(6,'(A,2F10.2)') '### TIMER(WRITE_ENS_MPI):',rtimer,rtimer-rtimer00
+    if (dortout) then !STEVE: force output to file
+      CLOSE(6)
+      OPEN(6,FILE=stdoutf,POSITION='APPEND',STATUS = 'OLD')
+    endif
+    rtimer00=rtimer
   endif
-  rtimer00=rtimer
 
   !----------------------------------------------------------------------------
   ! Write analysis ensemble mean and spread for evaluation

diff --git a/src/letkf/params_letkf.f90 b/src/letkf/params_letkf.f90
@@ -29,6 +29,7 @@ MODULE params_letkf
 !-------------------------------------------------------------------------------
 LOGICAL :: DO_WRITE_ENS_MEAN_SPRD = .false.  ! Indicates whether to write the ensemble mean and spread (fcst and analysis)
 LOGICAL :: DO_WRITE_OMB_MEAN = .false. ! whether write out mean[yo-h(x_i)] and other obs diag info
+LOGICAL :: DO_WRITE_ANAL_MEMBERS = .true. ! whether write out analysis into netcdf files.
 
 !-------------------------------------------------------------------------------
 ! From letkf_local.f90

diff --git a/src/model_specific/mom6/input_nml_mom6.f90 b/src/model_specific/mom6/input_nml_mom6.f90
@@ -77,6 +77,7 @@ MODULE input_nml_oceanmodel
                               DO_OBSERR_AOEI, &      ! logical flag to enable adaptive obs error inflation
                               DO_WRITE_ENS_MEAN_SPRD, &  ! logical flag to write the ensemble mean and spread for the forecast and analysis fields
                               DO_WRITE_OMB_MEAN,      &  ! logical flag to write the obs info (omb,qc,...)
+                              DO_WRITE_ANAL_MEMBERS,  &  ! logical flag to write out the analysis into netcdf files (default=.true.)
                               DO_DRIFTERS, &         ! logical flag to do lagrangian drifters assimilation
                               DO_ALTIMETRY, &        ! logical flag to do altimetry data assimilation
                               DO_SLA, &              ! logical flag to use SLA for altimetry