/[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.2 by jmc, Tue Oct 30 22:55:04 2007 UTC revision 1.5 by jmc, Mon Nov 26 23:46:21 2007 UTC
# Line 3  Line 3 
3  # $Header$  # $Header$
4  # $Name$  # $Name$
5    
6    rnp_loc()
7    {
8    #   rnp_loc arg_1 arg_2
9    #   rename files with prefix = 'arg_1' to files with prefix = 'arg_2'
10        echo 'rnp_loc:' $1 $2
11    #   rnp -s $1 $2
12        listY=`ls -1 $1*`
13        for yy in $listY
14        do
15          zz=`echo $yy | sed "s/^$1/$2/"`
16          mv $yy $zz
17        done
18    }
19    
20    dif_tiles()
21    {
22    #   dif_data_loc dd1 dd2 file
23    #   for all tiles  ==> do diff dd1/file dd2/file'
24        dd1=$1 ; dd2=$2 ; zz=$3
25        chkmeta='Y'
26        listY=`(cd $dd1 ; ls $zz.*.data | sed "s/$zz\.//" | sed "s/\.data//")`
27        echo '--> file='$zz', listY='$listY
28        for yy in $listY
29        do
30          echo diff $dd1/$zz.$yy.data $dd2
31          diff $dd1/$zz.$yy.data $dd2
32          out=$?
33          if test $out != 0
34          then echo 'Diff outp=' $out ' ==> stop'
35           exit ; fi
36          if test $chkmeta = 'Y'
37          then
38           echo diff $dd1/$zz.$yy.meta $dd2
39           diff $dd1/$zz.$yy.meta $dd2
40           out=$?
41           if test $out != 0
42           then echo 'Diff outp=' $out ' ==> stop'
43            exit ; fi
44          fi
45        done
46    }
47    
48  if test -f input_ocn/data.tst  if test -f input_ocn/data.tst
49  then  then
50   noc0=`sed -n 's/nIter0=//p' input_ocn/data.tst | sed 's/,//g' | sed 's/ //g'`   noc0=`sed -n 's/nIter0=//p' input_ocn/data.tst | sed 's/,//g' | sed 's/ //g'`
# Line 47  ksel=$1 Line 89  ksel=$1
89    
90  #make sure that local dir & ~jmc/bin is in the path:  #make sure that local dir & ~jmc/bin is in the path:
91  #export PATH=${PATH}:.  #export PATH=${PATH}:.
 RNPref=`ls ~jmc/bin/rnp`" -s"  
 Prt10c=`ls ~jmc/bin/add0upto10c`  
  difex=`ls ~jmc/bin/dif_data`  
 echo -n " use script: "  
 echo    "   Prt10c='$Prt10c' ,"  
 echo -n "   RNPref='$RNPref'"  
 echo " and difex='$difex'"  
92    
93  Nit=$NiOc  Nit=$NiOc
94  Dbl=`expr $Nit + $Nit`  Dbl=`expr $Nit + $Nit`
95  noc1=`expr $noc0 + $Nit`  noc1=`expr $noc0 + $Nit`
96  noc2=`expr $noc0 + $Dbl`  noc2=`expr $noc0 + $Dbl`
97  noc0c=`$Prt10c $noc0`  noc0c=`printf "%10.10i\n" $noc0`
98  noc1c=`$Prt10c $noc1`  noc1c=`printf "%10.10i\n" $noc1`
99  noc2c=`$Prt10c $noc2`  noc2c=`printf "%10.10i\n" $noc2`
100  echo 'Oce (rank_1):' $noc0  $noc1  $noc2  echo 'Oce (rank_1):' $noc0  $noc1  $noc2
101  echo $noc0c $noc1c $noc2c  echo $noc0c $noc1c $noc2c
102    
# Line 69  Nit=$NiAt Line 104  Nit=$NiAt
104  Dbl=`expr $Nit + $Nit`  Dbl=`expr $Nit + $Nit`
105  nat1=`expr $nat0 + $Nit`  nat1=`expr $nat0 + $Nit`
106  nat2=`expr $nat0 + $Dbl`  nat2=`expr $nat0 + $Dbl`
107  nat0c=`$Prt10c $nat0`  nat0c=`printf "%10.10i\n" $nat0`
108  nat1c=`$Prt10c $nat1`  nat1c=`printf "%10.10i\n" $nat1`
109  nat2c=`$Prt10c $nat2`  nat2c=`printf "%10.10i\n" $nat2`
110  echo 'Atm (rank_2):' $nat0  $nat1  $nat2  echo 'Atm (rank_2):' $nat0  $nat1  $nat2
111  echo $nat0c $nat1c $nat2c  echo $nat0c $nat1c $nat2c
112    
# Line 89  then Line 124  then
124   echo ' '   echo ' '
125   cd rank_1   cd rank_1
126   for xx in $listOc   for xx in $listOc
127   do   do rnp_loc $xx.ckptA $xx.$noc0c
   echo $RNPref $xx.ckptA $xx.$noc0c  
   $RNPref $xx.ckptA $xx.$noc0c  
128   done   done
129   cd ../rank_2   cd ../rank_2
130   for xx in $listAt   for xx in $listAt
131   do   do rnp_loc $xx.ckptA $xx.$nat0c
   echo $RNPref $xx.ckptA $xx.$nat0c  
   $RNPref $xx.ckptA $xx.$nat0c  
132   done   done
133   cd ..   cd ..
134  #-- prepare for running Dbl  #-- prepare for running Dbl
# Line 123  then Line 154  then
154   cd rank_1   cd rank_1
155   for xx in $listOc   for xx in $listOc
156   do   do
157    echo $RNPref $xx.ckptA $xx.$noc2c    rnp_loc $xx.ckptA $xx.$noc2c
   $RNPref $xx.ckptA $xx.$noc2c  
158   done   done
159   echo 'move pickups to' $dir1/ocn   echo 'move pickups to' $dir1/ocn
160   mv pickup*.$noc2c.* ../$dir1/ocn   mv pickup*.$noc2c.* ../$dir1/ocn
# Line 134  then Line 164  then
164   cd rank_2   cd rank_2
165   for xx in $listAt   for xx in $listAt
166   do   do
167    echo $RNPref $xx.ckptA $xx.$nat2c    rnp_loc $xx.ckptA $xx.$nat2c
   $RNPref $xx.ckptA $xx.$nat2c  
168   done   done
169   echo 'move pickups to' $dir1/atm   echo 'move pickups to' $dir1/atm
170   mv pickup*.$nat2c.* ../$dir1/atm   mv pickup*.$nat2c.* ../$dir1/atm
# Line 167  then Line 196  then
196   cd rank_1   cd rank_1
197   for xx in $listOc   for xx in $listOc
198   do   do
199    echo $RNPref $xx.ckptA $xx.$noc1c    rnp_loc $xx.ckptA $xx.$noc1c
   $RNPref $xx.ckptA $xx.$noc1c  
200   done   done
201   echo 'move pickups to' $dir2/ocn   echo 'move pickups to' $dir2/ocn
202   mv pickup*.$noc1c.* ../$dir2/ocn   mv pickup*.$noc1c.* ../$dir2/ocn
# Line 178  then Line 206  then
206   cd rank_2   cd rank_2
207   for xx in $listAt   for xx in $listAt
208   do   do
209    echo $RNPref $xx.ckptA $xx.$nat1c    rnp_loc $xx.ckptA $xx.$nat1c
   $RNPref $xx.ckptA $xx.$nat1c  
210   done   done
211   echo 'move pickups to' $dir2/atm   echo 'move pickups to' $dir2/atm
212   mv pickup*.$nat1c.* ../$dir2/atm   mv pickup*.$nat1c.* ../$dir2/atm
# Line 212  then Line 239  then
239   mkdir $dir3/ocn   mkdir $dir3/ocn
240   cd rank_1   cd rank_1
241   for xx in $listOc   for xx in $listOc
242   do   do rnp_loc $xx.ckptA $xx.$noc2c
   echo $RNPref $xx.ckptA $xx.$noc2c  
   $RNPref $xx.ckptA $xx.$noc2c  
243   done   done
244   echo 'move pickups to' $dir3/ocn   echo 'move pickups to' $dir3/ocn
245   mv pickup*.$noc2c.* ../$dir3/ocn   mv pickup*.$noc2c.* ../$dir3/ocn
246     rm -f pickup*.$noc1c.*
247   cd ..   cd ..
248    
249   mkdir $dir3/atm   mkdir $dir3/atm
250   cd rank_2   cd rank_2
251   for xx in $listAt   for xx in $listAt
252   do   do rnp_loc $xx.ckptA $xx.$nat2c
   echo $RNPref $xx.ckptA $xx.$nat2c  
   $RNPref $xx.ckptA $xx.$nat2c  
253   done   done
254   echo 'move pickups to' $dir3/atm   echo 'move pickups to' $dir3/atm
255   mv pickup*.$nat2c.* ../$dir3/atm   mv pickup*.$nat2c.* ../$dir3/atm
256     rm -f pickup*.$nat1c.*
257   cd ..   cd ..
258   echo 'restore std data files'   echo 'restore std data files'
259   cp -p input_cpl/data rank_0   cp -p input_cpl/data rank_0
# Line 291  then Line 316  then
316   echo ' ' ; echo '=== compare Ocn pickups: =============='   echo ' ' ; echo '=== compare Ocn pickups: =============='
317   for xx in $listOc   for xx in $listOc
318   do   do
319     $difex $dir1/ocn $dir3/ocn $xx.$noc2c     dif_tiles $dir1/ocn $dir3/ocn $xx.$noc2c
320   done   done
321   echo ' ' ; echo '=== compare Atm pickups: =============='   echo ' ' ; echo '=== compare Atm pickups: =============='
322   for xx in $listAt   for xx in $listAt
323   do   do
324     $difex $dir1/atm $dir3/atm $xx.$nat2c     dif_tiles $dir1/atm $dir3/atm $xx.$nat2c
325   done   done
326   exit   exit
327  fi  fi
   

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

  ViewVC Help
Powered by ViewVC 1.1.22