[Dart-dev] svn:mime-type/ Revision: 11876

dart at ucar.edu dart at ucar.edu
Thu Aug 3 15:33:47 MDT 2017


thoar at ucar.edu
2017-08-03 15:33:46 -0600 (Thu, 03 Aug 2017)
387
Added ability to geographically subset at the 'makedaily.sh' stage.
Improved error handling in the python scripts thanks to Romain Escudier.
Avoid recreating existing AVISO files, rename the output files only
if the output file was succesfully created.

The run_convert_aviso.csh script had an implicit need for an existing
input.nml, that need is now explicit.

Documentation updated.




Modified: DART/trunk/observations/AVISO/AVISO.html
===================================================================
--- DART/trunk/observations/AVISO/AVISO.html	2017-08-03 21:13:36 UTC (rev 11875)
+++ DART/trunk/observations/AVISO/AVISO.html	2017-08-03 21:33:46 UTC (rev 11876)
@@ -85,6 +85,18 @@
 have the output file have the name <em class=file>obs_seq.aviso</em>. 
 As such, there is no input namelist specifically for these parameters, but other 
 DART modules still require run-time crontrol specified by <em class=file>input.nml</em>.
+<br />
+<br />
+Fred wrote another python script (<em class=program>shell_scripts/convert_aviso_2.py</em>) 
+to repeatedly call <em class=program>convert_aviso</em> which was subsequently modified by
+Romain Escudier of Rutgers to be able to select which year you want to convert with a little
+more error-checking. Both <em class=program>convert_aviso.py</em> and 
+<em class=program>convert_aviso_2.py</em> must be edited to reflect the location of the data
+files on your system.
+<br />
+<br />
+DART provides <em class=program>shell_scripts/run_convert_aviso.csh</em>
+to prepare the input required and run <em class=program>convert_aviso_2.py</em> as a batch job.
 </P>
 
 <P>
@@ -143,15 +155,28 @@
 </tr>
 <tr><td valign=top><em class=program>shell_scripts/convert_aviso.py</em></td>
     <td valign=top>python script to convert a series of input files and datestamp
-        the output files. 
+        the output files.  This script works with a single satellite and issues
+        a warning if a particular date has no input data.
         </td>
 </tr>
+<tr><td valign=top><em class=program>shell_scripts/convert_aviso_2.py</em></td>
+    <td valign=top>python script to convert a series of input files and datestamp
+        the output files. This version works on many satellite platforms.
+        </td>
+</tr>
+<tr><td valign=top><em class=program>shell_scripts/run_convert_aviso.csh</em></td>
+    <td valign=top>shell script that can be run interactively or as a batch job
+        to run <em class=program>convert_aviso_2.py</em>.
+        </td>
+</tr>
 <tr><td valign=top><em class=program>shell_scripts/makedaily.sh</em></td>
     <td valign=top>shell script to repeatedly call <em class=program>obs_sequence_tool</em>
         to consolidate multiple observation sequence files into an observation
         sequence file that has ALL the observations from ALL platforms in a
-        single file. <em class=program>makedaily.sh</em> is capable of looping
-        over time ranges and creating observation sequences for each time range.
+        single file, as well as any other observations from pre-existing observation
+        sequence files (from the WOD, for example). <em class=program>makedaily.sh</em> 
+        is capable of looping over time ranges and creating observation sequences for 
+        each time range. This is also the place to do any geographic subsetting.
         </td>
 </tr>
 </table>

Modified: DART/trunk/observations/AVISO/convert_aviso.f90
===================================================================
--- DART/trunk/observations/AVISO/convert_aviso.f90	2017-08-03 21:13:36 UTC (rev 11875)
+++ DART/trunk/observations/AVISO/convert_aviso.f90	2017-08-03 21:33:46 UTC (rev 11876)
@@ -105,16 +105,16 @@
 
 call initialize_utilities('convert_aviso')
 
-! command line argument
+! read command line argument(s) - the second is optional.
 call getarg(1, input_file)
 call getarg(2, RE_file)
 
-! Read the DART namelist
+! Read the namelist.
 call find_namelist_in_file('input.nml', 'convert_aviso_nml', iunit)
 read(iunit, nml = convert_aviso_nml, iostat = io)
 call check_namelist_read(iunit, io, 'convert_aviso_nml')
 
-! Record the namelist values used for the run
+! Record the namelist values used for the run.
 call error_handler(E_MSG,'convert_aviso:','convert_aviso_nml values are',' ',' ',' ')
 if (do_output()) write(logfileunit, nml=convert_aviso_nml)
 if (do_output()) write(     *     , nml=convert_aviso_nml)

Added: DART/trunk/observations/AVISO/convert_aviso.nml
===================================================================
--- DART/trunk/observations/AVISO/convert_aviso.nml	                        (rev 0)
+++ DART/trunk/observations/AVISO/convert_aviso.nml	2017-08-03 21:33:46 UTC (rev 11876)
@@ -0,0 +1,4 @@
+&convert_aviso_nml
+   observation_error_std = 0.03  
+  /
+


Property changes on: DART/trunk/observations/AVISO/convert_aviso.nml
___________________________________________________________________
Added: svn:eol-style
## -0,0 +1 ##
+native
\ No newline at end of property
Added: svn:mime-type
## -0,0 +1 ##
+text/plain


More information about the Dart-dev mailing list