Procedures

Error handling for netCDF commands

Return number of observations from file

Return data assimilation interval from file

Deallocation of observation arrays

Deallocation of observation index arrays

CLM3.5 Initialization routine for TSMP-PDAF

TSMP PDAF comment out beginning

Set indices of grid cells with lon/lat smaller than observation locations

Set indices of grid cells with lon/lat smaller than observation locations

Read observation index arrays for C-code

Return the filename for output of CLM statistics

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 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Read NetCDF observation file