/[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.4 by jmc, Thu Mar 31 14:48:15 2005 UTC revision 1.6 by jmc, Sun Mar 26 23:21:50 2006 UTC
# Line 1  Line 1 
1  #!/bin/sh  #!/bin/sh
2    
3    # $Header$
4    # $Name$
5    
6  # 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)
7  Dbl=`sed -n 's/nTimeSteps=//p' data.tst | sed 's/,//g'`  Dbl=`sed -n 's/nTimeSteps=//p' data.tst | sed 's/,//g'`
8  Nit=`expr $Dbl / 2`  Nit=`expr $Dbl / 2`
# Line 17  case $# in Line 20  case $# in
20  esac  esac
21  flag=$1  flag=$1
22    
23    #- input :
 #- input :  
24  # 1) data.tst = parameter file "data" for 2.Nit iter run  # 1) data.tst = parameter file "data" for 2.Nit iter run
25  # 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data  # 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data
26    
27  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'`
28  #num0="72000"  #num0="72000"
29    
 gcmExc=mitgcmuv  
30  # list of pickup(s) that are needed for a restart :  # list of pickup(s) that are needed for a restart :
31  listP='pickup'  listP='pickup'
32  #listP='pickup pickup_cd'  #listP='pickup pickup_cd pickup_nh pickup_ic pickup_land'
33    
34    mpi=0
35    if [ $mpi -ge 1 ]
36    then gcmExc=run_ifc_mpi
37    else gcmExc=mitgcmuv
38    fi
39    
40  #- dir where to put the results :  #- dir where to put the results :
41  dir1=res_2it  dir1=res_2it
42  dir2=res_1iA  dir2=res_1iA
43  dir3=res_1iB  dir3=res_1iB
44    
45    #make sure that local dir is in the path:
46    export PATH=${PATH}:.
47    
48  #---------------------------  #---------------------------
49  Dbl=`expr $Nit + $Nit`  Dbl=`expr $Nit + $Nit`
50  num1=`expr $num0 + $Nit`  num1=`expr $num0 + $Nit`
# Line 42  num0c=`add0upto10c $num0` Line 53  num0c=`add0upto10c $num0`
53  num1c=`add0upto10c $num1`  num1c=`add0upto10c $num1`
54  num2c=`add0upto10c $num2`  num2c=`add0upto10c $num2`
55  echo $num0 $num1 $num2  echo $num0 $num1 $num2
56  echo $num0c $num1c $num2c  echo $num0c $num1c $num2c
57  echo 'gcmExc='$gcmExc  echo 'gcmExc='$gcmExc
58  # exit  # exit
59  #------------------------  #------------------------
# Line 55  mkdir $dir1 $dir2 $dir3 Line 66  mkdir $dir1 $dir2 $dir3
66  cp -p data.tst data  cp -p data.tst data
67  ln -s temp/pickup*.${num0c}* .  ln -s temp/pickup*.${num0c}* .
68  echo '==> START RUN 2 x' $Nit 'it'  echo '==> START RUN 2 x' $Nit 'it'
69  ./$gcmExc > std_outp.2it 2>&1  if [ $mpi -ge 1 ]
70  if test -f pickup.ckptA.001.001.data  then ./$gcmExc std_outp.2it > TTT.out.2it 2>&1
71    else ./$gcmExc > std_outp.2it 2>&1
72    fi
73    if test -f pickup.ckptA.001.001.data
74  then out=`egrep -c 'STOP ABNORMAL END' std_outp.2it`  then out=`egrep -c 'STOP ABNORMAL END' std_outp.2it`
75  else out=9 ; fi  else out=9 ; fi
76  if test $out != 0 ; then echo 'RUN 2it STOP before END => exit' ; exit ; fi  if test $out != 0 ; then echo 'RUN 2it STOP before END => exit' ; exit ; fi
# Line 66  do Line 80  do
80   rnp -s $xx.ckptA $xx.$num2c   rnp -s $xx.ckptA $xx.$num2c
81  done  done
82  mv *.0??.0??.* $dir1  mv *.0??.0??.* $dir1
83    mv R?.data R?.meta DR?.data DR?.meta $dir1
84  mv $dir1/pickup*.${num0c}* .  mv $dir1/pickup*.${num0c}* .
85  #--  #--
86  sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" data.tst > data  sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" data.tst > data
87  echo '==> START RUN 1iA'  echo '==> START RUN 1iA'
88  ./$gcmExc > std_outp.1iA 2>&1  if [ $mpi -ge 1 ]
89  if test -f pickup.ckptA.001.001.data  then ./$gcmExc std_outp.1iA > TTT.out.1iA 2>&1
90    else ./$gcmExc > std_outp.1iA 2>&1
91    fi
92    if test -f pickup.ckptA.001.001.data
93  then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iA`  then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iA`
94  else out=9 ; fi  else out=9 ; fi
95  if test $out != 0 ; then echo 'RUN 1iA STOP before END => exit' ; exit ; fi  if test $out != 0 ; then echo 'RUN 1iA STOP before END => exit' ; exit ; fi
# Line 81  do Line 99  do
99   rnp -s $xx.ckptA $xx.$num1c   rnp -s $xx.ckptA $xx.$num1c
100  done  done
101  mv *.0??.0??.* $dir2  mv *.0??.0??.* $dir2
102    mv R?.data R?.meta DR?.data DR?.meta $dir2
103  rm $dir2/pickup*.${num0c}*  rm $dir2/pickup*.${num0c}*
104  #--  #--
105  ln -s $dir2/pickup*.$num1c.* .  ln -s $dir2/pickup*.$num1c.* .
106  mv data data.tmp  mv data data.tmp
107  sed "s/^ nIter0=$num0/ nIter0=$num1/g" data.tmp > data  sed "s/^ nIter0=$num0/ nIter0=$num1/g" data.tmp > data
108  echo '==> START RUN 1iB'  echo '==> START RUN 1iB'
109  ./$gcmExc > std_outp.1iB 2>&1  if [ $mpi -ge 1 ]
110  if test -f pickup.ckptA.001.001.data  then ./$gcmExc std_outp.1iB > TTT.out.1iB 2>&1
111    else ./$gcmExc > std_outp.1iB 2>&1
112    fi
113    if test -f pickup.ckptA.001.001.data
114  then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iB`  then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iB`
115  else out=9 ; fi  else out=9 ; fi
116  if test $out != 0 ; then echo 'RUN 1iB STOP before END => exit' ; exit ; fi  if test $out != 0 ; then echo 'RUN 1iB STOP before END => exit' ; exit ; fi
# Line 98  do Line 120  do
120   rnp -s $xx.ckptA $xx.$num2c   rnp -s $xx.ckptA $xx.$num2c
121  done  done
122  mv *.0??.0??.* $dir3  mv *.0??.0??.* $dir3
123    mv R?.data R?.meta DR?.data DR?.meta $dir3
124  rm $dir3/pickup*.$num1c.*  rm $dir3/pickup*.$num1c.*
125  exit  exit
126  #------------  #------------
# Line 114  echo ' run 1iB:' Line 137  echo ' run 1iB:'
137        | sed 's/.* cg2d_init_res =//'        | sed 's/.* cg2d_init_res =//'
138  echo ' run 2it:'  echo ' run 2it:'
139   grep "cg2d_init_res" std_outp.2it \   grep "cg2d_init_res" std_outp.2it \
140        | sed 's/.* cg2d_init_res =//'        | sed 's/.* cg2d_init_res =//'
141  exit  exit
142  #------------  #------------
143  fi  fi
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22