[Dart-dev] DART/branches Revision: 11132

[dart at ucar.edu]

mizzi at ucar.edu
2017-02-22 12:55:59 -0700 (Wed, 22 Feb 2017)
34
Update WRF-Chem/DART repository.




Modified: DART/branches/mizzi/observations/MODIS/modis_ascii_to_obs.f90
===================================================================
--- DART/branches/mizzi/observations/MODIS/modis_ascii_to_obs.f90	2017-02-22 18:56:31 UTC (rev 11131)
+++ DART/branches/mizzi/observations/MODIS/modis_ascii_to_obs.f90	2017-02-22 19:55:59 UTC (rev 11132)
@@ -1,11 +1,15 @@
-! DART software - Copyright 2004 - 2013 UCAR. This open source software is
+! DART software - Copyright © 2004 - 2010 UCAR. This open source software is
 ! provided by UCAR, "as is", without charge, subject to all terms of use at
 ! http://www.image.ucar.edu/DAReS/DART/DART_download
-!
-! $Id$
 
 program modis_ascii_to_obs
 
+! <next few lines under version control, do not edit>
+! $URL$
+! $Id$
+! $Revision$
+! $Date$
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 !   modis_ascii_to_obs - a program that only needs minor customization to read
@@ -94,7 +98,7 @@
 
 ! the first one needs to contain the string 'observation' and the
 ! second needs the string 'QC'.
-call set_copy_meta_data(obs_seq, 1, 'observation')
+call set_copy_meta_data(obs_seq, 1, 'MODIS observation')
 call set_qc_meta_data(obs_seq, 1, 'MODIS QC index')
 
 ! if you want to append to existing files (e.g. you have a lot of
@@ -159,12 +163,19 @@
    if (debug) print *, 'next observation located at lat, lon = ', lat, lon
 
    ! check the lat/lon values to see if they are ok
-   if ( lat >  90.0_r8 .or. lat <  -90.0_r8 ) cycle obsloop
+   !if ( lat >  90.0_r8 .or. lat <  -90.0_r8 ) cycle obsloop
+   !if ( lon <   0.0_r8 .or. lon >  360.0_r8 ) cycle obsloop
 
+
    ! if lon comes in between -180 and 180, use these lines instead:
+   if ( lat >  90.0_r8 .or. lat <  -90.0_r8 ) cycle obsloop
    if ( lon > 180.0_r8 .or. lon < -180.0_r8 ) cycle obsloop
    if ( lon < 0.0_r8 )  lon = lon + 360.0_r8 ! changes into 0-360
 
+!
+! APM: skip for single obs cluster
+!   if ( lon < 263.0_r8 .or. lon > 267.0_r8 .or. lat < 38.0_r8 .or. lat > 42.0_r8 ) cycle obsloop
+!
    ! put date into a dart time format
    time_obs = set_date(year, month, day, hour, minute, second)
 
@@ -313,10 +324,3 @@
 end subroutine add_obs_to_seq
 
 end program modis_ascii_to_obs
-
-! <next few lines under version control, do not edit>
-! $URL$
-! $Id$
-! $Revision$
-! $Date$
-

Deleted: DART/branches/mizzi/observations/MODIS/read_modis.m
===================================================================
--- DART/branches/mizzi/observations/MODIS/read_modis.m	2017-02-22 18:56:31 UTC (rev 11131)
+++ DART/branches/mizzi/observations/MODIS/read_modis.m	2017-02-22 19:55:59 UTC (rev 11132)
@@ -1,196 +0,0 @@
-%% read_modis reads MODIS HDF files; process the data and dump the data as ascii file
-% the output of this code is the input of modis_ascii_to_obs in DART
-
-%% DART software - Copyright 2004 - 2013 UCAR. This open source software is
-% provided by UCAR, "as is", without charge, subject to all terms of use at
-% http://www.image.ucar.edu/DAReS/DART/DART_download
-%
-% DART $Id$
-
-   clear all;
-   iday_start=29;
-   iday_end=30;
-   ifile_start=1;
-%
-% where's the data?
-   rootdir='/glade/p/work/mizzi/TRUNK/DART_CHEM/observations/MODIS';
-   rundir='/glade/scratch/mizzi/MODIS_ASCII';
-   date_dir=2008060100:100:2008063000;
-   indx=1;
-%   for indx=1:30
-   for indx=iday_start:iday_end
-      cd(rundir);
-      date_str=int2str(date_dir(indx));
-      datadir=strcat('/glade/p/acd/mizzi/AVE_TEST_DATA/MODIS/',date_str);
-%
-% directory listing
-      list=dir(datadir);
-%
-% cut other files not HDF
-      list=list(3:end);
-%      list(1).name