/[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.1 by jmc, Tue Oct 30 22:41:33 2007 UTC revision 1.11 by jmc, Thu Aug 16 23:08:56 2012 UTC
# Line 1  Line 1 
1  #!/bin/sh  #! /usr/bin/env bash
2    
3  noc0=`sed -n 's/nIter0=//p' input_ocn/data.tst | sed 's/,//g' | sed 's/ //g'`  # $Header$
4  Dbl=`sed -n 's/nTimeSteps=//p' input_ocn/data.tst | sed 's/,//g'`  # $Name$
 NiOc=`expr $Dbl / 2`  
 nat0=`sed -n 's/nIter0=//p' input_atm/data.tst | sed 's/,//g' | sed 's/ //g'`  
 Dbl=`sed -n 's/nTimeSteps=//p' input_atm/data.tst | sed 's/,//g'`  
 NiAt=`expr $Dbl / 2`  
   
 Nit=$NiOc  
 case $# in  
   1) ;;  
    *) echo "Usage: `basename $0` flag"  
       echo " Check restart of coupled set-up: compare 1 run of 2 x $Nit it long"  
       echo "                       with 2 consecutive runs of $Nit it long each"  
       echo " where: flag = 0 -> prepare for 1rst run"  
       echo "        flag = 1 -> move res. after 1rst run & prepare for 2nd"  
       echo "        flag = 2 -> move res. after 2nd  run & prepare for 3rd"  
       echo "        flag = 3 -> move res. after 3rd  run"  
       echo "        flag = 4 -> compare Ocean results"  
       echo "        flag = 5 -> compare Atmos results"  
       echo "        flag = 6 -> diff pickup files"  
       echo " NB: need 2 data files: 'input_ocn/data.tst' & 'input_atm/data.tst'"  
       echo "                 corresponding to 1rst run"  
       exit ;;  
 esac  
 ksel=$1  
5    
6  #make sure that local dir & ~jmc/bin is in the path:  rnp_loc()
7  #export PATH=${PATH}:.  {
8  RNPref=`ls ~jmc/bin/rnp`" -s"  #   rnp_loc arg_1 arg_2
9  Prt10c=`ls ~jmc/bin/add0upto10c`  #   rename files with prefix = 'arg_1' to files with prefix = 'arg_2'
10   difex=`ls ~jmc/bin/dif_data`      if [ $prt -ge 1 ] ; then echo '  rnp_loc:' $1 $2 ; fi
11  echo -n " use script: "  #   rnp -s $1 $2
12  echo    "   Prt10c='$Prt10c' ,"      listY=`ls -1 $1*`
13  echo -n "   RNPref='$RNPref'"      for yy in $listY
14  echo " and difex='$difex'"      do
15          zz=`echo $yy | sed "s/^$1/$2/"`
16  Nit=$NiOc        if [ $prt -ge 2 ] ; then echo '   mv' $yy $zz ; fi
17  Dbl=`expr $Nit + $Nit`        if test $action = 1 ; then    mv $yy $zz ; fi
18  noc1=`expr $noc0 + $Nit`      done
19  noc2=`expr $noc0 + $Dbl`  }
20  noc0c=`$Prt10c $noc0`  
21  noc1c=`$Prt10c $noc1`  dif_tiles()
22  noc2c=`$Prt10c $noc2`  {
23  echo 'Oce (rank_1):' $noc0  $noc1  $noc2  #   dif_data_loc dd1 dd2 file
24  echo $noc0c $noc1c $noc2c  #   for all tiles  ==> do diff dd1/file dd2/file'
25        dd1=$1 ; dd2=$2 ; zz=$3
26  Nit=$NiAt      chkmeta='N'
27  Dbl=`expr $Nit + $Nit`      listY=`(cd $dd1 ; ls $zz.*.data | sed "s/$zz\.//" | sed "s/\.data//")`
28  nat1=`expr $nat0 + $Nit`      echo '--> file='$zz', listY='$listY
29  nat2=`expr $nat0 + $Dbl`      for yy in $listY
30  nat0c=`$Prt10c $nat0`      do
31  nat1c=`$Prt10c $nat1`        echo diff $dd1/$zz.$yy.data $dd2
32  nat2c=`$Prt10c $nat2`        diff $dd1/$zz.$yy.data $dd2
33  echo 'Atm (rank_2):' $nat0  $nat1  $nat2        out=$?
34  echo $nat0c $nat1c $nat2c        if test $out != 0
35          then echo 'Diff outp=' $out ' ==> stop'
36           exit ; fi
37          if test $chkmeta = 'Y' ; then
38           echo diff $dd1/$zz.$yy.meta $dd2
39           diff $dd1/$zz.$yy.meta $dd2
40           out=$?
41           if test $out != 0 ; then echo 'Diff outp=' $out ' ==> stop' ; exit ; fi
42          fi
43        done
44    }
45    
46    move_outp()
47    {
48    #   move_outp arg_1
49    #   move mitgcm output files to directory 'arg_1'
50        if [ $prt -ge 1 ] ; then echo ' move_outp:' $1 ; fi
51        mv *.data *.meta $1
52        listF=`ls -1 *.txt STD???.0000 *.log 2> /dev/null`
53        if test "x$listF" != x ; then mv $listF $1 ; fi
54        #- move back sym link:
55        listL=`find $1 -type l`
56        ttd=`echo $listL | wc -w`
57        #echo "ttd=$ttd"
58        if test $ttd != 0 ; then
59          #echo "listL='$listL'"
60          shortList=`ls $listL | sed "s|$1/||"`
61          if [ $prt -ge 2 ] ; then echo "  mv back:" $shortList ; fi
62          mv $listL .
63        fi
64    }
65    
66    #----------------------------------------------------------------------
67    
68    prepare_first()
69    {
70    #   prepare_first arg_1 arg_2 arg_3
71    #-----------------------------------------------------------
72    #   from a previous run in "arg_2" dir: collect info for
73    #     test: arg_1 (iters) + arg_1 (iter) = 2x arg_1 (iter) ;
74    #     move output files from "arg_2" dir to 'arg_3" dir ;
75    #     create "data.tst" in arg_2 dir ;
76    #     and prepare pickup and params file for 1rst run.
77    #-----------------------------------------------------------
78    
79      nIt=$1 ; dFr=$2 ; dTo=$3
80      if [ $prt -ge 1 ] ; then
81          echo " run 'prepare_first' : nIt=$nIt, from '$dFr' to '$dTo'" ; fi
82      outpRun="$dFr/STDOUT.0000"
83    #- last pickup written (from previous run):
84      if test -r $outpRun
85      then
86    #   grep ' %CHECKPOINT ' $outpRun | tail -1
87        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
186    #----------------------------------------------------------------------
187    }
188    
189  listOc="pickup pickup_cpl"  #===============================================================================
190  listAt="pickup pickup_cpl pickup_ic pickup_land"  #-- sequential part start here
191    #===============================================================================
192    
193    #- default:
194    Nit=2
195    prt=1    # 1 = print some informations ; 2 = print more informations
196    action=1
197    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 :  #- dir where to put the results :
204  dir1=res_2it  dir1=res_2it
205  dir2=res_1iA  dir2=res_1iA
206  dir3=res_1iB  dir3=res_1iB
207    
208  #-- after a run, rename pickup file to restart  #- get number of procs from file "pr_group":
209  if test $ksel = 0  if test -f pr_group ; then
210  then   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:
217    CMD='run_cpl_test'
218    MTH=
219    argList=$*
220    for xx in $argList
221    do
222     if test $xx = '-mth' ; then MTH=$xx ; shift ; continue ; fi
223     if test -x ${CMD}.$xx ; then CMD="${CMD}.$xx" ; shift ; continue ; fi
224    done
225    
226    if [ $# -ne 1 ] ; then
227       echo "Usage: `basename $0` [opt] flag"
228       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"
230       echo "opt = -mth : run script '$CMD' with option '-mth'"
231       echo "opt = sufx : run script 'run_cpl_test.{sufx}'"
232       echo " where: flag = 0  -> prepare for 1rst run"
233       echo "        flag = r1 -> do 1rst run and step 1"
234       echo "        flag =  1 -> move res. after 1rst run & prepare for 2nd"
235       echo "        flag = r2 -> do 2nd  run and step 2"
236       echo "        flag =  2 -> move res. after 2nd  run & prepare for 3rd"
237       echo "        flag = r3 -> do 3rd  run and step 3"
238       echo "        flag =  3 -> move res. after 3rd  run"
239       echo "        flag=(r)N+ , N=1,2,3 => do steps (+run) from: (r)N to 3"
240       echo "        flag = 4  -> compare Ocean results"
241       echo "        flag = 5  -> compare Atmos results"
242       echo "        flag = 6  -> diff pickup files"
243       echo "        flag = 7  -> clean-up output files"
244       exit
245    fi
246    arg1=$1
247    xx=`echo $arg1 | sed 's/\+$//'`
248    if test $xx != $arg1
249    then kUp=`echo $xx | sed 's/^r//'` ; arg1=$xx
250    else kUp=9 ; fi
251    xx=`echo $arg1 | sed 's/^r//'`
252    if test $xx != $arg1 ; then doRun=1 ; arg1=$xx ; else doRun=0 ; fi
253    ksel=$arg1
254    
255    CMD="./${CMD} $MTH 3"
256    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:
265    #export PATH=${PATH}:.
266    
267    #----------------------------------------------------------------------
268    
269    if test $ksel = 0 ; then
270    #-- 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 rank_0/data $tmpDir.0
285     mv rank_$rnkO/*.0001.clog $tmpDir.0
286     mv rank_$rnkA/*.0001.clog $tmpDir.0
287   echo ' '   echo ' '
288   cd rank_1   echo "=> prepare OCN dir:"
289   for xx in $listOc   prepare_first $NiOc rank_$rnkO $tmpDir.1
290   do   echo "=> prepare ATM dir:"
291    echo $RNPref $xx.ckptA $xx.$noc0c   prepare_first $NiAt rank_$rnkA $tmpDir.2
   $RNPref $xx.ckptA $xx.$noc0c  
  done  
  cd ../rank_2  
  for xx in $listAt  
  do  
   echo $RNPref $xx.ckptA $xx.$nat0c  
   $RNPref $xx.ckptA $xx.$nat0c  
  done  
  cd ..  
292  #-- prepare for running Dbl  #-- prepare for running Dbl
293   echo ' prepare for running Dbl:'   echo '=> prepare for running Dbl:'
294   Dbl=`expr $NiOc + $NiOc`   Dbl=`expr $NiOc + $NiOc`
295   sed "1 s/[0-9]*/$Dbl/" input_cpl/data > rank_0/data   sed "1 s/[0-9]*/$Dbl/" input_cpl/data > rank_0/data
296   cp -p input_ocn/data.tst rank_1/data   cp -p rank_$rnkO/data.tst rank_$rnkO/data
297   cp -p input_atm/data.tst rank_2/data   cp -p rank_$rnkA/data.tst rank_$rnkA/data
298   echo 'rank_0/data:' `head -1 rank_0/data`   echo 'rank_0/data:' `head -1 rank_0/data`
299   egrep 'nIter0|nTimeSteps' rank_?/data   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
300    fi
301    
302    #----------------------------------------------------------------------
303    
304    #- extract length of runs & starting iter number from file "data.tst"
305    #  Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)
306    pFile="rank_$rnkO/data.tst"
307    if test -f $pFile ; then
308      noc0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
309      Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
310      NiOc=`expr $Dbl / 2`
311      Nit=$NiOc
312    else echo " file: '$pFile' not found"
313      NiOc=0 ; Nit=2
314    fi
315    pFile="rank_$rnkA/data.tst"
316    if test -f $pFile ; then
317      nat0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
318      Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
319      NiAt=`expr $Dbl / 2`
320    else echo " file: '$pFile' not found"
321      NiAt=0
322    fi
323    if test $NiAt = 0 -o $NiOc = 0 ; then
324        echo " needs 2 data files: 'rank_$rnkO/data.tst' & 'rank_$rnkA/data.tst'"
325        echo "  (corresponding to 1rst run) to continue"
326      exit
327  fi  fi
328    
329    nit=$NiOc
330    Dbl=`expr $nit + $nit`
331    noc1=`expr $noc0 + $nit`
332    noc2=`expr $noc0 + $Dbl`
333    noc0c=`printf "%10.10i\n" $noc0`
334    noc1c=`printf "%10.10i\n" $noc1`
335    noc2c=`printf "%10.10i\n" $noc2`
336    echo " Oce (rank_$rnkO):" $noc0  $noc1  $noc2
337    echo "     " $noc0c $noc1c $noc2c
338    
339    nit=$NiAt
340    Dbl=`expr $nit + $nit`
341    nat1=`expr $nat0 + $nit`
342    nat2=`expr $nat0 + $Dbl`
343    nat0c=`printf "%10.10i\n" $nat0`
344    nat1c=`printf "%10.10i\n" $nat1`
345    nat2c=`printf "%10.10i\n" $nat2`
346    echo " Atm (rank_$rnkA):" $nat0  $nat1  $nat2
347    echo "     " $nat0c $nat1c $nat2c
348    
349    #----------------------------------------------------------------------
350    
351    #-- run & post-process Dbl:
352    if [ $ksel -eq 1 -o $kUp -le 1 ] ; then
353     echo ' '
354    #-- run coupled test (Dbl):
355     if [ $doRun -eq 1 -o $kUp -lt 1 ] ; then
356       echo "=> run (Dbl): $CMD"
357       $CMD
358     fi
359  #-- after running Dbl:  #-- after running Dbl:
360  if test $ksel = 1  #-  check for pickup
361  then   listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
362     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
363     if test "x$listOc" = x -o "x$listAt" = x
364     then echo "==> RUN 2 x $NiOc it STOP without writing $pickRef => exit" ; exit 1
365     else echo "==>  END  RUN 2 x $NiOc it" ; fi
366    #-  move output
367   rm -rf $dir1   rm -rf $dir1
368   mkdir $dir1   mkdir $dir1
369   echo ' ' ; echo 'move STDOUT to' $dir1   echo '=> move STDOUT to' $dir1
370   mv rank_1/STDOUT* $dir1/ocnSTDOUT.2it   mv std_outp rank_0/*.0000.clog rank_$rnkO/*.0001.clog rank_$rnkA/*.0001.clog $dir1
371   mv rank_2/STDOUT* $dir1/atmSTDOUT.2it   mv rank_$rnkO/STDOUT.0000 $dir1/ocnSTDOUT.0000
372     mv rank_$rnkA/STDOUT.0000 $dir1/atmSTDOUT.0000
373    
374   mkdir $dir1/ocn   mkdir $dir1/ocn
375   cd rank_1   cd rank_$rnkO
376   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
377   do   echo ' listP=' $listP
378    echo $RNPref $xx.ckptA $xx.$noc2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
379    $RNPref $xx.ckptA $xx.$noc2c   echo 'move output to' $dir1/ocn
380   done   move_outp  ../$dir1/ocn
  echo 'move pickups to' $dir1/ocn  
  mv pickup*.$noc2c.* ../$dir1/ocn  
381   cd ..   cd ..
382    
383   mkdir $dir1/atm   mkdir $dir1/atm
384   cd rank_2   cd rank_$rnkA
385   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
386   do   echo ' listP=' $listP
387    echo $RNPref $xx.ckptA $xx.$nat2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
388    $RNPref $xx.ckptA $xx.$nat2c   echo 'move output to' $dir1/atm
389   done   move_outp  ../$dir1/atm
  echo 'move pickups to' $dir1/atm  
  mv pickup*.$nat2c.* ../$dir1/atm  
390   cd ..   cd ..
391    
392  #-- prepare for running 1iA  #-- prepare for running 1iA
393   echo ' prepare for running 1iA:'   echo '=> prepare for running 1iA:'
394   Nit=$NiOc   nit=$NiOc ; Dbl=`expr $nit + $nit`
395   Dbl=`expr $Nit + $Nit`   sed "1 s/[0-9]*/$nit/" input_cpl/data > rank_0/data
396   sed "1 s/[0-9]*/$Nit/" input_cpl/data > rank_0/data   pFile=rank_$rnkO/data
397   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" input_ocn/data.tst > rank_1/data   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
398   Nit=$NiAt   nit=$NiAt ; Dbl=`expr $nit + $nit`
399   Dbl=`expr $Nit + $Nit`   pFile=rank_$rnkA/data
400   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" input_atm/data.tst > rank_2/data   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
401   echo 'rank_0/data:' `head -1 rank_0/data`   echo 'rank_0/data:' `head -1 rank_0/data`
402   egrep 'nIter0|nTimeSteps' rank_?/data   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
403  fi  fi
404    
405    #-- run & post-process 1iA:
406    if [ $ksel -eq 2 -o $kUp -le 2 ] ; then
407     echo ' '
408    #-- run coupled test (1iA):
409     if [ $doRun -eq 1 -o $kUp -lt 2 ] ; then
410       echo "=> run (1iA): $CMD"
411       $CMD
412     fi
413  #-- after running 1iA  #-- after running 1iA
414  if test $ksel = 2  #-  check for pickup
415  then   listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
416     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
417     if test "x$listOc" = x -o "x$listAt" = x
418     then echo "==> RUN 1iA STOP without writing $pickRef => exit" ; exit 2
419     else echo "==>  END  RUN 1iA" ; fi
420    #-  move output
421   rm -rf $dir2   rm -rf $dir2
422   mkdir $dir2   mkdir $dir2
423   echo ' ' ; echo 'move STDOUT to' $dir2   echo '=> move STDOUT to' $dir2
424   mv rank_1/STDOUT* $dir2/ocnSTDOUT.1iA   mv std_outp rank_0/*.clog rank_$rnkO/*.clog rank_$rnkA/*.clog $dir2
425   mv rank_2/STDOUT* $dir2/atmSTDOUT.1iA   mv rank_$rnkO/STDOUT.0000 $dir2/ocnSTDOUT.0000
426     mv rank_$rnkA/STDOUT.0000 $dir2/atmSTDOUT.0000
427    
428   mkdir $dir2/ocn   mkdir $dir2/ocn
429   cd rank_1   cd rank_$rnkO
430   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
431   do   echo ' listP=' $listP
432    echo $RNPref $xx.ckptA $xx.$noc1c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc1c ; done
433    $RNPref $xx.ckptA $xx.$noc1c  # remove pickup*.$noc0c* sym-link:
434   done   find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*' -exec rm {} \;
435   echo 'move pickups to' $dir2/ocn   echo 'move output to' $dir2/ocn
436   mv pickup*.$noc1c.* ../$dir2/ocn   move_outp  ../$dir2/ocn
437   cd ..   cd ..
438    
439   mkdir $dir2/atm   mkdir $dir2/atm
440   cd rank_2   cd rank_$rnkA
441   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
442   do   echo ' listP=' $listP
443    echo $RNPref $xx.ckptA $xx.$nat1c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat1c ; done
444    $RNPref $xx.ckptA $xx.$nat1c  # remove pickup*.$nat0c* sym-link:
445   done   find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*' -exec rm {} \;
446   echo 'move pickups to' $dir2/atm   echo 'move output to' $dir2/atm
447   mv pickup*.$nat1c.* ../$dir2/atm   move_outp  ../$dir2/atm
448   cd ..   cd ..
449    
450  #-- prepare for running 1iB  #-- prepare for running 1iB
451   echo ' prepare for running 1iB:'   echo '=> prepare for running 1iB:'
452   cd rank_1   cd rank_$rnkO
453   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp   ln -s ../$dir2/ocn/pickup*.$noc1c.* .
454   mv -f data.tmp data   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp_$$
455   ln -s ../$dir2/ocn/* .   mv -f data.tmp_$$ data
456   cd ../rank_2   cd ../rank_$rnkA
457   sed "s/^ nIter0=$nat0/ nIter0=$nat1/g" data > data.tmp   ln -s ../$dir2/atm/pickup*.$nat1c.* .
458   mv -f data.tmp data   sed "s/^ nIter0=$nat0/ nIter0=$nat1/g" data > data.tmp_$$
459   ln -s ../$dir2/atm/* .   mv -f data.tmp_$$ data
460   cd ..   cd ..
461   echo 'rank_0/data:' `head -1 rank_0/data`   echo 'rank_0/data:' `head -1 rank_0/data`
462   egrep 'nIter0|nTimeSteps' rank_?/data   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
463  fi  fi
464    
465    #-- run & post-process 1iB:
466    if [ $ksel -eq 3 -o $kUp -le 3 ] ; then
467     echo ' '
468    #-- run coupled test (1iB):
469     if [ $doRun -eq 1 -o $kUp -lt 3 ] ; then
470       echo "=> run (1iB): $CMD"
471       $CMD
472     fi
473  #-- after running 1iB  #-- after running 1iB
474  if test $ksel = 3  #-  check for pickup
475  then   listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
476     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
477     if test "x$listOc" = x -o "x$listAt" = x
478     then echo "==> RUN 1iB STOP without writing $pickRef => exit" ; exit 3
479     else echo "==>  END  RUN 1iB" ; fi
480    #-  move output
481   rm -rf $dir3   rm -rf $dir3
482   mkdir $dir3   mkdir $dir3
483   echo ' ' ; echo 'move STDOUT to' $dir3   echo '=> move STDOUT to' $dir3
484   mv rank_1/STDOUT* $dir3/ocnSTDOUT.1iB   mv std_outp rank_0/*.clog rank_$rnkO/*.clog rank_$rnkA/*.clog $dir3
485   mv rank_2/STDOUT* $dir3/atmSTDOUT.1iB   mv rank_$rnkO/STDOUT.0000 $dir3/ocnSTDOUT.0000
486     mv rank_$rnkA/STDOUT.0000 $dir3/atmSTDOUT.0000
487    
488   mkdir $dir3/ocn   mkdir $dir3/ocn
489   cd rank_1   cd rank_$rnkO
490   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
491   do   echo ' listP=' $listP
492    echo $RNPref $xx.ckptA $xx.$noc2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
493    $RNPref $xx.ckptA $xx.$noc2c  # remove pickup*.$noc1c* sym-link:
494   done   find . ! -name . -prune -type l -name 'pickup*.'${noc1c}'*' -exec rm {} \;
495   echo 'move pickups to' $dir3/ocn   echo 'move output to' $dir3/ocn
496   mv pickup*.$noc2c.* ../$dir3/ocn   move_outp  ../$dir3/ocn
497   cd ..   cd ..
498    
499   mkdir $dir3/atm   mkdir $dir3/atm
500   cd rank_2   cd rank_$rnkA
501   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
502   do   echo ' listP=' $listP
503    echo $RNPref $xx.ckptA $xx.$nat2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
504    $RNPref $xx.ckptA $xx.$nat2c  # remove pickup*.$nat1c* sym-link:
505   done   find . ! -name . -prune -type l -name 'pickup*.'${nat1c}'*' -exec rm {} \;
506   echo 'move pickups to' $dir3/atm   echo 'move output to' $dir3/atm
507   mv pickup*.$nat2c.* ../$dir3/atm   move_outp  ../$dir3/atm
508   cd ..   cd ..
509   echo 'restore std data files'  
  cp -p input_cpl/data rank_0  
  cp -p input_ocn/data rank_1  
  cp -p input_atm/data rank_2  
510  fi  fi
511    
512  #-- compare output:  #-- compare output:
513  if test $ksel = 4  if test $ksel = 4 ; then
 then  
514   rm -f ocn.txt   rm -f ocn.txt
515   echo ' run 2it: ---------- ' > tmp1.txt   echo ' run 2it: ---------- ' > tmp1.txt
516   grep "cg2d_init_res" $dir1/ocnSTDOUT.2it \   grep "cg2d_init_res" $dir1/ocnSTDOUT.0000 \
517        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
518   Nit=`expr $NiOc + 2`   nit=`expr $NiOc + 2`
519   sed "$Nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt   #echo "NiOc= $NiOc ; nit= $nit"
520     sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
521   rm -f tmp1.txt   rm -f tmp1.txt
522   echo ' run 1iA: ---------- ' > tmp1.txt   echo ' run 1iA: ---------- ' > tmp1.txt
523   grep "cg2d_init_res" $dir2/ocnSTDOUT.1iA \   grep "cg2d_init_res" $dir2/ocnSTDOUT.0000 \
524        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
525   echo ' run 1iB: ---------- ' >> tmp1.txt   echo ' run 1iB: ---------- ' >> tmp1.txt
526   grep "cg2d_init_res" $dir3/ocnSTDOUT.1iB \   grep "cg2d_init_res" $dir3/ocnSTDOUT.0000 \
527        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
528   nl tmp1.txt > tmp1.ttt   nl tmp1.txt > tmp1.ttt
529   nl tmp2.txt > tmp2.ttt   nl tmp2.txt > tmp2.ttt
# Line 242  then Line 534  then
534   exit   exit
535  fi  fi
536    
537  if test $ksel = 5  if test $ksel = 5 ; then
 then  
538   rm -f atm.txt   rm -f atm.txt
539   echo ' run 2it: ---------- ' > tmp1.txt   echo ' run 2it: ---------- ' > tmp1.txt
540   grep "cg2d_init_res" $dir1/atmSTDOUT.2it \   grep "cg2d_init_res" $dir1/atmSTDOUT.0000 \
541        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
542   Nit=`expr $NiAt + 2`   nit=`expr $NiAt + 2`
543   sed "$Nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt   sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
544   rm -f tmp1.txt   rm -f tmp1.txt
545   echo ' run 1iA: ---------- ' > tmp1.txt   echo ' run 1iA: ---------- ' > tmp1.txt
546   grep "cg2d_init_res" $dir2/atmSTDOUT.1iA \   grep "cg2d_init_res" $dir2/atmSTDOUT.0000 \
547        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
548   echo ' run 1iB: ---------- ' >> tmp1.txt   echo ' run 1iB: ---------- ' >> tmp1.txt
549   grep "cg2d_init_res" $dir3/atmSTDOUT.1iB \   grep "cg2d_init_res" $dir3/atmSTDOUT.0000 \
550        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
551   nl tmp1.txt > tmp1.ttt   nl tmp1.txt > tmp1.ttt
552   nl tmp2.txt > tmp2.ttt   nl tmp2.txt > tmp2.ttt
# Line 267  then Line 558  then
558  fi  fi
559    
560  #- compare pickup file:  #- compare pickup file:
561  if test $ksel = 6  if test $ksel = 6 ; then
 then  
562   echo ' ' ; echo '=== compare Ocn pickups: =============='   echo ' ' ; echo '=== compare Ocn pickups: =============='
563   for xx in $listOc   listP=`( cd $dir1/ocn ; ls pickup*.$noc2c* | sed 's/\..*//g' | uniq )`
564     for xx in $listP
565   do   do
566     $difex $dir1/ocn $dir3/ocn $xx.$noc2c     dif_tiles $dir1/ocn $dir3/ocn $xx.$noc2c
567   done   done
568   echo ' ' ; echo '=== compare Atm pickups: =============='   echo ' ' ; echo '=== compare Atm pickups: =============='
569   for xx in $listAt   listP=`( cd $dir1/atm ; ls pickup*.$nat2c* | sed 's/\..*//g' | uniq )`
570     for xx in $listP
571   do   do
572     $difex $dir1/atm $dir3/atm $xx.$nat2c     dif_tiles $dir1/atm $dir3/atm $xx.$nat2c
573   done   done
574   exit   exit
575  fi  fi
576    
577    #- clean-up and reset:
578    if test $ksel = 7 ; then
579     echo "== clean-up and reset :"
580     cd rank_$rnkO
581     listS=`find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*'`
582     if test "x$listS" != x ; then
583       echo '--> remove pickup*.'$noc0c'* sym-links:'
584       if [ $prt -ge 2 ] ; then echo " in rank_$rnkO: rm" $listS ; fi
585       rm -f $listS
586     fi
587     cd ../rank_$rnkA
588     listS=`find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*'`
589     if test "x$listS" != x ; then
590       echo '--> remove pickup*.'$nat0c'* sym-links:'
591       if [ $prt -ge 2 ] ; then echo " in rank_$rnkA: rm" $listS ; fi
592       rm -f $listS
593     fi
594     cd ..
595    #- move back files from temp_tst dir:
596     echo '=> restore std data files'
597     rm -f rank_{0,$rnkO,$rnkA}/data
598     rm -f rank_{$rnkO,$rnkA}/data.tst
599     rm -f rank_{$rnkO,$rnkA}/data.pkg
600     mv $tmpDir.0/*Ocean.0001.clog $tmpDir.1/* rank_$rnkO
601     mv $tmpDir.0/*Atmos.0001.clog $tmpDir.2/* rank_$rnkA
602     mv $tmpDir.0/std_outp .
603     mv $tmpDir.0/* rank_0
604     rmdir $tmpDir.0 $tmpDir.1 $tmpDir.2
605    #echo '--> remove output dir and output files:'
606     echo ' rm -rf' $dir1 $dir2 $dir3
607     rm -rf $dir1 $dir2 $dir3
608     rm -f atm.txt ocn.txt
609     exit
610    fi

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.11

  ViewVC Help
Powered by ViewVC 1.1.22