/[MITgcm]/MITgcm_contrib/jmc_script/tst_2+2_cpl
ViewVC logotype

Contents of /MITgcm_contrib/jmc_script/tst_2+2_cpl

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


Revision 1.14 - (show annotations) (download)
Thu Dec 24 21:29:23 2015 UTC (8 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: testing_tag01
Changes since 1.13: +18 -1 lines
adjust for global-file option

1 #! /usr/bin/env bash
2
3 # $Header: /u/gcmpack/MITgcm_contrib/jmc_script/tst_2+2_cpl,v 1.13 2015/10/25 18:02:14 jmc Exp $
4 # $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 if [ $prt -ge 1 ] ; then echo ' rnp_loc:' $1 $2 ; fi
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 if [ $prt -ge 2 ] ; then echo ' mv' $yy $zz ; fi
17 if test $action = 1 ; then mv $yy $zz ; fi
18 done
19 }
20
21 dif_tiles()
22 {
23 # dif_data_loc dd1 dd2 file
24 # for all tiles ==> do diff dd1/file dd2/file'
25 dd1=$1 ; dd2=$2 ; zz=$3
26 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//")`
45 echo '--> file='$zz', listY='$listY
46 for yy in $listY
47 do
48 echo diff $dd1/$zz.$yy.data $dd2
49 diff $dd1/$zz.$yy.data $dd2
50 out=$?
51 if test $out != 0
52 then echo 'Diff outp=' $out ' ==> stop'
53 exit ; fi
54 if test $chkmeta = 'Y' ; then
55 echo diff $dd1/$zz.$yy.meta $dd2
56 diff $dd1/$zz.$yy.meta $dd2
57 out=$?
58 if test $out != 0 ; then echo 'Diff outp=' $out ' ==> stop' ; exit ; fi
59 fi
60 done
61 }
62
63 move_outp()
64 {
65 # move_outp arg_1
66 # move mitgcm output files to directory 'arg_1'
67 if [ $prt -ge 1 ] ; then echo ' move_outp:' $1 ; fi
68 mv *.data *.meta $1
69 listF=`ls -1 *.txt STD???.0000 *.log 2> /dev/null`
70 if test "x$listF" != x ; then mv $listF $1 ; fi
71 #- 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
203 #----------------------------------------------------------------------
204 }
205
206 #===============================================================================
207 #-- 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 CMD="./${CMD} $MTH 3"
273 echo " doRun='$doRun' ; ksel='$ksel' ; kUp='$kUp' ; command='$CMD'"
274
275 #- set dir of Master-Proc for 3 components (coupler == rank_0):
276 if test $Npr = 0 ; then
277 echo " needs file: 'pr_group' to continue" ; exit
278 fi
279 rnkO=1 ; rnkA=`expr $rnkO + $NpOc`
280
281 #make sure that local dir & ~jmc/bin is in the path:
282 #export PATH=${PATH}:.
283
284 #----------------------------------------------------------------------
285
286 if test $ksel = 0 ; then
287 #-- from STDOUT of previous run, get Nb of iters for each component,
288 yy='(PID.TID 0000.0001) nTimeSteps = '
289 outpRun="rank_$rnkO/STDOUT.0000"
290 nTsO=`grep -A1 "$yy" $outpRun | tail -1 | sed 's/(PID.TID 0000.0001) *//'`
291 outpRun="rank_$rnkA/STDOUT.0000"
292 nTsA=`grep -A1 "$yy" $outpRun | tail -1 | sed 's/(PID.TID 0000.0001) *//'`
293 #- then calculate the ratio and set Nb of time-step for Atmos comp.
294 RatioTs=`expr $nTsA / $nTsO`
295 if [ $prt -ge 1 ] ; then echo " nTsA= $nTsA , nTsO= $nTsO , Ratio= $RatioTs"; fi
296 NiOc=$Nit ;
297 NiAt=`expr $NiOc \* $RatioTs`
298 #---- clean up current dir: rename pickup files and move output to temp_tst dir:
299 /bin/rm -rf $tmpDir.0 $tmpDir.1 $tmpDir.2
300 mkdir $tmpDir.0 $tmpDir.1 $tmpDir.2
301 mv std_outp rank_0/*.clog $tmpDir.0
302 mv rank_$rnkO/*.0001.clog $tmpDir.0
303 mv rank_$rnkA/*.0001.clog $tmpDir.0
304 echo ' '
305 echo "=> prepare OCN dir:"
306 prepare_first $NiOc rank_$rnkO $tmpDir.1
307 echo "=> prepare ATM dir:"
308 prepare_first $NiAt rank_$rnkA $tmpDir.2
309 #-- prepare for running Dbl
310 echo '=> prepare for running Dbl:'
311 Dbl=`expr $NiOc + $NiOc`
312 cp -p rank_$rnkO/data.tst rank_$rnkO/data
313 cp -p rank_$rnkA/data.tst rank_$rnkA/data
314 egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
315 fi
316
317 #----------------------------------------------------------------------
318
319 #- extract length of runs & starting iter number from file "data.tst"
320 # Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)
321 pFile="rank_$rnkO/data.tst"
322 if test -f $pFile ; then
323 noc0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
324 Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
325 NiOc=`expr $Dbl / 2`
326 Nit=$NiOc
327 else echo " file: '$pFile' not found"
328 NiOc=0 ; Nit=2
329 fi
330 pFile="rank_$rnkA/data.tst"
331 if test -f $pFile ; then
332 nat0=`sed -n 's/^ nIter0=//p' $pFile | sed 's/,//g' | sed 's/ //g'`
333 Dbl=`sed -n 's/^ nTimeSteps=//p' $pFile | sed 's/,//g'`
334 NiAt=`expr $Dbl / 2`
335 else echo " file: '$pFile' not found"
336 NiAt=0
337 fi
338 if test $NiAt = 0 -o $NiOc = 0 ; then
339 echo " needs 2 data files: 'rank_$rnkO/data.tst' & 'rank_$rnkA/data.tst'"
340 echo " (corresponding to 1rst run) to continue"
341 exit
342 fi
343
344 nit=$NiOc
345 Dbl=`expr $nit + $nit`
346 noc1=`expr $noc0 + $nit`
347 noc2=`expr $noc0 + $Dbl`
348 noc0c=`printf "%10.10i\n" $noc0`
349 noc1c=`printf "%10.10i\n" $noc1`
350 noc2c=`printf "%10.10i\n" $noc2`
351 echo " Oce (rank_$rnkO):" $noc0 $noc1 $noc2
352 echo " " $noc0c $noc1c $noc2c
353
354 nit=$NiAt
355 Dbl=`expr $nit + $nit`
356 nat1=`expr $nat0 + $nit`
357 nat2=`expr $nat0 + $Dbl`
358 nat0c=`printf "%10.10i\n" $nat0`
359 nat1c=`printf "%10.10i\n" $nat1`
360 nat2c=`printf "%10.10i\n" $nat2`
361 echo " Atm (rank_$rnkA):" $nat0 $nat1 $nat2
362 echo " " $nat0c $nat1c $nat2c
363
364 #----------------------------------------------------------------------
365
366 #-- run & post-process Dbl:
367 if [ $ksel -eq 1 -o $kUp -le 1 ] ; then
368 echo ' '
369 #-- run coupled test (Dbl):
370 if [ $doRun -eq 1 -o $kUp -lt 1 ] ; then
371 echo "=> run (Dbl): $CMD"
372 $CMD
373 fi
374 #-- after running Dbl:
375 #- check for pickup
376 listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
377 listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
378 if test "x$listOc" = x -o "x$listAt" = x
379 then echo "==> RUN 2 x $NiOc it STOP without writing $pickRef => exit" ; exit 1
380 else echo "==> END RUN 2 x $NiOc it" ; fi
381 #- move output
382 rm -rf $dir1
383 mkdir $dir1
384 echo '=> move STDOUT to' $dir1
385 mv std_outp rank_0/*.0000.clog rank_$rnkO/*.0001.clog rank_$rnkA/*.0001.clog $dir1
386 mv rank_$rnkO/STDOUT.0000 $dir1/ocnSTDOUT.0000
387 mv rank_$rnkA/STDOUT.0000 $dir1/atmSTDOUT.0000
388
389 mkdir $dir1/ocn
390 cd rank_$rnkO
391 listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
392 echo ' listP=' $listP
393 for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
394 echo 'move output to' $dir1/ocn
395 move_outp ../$dir1/ocn
396 cd ..
397
398 mkdir $dir1/atm
399 cd rank_$rnkA
400 listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
401 echo ' listP=' $listP
402 for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
403 echo 'move output to' $dir1/atm
404 move_outp ../$dir1/atm
405 cd ..
406
407 #-- prepare for running 1iA
408 echo '=> prepare for running 1iA:'
409 nit=$NiOc ; Dbl=`expr $nit + $nit`
410 pFile=rank_$rnkO/data
411 sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
412 nit=$NiAt ; Dbl=`expr $nit + $nit`
413 pFile=rank_$rnkA/data
414 sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$nit/g" $pFile.tst > $pFile
415 egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
416 fi
417
418 #-- run & post-process 1iA:
419 if [ $ksel -eq 2 -o $kUp -le 2 ] ; then
420 echo ' '
421 #-- run coupled test (1iA):
422 if [ $doRun -eq 1 -o $kUp -lt 2 ] ; then
423 echo "=> run (1iA): $CMD"
424 $CMD
425 fi
426 #-- after running 1iA
427 #- check for pickup
428 listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
429 listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
430 if test "x$listOc" = x -o "x$listAt" = x
431 then echo "==> RUN 1iA STOP without writing $pickRef => exit" ; exit 2
432 else echo "==> END RUN 1iA" ; fi
433 #- move output
434 rm -rf $dir2
435 mkdir $dir2
436 echo '=> move STDOUT to' $dir2
437 mv std_outp rank_0/*.clog rank_$rnkO/*.clog rank_$rnkA/*.clog $dir2
438 mv rank_$rnkO/STDOUT.0000 $dir2/ocnSTDOUT.0000
439 mv rank_$rnkA/STDOUT.0000 $dir2/atmSTDOUT.0000
440
441 mkdir $dir2/ocn
442 cd rank_$rnkO
443 listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
444 echo ' listP=' $listP
445 for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc1c ; done
446 # remove pickup*.$noc0c* sym-link:
447 find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*' -exec rm {} \;
448 echo 'move output to' $dir2/ocn
449 move_outp ../$dir2/ocn
450 cd ..
451
452 mkdir $dir2/atm
453 cd rank_$rnkA
454 listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
455 echo ' listP=' $listP
456 for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat1c ; done
457 # remove pickup*.$nat0c* sym-link:
458 find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*' -exec rm {} \;
459 echo 'move output to' $dir2/atm
460 move_outp ../$dir2/atm
461 cd ..
462
463 #-- prepare for running 1iB
464 echo '=> prepare for running 1iB:'
465 cd rank_$rnkO
466 ln -s ../$dir2/ocn/pickup*.$noc1c.* .
467 sed "s/^ nIter0=$noc0/ nIter0=$noc1/g" data > data.tmp_$$
468 mv -f data.tmp_$$ data
469 cd ../rank_$rnkA
470 ln -s ../$dir2/atm/pickup*.$nat1c.* .
471 sed "s/^ nIter0=$nat0/ nIter0=$nat1/g" data > data.tmp_$$
472 mv -f data.tmp_$$ data
473 cd ..
474 egrep 'nIter0|nTimeSteps' rank_$rnkO/data rank_$rnkA/data
475 fi
476
477 #-- run & post-process 1iB:
478 if [ $ksel -eq 3 -o $kUp -le 3 ] ; then
479 echo ' '
480 #-- run coupled test (1iB):
481 if [ $doRun -eq 1 -o $kUp -lt 3 ] ; then
482 echo "=> run (1iB): $CMD"
483 $CMD
484 fi
485 #-- after running 1iB
486 #- check for pickup
487 listOc=`ls -1 rank_$rnkO/$pickRef.ckptA*.data 2> /dev/null`
488 listAt=`ls -1 rank_$rnkA/$pickRef.ckptA*.data 2> /dev/null`
489 if test "x$listOc" = x -o "x$listAt" = x
490 then echo "==> RUN 1iB STOP without writing $pickRef => exit" ; exit 3
491 else echo "==> END RUN 1iB" ; fi
492 #- move output
493 rm -rf $dir3
494 mkdir $dir3
495 echo '=> move STDOUT to' $dir3
496 mv std_outp rank_0/*.clog rank_$rnkO/*.clog rank_$rnkA/*.clog $dir3
497 mv rank_$rnkO/STDOUT.0000 $dir3/ocnSTDOUT.0000
498 mv rank_$rnkA/STDOUT.0000 $dir3/atmSTDOUT.0000
499
500 mkdir $dir3/ocn
501 cd rank_$rnkO
502 listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
503 echo ' listP=' $listP
504 for xx in $listP ; do rnp_loc $xx.ckptA $xx.$noc2c ; done
505 # remove pickup*.$noc1c* sym-link:
506 find . ! -name . -prune -type l -name 'pickup*.'${noc1c}'*' -exec rm {} \;
507 echo 'move output to' $dir3/ocn
508 move_outp ../$dir3/ocn
509 cd ..
510
511 mkdir $dir3/atm
512 cd rank_$rnkA
513 listP=`ls pickup*.ckptA* | sed 's/\..*//g' | uniq`
514 echo ' listP=' $listP
515 for xx in $listP ; do rnp_loc $xx.ckptA $xx.$nat2c ; done
516 # remove pickup*.$nat1c* sym-link:
517 find . ! -name . -prune -type l -name 'pickup*.'${nat1c}'*' -exec rm {} \;
518 echo 'move output to' $dir3/atm
519 move_outp ../$dir3/atm
520 cd ..
521
522 fi
523
524 #-- compare output:
525 if test $ksel = 4 ; then
526 rm -f ocn.txt
527 echo ' run 2it: ---------- ' > tmp1.txt
528 grep "cg2d_init_res" $dir1/ocnSTDOUT.0000 \
529 | sed 's/.* cg2d_init_res =//' >> tmp1.txt
530 nit=`expr $NiOc + 2`
531 #echo "NiOc= $NiOc ; nit= $nit"
532 sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
533 rm -f tmp1.txt
534 echo ' run 1iA: ---------- ' > tmp1.txt
535 grep "cg2d_init_res" $dir2/ocnSTDOUT.0000 \
536 | sed 's/.* cg2d_init_res =//' >> tmp1.txt
537 echo ' run 1iB: ---------- ' >> tmp1.txt
538 grep "cg2d_init_res" $dir3/ocnSTDOUT.0000 \
539 | sed 's/.* cg2d_init_res =//' >> tmp1.txt
540 nl tmp1.txt > tmp1.ttt
541 nl tmp2.txt > tmp2.ttt
542 join tmp1.ttt tmp2.ttt | sed 's/[0-9]*//' > ocn.txt
543 rm -f tmp1.ttt tmp2.ttt tmp1.txt tmp2.txt
544 echo ' ' ; echo '-- compare cg2d_init_res (Ocean):'
545 cat ocn.txt
546 exit
547 fi
548
549 if test $ksel = 5 ; then
550 rm -f atm.txt
551 echo ' run 2it: ---------- ' > tmp1.txt
552 grep "cg2d_init_res" $dir1/atmSTDOUT.0000 \
553 | sed 's/.* cg2d_init_res =//' >> tmp1.txt
554 nit=`expr $NiAt + 2`
555 sed "$nit i \ ... 2it continue ..." tmp1.txt > tmp2.txt
556 rm -f tmp1.txt
557 echo ' run 1iA: ---------- ' > tmp1.txt
558 grep "cg2d_init_res" $dir2/atmSTDOUT.0000 \
559 | sed 's/.* cg2d_init_res =//' >> tmp1.txt
560 echo ' run 1iB: ---------- ' >> tmp1.txt
561 grep "cg2d_init_res" $dir3/atmSTDOUT.0000 \
562 | sed 's/.* cg2d_init_res =//' >> tmp1.txt
563 nl tmp1.txt > tmp1.ttt
564 nl tmp2.txt > tmp2.ttt
565 join tmp1.ttt tmp2.ttt | sed 's/[0-9]*//' > atm.txt
566 rm -f tmp1.ttt tmp2.ttt tmp1.txt tmp2.txt
567 echo ' ' ; echo '-- compare cg2d_init_res (Atmos):'
568 cat atm.txt
569 exit
570 fi
571
572 #- compare pickup file:
573 if test $ksel = 6 ; then
574 echo ' ' ; echo '=== compare Ocn pickups: =============='
575 listP=`( cd $dir1/ocn ; ls pickup*.$noc2c* | sed 's/\..*//g' | uniq )`
576 for xx in $listP
577 do
578 dif_tiles $dir1/ocn $dir3/ocn $xx.$noc2c
579 done
580 echo ' ' ; echo '=== compare Atm pickups: =============='
581 listP=`( cd $dir1/atm ; ls pickup*.$nat2c* | sed 's/\..*//g' | uniq )`
582 for xx in $listP
583 do
584 dif_tiles $dir1/atm $dir3/atm $xx.$nat2c
585 done
586 exit
587 fi
588
589 #- clean-up and reset:
590 if test $ksel = 7 ; then
591 echo "== clean-up and reset :"
592 cd rank_$rnkO
593 listS=`find . ! -name . -prune -type l -name 'pickup*.'${noc0c}'*'`
594 if test "x$listS" != x ; then
595 echo '--> remove pickup*.'$noc0c'* sym-links:'
596 if [ $prt -ge 2 ] ; then echo " in rank_$rnkO: rm" $listS ; fi
597 rm -f $listS
598 fi
599 cd ../rank_$rnkA
600 listS=`find . ! -name . -prune -type l -name 'pickup*.'${nat0c}'*'`
601 if test "x$listS" != x ; then
602 echo '--> remove pickup*.'$nat0c'* sym-links:'
603 if [ $prt -ge 2 ] ; then echo " in rank_$rnkA: rm" $listS ; fi
604 rm -f $listS
605 fi
606 cd ..
607 #- move back files from temp_tst dir:
608 echo '=> restore std data files'
609 rm -f rank_{$rnkO,$rnkA}/data
610 rm -f rank_{$rnkO,$rnkA}/data.tst
611 rm -f rank_{$rnkO,$rnkA}/data.pkg
612 mv $tmpDir.0/*Ocean.0001.clog $tmpDir.1/* rank_$rnkO
613 mv $tmpDir.0/*Atmos.0001.clog $tmpDir.2/* rank_$rnkA
614 mv $tmpDir.0/std_outp .
615 mv $tmpDir.0/* rank_0
616 rmdir $tmpDir.0 $tmpDir.1 $tmpDir.2
617 #echo '--> remove output dir and output files:'
618 echo ' rm -rf' $dir1 $dir2 $dir3
619 rm -rf $dir1 $dir2 $dir3
620 rm -f atm.txt ocn.txt
621 exit
622 fi

  ViewVC Help
Powered by ViewVC 1.1.22