Procedures

@author Wolfgang Kurtz, Guowei He, Mukund Pondkule @date 03.03.2023 @brief Error handling for netCDF commands @param[in] status netCDF command status @details This subroutine checks the status of a netCDF command and prints an error message if necessary.

@author Wolfgang Kurtz, Guowei He @date 21.03.2022 @brief Return number of observations from file @param[in] fn Filename of the observation file @param[out] nn number of observations in fn @details Reads the content of the variable (!) named no_obs from NetCDF file fn.

@author Yorck Ewerdwalbesloh, Johannes Keller @date 11.09.2023 @brief Return data assimilation interval from file @param[in] fn Filename of the observation file @param[out] aa new da_interval (number of time steps until next assimilation time step) @details Reads the content of the variable name da_interval from NetCDF file fn using subroutines from the NetCDF module. The result is returned in aa.

@author Wolfgang Kurtz, Guowei He, Mukund Pondkule @date 03.03.2023 @brief Deallocation of observation arrays @details This subroutine deallocates the observation arrays used in subroutine read_obs_nc.

@author Wolfgang Kurtz, Guowei He, Mukund Pondkule @date 03.03.2023 @brief Deallocation of observation index arrays @details This subroutine deallocates the observation index arrays used in subroutine get_obsindex_currentobsfile.

@author Wolfgang Kurtz, Guowei He @date 25.02.2022 @brief CLM3.5 Initialization routine for TSMP-PDAF @param[in] finname CLM input file name @details Initialization routines for the land model analogous to clm3_5/src/main/program_off.F90

TSMP PDAF comment out beginning

TSMP PDAF comment out beginning

@author Mukund Pondkule, Johannes Keller @date 27.03.2023 @brief Set indices of grid cells with lon/lat smaller than observation locations @details This routine sets the indices of grid cells with lon/lat smaller than observation locations.

@author Mukund Pondkule, Johannes Keller @date 27.03.2023 @brief Set indices of grid cells with lon/lat smaller than observation locations @details This routine sets the indices of grid cells with lon/lat smaller than observation locations.

@author Mukund Pondkule, Johannes Keller @date 27.03.2023 @brief Set indices of grid cells with lon/lat smaller than observation locations @details This routine sets the indices of grid cells with lon/lat smaller than observation locations.

@author Wolfgang Kurtz, Guowei He, Mukund Pondkule @date 03.03.2023 @brief Read observation index arrays for C-code @param[out] no_obs Number of observations @details This subroutine reads the observation index arrays for usage in the enkf_parflow.c for groundwater masking.

@author Wolfgang Kurtz, Guowei He @date 21.03.2022 @brief Return the filename for output of CLM statistics @return character(len=256) get_statistic_filename Filename @details This functions has likely been adapted from set_hist_filename from histFileMod.F90 in CLM3.5

broadcast crns_flag

broadcast crns_flag

$ ! Due to domain decomposition in our case the mean variance is computed $ ! for the full domain using the function MPI_Allreduce $#ifndef CLMSA $ if (model .eq. tag_model_parflow) then $ meanvar_p = 0 $ sum_p = 0 $ count = 0 $ do i = 1, dim_obs_p $ if(pressure_obserr_p(i) /= 0) then $ sum_p = sum_p + pressure_obserr_p(i) $ count = count + 1 $ endif $ enddo $ ! averaging the sum of observation errors with total no of non-zero observations $ meanvar_p = sum_p/count $ ! summing the average of observation errors and communicating it back to each rank $ call MPI_Allreduce(meanvar_p, meanvar, 1, MPI_REAL8, MPI_SUM, COMM_filter, MPIerr) $ ! to get the mean dividing the mean observation error by size of processors $ meanvar = meanvar/npes_filter $ end if $#endif $ $#if defined CLMSA $ if(model .eq. tag_model_clm) then $ meanvar_p = 0 $ sum_p = 0 $ count = 0 $ do i = 1, dim_obs_p $ if(clm_obserr_p(i) /= 0) then $ sum_p = sum_p + clm_obserr_p(i) $ count = count + 1 $ endif $ enddo $ ! averaging the sum of observation errors with total no of non-zero observations $ meanvar_p = sum_p/count $ ! summing the average of observation errors and communicating it back to each rank $ call MPI_Allreduce(meanvar_p, meanvar, 1, MPI_REAL8, MPI_SUM, COMM_filter, MPIerr) $ ! to get the mean dividing the mean observation error by size of processors $ meanvar = meanvar/npes_filter $ end if $#endif

print *, "stepnow", stepnow

!!!!!!!!!!!!!!!!!!!!!! case below contains dummy CLM component !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

@author Wolfgang Kurtz, Guowei He, Mukund Pondkule @date 03.03.2023 @brief Read NetCDF observation file @param[in] current_observation_filename Name of observation file @details This subroutine reads the observation file and stores the data in the corresponding arrays.