/[MITgcm]/MITgcm_contrib/jmc_script/tst1+1
ViewVC logotype

Diff of /MITgcm_contrib/jmc_script/tst1+1

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

revision 1.10 by jmc, Mon Oct 29 23:18:18 2007 UTC revision 1.11 by jmc, Sun Nov 4 20:11:07 2007 UTC
# Line 3  Line 3 
3  # $Header$  # $Header$
4  # $Name$  # $Name$
5    
6    #- uses scripts in the path (1rst) or if not in the path, from dScript directory:
7    dScript=../../../../MITgcm_contrib/jmc_script
8    
9    #- input :
10    # 1) data.tst = parameter file "data" for 2.Nit iter run
11    # 2) here or in dir temp : pickup(_cd).0..0{num0}.001.001.data
12    
13    NbArg=$#
14    if test -r data.tst ; then
15  # Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)  # Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)
16  Dbl=`sed -n 's/nTimeSteps=//p' data.tst | sed 's/,//g'`   Dbl=`sed -n 's/nTimeSteps=//p' data.tst | sed 's/,//g'`
17  Nit=`expr $Dbl / 2`   Nit=`expr $Dbl / 2`
18  #Nit=2  else
19     echo '=> Need a readable "data.tst" file for the 2 x Nit iterations run'
20     Nit=2 ; NbArg=0;
21    fi
22    
23  case $# in  case $NbArg in
24    1) ;;    1) ;;
25     *) echo "Usage: `basename $0` flag"     *) echo "Usage: `basename $0` flag"
26        echo " Check restart: compare 1 run of 2 x $Nit it long"        echo " Check restart: compare 1 run of 2 x $Nit it long"
27        echo "     with 2 consecutive runs of $Nit it long each"        echo "     with 2 consecutive runs of $Nit it long each"
28        echo " where: flag = 1 -> do the 3 run"        echo " where: flag = 1 -> do the 3 runs"
29        echo "        flag = 2 -> compare std_outp"        echo "        flag = 2 -> compare std_outp"
30        echo "        flag = 3 -> compare pickup files"        echo "        flag = 3 -> compare pickup files"
31        exit ;;        exit ;;
32  esac  esac
33  flag=$1  flag=$1
34    
 #- input :  
 # 1) data.tst = parameter file "data" for 2.Nit iter run  
 # 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data  
   
35  num0=`sed -n 's/nIter0=//p' data.tst | sed 's/,//g' | sed 's/ //g'`  num0=`sed -n 's/nIter0=//p' data.tst | sed 's/,//g' | sed 's/ //g'`
36  #num0="72000"  #num0="72000"
37    
38  # list of pickup(s) that are needed for a restart :  # list of pickup(s) that are needed for a restart :
39  #listP='pickup pickup_cd pickup_nh pickup_ic pickup_land pickup_seaice'  #listP='pickup pickup_cd pickup_nh pickup_ic pickup_land pickup_seaice'
40    
41    #- executable:
42  mpi=0  mpi=0
43  if [ $mpi -ge 1 ]  if [ $mpi -ge 1 ]
44  then gcmExc=run_ifc_mpi  then gcmExc=run_ifc_mpi
# Line 41  dir1=res_2it Line 50  dir1=res_2it
50  dir2=res_1iA  dir2=res_1iA
51  dir3=res_1iB  dir3=res_1iB
52    
53  #make sure that local dir is in the path:  #---------------------------
54  export PATH=${PATH}:.  # find where the 3 little scripts are:
55    
56    #export PATH=${PATH}:$HOME/MITgcm_contrib/jmc_script
57    
58    if test -d $dScript ; then : ; else echo 'dScript='$dScript 'is not a directory' ; fi
59    
60    ff=rnp ;
61     RNPs=`which $ff 2> /dev/null`
62     if test "x$RNPs" = x -a -x $dScript/$ff ; then RNPs=$dScript/$ff ; fi
63     if test "x$RNPs" != x ; then RNPs=$RNPs" -s" ; fi
64    ff=add0upto10c ;
65     P10c=`which $ff 2> /dev/null`
66     if test "x$P10c" = x -a -x $dScript/$ff ; then P10c=$dScript/$ff ; fi
67    ff=dif_data ;
68     DIFd=`which $ff 2> /dev/null`
69     if test "x$DIFd" = x -a -x $dScript/$ff ; then DIFd=$dScript/$ff ; fi
70    echo " use script: P10c='$P10c'"
71    echo "             RNPs='$RNPs'"
72    echo "        and  DIFd='$DIFd'"
73    if test "x$RNPs" = x -o "x$P10c" = x -o "x$DIFd" = x
74    then echo 'missing some scripts' ; exit ; fi
75  #---------------------------  #---------------------------
76  Dbl=`expr $Nit + $Nit`  Dbl=`expr $Nit + $Nit`
77  num1=`expr $num0 + $Nit`  num1=`expr $num0 + $Nit`
78  num2=`expr $num0 + $Dbl`  num2=`expr $num0 + $Dbl`
79  num0c=`add0upto10c $num0`  num0c=`$P10c $num0`
80  num1c=`add0upto10c $num1`  num1c=`$P10c $num1`
81  num2c=`add0upto10c $num2`  num2c=`$P10c $num2`
82  echo $num0 $num1 $num2  echo $num0 $num1 $num2
83  echo $num0c $num1c $num2c  echo $num0c $num1c $num2c
84  echo 'gcmExc='$gcmExc  echo 'gcmExc='$gcmExc
85  # exit  # exit
86  #------------------------  #---------------------------
87    
88  if test $flag = 1 ; then  if test $flag = 1 ; then
89  #------------  #------------
# Line 71  rm -f STDERR.0000 Line 99  rm -f STDERR.0000
99    
100  cp -p data.tst data  cp -p data.tst data
101  numF=`ls temp/pickup*.${num0c}* 2> /dev/null | grep -c pickup`  numF=`ls temp/pickup*.${num0c}* 2> /dev/null | grep -c pickup`
102  if test $numF != 0 ; then ln -s temp/pickup*.${num0c}* . ; fi  if test $numF != 0 ; then ln -s temp/pickup*.${num0c}* . 2> /dev/null ; fi
103  echo '==> START RUN 2 x' $Nit 'it'  echo '==> START RUN 2 x' $Nit 'it'
104  if [ $mpi -ge 1 ]  if [ $mpi -ge 1 ]
105  then ./$gcmExc std_outp.2it > TTT.out.2it 2>&1  then ./$gcmExc std_outp.2it > TTT.out.2it 2>&1
# Line 87  listP=`ls pickup*.ckptA* | sed 's/\..*// Line 115  listP=`ls pickup*.ckptA* | sed 's/\..*//
115  echo ' listP=' $listP  echo ' listP=' $listP
116  for xx in $listP  for xx in $listP
117  do  do
118   rnp -s $xx.ckptA $xx.$num2c   $RNPs $xx.ckptA $xx.$num2c
119  done  done
120  mv *.0??.0??.* $dir1  mv *.0??.0??.* $dir1
121  txtF=`ls -1 *.txt 2> /dev/null`  txtF=`ls -1 *.txt 2> /dev/null`
# Line 109  if test $out != 0 ; then echo 'RUN 1iA S Line 137  if test $out != 0 ; then echo 'RUN 1iA S
137  echo '==>  END  RUN 1iA'  echo '==>  END  RUN 1iA'
138  for xx in $listP  for xx in $listP
139  do  do
140   rnp -s $xx.ckptA $xx.$num1c   $RNPs $xx.ckptA $xx.$num1c
141  done  done
142  mv *.0??.0??.* $dir2  mv *.0??.0??.* $dir2
143  txtF=`ls -1 *.txt 2> /dev/null`  txtF=`ls -1 *.txt 2> /dev/null`
# Line 133  if test $out != 0 ; then echo 'RUN 1iB S Line 161  if test $out != 0 ; then echo 'RUN 1iB S
161  echo '==>  END  RUN 1iB'  echo '==>  END  RUN 1iB'
162  for xx in $listP  for xx in $listP
163  do  do
164   rnp -s $xx.ckptA $xx.$num2c   $RNPs $xx.ckptA $xx.$num2c
165  done  done
166  mv *.0??.0??.* $dir3  mv *.0??.0??.* $dir3
167  txtF=`ls -1 *.txt 2> /dev/null`  txtF=`ls -1 *.txt 2> /dev/null`
# Line 161  fi Line 189  fi
189    
190  if test $flag = 3 ; then  if test $flag = 3 ; then
191  #------------  #------------
 difex=`ls ~jmc/bin/dif_data`  
192  listP=`( cd $dir1 ; ls pickup*.$num2c* | sed 's/\..*//g' | uniq )`  listP=`( cd $dir1 ; ls pickup*.$num2c* | sed 's/\..*//g' | uniq )`
193  for xx in $listP  for xx in $listP
194  do  do
195   $difex $dir1 $dir3 $xx.$num2c   $DIFd $dir1 $dir3 $xx.$num2c
196  done  done
197  exit  exit
198  #------------  #------------
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22