/[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.9 by jmc, Mon Dec 6 02:57:56 2010 UTC revision 1.12 by jmc, Tue Dec 3 00:03:58 2013 UTC
# Line 7  rnp_loc() Line 7  rnp_loc()
7  {  {
8  #   rnp_loc arg_1 arg_2  #   rnp_loc arg_1 arg_2
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      if [ $prt -ge 1 ] ; then echo '  rnp_loc:' $1 $2 ; fi
11  #   rnp -s $1 $2  #   rnp -s $1 $2
12      listY=`ls -1 $1*`      listY=`ls -1 $1*`
13      for yy in $listY      for yy in $listY
14      do      do
15        zz=`echo $yy | sed "s/^$1/$2/"`        zz=`echo $yy | sed "s/^$1/$2/"`
16        mv $yy $zz        if [ $prt -ge 2 ] ; then echo '   mv' $yy $zz ; fi
17          if test $action = 1 ; then    mv $yy $zz ; fi
18      done      done
19  }  }
20    
# Line 62  move_outp() Line 63  move_outp()
63      fi      fi
64  }  }
65    
66  #===============================================================================  #----------------------------------------------------------------------
67    
68  if test -f input_ocn/data.tst ; then  prepare_first()
69   noc0=`sed -n 's/nIter0=//p' input_ocn/data.tst | sed 's/,//g' | sed 's/ //g'`  {
70   Dbl=`sed -n 's/nTimeSteps=//p' input_ocn/data.tst | sed 's/,//g'`  #   prepare_first arg_1 arg_2 arg_3
71   NiOc=`expr $Dbl / 2`  #-----------------------------------------------------------
72   Nit=$NiOc  #   from a previous run in "arg_2" dir: collect info for
73  else echo " file: 'input_ocn/data.tst' not found"  #     test: arg_1 (iters) + arg_1 (iter) = 2x arg_1 (iter) ;
74   NiOc=0 ; Nit=2  #     move output files from "arg_2" dir to 'arg_3" dir ;
75  fi  #     create "data.tst" in arg_2 dir ;
76  if test -f input_atm/data.tst ; then  #     and prepare pickup and params file for 1rst run.
77   nat0=`sed -n 's/nIter0=//p' input_atm/data.tst | sed 's/,//g' | sed 's/ //g'`  #-----------------------------------------------------------
78   Dbl=`sed -n 's/nTimeSteps=//p' input_atm/data.tst | sed 's/,//g'`  
79   NiAt=`expr $Dbl / 2`    nIt=$1 ; dFr=$2 ; dTo=$3
80  else echo " file: 'input_atm/data.tst' not found"    if [ $prt -ge 1 ] ; then
81   NiAt=0        echo " run 'prepare_first' : nIt=$nIt, from '$dFr' to '$dTo'" ; fi
82  fi    outpRun="$dFr/STDOUT.0000"
83  if test -f pr_group ; then  #- last pickup written (from previous run):
84   Npr=`cat pr_group | wc -l`    if test -r $outpRun
85   NpOc=`grep -c '\/build_ocn\/mitgcmuv' pr_group`    then
86  else echo " file: 'pr_group' not found"  #   grep ' %CHECKPOINT ' $outpRun | tail -1
87   Npr=0      lastPick=`grep ' %CHECKPOINT ' $outpRun | sed 's/^.*%CHECKPOINT //' | tail -1`
88        if test "x$lastPick" = x
89        then
90         echo "no pickup information written in '$outpRun' from previous run"
91         exit 8
92        else
93         iter=`echo $lastPick | awk '{print $1}'`
94         sufx=`echo $lastPick | awk '{print $2}'`
95         echo " from '$outpRun' : lastPick='$lastPick' > iter='$iter', sufx='$sufx'"
96        fi
97      else
98        echo "no standard output file '$outpRun' from previous run"
99        exit 7
100      fi
101    #- create parameter file "data.tst":
102      echo " preparing parameter file 'data.tst' :"
103      pFile="$dFr/data.tst"
104      if test -h $pFile ; then rm -f $pFile ; fi
105      cp -p -f $dFr/data $pFile
106      listOut='startTime endTime nIter0 nTimeSteps nEndIter chkPtFreq pChkPtFreq'
107      listOut=$listOut' pickupSuff usePickupBeforeC54 mdsioLocalDir pickupStrictlyMatch'
108    #- comment out (in data.tst) all parameters from list: 'listOut'
109      for yy in $listOut
110      do
111        if [ $prt -ge 2 ] ; then echo -n '  commenting out:' $yy ; fi
112        nn=`grep -i -c "^ *$yy *=" $pFile`
113        if [ $nn -ge 1 ]
114        then
115          listN=`grep -i -n "^ *$yy *=" $pFile | sed 's/:.*$//'`
116          if [ $prt -ge 2 ] ; then echo " : found on line: $listN" ; fi
117          for ll in $listN
118          do
119            sed "$ll s/^./#&/" $pFile > data.tmp_$$
120           #diff data.tmp_$$ $pFile
121            mv -f data.tmp_$$ $pFile
122          done
123        else if [ $prt -ge 2 ] ; then echo ' : not found' ; fi
124        fi
125      done
126      if [ $prt -ge 2 ] ; then echo ' ' ; fi
127    # add nIter0 & nTimeSteps in namelist "PARM03":
128      Dbl=`expr $nIt \* 2`
129      sed "/^ *\&PARM03/a \ nTimeSteps=$Dbl," $pFile > data.tmp_$$
130      mv -f data.tmp_$$ $pFile
131      sed "/^ *\&PARM03/a \ nIter0=$iter," $pFile > data.tmp_$$
132      mv -f data.tmp_$$ $pFile
133      echo "preparing file '$pFile' : done"
134      if  [ $prt -ge 1 ] ; then
135        echo " diff $pFile $dFr/data" ; diff $pFile $dFr/data ; echo ' ' ;
136      fi
137    #---- clean up current dir: move output to temp_tst dir:
138    if test $action = 1 ; then
139      ( cd $dFr ; move_outp ../$dTo )
140      listF=`ls -1 $dFr/STD???.0[0-9][0-9][0-9] 2> /dev/null`
141      if test "x$listF" != x ; then mv $listF $dTo ; fi
142    #- move main parameter file "data":
143      mv $dFr/data $dTo
144    #- do not deal with MNC pickup at all:
145      cp -p -f $dFr/data.pkg data.tmp_$$
146      mv -f $dFr/data.pkg $dTo
147      mv -f data.tmp_$$ $dFr/data.pkg
148      yy='useMNC'
149      if [ $prt -ge 2 ] ; then echo -n "  commenting out: $yy in data.pkg" ; fi
150      nn=`grep -i -c "^ *$yy *=" $dFr/data.pkg`
151      if [ $nn -ge 1 ] ; then
152        listN=`grep -i -n "^ *$yy *=" $dFr/data.pkg | sed 's/:.*$//'`
153        if [ $prt -ge 2 ] ; then echo " : found on line: $listN" ; fi
154        for ll in $listN
155        do
156          sed "$ll s/^./#&/" $dFr/data.pkg > data.tmp_$$
157          if  [ $prt -ge 1 ] ; then echo ' changes in data.pkg'
158            diff data.tmp_$$ $dFr/data.pkg ; echo ' ' ; fi
159          mv -f data.tmp_$$ $dFr/data.pkg
160        done
161      else if [ $prt -ge 2 ] ; then echo ' : not found' ; echo ' ' ; fi
162      fi
163    #-- link back pickup files
164      listS=`ls $dTo/pickup*.${sufx}.* 2> /dev/null | head -1`
165      if test "x$listS" != x ; then
166        if [ $prt -ge 1 ] ; then echo ' link back:' ../$dTo'/pickup*.'${sufx}'*' ; fi
167        ( cd $dFr ; ln -s ../$dTo/pickup*.${sufx}* . )
168      fi
169    #- rename pickup (if needed)
170      nIt0c=`printf "%10.10i\n" $iter`
171      if test $sufx = $nIt0c
172      then
173        echo "no need to change pickup sufix"
174      else
175        listP=`(cd $dFr ; ls pickup*.${sufx}* 2> /dev/null | sed 's/\..*//g' | uniq )`
176        if test "x$listP" = x
177        then echo ' rename' $sufx  '->' $nIt0c ': no file to rename'
178        else
179         if [ $prt -ge 1 ] ; then echo ' rename' $sufx  '->' $nIt0c 'for all:' $listP ; fi
180         ( cd $dFr ; for yy in $listP
181           do rnp_loc $yy.${sufx} $yy.$nIt0c
182           done )
183        fi
184      fi
185  fi  fi
186    #----------------------------------------------------------------------
187    }
188    
189    #===============================================================================
190    #-- sequential part start here
191    #===============================================================================
192    
193  #- default:  #- default:
194    Nit=2
195  prt=1    # 1 = print some informations ; 2 = print more informations  prt=1    # 1 = print some informations ; 2 = print more informations
196  action=1  action=1
197  tmpDir='temp_tst'  tmpDir='temp_tst'
198    
199    pickRef='pickup'
200    # offLine=`grep -i -c '^ *useOffLine *= *\.TRUE\.' data.pkg`
201    # if [ $offLine -ge 1 ] ; then pickRef='pickup_ptracers' ; fi
202    
203    #- dir where to put the results :
204    dir1=res_2it
205    dir2=res_1iA
206    dir3=res_1iB
207    
208    #- get number of procs from file "pr_group":
209    if test -f pr_group ; then
210     Npr=`cat pr_group | wc -l`
211     NpOc=`grep -c '\/build_ocn\/mitgcmuv' pr_group`
212    else echo " file: 'pr_group' not found"
213     Npr=0
214    fi
215    
216  #- parse options:  #- parse options:
217  CMD='run_cpl_test'  CMD='run_cpl_test'
218  MTH=  MTH=
# Line 98  argList=$* Line 220  argList=$*
220  for xx in $argList  for xx in $argList
221  do  do
222   if test $xx = '-mth' ; then MTH=$xx ; shift ; continue ; fi   if test $xx = '-mth' ; then MTH=$xx ; shift ; continue ; fi
223   if test -x ${CMD}_$xx ; then CMD="${CMD}_$xx" ; shift ; continue ; fi   if test -x ${CMD}.$xx ; then CMD="${CMD}.$xx" ; shift ; continue ; fi
224  done  done
225    
226  if [ $# -ne 1 ] ; then  if [ $# -ne 1 ] ; then
# Line 106  if [ $# -ne 1 ] ; then Line 228  if [ $# -ne 1 ] ; then
228     echo " Check restart of coupled set-up: compare 1 run of 2 x $Nit it long"     echo " Check restart of coupled set-up: compare 1 run of 2 x $Nit it long"
229     echo "                       with 2 consecutive runs of $Nit it long each"     echo "                       with 2 consecutive runs of $Nit it long each"
230     echo "opt = -mth : run script '$CMD' with option '-mth'"     echo "opt = -mth : run script '$CMD' with option '-mth'"
231     echo "opt = sufx : run script 'run_cpl_test_{sufx}'"     echo "opt = sufx : run script 'run_cpl_test.{sufx}'"
232     echo " where: flag = 0  -> prepare for 1rst run"     echo " where: flag = 0  -> prepare for 1rst run"
233     echo "        flag = r1 -> do 1rst run and step 1"     echo "        flag = r1 -> do 1rst run and step 1"
234     echo "        flag =  1 -> move res. after 1rst run & prepare for 2nd"     echo "        flag =  1 -> move res. after 1rst run & prepare for 2nd"
# Line 121  if [ $# -ne 1 ] ; then Line 243  if [ $# -ne 1 ] ; then
243     echo "        flag = 7  -> clean-up output files"     echo "        flag = 7  -> clean-up output files"
244     exit     exit
245  fi  fi
 if test $NiAt = 0 -o $NiOc = 0 ; then  
     echo " needs 2 data files: 'input_ocn/data.tst' & 'input_atm/data.tst'"  
     echo "  (corresponding to 1rst run) to continue"  
   exit  
 fi  
 if test $Npr = 0 ; then  
     echo " needs file: 'pr_group' to continue" ; exit  
 fi  
246  arg1=$1  arg1=$1
247  xx=`echo $arg1 | sed 's/\+$//'`  xx=`echo $arg1 | sed 's/\+$//'`
248  if test $xx != $arg1  if test $xx != $arg1
# Line 141  ksel=$arg1 Line 255  ksel=$arg1
255  CMD="./${CMD} $MTH 3"  CMD="./${CMD} $MTH 3"
256  echo " doRun='$doRun' ; ksel='$ksel' ; kUp='$kUp' ; command='$CMD'"  echo " doRun='$doRun' ; ksel='$ksel' ; kUp='$kUp' ; command='$CMD'"
257    
258    #- set dir of Master-Proc for 3 components (coupler == rank_0):
259    if test $Npr = 0 ; then
260        echo " needs file: 'pr_group' to continue" ; exit
261    fi
262    rnkO=1 ; rnkA=`expr $rnkO + $NpOc`
263    
264  #make sure that local dir & ~jmc/bin is in the path:  #make sure that local dir & ~jmc/bin is in the path:
265  #export PATH=${PATH}:.  #export PATH=${PATH}:.
266    
267  rnkO=1 ; rnkA=`expr $rnkO + $NpOc`  #----------------------------------------------------------------------
268  Nit=$NiOc  
269  Dbl=`expr $Nit + $Nit`  if test $ksel = 0 ; then
270  noc1=`expr $noc0 + $Nit`  #-- from STDOUT of previous run, get Nb of iters for each component,
271      yy='(PID.TID 0000.0001) nTimeSteps = '
272      outpRun="rank_$rnkO/STDOUT.0000"
273      nTsO=`grep -A1 "$yy" $outpRun | tail -1 | sed 's/(PID.TID 0000.0001) *//'`
274      outpRun="rank_$rnkA/STDOUT.0000"
275      nTsA=`grep -A1 "$yy" $outpRun | tail -1 | sed 's/(PID.TID 0000.0001) *//'`
276    #-  then calculate the ratio and set Nb of time-step for Atmos comp.
277      RatioTs=`expr $nTsA / $nTsO`
278      if [ $prt -ge 1 ] ; then echo " nTsA= $nTsA , nTsO= $nTsO , Ratio= $RatioTs"; fi
279      NiOc=$Nit ;
280      NiAt=`expr $NiOc \* $RatioTs`
281    #---- clean up current dir: rename pickup files and move output to temp_tst dir:
282     /bin/rm -rf $tmpDir.0 $tmpDir.1 $tmpDir.2
283     mkdir $tmpDir.0 $tmpDir.1 $tmpDir.2
284     mv std_outp rank_0/*.clog $tmpDir.0
285     mv rank_$rnkO/*.0001.clog $tmpDir.0
286     mv rank_$rnkA/*.0001.clog $tmpDir.0
287     echo ' '
288     echo "=> prepare OCN dir:"
289     prepare_first $NiOc rank_$rnkO $tmpDir.1
290     echo "=> prepare ATM dir:"
291     prepare_first $NiAt rank_$rnkA $tmpDir.2
292    #-- prepare for running Dbl
293     echo '=> prepare for running Dbl:'
294     Dbl=`expr $NiOc + $NiOc`
295     cp -p rank_$rnkO/data.tst rank_$rnkO/data
296     cp -p rank_$rnkA/data.tst rank_$rnkA/data
297     egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
298    fi
299    
300    #----------------------------------------------------------------------
301    
302    #- extract length of runs & starting iter number from file "data.tst"
303    #  Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)
304    pFile="rank_$rnkO/data.tst"
305    if test -f $pFile ; then
306      noc0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
307      Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
308      NiOc=`expr $Dbl / 2`
309      Nit=$NiOc
310    else echo " file: '$pFile' not found"
311      NiOc=0 ; Nit=2
312    fi
313    pFile="rank_$rnkA/data.tst"
314    if test -f $pFile ; then
315      nat0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
316      Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
317      NiAt=`expr $Dbl / 2`
318    else echo " file: '$pFile' not found"
319      NiAt=0
320    fi
321    if test $NiAt = 0 -o $NiOc = 0 ; then
322        echo " needs 2 data files: 'rank_$rnkO/data.tst' & 'rank_$rnkA/data.tst'"
323        echo "  (corresponding to 1rst run) to continue"
324      exit
325    fi
326    
327    nit=$NiOc
328    Dbl=`expr $nit + $nit`
329    noc1=`expr $noc0 + $nit`
330  noc2=`expr $noc0 + $Dbl`  noc2=`expr $noc0 + $Dbl`
331  noc0c=`printf "%10.10i\n" $noc0`  noc0c=`printf "%10.10i\n" $noc0`
332  noc1c=`printf "%10.10i\n" $noc1`  noc1c=`printf "%10.10i\n" $noc1`
# Line 155  noc2c=`printf "%10.10i\n" $noc2` Line 334  noc2c=`printf "%10.10i\n" $noc2`
334  echo " Oce (rank_$rnkO):" $noc0  $noc1  $noc2  echo " Oce (rank_$rnkO):" $noc0  $noc1  $noc2
335  echo "     " $noc0c $noc1c $noc2c  echo "     " $noc0c $noc1c $noc2c
336    
337  Nit=$NiAt  nit=$NiAt
338  Dbl=`expr $Nit + $Nit`  Dbl=`expr $nit + $nit`
339  nat1=`expr $nat0 + $Nit`  nat1=`expr $nat0 + $nit`
340  nat2=`expr $nat0 + $Dbl`  nat2=`expr $nat0 + $Dbl`
341  nat0c=`printf "%10.10i\n" $nat0`  nat0c=`printf "%10.10i\n" $nat0`
342  nat1c=`printf "%10.10i\n" $nat1`  nat1c=`printf "%10.10i\n" $nat1`
# Line 165  nat2c=`printf "%10.10i\n" $nat2` Line 344  nat2c=`printf "%10.10i\n" $nat2`
344  echo " Atm (rank_$rnkA):" $nat0  $nat1  $nat2  echo " Atm (rank_$rnkA):" $nat0  $nat1  $nat2
345  echo "     " $nat0c $nat1c $nat2c  echo "     " $nat0c $nat1c $nat2c
346    
347  listOc="pickup"  #----------------------------------------------------------------------
 listAt="pickup pickup_cpl pickup_ic pickup_land"  
   
 #- dir where to put the results :  
 dir1=res_2it  
 dir2=res_1iA  
 dir3=res_1iB  
   
 #---- clean up current dir: rename pickup files and move output to temp_tst dir:  
 if test $ksel = 0 ; then  
  /bin/rm -rf $tmpDir.0 $tmpDir.1 $tmpDir.2  
  mkdir $tmpDir.0 $tmpDir.1 $tmpDir.2  
  mv std_outp rank_0/*.clog rank_0/data $tmpDir.0  
  echo ' '  
  cd rank_$rnkO  
  for xx in $listOc  
  do rnp_loc $xx.ckptA $xx.$noc0c  
  done  
  move_outp ../$tmpDir.1  
  mv data ../$tmpDir.1  
  mv *.0001.clog ../$tmpDir.0  
 #-- link back pickup files  
  sufx=$noc0c  
  listS=`ls ../$tmpDir.1/pickup*.${sufx}.* 2> /dev/null | head -1`  
  if test "x$listS" != x ; then  
    if [ $prt -ge 1 ] ; then echo ' link back:' ../$tmpDir'.1/pickup*.'${sufx}'*' ; fi  
    if test $action = 1 ; then ln -s ../$tmpDir.1/pickup*.${sufx}* . ; fi  
  fi  
  cd ..  
  cd rank_$rnkA  
  for xx in $listAt  
  do rnp_loc $xx.ckptA $xx.$nat0c  
  done  
  move_outp ../$tmpDir.2  
  mv data ../$tmpDir.2  
  mv *.0001.clog ../$tmpDir.0  
 #-- link back pickup files  
  sufx=$nat0c  
  listS=`ls ../$tmpDir.2/pickup*.${sufx}.* 2> /dev/null | head -1`  
  if test "x$listS" != x ; then  
    if [ $prt -ge 1 ] ; then echo ' link back:' ../$tmpDir'.2/pickup*.'${sufx}'*' ; fi  
    if test $action = 1 ; then ln -s ../$tmpDir.2/pickup*.${sufx}* . ; fi  
  fi  
  cd ..  
 #-- prepare for running Dbl  
  echo '=> prepare for running Dbl:'  
  Dbl=`expr $NiOc + $NiOc`  
  sed "1 s/[0-9]*/$Dbl/" input_cpl/data > rank_0/data  
  cp -p input_ocn/data.tst rank_$rnkO/data  
  cp -p input_atm/data.tst rank_$rnkA/data  
  echo 'rank_0/data:' `head -1 rank_0/data`  
  egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data  
 fi  
348    
349  #-- run & post-process Dbl:  #-- run & post-process Dbl:
350  if [ $ksel -eq 1 -o $kUp -le 1 ] ; then  if [ $ksel -eq 1 -o $kUp -le 1 ] ; then
# Line 228  if [ $ksel -eq 1 -o $kUp -le 1 ] ; then Line 355  if [ $ksel -eq 1 -o $kUp -le 1 ] ; then
355     $CMD     $CMD
356   fi   fi
357  #-- after running Dbl:  #-- after running Dbl:
358    #-  check for pickup
359     listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
360     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
361     if test "x$listOc" = x -o "x$listAt" = x
362     then echo "==> RUN 2 x $NiOc it STOP without writing $pickRef => exit" ; exit 1
363     else echo "==>  END  RUN 2 x $NiOc it" ; fi
364    #-  move output
365   rm -rf $dir1   rm -rf $dir1
366   mkdir $dir1   mkdir $dir1
367   echo '=> move STDOUT to' $dir1   echo '=> move STDOUT to' $dir1
# Line 237  if [ $ksel -eq 1 -o $kUp -le 1 ] ; then Line 371  if [ $ksel -eq 1 -o $kUp -le 1 ] ; then
371    
372   mkdir $dir1/ocn   mkdir $dir1/ocn
373   cd rank_$rnkO   cd rank_$rnkO
374   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
375   do   echo ' listP=' $listP
376    rnp_loc $xx.ckptA $xx.$noc2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
  done  
377   echo 'move output to' $dir1/ocn   echo 'move output to' $dir1/ocn
378   move_outp  ../$dir1/ocn   move_outp  ../$dir1/ocn
379   cd ..   cd ..
380    
381   mkdir $dir1/atm   mkdir $dir1/atm
382   cd rank_$rnkA   cd rank_$rnkA
383   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
384   do   echo ' listP=' $listP
385    rnp_loc $xx.ckptA $xx.$nat2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
  done  
386   echo 'move output to' $dir1/atm   echo 'move output to' $dir1/atm
387   move_outp  ../$dir1/atm   move_outp  ../$dir1/atm
388   cd ..   cd ..
389    
390  #-- prepare for running 1iA  #-- prepare for running 1iA
391   echo '=> prepare for running 1iA:'   echo '=> prepare for running 1iA:'
392   Nit=$NiOc   nit=$NiOc ; Dbl=`expr $nit + $nit`
393   Dbl=`expr $Nit + $Nit`   pFile=rank_$rnkO/data
394   sed "1 s/[0-9]*/$Nit/" input_cpl/data > rank_0/data   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
395   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" input_ocn/data.tst > rank_$rnkO/data   nit=$NiAt ; Dbl=`expr $nit + $nit`
396   Nit=$NiAt   pFile=rank_$rnkA/data
397   Dbl=`expr $Nit + $Nit`   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
  sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" input_atm/data.tst > rank_$rnkA/data  
  echo 'rank_0/data:' `head -1 rank_0/data`  
398   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
399  fi  fi
400    
# Line 277  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then Line 407  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then
407     $CMD     $CMD
408   fi   fi
409  #-- after running 1iA  #-- after running 1iA
410    #-  check for pickup
411     listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
412     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
413     if test "x$listOc" = x -o "x$listAt" = x
414     then echo "==> RUN 1iA STOP without writing $pickRef => exit" ; exit 2
415     else echo "==>  END  RUN 1iA" ; fi
416    #-  move output
417   rm -rf $dir2   rm -rf $dir2
418   mkdir $dir2   mkdir $dir2
419   echo '=> move STDOUT to' $dir2   echo '=> move STDOUT to' $dir2
# Line 286  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then Line 423  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then
423    
424   mkdir $dir2/ocn   mkdir $dir2/ocn
425   cd rank_$rnkO   cd rank_$rnkO
426   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
427   do   echo ' listP=' $listP
428    rnp_loc $xx.ckptA $xx.$noc1c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc1c ; done
  done  
429  # remove pickup*.$noc0c* sym-link:  # remove pickup*.$noc0c* sym-link:
430   find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*' -exec rm {} \;   find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*' -exec rm {} \;
431   echo 'move output to' $dir2/ocn   echo 'move output to' $dir2/ocn
# Line 298  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then Line 434  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then
434    
435   mkdir $dir2/atm   mkdir $dir2/atm
436   cd rank_$rnkA   cd rank_$rnkA
437   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
438   do   echo ' listP=' $listP
439    rnp_loc $xx.ckptA $xx.$nat1c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat1c ; done
  done  
440  # remove pickup*.$nat0c* sym-link:  # remove pickup*.$nat0c* sym-link:
441   find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*' -exec rm {} \;   find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*' -exec rm {} \;
442   echo 'move output to' $dir2/atm   echo 'move output to' $dir2/atm
# Line 312  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then Line 447  if [ $ksel -eq 2 -o $kUp -le 2 ] ; then
447   echo '=> prepare for running 1iB:'   echo '=> prepare for running 1iB:'
448   cd rank_$rnkO   cd rank_$rnkO
449   ln -s ../$dir2/ocn/pickup*.$noc1c.* .   ln -s ../$dir2/ocn/pickup*.$noc1c.* .
450   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp_$$
451   mv -f data.tmp data   mv -f data.tmp_$$ data
452   cd ../rank_$rnkA   cd ../rank_$rnkA
453   ln -s ../$dir2/atm/pickup*.$nat1c.* .   ln -s ../$dir2/atm/pickup*.$nat1c.* .
454   sed "s/^ nIter0=$nat0/ nIter0=$nat1/g" data > data.tmp   sed "s/^ nIter0=$nat0/ nIter0=$nat1/g" data > data.tmp_$$
455   mv -f data.tmp data   mv -f data.tmp_$$ data
456   cd ..   cd ..
  echo 'rank_0/data:' `head -1 rank_0/data`  
457   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
458  fi  fi
459    
# Line 332  if [ $ksel -eq 3 -o $kUp -le 3 ] ; then Line 466  if [ $ksel -eq 3 -o $kUp -le 3 ] ; then
466     $CMD     $CMD
467   fi   fi
468  #-- after running 1iB  #-- after running 1iB
469    #-  check for pickup
470     listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
471     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
472     if test "x$listOc" = x -o "x$listAt" = x
473     then echo "==> RUN 1iB STOP without writing $pickRef => exit" ; exit 3
474     else echo "==>  END  RUN 1iB" ; fi
475    #-  move output
476   rm -rf $dir3   rm -rf $dir3
477   mkdir $dir3   mkdir $dir3
478   echo '=> move STDOUT to' $dir3   echo '=> move STDOUT to' $dir3
# Line 341  if [ $ksel -eq 3 -o $kUp -le 3 ] ; then Line 482  if [ $ksel -eq 3 -o $kUp -le 3 ] ; then
482    
483   mkdir $dir3/ocn   mkdir $dir3/ocn
484   cd rank_$rnkO   cd rank_$rnkO
485   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
486   do rnp_loc $xx.ckptA $xx.$noc2c   echo ' listP=' $listP
487   done   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
488  # remove pickup*.$noc1c* sym-link:  # remove pickup*.$noc1c* sym-link:
489   find . ! -name . -prune -type l -name 'pickup*.'${noc1c}'*' -exec rm {} \;   find . ! -name . -prune -type l -name 'pickup*.'${noc1c}'*' -exec rm {} \;
490   echo 'move output to' $dir3/ocn   echo 'move output to' $dir3/ocn
# Line 352  if [ $ksel -eq 3 -o $kUp -le 3 ] ; then Line 493  if [ $ksel -eq 3 -o $kUp -le 3 ] ; then
493    
494   mkdir $dir3/atm   mkdir $dir3/atm
495   cd rank_$rnkA   cd rank_$rnkA
496   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
497   do rnp_loc $xx.ckptA $xx.$nat2c   echo ' listP=' $listP
498   done   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
499  # remove pickup*.$nat1c* sym-link:  # remove pickup*.$nat1c* sym-link:
500   find . ! -name . -prune -type l -name 'pickup*.'${nat1c}'*' -exec rm {} \;   find . ! -name . -prune -type l -name 'pickup*.'${nat1c}'*' -exec rm {} \;
501   echo 'move output to' $dir3/atm   echo 'move output to' $dir3/atm
502   move_outp  ../$dir3/atm   move_outp  ../$dir3/atm
503   cd ..   cd ..
504    
505  fi  fi
506    
507  #-- compare output:  #-- compare output:
# Line 368  if test $ksel = 4 ; then Line 510  if test $ksel = 4 ; then
510   echo ' run 2it: ---------- ' > tmp1.txt   echo ' run 2it: ---------- ' > tmp1.txt
511   grep "cg2d_init_res" $dir1/ocnSTDOUT.0000 \   grep "cg2d_init_res" $dir1/ocnSTDOUT.0000 \
512        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
513   Nit=`expr $NiOc + 2`   nit=`expr $NiOc + 2`
514   echo "NiOc= $NiOc ; Nit= $Nit"   #echo "NiOc= $NiOc ; nit= $nit"
515   sed "$Nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt   sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
516   rm -f tmp1.txt   rm -f tmp1.txt
517   echo ' run 1iA: ---------- ' > tmp1.txt   echo ' run 1iA: ---------- ' > tmp1.txt
518   grep "cg2d_init_res" $dir2/ocnSTDOUT.0000 \   grep "cg2d_init_res" $dir2/ocnSTDOUT.0000 \
# Line 392  if test $ksel = 5 ; then Line 534  if test $ksel = 5 ; then
534   echo ' run 2it: ---------- ' > tmp1.txt   echo ' run 2it: ---------- ' > tmp1.txt
535   grep "cg2d_init_res" $dir1/atmSTDOUT.0000 \   grep "cg2d_init_res" $dir1/atmSTDOUT.0000 \
536        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
537   Nit=`expr $NiAt + 2`   nit=`expr $NiAt + 2`
538   sed "$Nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt   sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
539   rm -f tmp1.txt   rm -f tmp1.txt
540   echo ' run 1iA: ---------- ' > tmp1.txt   echo ' run 1iA: ---------- ' > tmp1.txt
541   grep "cg2d_init_res" $dir2/atmSTDOUT.0000 \   grep "cg2d_init_res" $dir2/atmSTDOUT.0000 \
# Line 413  fi Line 555  fi
555  #- compare pickup file:  #- compare pickup file:
556  if test $ksel = 6 ; then  if test $ksel = 6 ; then
557   echo ' ' ; echo '=== compare Ocn pickups: =============='   echo ' ' ; echo '=== compare Ocn pickups: =============='
558   for xx in $listOc   listP=`( cd $dir1/ocn ; ls pickup*.$noc2c* | sed 's/\..*//g' | uniq )`
559     for xx in $listP
560   do   do
561     dif_tiles $dir1/ocn $dir3/ocn $xx.$noc2c     dif_tiles $dir1/ocn $dir3/ocn $xx.$noc2c
562   done   done
563   echo ' ' ; echo '=== compare Atm pickups: =============='   echo ' ' ; echo '=== compare Atm pickups: =============='
564   for xx in $listAt   listP=`( cd $dir1/atm ; ls pickup*.$nat2c* | sed 's/\..*//g' | uniq )`
565     for xx in $listP
566   do   do
567     dif_tiles $dir1/atm $dir3/atm $xx.$nat2c     dif_tiles $dir1/atm $dir3/atm $xx.$nat2c
568   done   done
# Line 428  fi Line 572  fi
572  #- clean-up and reset:  #- clean-up and reset:
573  if test $ksel = 7 ; then  if test $ksel = 7 ; then
574   echo "== clean-up and reset :"   echo "== clean-up and reset :"
575  #echo '--> remove pickup*.'$noc0c'* sym-links:'   cd rank_$rnkO
576  #cd rank_$rnkO   listS=`find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*'`
577  #  listS=`find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*'`   if test "x$listS" != x ; then
578  #  if [ $prt -ge 2 ] ; then echo " in rank_$rnkO: rm" $listS ; fi     echo '--> remove pickup*.'$noc0c'* sym-links:'
579  #  rm -f $listS     if [ $prt -ge 2 ] ; then echo " in rank_$rnkO: rm" $listS ; fi
580  #echo '--> remove pickup*.'$nat0c'* sym-links:'     rm -f $listS
581  #cd ../rank_$rnkA   fi
582  #  listS=`find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*'`   cd ../rank_$rnkA
583  #  if [ $prt -ge 2 ] ; then echo " in rank_$rnkA: rm" $listS ; fi   listS=`find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*'`
584  #  rm -f $listS   if test "x$listS" != x ; then
585  #cd ..     echo '--> remove pickup*.'$nat0c'* sym-links:'
586       if [ $prt -ge 2 ] ; then echo " in rank_$rnkA: rm" $listS ; fi
587       rm -f $listS
588     fi
589     cd ..
590  #- move back files from temp_tst dir:  #- move back files from temp_tst dir:
591   echo '=> restore std data files'   echo '=> restore std data files'
592   rm -f rank_0/data  rank_$rnkO/data rank_$rnkA/data   rm -f rank_{0,$rnkO,$rnkA}/data
593     rm -f rank_{$rnkO,$rnkA}/data.tst
594     rm -f rank_{$rnkO,$rnkA}/data.pkg
595   mv $tmpDir.0/*Ocean.0001.clog $tmpDir.1/* rank_$rnkO   mv $tmpDir.0/*Ocean.0001.clog $tmpDir.1/* rank_$rnkO
596   mv $tmpDir.0/*Atmos.0001.clog $tmpDir.2/* rank_$rnkA   mv $tmpDir.0/*Atmos.0001.clog $tmpDir.2/* rank_$rnkA
597   mv $tmpDir.0/std_outp .   mv $tmpDir.0/std_outp .
598   mv $tmpDir.0/* rank_0   mv $tmpDir.0/* rank_0
599   rmdir $tmpDir.0 $tmpDir.1 $tmpDir.2   rmdir $tmpDir.0 $tmpDir.1 $tmpDir.2
600  #echo '--> remove output dir and output files:'  #echo '--> remove output dir and output files:'
601  #echo ' rm -rf' $dir1 $dir2 $dir3   echo ' rm -rf' $dir1 $dir2 $dir3
602  #rm -rf $dir1 $dir2 $dir3   rm -rf $dir1 $dir2 $dir3
603  #rm -f atm.txt ocn.txt   rm -f atm.txt ocn.txt
604   exit   exit
605  fi  fi
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22