[Dart-dev] DART/branches Revision: 11700
dart at ucar.edu
dart at ucar.edu
Tue Jun 6 19:36:04 MDT 2017
thoar at ucar.edu
2017-06-06 19:36:03 -0600 (Tue, 06 Jun 2017)
432
Adding observation converter support for CICE.
Most of these read a CICE restart file that has had noise added to it.
That happens offline, and then those perturbed restart files are run
through the different converters. Since the observations are
from every gridpoint, Yongfei found it easier to do this than generating
an obs_seq.in for each gridpoint.
These converters were developed by Yonfei Zhang - University of Washington.
Modified: DART/branches/rma_trunk/assimilation_code/modules/observations/DEFAULT_obs_kind_mod.F90
===================================================================
--- DART/branches/rma_trunk/assimilation_code/modules/observations/DEFAULT_obs_kind_mod.F90 2017-06-05 19:29:30 UTC (rev 11699)
+++ DART/branches/rma_trunk/assimilation_code/modules/observations/DEFAULT_obs_kind_mod.F90 2017-06-07 01:36:03 UTC (rev 11700)
@@ -435,19 +435,20 @@
QTY_SEAICE_SNOWENTHALPY003 = 352, &
QTY_SOM_TEMPERATURE = 353, &
QTY_SEAICE_FY = 354, &
- QTY_SEAICE_AGREG_FY = 355
+ QTY_SEAICE_AGREG_FY = 355, &
+ QTY_SEAICE_AGREG_SURFACETEMP = 356
integer, parameter, public :: &
- QTY_SEA_SURFACE_ANOMALY = 356
+ QTY_SEA_SURFACE_ANOMALY = 360
integer, parameter, public :: &
- QTY_LARGE_SCALE_STATE = 357, &
- QTY_SMALL_SCALE_STATE = 358
+ QTY_LARGE_SCALE_STATE = 361, &
+ QTY_SMALL_SCALE_STATE = 362
! max_defined_quantities is private to this module. see comment below near the max_obs_specific
! declaration for more info about publics and private values.
-integer, parameter :: max_defined_quantities = 358
+integer, parameter :: max_defined_quantities = 362
!----------------------------------------------------------------------------
! This list is autogenerated by the 'preprocess' program. To add new
@@ -846,10 +847,12 @@
obs_kind_names(353) = obs_kind_type(QTY_SOM_TEMPERATURE, 'QTY_SOM_TEMPERATURE')
obs_kind_names(354) = obs_kind_type(QTY_SEAICE_FY, 'QTY_SEAICE_FY')
obs_kind_names(355) = obs_kind_type(QTY_SEAICE_AGREG_FY, 'QTY_SEAICE_AGREG_FY')
-obs_kind_names(356) = obs_kind_type(QTY_SEA_SURFACE_ANOMALY, 'QTY_SEA_SURFACE_ANOMALY')
-obs_kind_names(357) = obs_kind_type(QTY_LARGE_SCALE_STATE, 'QTY_LARGE_SCALE_STATE')
-obs_kind_names(358) = obs_kind_type(QTY_SMALL_SCALE_STATE, 'QTY_SMALL_SCALE_STATE')
+obs_kind_names(356) = obs_kind_type(QTY_SEAICE_AGREG_SURFACETEMP,'QTY_SEAICE_AGREG_SURFACETEMP')
+obs_kind_names(360) = obs_kind_type(QTY_SEA_SURFACE_ANOMALY, 'QTY_SEA_SURFACE_ANOMALY')
+obs_kind_names(361) = obs_kind_type(QTY_LARGE_SCALE_STATE, 'QTY_LARGE_SCALE_STATE')
+obs_kind_names(362) = obs_kind_type(QTY_SMALL_SCALE_STATE, 'QTY_SMALL_SCALE_STATE')
+
! count here, then output below
num_kind_assimilate = 0
Modified: DART/branches/rma_trunk/observations/forward_operators/obs_def_cice_mod.f90
===================================================================
--- DART/branches/rma_trunk/observations/forward_operators/obs_def_cice_mod.f90 2017-06-05 19:29:30 UTC (rev 11699)
+++ DART/branches/rma_trunk/observations/forward_operators/obs_def_cice_mod.f90 2017-06-07 01:36:03 UTC (rev 11700)
@@ -23,6 +23,7 @@
!SAT_SEAICE_VOLUME, QTY_SEAICE_VOLUME, COMMON_CODE
!SAT_SEAICE_SNOWVOLUME, QTY_SEAICE_SNOWVOLUME, COMMON_CODE
!SAT_SEAICE_AGREG_FY, QTY_SEAICE_AGREG_FY, COMMON_CODE
+!SAT_SEAICE_AGREG_SURFACETEMP, QTY_SEAICE_AGREG_SURFACETEMP, COMMON_CODE
! END DART PREPROCESS KIND LIST
! <next few lines under version control, do not edit>
Modified: DART/branches/rma_trunk/observations/obs_converters/cice/cice_to_obs.f90
===================================================================
--- DART/branches/rma_trunk/observations/obs_converters/cice/cice_to_obs.f90 2017-06-05 19:29:30 UTC (rev 11699)
+++ DART/branches/rma_trunk/observations/obs_converters/cice/cice_to_obs.f90 2017-06-07 01:36:03 UTC (rev 11700)
@@ -2,37 +2,31 @@
! by UCAR, "as is", without charge, subject to all terms of use at
! http://www.image.ucar.edu/DAReS/DART/DART_download
!
-! $Id$
+! DART $Id$
-!> convert snow and ice data center files from binary to obs_seq format.
+!> cice_to_obs - read binary files from the national snow and ice data
+!> center and convert percentages into sea ice concentration
+!> observation sequence files.
!>
-!> ;============data description from CC matlab script=================================
-!> ; Data are stored as flat two-byte integers representing sea ice concentration
-!> ; values. The sea ice concentration data values are packed into integer
-!> ; format by multiplying the original sea ice concentration values by
-!> ; 10. These values range from 0 to 1,000, with land registered as -800 and
-!> ; -100 being the northern hemisphere hole poleward of 85 degrees for the
-!> ; SMMR data and poleward of 87 degrees for the SSM/I data, which the
-!> ; satellite cannot cover due to orbit inclination.
-!> ;
-!> ; (nsc -- what are these? ...)
-!> ; little-endian byte order no mention of header
-!> ; One-byte (scaled, unsigned integer) flat binary arrays preceded by a 300-byte header;
-!> ; *uint8
-!> ;=====================================================================================
-!>
+!> nsc 28 june 2016 based on conversion programs from Yongfei Zhang,
+!> now at University of Washington working with Cecelia Bitz.
program cice_to_obs
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!
-! cice_to_obs - read binary files from the national snow and ice data
-! center and convert percentages into sea ice concentration observations.
-!
More information about the Dart-dev
mailing list