[Dart-dev] [4055] DART/trunk/models/POP: The advance_model. csh script now enforces that the 'advance_time_present' argument
nancy at ucar.edu
nancy at ucar.edu
Fri Sep 18 15:07:51 MDT 2009
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20090918/8519d604/attachment.html
-------------- next part --------------
Modified: DART/trunk/models/POP/dart_pop_mod.f90
===================================================================
--- DART/trunk/models/POP/dart_pop_mod.f90 2009-09-17 15:09:17 UTC (rev 4054)
+++ DART/trunk/models/POP/dart_pop_mod.f90 2009-09-18 21:07:51 UTC (rev 4055)
@@ -458,7 +458,7 @@
INQUIRE(iolength=reclength) ULAT
open(grid_unit, file=trim(horiz_grid_file), form='unformatted', &
- access='direct', recl=reclength, status='old' )
+ access='direct', recl=reclength, status='old', action='read' )
read(grid_unit, rec=1) ULAT
read(grid_unit, rec=2) ULON
!read(grid_unit, rec=3) HTN
@@ -628,7 +628,7 @@
INQUIRE(iolength=reclength) KMT
open( topo_unit, file=trim(topography_file), form='unformatted', &
- access='direct', recl=reclength, status='old' )
+ access='direct', recl=reclength, status='old', action='read' )
read( topo_unit, rec=1) KMT
close(topo_unit)
Modified: DART/trunk/models/POP/model_mod.f90
===================================================================
--- DART/trunk/models/POP/model_mod.f90 2009-09-17 15:09:17 UTC (rev 4054)
+++ DART/trunk/models/POP/model_mod.f90 2009-09-18 21:07:51 UTC (rev 4055)
@@ -347,7 +347,6 @@
do i = 1, nx
if(ulon(i, 1) /= ulon(i, ny)) then
dipole_grid = .true.
- write(*, *) 'dipole grid is ', dipole_grid
call init_dipole_interp()
return
endif
Modified: DART/trunk/models/POP/shell_scripts/MakeInitialEnsemble.csh
===================================================================
--- DART/trunk/models/POP/shell_scripts/MakeInitialEnsemble.csh 2009-09-17 15:09:17 UTC (rev 4054)
+++ DART/trunk/models/POP/shell_scripts/MakeInitialEnsemble.csh 2009-09-18 21:07:51 UTC (rev 4055)
@@ -33,6 +33,7 @@
# The call to restart_file_tool changes the time tag in the file
# to be whatever is in the namelist. In this case, we want the time
# tag to be 1996 01 01 00Z ... in DART-speak ... 144270 days 0 seconds.
+# tag to be 2000 01 01 00Z ... in DART-speak ... 145731 days 0 seconds.
#
# TJH - 28 July 2008
@@ -100,7 +101,7 @@
cp ${DARTNML} .
cp ${POPDIR}/pop_in .
cp ${POPDIR}/horiz_grid.gx3v5.r8ieee.le .
-cp ${POPDIR}/topography.gx3v5.r8ieee.le .
+cp ${POPDIR}/topography.gx3v5.i4ieee.le .
cp ${POPDIR}/vert_grid.gx3v5 .
#-------------------------------------------------------------------------
@@ -140,7 +141,7 @@
# single_restart_file_out = .true.,
# write_binary_restart_files = .true.,
# overwrite_data_time = .true.,
-# new_data_days = 144270,
+# new_data_days = 145731,
# new_data_secs = 0,
# input_is_model_advance_file = .false.,
# output_is_model_advance_file = .false.,
@@ -157,7 +158,7 @@
echo '/overwrite_data_time' >> vi_script
echo ':s/.false./.true./' >> vi_script
echo '/new_data_days' >> vi_script
-echo ':s/-1/144270/' >> vi_script
+echo ':s/-1/145731/' >> vi_script
echo '/new_data_secs' >> vi_script
echo ':s/-1/0/' >> vi_script
echo ':wq' >> vi_script
@@ -170,7 +171,7 @@
# Loop over all the restart files.
#-------------------------------------------------------------------------
-@ memcount = 0
+@ member = 0
foreach FILE ( ${RESTARTDIR}/*/cx3.dart.001.pop.r.*.nc )
@@ -180,8 +181,8 @@
ln -sf $FILE pop.r.nc
- echo "pop.r.nc" >! pop_pointer.restart
- echo "RESTART_FMT=nc" >> pop_pointer.restart
+ echo "pop.r.nc" >! rpointer.ocn.restart
+ echo "RESTART_FMT=nc" >> rpointer.ocn.restart
echo "Converting restart $FILE ..."
@@ -211,13 +212,13 @@
# we could probably use just one, but we've got them, and we're just
# creating links ... it's easy enough.
#----------------------------------------------------------------------
- @ memcount = $memcount + 1
- set OFNAME = `printf ens_mem.%03d $memcount`
+ @ member = $member + 1
+ set OFNAME = `printf ens_mem.%03d $member`
mv filter_updated_restart $OFNAME
- echo $FILE:t >! rpointer.ocn.${memcount}.restart
- echo "RESTART_FMT=nc" >> rpointer.ocn.${memcount}.restart
- ln -s $FILE $FILE:t
+ echo $FILE:t >! rpointer.ocn.${member}.restart
+ echo "RESTART_FMT=nc" >> rpointer.ocn.${member}.restart
+ ln -sf $FILE $FILE:t
#----------------------------------------------------------------------
# remove input files to prep for next iteration
Modified: DART/trunk/models/POP/shell_scripts/advance_model.csh
===================================================================
--- DART/trunk/models/POP/shell_scripts/advance_model.csh 2009-09-17 15:09:17 UTC (rev 4054)
+++ DART/trunk/models/POP/shell_scripts/advance_model.csh 2009-09-18 21:07:51 UTC (rev 4055)
@@ -52,6 +52,18 @@
cp -pv $FILE . || exit 1
end
+# Try to ensure that the input.nml has the required value for
+# dart_to_pop_nml:advance_time_present for this context.
+
+echo '1' >! ex_commands
+echo '/dart_to_pop_nml' >> ex_commands
+echo '/advance_time_present' >> ex_commands
+echo ':s/\.false\./\.true\./' >> ex_commands
+echo ':wq' >> ex_commands
+
+( ex input.nml < ex_commands )
+\rm -f ex_commands
+
# copy the files used by
foreach FILE ( ../horiz_grid.gx3v5.* \
../topography.gx3v5.* \
Modified: DART/trunk/models/POP/shell_scripts/run_perfect_model_obs.batch
===================================================================
--- DART/trunk/models/POP/shell_scripts/run_perfect_model_obs.batch 2009-09-17 15:09:17 UTC (rev 4054)
+++ DART/trunk/models/POP/shell_scripts/run_perfect_model_obs.batch 2009-09-18 21:07:51 UTC (rev 4055)
@@ -163,7 +163,7 @@
${COPY} ${POPDIR}/vert_grid.gx3v5 .
${COPY} ${POPDIR}/horiz_grid.gx3v5.r8ieee.le .
- ${COPY} ${POPDIR}/topography.gx3v5.r8ieee.le .
+ ${COPY} ${POPDIR}/topography.gx3v5.i4ieee.le .
#${COPY} ${POPDIR}/chl_mm_SeaWiFs97-01_20031205.ieeer8 .
#${COPY} ${POPDIR}/sfwf_20040517.ieeer8 .
More information about the Dart-dev
mailing list