ESCOMP / SimpleLand

Simple Land Model for CESM --- *** IN DEVELOPMENT *** --- please contact for more info. See supplemental information of https://journals.ametsoc.org/doi/abs/10.1175/JCLI-D-18-0812.1 for a description of SLIM physics. Implementation of SLIM into the main CESM trunk is ongoing. SLIM currently works with the CESM2.1 release, but must be downloaded from this repository until we finish implementing it properly into the main CESM code.
Other
12 stars 7 forks source link

Data assimilation pause-resume needs some updates from release-clm5.0.35-ctsm5.1.dev111 #47

Closed ekluzek closed 1 year ago

ekluzek commented 1 year ago

I don't think this is a priority to do, especially right now.

But, there were some changes/fixes to the data assimilation pause-resume mode for CLM that would need to go into SLIM to get it working.

Here's some example changes in lnd_comp_mct.F90 (because of directory changes the best way to do the comparison is to compare the release version in src/cpl to the latest in src/cpl/mct.

         ! Fill in infodata settings
         call seq_infodata_PutData(infodata, lnd_prognostic=.true.)
         call seq_infodata_PutData(infodata, lnd_nx=ldomain%ni, lnd_ny=ldomain%nj)
!       call lnd_handle_resume( cdata_l )
--- 607,630 ----
    end subroutine lnd_domain_mct

    !====================================================================================
!   subroutine lnd_handle_resume( cdata_l )
      !
      ! !DESCRIPTION:
      ! Handle resume signals for Data Assimilation (DA)
      !
      ! !USES:
      use clm_time_manager , only : update_DA_nstep
!     use seq_cdata_mod    , only : seq_cdata, seq_cdata_setptrs
      implicit none
      ! !ARGUMENTS:
!     type(seq_cdata),            intent(inout) :: cdata_l          ! Input land-model driver data
      ! !LOCAL VARIABLES:
      logical :: resume_from_data_assim                      ! flag if we are resuming after data assimulation was done
      !---------------------------------------------------------------------------

!     ! Check to see if restart was modified and we are resuming from data
!     ! assimilation
!     call seq_cdata_setptrs(cdata_l, post_assimilation=resume_from_data_assim)
      if ( resume_from_data_assim ) call update_DA_nstep()

    end subroutine lnd_handle_resume
ekluzek commented 1 year ago

Closing this as something we aren't going to do.