[Dart-dev] DART/branches Revision: 12550

[dart at ucar.edu]

thoar at ucar.edu
2018-04-20 16:32:30 -0600 (Fri, 20 Apr 2018)
283
Do not remove the from_dart.semaphore in the qsub script.
It is deleted in the Fortran code ... there appears to be
a timing issue.

Now the model is not actually using the time to advance from the semaphore file,
so the DART script is linking to non-existent observation files ...




Modified: DART/branches/openggcm/models/openggcm/shell_scripts/run_filter.qsub
===================================================================
--- DART/branches/openggcm/models/openggcm/shell_scripts/run_filter.qsub	2018-04-20 21:16:58 UTC (rev 12549)
+++ DART/branches/openggcm/models/openggcm/shell_scripts/run_filter.qsub	2018-04-20 22:32:30 UTC (rev 12550)
@@ -238,7 +238,7 @@
       foreach FILE ( preassim_*.nc output_*.nc obs_seq.final )
          set BASE = $FILE:r
          set EXT  = $FILE:e
-         \mv -v $FILE ${BASE}_${ModelTime}.${EXT}
+         \mv -v $FILE ${BASE}.${ModelTime}.${EXT}
       end
 
       #------------------------------------------------------------
@@ -268,17 +268,23 @@
 
       @ ncycle ++  # debug
 
+      echo "600"                          >! temp_file.txt
+      echo "Assimilation for $TimeString" >> temp_file.txt
+
       foreach FILE ( `cat semaphore_files.txt` )
          set NEWFILE = `echo ${FILE} | sed -e "s/to_dart/from_dart/"`
          echo "making ${NEWFILE} to signal openggcm at cycle $ncycle"
-         ${REMOVE} ${FILE} ${NEWFILE}
-         echo "600"                          >! temp_file.txt
-         echo "Assimilation for $TimeString" >> temp_file.txt
+
+         ${REMOVE} ${FILE}
+#         ${REMOVE} ${FILE} ${NEWFILE}
+
          ${COPY} temp_file.txt ${NEWFILE}
          ${COPY} temp_file.txt ${NEWFILE}.safety
 
-         if ( $ncycle > 3 ) echo "DONE" >! ${NEWFILE}
-         if ( $ncycle > 3 ) echo "DONE" >! ${NEWFILE}.safety
+         if ( $ncycle > 3 ) echo "DONE"
+#        if ( $ncycle > 3 ) echo "DONE" >! ${NEWFILE}
+#        if ( $ncycle > 3 ) echo "DONE" >! ${NEWFILE}.safety
+
       end
       if ( $ncycle > 3 ) then
          echo "${JOBNAME} ($JOBID) finished at "`date`
@@ -285,6 +291,8 @@
          exit
       endif
 
+      ${REMOVE} semaphore_files.txt
+
    endif
 
    sleep 1

Modified: DART/branches/openggcm/models/openggcm/work/input.nml
===================================================================
--- DART/branches/openggcm/models/openggcm/work/input.nml	2018-04-20 21:16:58 UTC (rev 12549)
+++ DART/branches/openggcm/models/openggcm/work/input.nml	2018-04-20 22:32:30 UTC (rev 12550)
@@ -38,16 +38,16 @@
    last_obs_days            = -1
    last_obs_seconds         = -1
 
-   stages_to_write          = 'preassim', 'analysis'
-   output_members           = .false.
+   stages_to_write          = 'preassim', 'output'
+   output_members           = .true.
    output_mean              = .true.
    output_sd                = .false.
-   write_all_stages_at_end  = .true.
+   write_all_stages_at_end  = .false.
    output_state_file_list   = "output_list.txt"
 
    inf_flavor                  = 2,             0
-   inf_initial_from_restart    = .false.,       .false.
-   inf_sd_initial_from_restart = .false.,       .false.
+   inf_initial_from_restart    = .true.,       .false.
+   inf_sd_initial_from_restart = .true.,       .false.
    inf_initial                 = 1.0,           1.0
    inf_sd_initial              = 0.6,           0.6
    inf_damping                 = 0.9,           0.9
@@ -132,7 +132,7 @@
 
 # maybe later   assimilate_these_obs_types = 'SUPER_DARN_ELECTRIC_POTENTIAL',
 &obs_kind_nml
-   assimilate_these_obs_types = 'GPS_PROFILE'
+   assimilate_these_obs_types = 'COSMIC_ELECTRON_DENSITY'
     evaluate_these_obs_types  = ''
   /
 
@@ -236,15 +236,15 @@
    loc_of_interest       = 54.0, -88.0, 425835.8  ! (40,2,4)   516815258820.758
 &model_mod_check_nml
    verbose               = .true.
-   test1thru             = 4
+   test1thru             = 5
    x_ind                 = 3951
    loc_of_interest       = 54.0, -88.0, 425835.8
    quantity_of_interest  = 'QTY_ELECTRON_DENSITY'
    interp_test_lonrange  = 0.0, 360.0
-   interp_test_dlon      = 1.0
+   interp_test_dlon      = 0.5
    interp_test_latrange  = -90.0, 90.0
-   interp_test_dlat      = 1.0