/[MITgcm]/MITgcm_contrib/jmc_script/tst_2+2_cpl
ViewVC logotype

Diff of /MITgcm_contrib/jmc_script/tst_2+2_cpl

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by jmc, Wed Nov 7 01:40:47 2007 UTC revision 1.7 by jmc, Sat Dec 26 23:43:22 2009 UTC
# Line 9  rnp_loc() Line 9  rnp_loc()
9  #   rename files with prefix = 'arg_1' to files with prefix = 'arg_2'  #   rename files with prefix = 'arg_1' to files with prefix = 'arg_2'
10      echo 'rnp_loc:' $1 $2      echo 'rnp_loc:' $1 $2
11  #   rnp -s $1 $2  #   rnp -s $1 $2
12      yy=$1 ; zz=$2      listY=`ls -1 $1*`
13      listF=`ls -1 ${yy}*`      for yy in $listY
     for xx in $listF  
14      do      do
15        ff=`echo $xx | sed "s/^$yy/$zz/"`        zz=`echo $yy | sed "s/^$1/$2/"`
16        mv $xx $ff        mv $yy $zz
17      done      done
18  }  }
19    
20  dif_tiles()  dif_tiles()
21  {  {
22  #   dif_data_loc dir1 dir2 file  #   dif_data_loc dd1 dd2 file
23  #   for all tiles  ==> do diff dir1/file dir2/file'  #   for all tiles  ==> do diff dd1/file dd2/file'
24      dir1=$1 ; dir2=$2 ; zz=$3      dd1=$1 ; dd2=$2 ; zz=$3
25      chkmeta='Y'      chkmeta='Y'
26      listX=`(cd $dir1 ; ls $zz.*.data | sed "s/$zz\.//" | sed "s/\.data//")`      listY=`(cd $dd1 ; ls $zz.*.data | sed "s/$zz\.//" | sed "s/\.data//")`
27      echo '--> file=' $zz ', listX=' $listX      echo '--> file='$zz', listY='$listY
28      for xx in $listX      for yy in $listY
29      do      do
30        echo diff $dir1/$zz.$xx.data $dir2        echo diff $dd1/$zz.$yy.data $dd2
31        diff $dir1/$zz.$xx.data $dir2        diff $dd1/$zz.$yy.data $dd2
32        out=$?        out=$?
33        if test $out != 0        if test $out != 0
34        then echo 'Diff outp=' $out ' ==> stop'        then echo 'Diff outp=' $out ' ==> stop'
35         exit ; fi         exit ; fi
36        if test $chkmeta = 'Y'        if test $chkmeta = 'Y'
37        then        then
38         echo diff $dir1/$zz.$xx.meta $dir3         echo diff $dd1/$zz.$yy.meta $dd2
39         diff $dir1/$zz.$xx.meta $dir3         diff $dd1/$zz.$yy.meta $dd2
40         out=$?         out=$?
41         if test $out != 0         if test $out != 0
42         then echo 'Diff outp=' $out ' ==> stop'         then echo 'Diff outp=' $out ' ==> stop'
# Line 66  else Line 65  else
65   NiAt=0   NiAt=0
66  fi  fi
67    
68  case $# in  #- parse options:
69    1) ;;  CMD='run_cpl_test'
70     *) echo "Usage: `basename $0` flag"  MTH=
71        echo " Check restart of coupled set-up: compare 1 run of 2 x $Nit it long"  argList=$*
72        echo "                       with 2 consecutive runs of $Nit it long each"  for xx in $argList
73        echo " where: flag = 0 -> prepare for 1rst run"  do
74        echo "        flag = 1 -> move res. after 1rst run & prepare for 2nd"   if test $xx = '-mth' ; then MTH=$xx ; shift ; continue ; fi
75        echo "        flag = 2 -> move res. after 2nd  run & prepare for 3rd"   if test -x ${CMD}_$xx ; then CMD="${CMD}_$xx" ; shift ; continue ; fi
76        echo "        flag = 3 -> move res. after 3rd  run"  done
77        echo "        flag = 4 -> compare Ocean results"  
78        echo "        flag = 5 -> compare Atmos results"  if [ $# -ne 1 ]
79        echo "        flag = 6 -> diff pickup files"  then
80        exit ;;     echo "Usage: `basename $0` [opt] flag"
81  esac     echo " Check restart of coupled set-up: compare 1 run of 2 x $Nit it long"
82       echo "                       with 2 consecutive runs of $Nit it long each"
83       echo "opt = -mth : run script '$CMD' with option '-mth'"
84       echo "opt = sufx : run script 'run_cpl_test_{sufx}'"
85       echo " where: flag = 0  -> prepare for 1rst run"
86       echo "        flag = r1 -> do 1rst run and step 1"
87       echo "        flag =  1 -> move res. after 1rst run & prepare for 2nd"
88       echo "        flag = r2 -> do 2nd  run and step 2"
89       echo "        flag =  2 -> move res. after 2nd  run & prepare for 3rd"
90       echo "        flag = r3 -> do 3rd  run and step 3"
91       echo "        flag =  3 -> move res. after 3rd  run"
92       echo "        flag=(r)N+ , N=1,2,3 => do steps (+run) from: (r)N to 3"
93       echo "        flag = 4  -> compare Ocean results"
94       echo "        flag = 5  -> compare Atmos results"
95       echo "        flag = 6  -> diff pickup files"
96       exit
97    fi
98  if test $NiAt = 0 -o $NiOc = 0  if test $NiAt = 0 -o $NiOc = 0
99  then  then
100      echo " needs 2 data files: 'input_ocn/data.tst' & 'input_atm/data.tst'"      echo " needs 2 data files: 'input_ocn/data.tst' & 'input_atm/data.tst'"
101      echo "  (corresponding to 1rst run) to continue"      echo "  (corresponding to 1rst run) to continue"
102    exit    exit
103  fi  fi
104  ksel=$1  arg1=$1
105    xx=`echo $arg1 | sed 's/\+$//'`
106    if test $xx != $arg1
107    then kUp=`echo $xx | sed 's/^r//'` ; arg1=$xx
108    else kUp=9 ; fi
109    xx=`echo $arg1 | sed 's/^r//'`
110    if test $xx != $arg1 ; then doRun=1 ; arg1=$xx ; else doRun=0 ; fi
111    ksel=$arg1
112    
113    CMD="./${CMD} $MTH 3"
114    echo " doRun='$doRun' ; ksel='$ksel' ; kUp='$kUp' ; command='$CMD'"
115    
116  #make sure that local dir & ~jmc/bin is in the path:  #make sure that local dir & ~jmc/bin is in the path:
117  #export PATH=${PATH}:.  #export PATH=${PATH}:.
# Line 98  noc2=`expr $noc0 + $Dbl` Line 123  noc2=`expr $noc0 + $Dbl`
123  noc0c=`printf "%10.10i\n" $noc0`  noc0c=`printf "%10.10i\n" $noc0`
124  noc1c=`printf "%10.10i\n" $noc1`  noc1c=`printf "%10.10i\n" $noc1`
125  noc2c=`printf "%10.10i\n" $noc2`  noc2c=`printf "%10.10i\n" $noc2`
126  echo 'Oce (rank_1):' $noc0  $noc1  $noc2  echo ' Oce (rank_1):' $noc0  $noc1  $noc2
127  echo $noc0c $noc1c $noc2c  echo ' ' $noc0c $noc1c $noc2c
128    
129  Nit=$NiAt  Nit=$NiAt
130  Dbl=`expr $Nit + $Nit`  Dbl=`expr $Nit + $Nit`
# Line 108  nat2=`expr $nat0 + $Dbl` Line 133  nat2=`expr $nat0 + $Dbl`
133  nat0c=`printf "%10.10i\n" $nat0`  nat0c=`printf "%10.10i\n" $nat0`
134  nat1c=`printf "%10.10i\n" $nat1`  nat1c=`printf "%10.10i\n" $nat1`
135  nat2c=`printf "%10.10i\n" $nat2`  nat2c=`printf "%10.10i\n" $nat2`
136  echo 'Atm (rank_2):' $nat0  $nat1  $nat2  echo ' Atm (rank_2):' $nat0  $nat1  $nat2
137  echo $nat0c $nat1c $nat2c  echo ' ' $nat0c $nat1c $nat2c
138    
139  listOc="pickup pickup_cpl"  listOc="pickup"
140  listAt="pickup pickup_cpl pickup_ic pickup_land"  listAt="pickup pickup_cpl pickup_ic pickup_land"
141    
142  #- dir where to put the results :  #- dir where to put the results :
# Line 133  then Line 158  then
158   done   done
159   cd ..   cd ..
160  #-- prepare for running Dbl  #-- prepare for running Dbl
161   echo ' prepare for running Dbl:'   echo '=> prepare for running Dbl:'
162   Dbl=`expr $NiOc + $NiOc`   Dbl=`expr $NiOc + $NiOc`
163   sed "1 s/[0-9]*/$Dbl/" input_cpl/data > rank_0/data   sed "1 s/[0-9]*/$Dbl/" input_cpl/data > rank_0/data
164   cp -p input_ocn/data.tst rank_1/data   cp -p input_ocn/data.tst rank_1/data
# Line 142  then Line 167  then
167   egrep 'nIter0|nTimeSteps' rank_?/data   egrep 'nIter0|nTimeSteps' rank_?/data
168  fi  fi
169    
170  #-- after running Dbl:  #-- run & post-process Dbl:
171  if test $ksel = 1  if [ $ksel -eq 1 -o $kUp -le 1 ]
172  then  then
173     echo ' '
174    #-- run coupled test (Dbl):
175     if [ $doRun -eq 1 -o $kUp -lt 1 ]
176     then
177       echo "=> run (Dbl): $CMD"
178       $CMD
179     fi
180    #-- after running Dbl:
181   rm -rf $dir1   rm -rf $dir1
182   mkdir $dir1   mkdir $dir1
183   echo ' ' ; echo 'move STDOUT to' $dir1   echo '=> move STDOUT to' $dir1
184   mv rank_1/STDOUT* $dir1/ocnSTDOUT.2it   mv rank_1/STDOUT* $dir1/ocnSTDOUT.2it
185   mv rank_2/STDOUT* $dir1/atmSTDOUT.2it   mv rank_2/STDOUT* $dir1/atmSTDOUT.2it
186    
# Line 172  then Line 205  then
205   cd ..   cd ..
206    
207  #-- prepare for running 1iA  #-- prepare for running 1iA
208   echo ' prepare for running 1iA:'   echo '=> prepare for running 1iA:'
209   Nit=$NiOc   Nit=$NiOc
210   Dbl=`expr $Nit + $Nit`   Dbl=`expr $Nit + $Nit`
211   sed "1 s/[0-9]*/$Nit/" input_cpl/data > rank_0/data   sed "1 s/[0-9]*/$Nit/" input_cpl/data > rank_0/data
# Line 184  then Line 217  then
217   egrep 'nIter0|nTimeSteps' rank_?/data   egrep 'nIter0|nTimeSteps' rank_?/data
218  fi  fi
219    
220  #-- after running 1iA  #-- run & post-process 1iA:
221  if test $ksel = 2  if [ $ksel -eq 2 -o $kUp -le 2 ]
222  then  then
223     echo ' '
224    #-- run coupled test (1iA):
225     if [ $doRun -eq 1 -o $kUp -lt 2 ]
226     then
227       echo "=> run (1iA): $CMD"
228       $CMD
229     fi
230    #-- after running 1iA
231   rm -rf $dir2   rm -rf $dir2
232   mkdir $dir2   mkdir $dir2
233   echo ' ' ; echo 'move STDOUT to' $dir2   echo '=> move STDOUT to' $dir2
234   mv rank_1/STDOUT* $dir2/ocnSTDOUT.1iA   mv rank_1/STDOUT* $dir2/ocnSTDOUT.1iA
235   mv rank_2/STDOUT* $dir2/atmSTDOUT.1iA   mv rank_2/STDOUT* $dir2/atmSTDOUT.1iA
236    
# Line 214  then Line 255  then
255   cd ..   cd ..
256    
257  #-- prepare for running 1iB  #-- prepare for running 1iB
258   echo ' prepare for running 1iB:'   echo '=> prepare for running 1iB:'
259   cd rank_1   cd rank_1
260   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp
261   mv -f data.tmp data   mv -f data.tmp data
# Line 228  then Line 269  then
269   egrep 'nIter0|nTimeSteps' rank_?/data   egrep 'nIter0|nTimeSteps' rank_?/data
270  fi  fi
271    
272  #-- after running 1iB  #-- run & post-process 1iB:
273  if test $ksel = 3  if [ $ksel -eq 3 -o $kUp -le 3 ]
274  then  then
275     echo ' '
276    #-- run coupled test (1iB):
277     if [ $doRun -eq 1 -o $kUp -lt 3 ]
278     then
279       echo "=> run (1iB): $CMD"
280       $CMD
281     fi
282    #-- after running 1iB
283   rm -rf $dir3   rm -rf $dir3
284   mkdir $dir3   mkdir $dir3
285   echo ' ' ; echo 'move STDOUT to' $dir3   echo '=> move STDOUT to' $dir3
286   mv rank_1/STDOUT* $dir3/ocnSTDOUT.1iB   mv rank_1/STDOUT* $dir3/ocnSTDOUT.1iB
287   mv rank_2/STDOUT* $dir3/atmSTDOUT.1iB   mv rank_2/STDOUT* $dir3/atmSTDOUT.1iB
288    
# Line 244  then Line 293  then
293   done   done
294   echo 'move pickups to' $dir3/ocn   echo 'move pickups to' $dir3/ocn
295   mv pickup*.$noc2c.* ../$dir3/ocn   mv pickup*.$noc2c.* ../$dir3/ocn
296     rm -f pickup*.$noc1c.*
297   cd ..   cd ..
298    
299   mkdir $dir3/atm   mkdir $dir3/atm
# Line 253  then Line 303  then
303   done   done
304   echo 'move pickups to' $dir3/atm   echo 'move pickups to' $dir3/atm
305   mv pickup*.$nat2c.* ../$dir3/atm   mv pickup*.$nat2c.* ../$dir3/atm
306     rm -f pickup*.$nat1c.*
307   cd ..   cd ..
308   echo 'restore std data files'   echo '=> restore std data files'
309   cp -p input_cpl/data rank_0   cp -p input_cpl/data rank_0
310   cp -p input_ocn/data rank_1   cp -p input_ocn/data rank_1
311   cp -p input_atm/data rank_2   cp -p input_atm/data rank_2
# Line 324  then Line 375  then
375   done   done
376   exit   exit
377  fi  fi
   
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22