[Dart-dev] DART/branches Revision: 12910

[dart at ucar.edu]

thoar at ucar.edu
2018-10-24 10:40:07 -0600 (Wed, 24 Oct 2018)
163
The forcing files have had the outliers removed; CLM ran through the end of 2010.
The CleanForcing.csh script was used for testing but not for the final cleanse.




Modified: DART/branches/cesm_clm/assimilation_code/programs/forcing_check/shell_scripts/CleanForcing.csh
===================================================================
--- DART/branches/cesm_clm/assimilation_code/programs/forcing_check/shell_scripts/CleanForcing.csh	2018-10-24 15:37:44 UTC (rev 12909)
+++ DART/branches/cesm_clm/assimilation_code/programs/forcing_check/shell_scripts/CleanForcing.csh	2018-10-24 16:40:07 UTC (rev 12910)
@@ -63,8 +63,9 @@
 
 # set DSSBASE = /glade/p/rda/data/ds199.1/CAM_DATM.cpl
 set OURBASE = /glade/p/image/thoar/CAM_DATM/4xdaily/CAM_DATM.cpl
+set OURBASE = /glade/p/image/thoar/ds199.1/CAM_DATM.cpl
 
-@ YEAR = 2006
+@ YEAR = 1997
 while ($YEAR <= 2010)
 
    @ MEMBER = 1

Modified: DART/branches/cesm_clm/models/clm/clean_forcing.f90
===================================================================
--- DART/branches/cesm_clm/models/clm/clean_forcing.f90	2018-10-24 15:37:44 UTC (rev 12909)
+++ DART/branches/cesm_clm/models/clm/clean_forcing.f90	2018-10-24 16:40:07 UTC (rev 12910)
@@ -38,9 +38,12 @@
 !           a2x6h_Faxa_lwdn:internal_dname = "a2x6h" ;
 !           a2x6h_Faxa_lwdn:cell_methods = "time: mean" ;
 !
+! swvdf timestep  352 max variance 25982340.000 144  76 max/median    24584.617 144  76
+! member 64 is 45594.96 ... others are x.y
+!
 !-------------------------------------------------------------------------------
 
-use            types_mod, only : r8,i8
+use            types_mod, only : r8, i8, MISSING_R8
 
 use             sort_mod, only : index_sort
 
@@ -52,6 +55,7 @@
                                  finalize_utilities, &
                                  find_namelist_in_file, &
                                  check_namelist_read, &
+                                 open_file, close_file, &
                                  nmlfileunit, &
                                  do_nml_file, &
                                  do_nml_term, &
@@ -102,31 +106,29 @@
 character(len=256) :: input_file(ensemble_size)
 
 real(r8), allocatable :: tensor(:,:,:) ! nx,ny,ensemble_size
-real(r8) :: minvalue, maxvalue, q1, q2, q3, iqr
-real(r8) :: mean, stddev, noise
-real(r8) :: newhigh, newlow, variance, new_variance
-logical  :: suspicious, newmin, newmax
+
 type(random_seq_type) :: r
 
-integer :: imember, itime, iy, ix, nT, ny, nx
+integer :: imember, nT, ny, nx
 integer :: ncstart(3)
 integer :: nccount(3)
 integer :: numdims, dimlens(NF90_MAX_VAR_DIMS)
-integer :: varid
+integer :: ivar, varid
 
 integer            :: iunit, io
 character(len=512) :: string1, string2, string3
 
-integer(i8) :: num_outliers
-
 !-------------------------------------------------------------------------------
 ! namelist
 !-------------------------------------------------------------------------------
 
-integer :: year = 2008
+integer  :: year = 2008
 real(r8) :: criterion = 100.0_r8
+real(r8) :: lat = MISSING_R8
+real(r8) :: lon = MISSING_R8
+integer  :: varnum = 1
 
-namelist /clean_forcing_nml/ year, criterion
+namelist /clean_forcing_nml/ year, criterion, lat, lon, varnum
 
 !===============================================================================
 
@@ -151,18 +153,107 @@
    write(*,*)'Opened '//trim(input_file(imember))
 enddo
 
-varname = variables(1)
+do ivar = 1,5
 
-call nc_get_variable_num_dimensions(ncid(1), varname, numdims)
-call nc_get_variable_size(ncid(1), varname, dimlens)
+   varname = variables(ivar)
 
-if ( .true. ) write(*,*)'dimlens are ',dimlens(1:numdims)
-nx = dimlens(1)
-ny = dimlens(2)
-nT = dimlens(3)
+   call nc_get_variable_num_dimensions(ncid(1), varname, numdims)
+   call nc_get_variable_size(ncid(1), varname, dimlens)
 
-allocate(tensor(ensemble_size,nx,ny))