/[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.5 by jmc, Mon Nov 26 23:46:21 2007 UTC revision 1.15 by jmc, Thu Mar 17 20:01:59 2022 UTC
# Line 1  Line 1 
1  #!/bin/sh  #! /usr/bin/env bash
2    
3  # $Header$  # $Header$
4  # $Name$  # $Name$
# 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 22  dif_tiles() Line 23  dif_tiles()
23  #   dif_data_loc dd1 dd2 file  #   dif_data_loc dd1 dd2 file
24  #   for all tiles  ==> do diff dd1/file dd2/file'  #   for all tiles  ==> do diff dd1/file dd2/file'
25      dd1=$1 ; dd2=$2 ; zz=$3      dd1=$1 ; dd2=$2 ; zz=$3
26      chkmeta='Y'      chkmeta='N'
27        listY=`(cd $dd1 ; ls $zz.data | sed "s/\.data//") 2> /dev/null`
28        if test "x$listY" = "x$zz" ; then
29          echo '--> file='$zz' (global)'
30          echo diff $dd1/$zz.data $dd2
31          diff $dd1/$zz.data $dd2
32          out=$?
33          if test $out != 0
34          then echo 'Diff outp=' $out ' ==> stop'
35           exit ; fi
36          if test $chkmeta = 'Y' ; then
37           echo diff $dd1/$zz.meta $dd2
38           diff $dd1/$zz.meta $dd2
39           out=$?
40           if test $out != 0 ; then echo 'Diff outp=' $out ' ==> stop' ; exit ; fi
41          fi
42          return
43        fi
44      listY=`(cd $dd1 ; ls $zz.*.data | sed "s/$zz\.//" | sed "s/\.data//")`      listY=`(cd $dd1 ; ls $zz.*.data | sed "s/$zz\.//" | sed "s/\.data//")`
45      echo '--> file='$zz', listY='$listY      echo '--> file='$zz', listY='$listY
46      for yy in $listY      for yy in $listY
# Line 33  dif_tiles() Line 51  dif_tiles()
51        if test $out != 0        if test $out != 0
52        then echo 'Diff outp=' $out ' ==> stop'        then echo 'Diff outp=' $out ' ==> stop'
53         exit ; fi         exit ; fi
54        if test $chkmeta = 'Y'        if test $chkmeta = 'Y' ; then
       then  
55         echo diff $dd1/$zz.$yy.meta $dd2         echo diff $dd1/$zz.$yy.meta $dd2
56         diff $dd1/$zz.$yy.meta $dd2         diff $dd1/$zz.$yy.meta $dd2
57         out=$?         out=$?
58         if test $out != 0         if test $out != 0 ; then echo 'Diff outp=' $out ' ==> stop' ; exit ; fi
        then echo 'Diff outp=' $out ' ==> stop'  
         exit ; fi  
59        fi        fi
60      done      done
61  }  }
62    
63  if test -f input_ocn/data.tst  move_outp()
64  then  {
65   noc0=`sed -n 's/nIter0=//p' input_ocn/data.tst | sed 's/,//g' | sed 's/ //g'`  #   move_outp arg_1
66   Dbl=`sed -n 's/nTimeSteps=//p' input_ocn/data.tst | sed 's/,//g'`  #   move mitgcm output files to directory 'arg_1'
67   NiOc=`expr $Dbl / 2`      if [ $prt -ge 1 ] ; then echo ' move_outp:' $1 ; fi
68   Nit=$NiOc      mv *.data *.meta $1
69  else      listF=`ls -1 *.txt STD???.0000 *.log 2> /dev/null`
70   echo " file: 'input_ocn/data.tst' not found"      if test "x$listF" != x ; then mv $listF $1 ; fi
71   NiOc=0 ; Nit=2      #- move back sym link:
72        listL=`find $1 -type l`
73        ttd=`echo $listL | wc -w`
74        #echo "ttd=$ttd"
75        if test $ttd != 0 ; then
76          #echo "listL='$listL'"
77          shortList=`ls $listL | sed "s|$1/||"`
78          if [ $prt -ge 2 ] ; then echo "  mv back:" $shortList ; fi
79          mv $listL .
80        fi
81    }
82    
83    #----------------------------------------------------------------------
84    
85    prepare_first()
86    {
87    #   prepare_first arg_1 arg_2 arg_3
88    #-----------------------------------------------------------
89    #   from a previous run in "arg_2" dir: collect info for
90    #     test: arg_1 (iters) + arg_1 (iter) = 2x arg_1 (iter) ;
91    #     move output files from "arg_2" dir to 'arg_3" dir ;
92    #     create "data.tst" in arg_2 dir ;
93    #     and prepare pickup and params file for 1rst run.
94    #-----------------------------------------------------------
95    
96      nIt=$1 ; dFr=$2 ; dTo=$3
97      if [ $prt -ge 1 ] ; then
98          echo " run 'prepare_first' : nIt=$nIt, from '$dFr' to '$dTo'" ; fi
99      outpRun="$dFr/STDOUT.0000"
100    #- last pickup written (from previous run):
101      if test -r $outpRun
102      then
103    #   grep ' %CHECKPOINT ' $outpRun | tail -1
104        lastPick=`grep ' %CHECKPOINT ' $outpRun | sed 's/^.*%CHECKPOINT //' | tail -1`
105        if test "x$lastPick" = x
106        then
107         echo "no pickup information written in '$outpRun' from previous run"
108         exit 8
109        else
110         iter=`echo $lastPick | awk '{print $1}'`
111         sufx=`echo $lastPick | awk '{print $2}'`
112         echo " from '$outpRun' : lastPick='$lastPick' > iter='$iter', sufx='$sufx'"
113        fi
114      else
115        echo "no standard output file '$outpRun' from previous run"
116        exit 7
117      fi
118    #- create parameter file "data.tst":
119      echo " preparing parameter file 'data.tst' :"
120      pFile="$dFr/data.tst"
121      if test -h $pFile ; then rm -f $pFile ; fi
122      cp -p -f $dFr/data $pFile
123      listOut='startTime endTime nIter0 nTimeSteps nEndIter chkPtFreq pChkPtFreq'
124      listOut=$listOut' pickupSuff usePickupBeforeC54 mdsioLocalDir pickupStrictlyMatch'
125    #- comment out (in data.tst) all parameters from list: 'listOut'
126      for yy in $listOut
127      do
128        if [ $prt -ge 2 ] ; then echo -n '  commenting out:' $yy ; fi
129        nn=`grep -i -c "^ *$yy *=" $pFile`
130        if [ $nn -ge 1 ]
131        then
132          listN=`grep -i -n "^ *$yy *=" $pFile | sed 's/:.*$//'`
133          if [ $prt -ge 2 ] ; then echo " : found on line: $listN" ; fi
134          for ll in $listN
135          do
136            sed "$ll s/^./#&/" $pFile > data.tmp_$$
137           #diff data.tmp_$$ $pFile
138            mv -f data.tmp_$$ $pFile
139          done
140        else if [ $prt -ge 2 ] ; then echo ' : not found' ; fi
141        fi
142      done
143      if [ $prt -ge 2 ] ; then echo ' ' ; fi
144    # add nIter0 & nTimeSteps in namelist "PARM03":
145      Dbl=`expr $nIt \* 2`
146      sed "/^ *\&PARM03/a \ nTimeSteps=$Dbl," $pFile > data.tmp_$$
147      mv -f data.tmp_$$ $pFile
148      sed "/^ *\&PARM03/a \ nIter0=$iter," $pFile > data.tmp_$$
149      mv -f data.tmp_$$ $pFile
150      echo "preparing file '$pFile' : done"
151      if  [ $prt -ge 1 ] ; then
152        echo " diff $pFile $dFr/data" ; diff $pFile $dFr/data ; echo ' ' ;
153      fi
154    #---- clean up current dir: move output to temp_tst dir:
155    if test $action = 1 ; then
156      ( cd $dFr ; move_outp ../$dTo )
157      listF=`ls -1 $dFr/STD???.0[0-9][0-9][0-9] 2> /dev/null`
158      if test "x$listF" != x ; then mv $listF $dTo ; fi
159    #- move main parameter file "data":
160      mv $dFr/data $dTo
161    #- do not deal with MNC pickup at all:
162      cp -p -f $dFr/data.pkg data.tmp_$$
163      mv -f $dFr/data.pkg $dTo
164      mv -f data.tmp_$$ $dFr/data.pkg
165      yy='useMNC'
166      if [ $prt -ge 2 ] ; then echo -n "  commenting out: $yy in data.pkg" ; fi
167      nn=`grep -i -c "^ *$yy *=" $dFr/data.pkg`
168      if [ $nn -ge 1 ] ; then
169        listN=`grep -i -n "^ *$yy *=" $dFr/data.pkg | sed 's/:.*$//'`
170        if [ $prt -ge 2 ] ; then echo " : found on line: $listN" ; fi
171        for ll in $listN
172        do
173          sed "$ll s/^./#&/" $dFr/data.pkg > data.tmp_$$
174          if  [ $prt -ge 1 ] ; then echo ' changes in data.pkg'
175            diff data.tmp_$$ $dFr/data.pkg ; echo ' ' ; fi
176          mv -f data.tmp_$$ $dFr/data.pkg
177        done
178      else if [ $prt -ge 2 ] ; then echo ' : not found' ; echo ' ' ; fi
179      fi
180    #-- link back pickup files
181      listS=`ls $dTo/pickup*.${sufx}.* 2> /dev/null | head -1`
182      if test "x$listS" != x ; then
183        if [ $prt -ge 1 ] ; then echo ' link back:' ../$dTo'/pickup*.'${sufx}'*' ; fi
184        ( cd $dFr ; ln -s ../$dTo/pickup*.${sufx}* . )
185      fi
186    #- rename pickup (if needed)
187      nIt0c=`printf "%10.10i\n" $iter`
188      if test $sufx = $nIt0c
189      then
190        echo "no need to change pickup sufix"
191      else
192        listP=`(cd $dFr ; ls pickup*.${sufx}* 2> /dev/null | sed 's/\..*//g' | uniq )`
193        if test "x$listP" = x
194        then echo ' rename' $sufx  '->' $nIt0c ': no file to rename'
195        else
196         if [ $prt -ge 1 ] ; then echo ' rename' $sufx  '->' $nIt0c 'for all:' $listP ; fi
197         ( cd $dFr ; for yy in $listP
198           do rnp_loc $yy.${sufx} $yy.$nIt0c
199           done )
200        fi
201      fi
202  fi  fi
203  if test -f input_atm/data.tst  #----------------------------------------------------------------------
204  then  }
205   nat0=`sed -n 's/nIter0=//p' input_atm/data.tst | sed 's/,//g' | sed 's/ //g'`  
206   Dbl=`sed -n 's/nTimeSteps=//p' input_atm/data.tst | sed 's/,//g'`  #===============================================================================
207   NiAt=`expr $Dbl / 2`  #-- sequential part start here
208    #===============================================================================
209    
210    #- default:
211    Nit=2
212    prt=1    # 1 = print some informations ; 2 = print more informations
213    action=1
214    tmpDir='temp_tst'
215    
216    pickRef='pickup'
217    # offLine=`grep -i -c '^ *useOffLine *= *\.TRUE\.' data.pkg`
218    # if [ $offLine -ge 1 ] ; then pickRef='pickup_ptracers' ; fi
219    
220    #- dir where to put the results :
221    dir1=res_2it
222    dir2=res_1iA
223    dir3=res_1iB
224    
225    #- get number of procs from file "pr_group":
226    if test -f pr_group ; then
227     Npr=`cat pr_group | wc -l`
228     NpOc=`grep -c '\/build_ocn\/mitgcmuv' pr_group`
229    else echo " file: 'pr_group' not found"
230     Npr=0
231    fi
232    
233    #- parse options:
234    CMD='run_cpl_test'
235    MTH=
236    argList=$*
237    for xx in $argList
238    do
239     if test $xx = '-mth' ; then MTH=$xx ; shift ; continue ; fi
240     if test -x ${CMD}.$xx ; then CMD="${CMD}.$xx" ; shift ; continue ; fi
241    done
242    
243    if [ $# -ne 1 ] ; then
244       echo "Usage: `basename $0` [opt] flag"
245       echo " Check restart of coupled set-up: compare 1 run of 2 x $Nit it long"
246       echo "                       with 2 consecutive runs of $Nit it long each"
247       echo "opt = -mth : run script '$CMD' with option '-mth'"
248       echo "opt = sufx : run script 'run_cpl_test.{sufx}'"
249       echo " where: flag = 0  -> prepare for 1rst run"
250       echo "        flag = r1 -> do 1rst run and step 1"
251       echo "        flag =  1 -> move res. after 1rst run & prepare for 2nd"
252       echo "        flag = r2 -> do 2nd  run and step 2"
253       echo "        flag =  2 -> move res. after 2nd  run & prepare for 3rd"
254       echo "        flag = r3 -> do 3rd  run and step 3"
255       echo "        flag =  3 -> move res. after 3rd  run"
256       echo "        flag=(r)N+ , N=1,2,3 => do steps (+run) from: (r)N to 3"
257       echo "        flag = 4  -> compare Ocean results"
258       echo "        flag = 5  -> compare Atmos results"
259       echo "        flag = 6  -> diff pickup files"
260       echo "        flag = 7  -> clean-up output files"
261       exit
262    fi
263    arg1=$1
264    xx=`echo $arg1 | sed 's/\+$//'`
265    if test $xx != $arg1
266    then kUp=`echo $xx | sed 's/^r//'` ; arg1=$xx
267    else kUp=9 ; fi
268    xx=`echo $arg1 | sed 's/^r//'`
269    if test $xx != $arg1 ; then doRun=1 ; arg1=$xx ; else doRun=0 ; fi
270    ksel=$arg1
271    
272    if test -x ${CMD} ; then
273      CMD="./${CMD} $MTH 3"
274  else  else
275   echo " file: 'input_atm/data.tst' not found"    if test -x ../../tools/${CMD} ; then
276   NiAt=0      CMD="../../tools/${CMD} $MTH 3"
277      else
278        echo " no test-script '${CMD}' found neither here nor in '../../tools/'"
279        exit
280      fi
281  fi  fi
282    echo " doRun='$doRun' ; ksel='$ksel' ; kUp='$kUp' ; command='$CMD'"
283    
284  case $# in  #- set dir of Master-Proc for 3 components (coupler == rank_0):
285    1) ;;  if test $Npr = 0 ; then
286     *) echo "Usage: `basename $0` flag"      echo " needs file: 'pr_group' to continue" ; exit
       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"  
       exit ;;  
 esac  
 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  
287  fi  fi
288  ksel=$1  rnkO=1 ; rnkA=`expr $rnkO + $NpOc`
289    
290  #make sure that local dir & ~jmc/bin is in the path:  #make sure that local dir & ~jmc/bin is in the path:
291  #export PATH=${PATH}:.  #export PATH=${PATH}:.
292    
293  Nit=$NiOc  #----------------------------------------------------------------------
294  Dbl=`expr $Nit + $Nit`  
295  noc1=`expr $noc0 + $Nit`  if test $ksel = 0 ; then
296    #-- from STDOUT of previous run, get Nb of iters for each component,
297      yy='(PID.TID 0000.0001) nTimeSteps = '
298      outpRun="rank_$rnkO/STDOUT.0000"
299      nTsO=`grep -A1 "$yy" $outpRun | tail -1 | sed 's/(PID.TID 0000.0001) *//'`
300      outpRun="rank_$rnkA/STDOUT.0000"
301      nTsA=`grep -A1 "$yy" $outpRun | tail -1 | sed 's/(PID.TID 0000.0001) *//'`
302    #-  then calculate the ratio and set Nb of time-step for Atmos comp.
303      RatioTs=`expr $nTsA / $nTsO`
304      if [ $prt -ge 1 ] ; then echo " nTsA= $nTsA , nTsO= $nTsO , Ratio= $RatioTs"; fi
305      NiOc=$Nit ;
306      NiAt=`expr $NiOc \* $RatioTs`
307    #---- clean up current dir: rename pickup files and move output to temp_tst dir:
308     /bin/rm -rf $tmpDir.0 $tmpDir.1 $tmpDir.2
309     mkdir $tmpDir.0 $tmpDir.1 $tmpDir.2
310     mv std_outp rank_0/*.clog $tmpDir.0
311     mv rank_$rnkO/*.0001.clog $tmpDir.0
312     mv rank_$rnkA/*.0001.clog $tmpDir.0
313     echo ' '
314     echo "=> prepare OCN dir:"
315     prepare_first $NiOc rank_$rnkO $tmpDir.1
316     echo "=> prepare ATM dir:"
317     prepare_first $NiAt rank_$rnkA $tmpDir.2
318    #-- prepare for running Dbl
319     echo '=> prepare for running Dbl:'
320     Dbl=`expr $NiOc + $NiOc`
321     cp -p rank_$rnkO/data.tst rank_$rnkO/data
322     cp -p rank_$rnkA/data.tst rank_$rnkA/data
323     egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
324    fi
325    
326    #----------------------------------------------------------------------
327    
328    #- extract length of runs & starting iter number from file "data.tst"
329    #  Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)
330    pFile="rank_$rnkO/data.tst"
331    if test -f $pFile ; then
332      noc0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
333      Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
334      NiOc=`expr $Dbl / 2`
335      Nit=$NiOc
336    else echo " file: '$pFile' not found"
337      NiOc=0 ; Nit=2
338    fi
339    pFile="rank_$rnkA/data.tst"
340    if test -f $pFile ; then
341      nat0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
342      Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
343      NiAt=`expr $Dbl / 2`
344    else echo " file: '$pFile' not found"
345      NiAt=0
346    fi
347    if test $NiAt = 0 -o $NiOc = 0 ; then
348        echo " needs 2 data files: 'rank_$rnkO/data.tst' & 'rank_$rnkA/data.tst'"
349        echo "  (corresponding to 1rst run) to continue"
350      exit
351    fi
352    
353    nit=$NiOc
354    Dbl=`expr $nit + $nit`
355    noc1=`expr $noc0 + $nit`
356  noc2=`expr $noc0 + $Dbl`  noc2=`expr $noc0 + $Dbl`
357  noc0c=`printf "%10.10i\n" $noc0`  noc0c=`printf "%10.10i\n" $noc0`
358  noc1c=`printf "%10.10i\n" $noc1`  noc1c=`printf "%10.10i\n" $noc1`
359  noc2c=`printf "%10.10i\n" $noc2`  noc2c=`printf "%10.10i\n" $noc2`
360  echo 'Oce (rank_1):' $noc0  $noc1  $noc2  echo " Oce (rank_$rnkO):" $noc0  $noc1  $noc2
361  echo $noc0c $noc1c $noc2c  echo "     " $noc0c $noc1c $noc2c
362    
363  Nit=$NiAt  nit=$NiAt
364  Dbl=`expr $Nit + $Nit`  Dbl=`expr $nit + $nit`
365  nat1=`expr $nat0 + $Nit`  nat1=`expr $nat0 + $nit`
366  nat2=`expr $nat0 + $Dbl`  nat2=`expr $nat0 + $Dbl`
367  nat0c=`printf "%10.10i\n" $nat0`  nat0c=`printf "%10.10i\n" $nat0`
368  nat1c=`printf "%10.10i\n" $nat1`  nat1c=`printf "%10.10i\n" $nat1`
369  nat2c=`printf "%10.10i\n" $nat2`  nat2c=`printf "%10.10i\n" $nat2`
370  echo 'Atm (rank_2):' $nat0  $nat1  $nat2  echo " Atm (rank_$rnkA):" $nat0  $nat1  $nat2
371  echo $nat0c $nat1c $nat2c  echo "     " $nat0c $nat1c $nat2c
   
 listOc="pickup pickup_cpl"  
 listAt="pickup pickup_cpl pickup_ic pickup_land"  
372    
373  #- dir where to put the results :  #----------------------------------------------------------------------
 dir1=res_2it  
 dir2=res_1iA  
 dir3=res_1iB  
374    
375  #-- after a run, rename pickup file to restart  #-- run & post-process Dbl:
376  if test $ksel = 0  if [ $ksel -eq 1 -o $kUp -le 1 ] ; then
 then  
377   echo ' '   echo ' '
378   cd rank_1  #-- run coupled test (Dbl):
379   for xx in $listOc   if [ $doRun -eq 1 -o $kUp -lt 1 ] ; then
380   do rnp_loc $xx.ckptA $xx.$noc0c     echo "=> run (Dbl): $CMD"
381   done     $CMD
382   cd ../rank_2   fi
  for xx in $listAt  
  do rnp_loc $xx.ckptA $xx.$nat0c  
  done  
  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_1/data  
  cp -p input_atm/data.tst rank_2/data  
  echo 'rank_0/data:' `head -1 rank_0/data`  
  egrep 'nIter0|nTimeSteps' rank_?/data  
 fi  
   
383  #-- after running Dbl:  #-- after running Dbl:
384  if test $ksel = 1  #-  check for pickup
385  then   listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
386     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
387     if test "x$listOc" = x -o "x$listAt" = x
388     then echo "==> RUN 2 x $NiOc it STOP without writing $pickRef => exit" ; exit 1
389     else echo "==>  END  RUN 2 x $NiOc it" ; fi
390    #-  move output
391   rm -rf $dir1   rm -rf $dir1
392   mkdir $dir1   mkdir $dir1
393   echo ' ' ; echo 'move STDOUT to' $dir1   echo '=> move STDOUT to' $dir1
394   mv rank_1/STDOUT* $dir1/ocnSTDOUT.2it   mv std_outp rank_0/*.0000.clog rank_$rnkO/*.0001.clog rank_$rnkA/*.0001.clog $dir1
395   mv rank_2/STDOUT* $dir1/atmSTDOUT.2it   mv rank_$rnkO/STDOUT.0000 $dir1/ocnSTDOUT.0000
396     mv rank_$rnkA/STDOUT.0000 $dir1/atmSTDOUT.0000
397    
398   mkdir $dir1/ocn   mkdir $dir1/ocn
399   cd rank_1   cd rank_$rnkO
400   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
401   do   echo ' listP=' $listP
402    rnp_loc $xx.ckptA $xx.$noc2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
403   done   echo 'move output to' $dir1/ocn
404   echo 'move pickups to' $dir1/ocn   move_outp  ../$dir1/ocn
  mv pickup*.$noc2c.* ../$dir1/ocn  
405   cd ..   cd ..
406    
407   mkdir $dir1/atm   mkdir $dir1/atm
408   cd rank_2   cd rank_$rnkA
409   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
410   do   echo ' listP=' $listP
411    rnp_loc $xx.ckptA $xx.$nat2c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
412   done   echo 'move output to' $dir1/atm
413   echo 'move pickups to' $dir1/atm   move_outp  ../$dir1/atm
  mv pickup*.$nat2c.* ../$dir1/atm  
414   cd ..   cd ..
415    
416  #-- prepare for running 1iA  #-- prepare for running 1iA
417   echo ' prepare for running 1iA:'   echo '=> prepare for running 1iA:'
418   Nit=$NiOc   nit=$NiOc ; Dbl=`expr $nit + $nit`
419   Dbl=`expr $Nit + $Nit`   pFile=rank_$rnkO/data
420   sed "1 s/[0-9]*/$Nit/" input_cpl/data > rank_0/data   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
421   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" input_ocn/data.tst > rank_1/data   nit=$NiAt ; Dbl=`expr $nit + $nit`
422   Nit=$NiAt   pFile=rank_$rnkA/data
423   Dbl=`expr $Nit + $Nit`   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
424   sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" input_atm/data.tst > rank_2/data   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
  echo 'rank_0/data:' `head -1 rank_0/data`  
  egrep 'nIter0|nTimeSteps' rank_?/data  
425  fi  fi
426    
427    #-- run & post-process 1iA:
428    if [ $ksel -eq 2 -o $kUp -le 2 ] ; then
429     echo ' '
430    #-- run coupled test (1iA):
431     if [ $doRun -eq 1 -o $kUp -lt 2 ] ; then
432       echo "=> run (1iA): $CMD"
433       $CMD
434     fi
435  #-- after running 1iA  #-- after running 1iA
436  if test $ksel = 2  #-  check for pickup
437  then   listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
438     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
439     if test "x$listOc" = x -o "x$listAt" = x
440     then echo "==> RUN 1iA STOP without writing $pickRef => exit" ; exit 2
441     else echo "==>  END  RUN 1iA" ; fi
442    #-  move output
443   rm -rf $dir2   rm -rf $dir2
444   mkdir $dir2   mkdir $dir2
445   echo ' ' ; echo 'move STDOUT to' $dir2   echo '=> move STDOUT to' $dir2
446   mv rank_1/STDOUT* $dir2/ocnSTDOUT.1iA   mv std_outp rank_0/*.clog rank_$rnkO/*.clog rank_$rnkA/*.clog $dir2
447   mv rank_2/STDOUT* $dir2/atmSTDOUT.1iA   mv rank_$rnkO/STDOUT.0000 $dir2/ocnSTDOUT.0000
448     mv rank_$rnkA/STDOUT.0000 $dir2/atmSTDOUT.0000
449    
450   mkdir $dir2/ocn   mkdir $dir2/ocn
451   cd rank_1   cd rank_$rnkO
452   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
453   do   echo ' listP=' $listP
454    rnp_loc $xx.ckptA $xx.$noc1c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc1c ; done
455   done  # remove pickup*.$noc0c* sym-link:
456   echo 'move pickups to' $dir2/ocn   find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*' -exec rm {} \;
457   mv pickup*.$noc1c.* ../$dir2/ocn   echo 'move output to' $dir2/ocn
458     move_outp  ../$dir2/ocn
459   cd ..   cd ..
460    
461   mkdir $dir2/atm   mkdir $dir2/atm
462   cd rank_2   cd rank_$rnkA
463   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
464   do   echo ' listP=' $listP
465    rnp_loc $xx.ckptA $xx.$nat1c   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat1c ; done
466   done  # remove pickup*.$nat0c* sym-link:
467   echo 'move pickups to' $dir2/atm   find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*' -exec rm {} \;
468   mv pickup*.$nat1c.* ../$dir2/atm   echo 'move output to' $dir2/atm
469     move_outp  ../$dir2/atm
470   cd ..   cd ..
471    
472  #-- prepare for running 1iB  #-- prepare for running 1iB
473   echo ' prepare for running 1iB:'   echo '=> prepare for running 1iB:'
474   cd rank_1   cd rank_$rnkO
475   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp   ln -s ../$dir2/ocn/pickup*.$noc1c.* .
476   mv -f data.tmp data   sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp_$$
477   ln -s ../$dir2/ocn/* .   mv -f data.tmp_$$ data
478   cd ../rank_2   cd ../rank_$rnkA
479   sed "s/^ nIter0=$nat0/ nIter0=$nat1/g" data > data.tmp   ln -s ../$dir2/atm/pickup*.$nat1c.* .
480   mv -f data.tmp data   sed "s/^ nIter0=$nat0/ nIter0=$nat1/g" data > data.tmp_$$
481   ln -s ../$dir2/atm/* .   mv -f data.tmp_$$ data
482   cd ..   cd ..
483   echo 'rank_0/data:' `head -1 rank_0/data`   egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
  egrep 'nIter0|nTimeSteps' rank_?/data  
484  fi  fi
485    
486    #-- run & post-process 1iB:
487    if [ $ksel -eq 3 -o $kUp -le 3 ] ; then
488     echo ' '
489    #-- run coupled test (1iB):
490     if [ $doRun -eq 1 -o $kUp -lt 3 ] ; then
491       echo "=> run (1iB): $CMD"
492       $CMD
493     fi
494  #-- after running 1iB  #-- after running 1iB
495  if test $ksel = 3  #-  check for pickup
496  then   listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
497     listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
498     if test "x$listOc" = x -o "x$listAt" = x
499     then echo "==> RUN 1iB STOP without writing $pickRef => exit" ; exit 3
500     else echo "==>  END  RUN 1iB" ; fi
501    #-  move output
502   rm -rf $dir3   rm -rf $dir3
503   mkdir $dir3   mkdir $dir3
504   echo ' ' ; echo 'move STDOUT to' $dir3   echo '=> move STDOUT to' $dir3
505   mv rank_1/STDOUT* $dir3/ocnSTDOUT.1iB   mv std_outp rank_0/*.clog rank_$rnkO/*.clog rank_$rnkA/*.clog $dir3
506   mv rank_2/STDOUT* $dir3/atmSTDOUT.1iB   mv rank_$rnkO/STDOUT.0000 $dir3/ocnSTDOUT.0000
507     mv rank_$rnkA/STDOUT.0000 $dir3/atmSTDOUT.0000
508    
509   mkdir $dir3/ocn   mkdir $dir3/ocn
510   cd rank_1   cd rank_$rnkO
511   for xx in $listOc   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
512   do rnp_loc $xx.ckptA $xx.$noc2c   echo ' listP=' $listP
513   done   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
514   echo 'move pickups to' $dir3/ocn  # remove pickup*.$noc1c* sym-link:
515   mv pickup*.$noc2c.* ../$dir3/ocn   find . ! -name . -prune -type l -name 'pickup*.'${noc1c}'*' -exec rm {} \;
516   rm -f pickup*.$noc1c.*   echo 'move output to' $dir3/ocn
517     move_outp  ../$dir3/ocn
518   cd ..   cd ..
519    
520   mkdir $dir3/atm   mkdir $dir3/atm
521   cd rank_2   cd rank_$rnkA
522   for xx in $listAt   listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
523   do rnp_loc $xx.ckptA $xx.$nat2c   echo ' listP=' $listP
524   done   for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
525   echo 'move pickups to' $dir3/atm  # remove pickup*.$nat1c* sym-link:
526   mv pickup*.$nat2c.* ../$dir3/atm   find . ! -name . -prune -type l -name 'pickup*.'${nat1c}'*' -exec rm {} \;
527   rm -f pickup*.$nat1c.*   echo 'move output to' $dir3/atm
528     move_outp  ../$dir3/atm
529   cd ..   cd ..
530   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  
531  fi  fi
532    
533  #-- compare output:  #-- compare output:
534  if test $ksel = 4  if test $ksel = 4 ; then
 then  
535   rm -f ocn.txt   rm -f ocn.txt
536   echo ' run 2it: ---------- ' > tmp1.txt   echo ' run 2it: ---------- ' > tmp1.txt
537   grep "cg2d_init_res" $dir1/ocnSTDOUT.2it \   grep "cg2d_init_res" $dir1/ocnSTDOUT.0000 \
538        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
539   Nit=`expr $NiOc + 2`   nit=`expr $NiOc + 2`
540   sed "$Nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt   #echo "NiOc= $NiOc ; nit= $nit"
541     sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
542   rm -f tmp1.txt   rm -f tmp1.txt
543   echo ' run 1iA: ---------- ' > tmp1.txt   echo ' run 1iA: ---------- ' > tmp1.txt
544   grep "cg2d_init_res" $dir2/ocnSTDOUT.1iA \   grep "cg2d_init_res" $dir2/ocnSTDOUT.0000 \
545        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
546   echo ' run 1iB: ---------- ' >> tmp1.txt   echo ' run 1iB: ---------- ' >> tmp1.txt
547   grep "cg2d_init_res" $dir3/ocnSTDOUT.1iB \   grep "cg2d_init_res" $dir3/ocnSTDOUT.0000 \
548        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
549   nl tmp1.txt > tmp1.ttt   nl tmp1.txt > tmp1.ttt
550   nl tmp2.txt > tmp2.ttt   nl tmp2.txt > tmp2.ttt
# Line 286  then Line 555  then
555   exit   exit
556  fi  fi
557    
558  if test $ksel = 5  if test $ksel = 5 ; then
 then  
559   rm -f atm.txt   rm -f atm.txt
560   echo ' run 2it: ---------- ' > tmp1.txt   echo ' run 2it: ---------- ' > tmp1.txt
561   grep "cg2d_init_res" $dir1/atmSTDOUT.2it \   grep "cg2d_init_res" $dir1/atmSTDOUT.0000 \
562        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
563   Nit=`expr $NiAt + 2`   nit=`expr $NiAt + 2`
564   sed "$Nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt   sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
565   rm -f tmp1.txt   rm -f tmp1.txt
566   echo ' run 1iA: ---------- ' > tmp1.txt   echo ' run 1iA: ---------- ' > tmp1.txt
567   grep "cg2d_init_res" $dir2/atmSTDOUT.1iA \   grep "cg2d_init_res" $dir2/atmSTDOUT.0000 \
568        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
569   echo ' run 1iB: ---------- ' >> tmp1.txt   echo ' run 1iB: ---------- ' >> tmp1.txt
570   grep "cg2d_init_res" $dir3/atmSTDOUT.1iB \   grep "cg2d_init_res" $dir3/atmSTDOUT.0000 \
571        | sed 's/.* cg2d_init_res =//' >> tmp1.txt        | sed 's/.* cg2d_init_res =//' >> tmp1.txt
572   nl tmp1.txt > tmp1.ttt   nl tmp1.txt > tmp1.ttt
573   nl tmp2.txt > tmp2.ttt   nl tmp2.txt > tmp2.ttt
# Line 311  then Line 579  then
579  fi  fi
580    
581  #- compare pickup file:  #- compare pickup file:
582  if test $ksel = 6  if test $ksel = 6 ; then
 then  
583   echo ' ' ; echo '=== compare Ocn pickups: =============='   echo ' ' ; echo '=== compare Ocn pickups: =============='
584   for xx in $listOc   listP=`( cd $dir1/ocn ; ls pickup*.$noc2c* | sed 's/\..*//g' | uniq )`
585     for xx in $listP
586   do   do
587     dif_tiles $dir1/ocn $dir3/ocn $xx.$noc2c     dif_tiles $dir1/ocn $dir3/ocn $xx.$noc2c
588   done   done
589   echo ' ' ; echo '=== compare Atm pickups: =============='   echo ' ' ; echo '=== compare Atm pickups: =============='
590   for xx in $listAt   listP=`( cd $dir1/atm ; ls pickup*.$nat2c* | sed 's/\..*//g' | uniq )`
591     for xx in $listP
592   do   do
593     dif_tiles $dir1/atm $dir3/atm $xx.$nat2c     dif_tiles $dir1/atm $dir3/atm $xx.$nat2c
594   done   done
595   exit   exit
596  fi  fi
597    
598    #- clean-up and reset:
599    if test $ksel = 7 ; then
600     echo "== clean-up and reset :"
601     cd rank_$rnkO
602     listS=`find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*'`
603     if test "x$listS" != x ; then
604       echo '--> remove pickup*.'$noc0c'* sym-links:'
605       if [ $prt -ge 2 ] ; then echo " in rank_$rnkO: rm" $listS ; fi
606       rm -f $listS
607     fi
608     cd ../rank_$rnkA
609     listS=`find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*'`
610     if test "x$listS" != x ; then
611       echo '--> remove pickup*.'$nat0c'* sym-links:'
612       if [ $prt -ge 2 ] ; then echo " in rank_$rnkA: rm" $listS ; fi
613       rm -f $listS
614     fi
615     cd ..
616    #- move back files from temp_tst dir:
617     echo '=> restore std data files'
618     rm -f rank_{$rnkO,$rnkA}/data
619     rm -f rank_{$rnkO,$rnkA}/data.tst
620     rm -f rank_{$rnkO,$rnkA}/data.pkg
621     mv $tmpDir.0/*Ocean.0001.clog $tmpDir.1/* rank_$rnkO
622     mv $tmpDir.0/*Atmos.0001.clog $tmpDir.2/* rank_$rnkA
623     mv $tmpDir.0/std_outp .
624     mv $tmpDir.0/* rank_0
625     rmdir $tmpDir.0 $tmpDir.1 $tmpDir.2
626    #echo '--> remove output dir and output files:'
627     echo ' rm -rf' $dir1 $dir2 $dir3
628     rm -rf $dir1 $dir2 $dir3
629     rm -f atm.txt ocn.txt
630     exit
631    fi

Legend:
Removed from v.1.5  
changed lines
  Added in v.1.15

  ViewVC Help
Powered by ViewVC 1.1.22