Subversion Repositories lagranto.um

Rev

Rev 3 | Rev 13 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 3 Rev 5
Line 254... Line 254...
254 
      @ timediff_mm = 60 + $min
 254 
      @ timediff_mm = 60 + $min
255 
    else
 255 
    else
256 
      @ timediff_mm = 0 - $min
 256 
      @ timediff_mm = 0 - $min
257 
    endif
 257 
    endif
258 
  endif
 258 
  endif
- 
 
 259 
else
- 
 
 260 
  set timediff_mm = 0
259 
endif
 261 
endif
260 
 
 262 
 
261 
# Write status information
 263 
# Write status information
262 
echo "Time range      : ${startdate} -> ${enddate}"
 264 
echo "Time range      : ${startdate} -> ${enddate}"
263 
if ( ${timediff_mm} != "" ) then
 265 
if ( ${timediff_mm} != "" ) then
Line 290... Line 292...
290 
endif
 292 
endif
291 
 
 293 
 
292 
# Determine timeinc (the time difference in hours between two data file)
 294 
# Determine timeinc (the time difference in hours between two data file)
293 
# if not already defined with option -i
 295 
# if not already defined with option -i
294 
if ( ${timeinc} == "nil" ) then
 296 
if ( ${timeinc} == "nil" ) then
295 
  set date1=`echo $file | cut -c 2-12`
 297 
  set date1=`echo $file | cut -c 2-`
296 
  set n=`ls ${charp}[0-9_]*[0-9] | grep -n $date1 | awk -F: '{print $1}'`
 298 
  set n=`ls ${charp}[0-9_]*[0-9] | grep -n $date1 | awk -F: '{print $1}'`
297 
  @ n ++
 299 
  @ n ++
298 
  set date2=`ls ${charp}[0-9_]*[0-9] | head -$n | tail -1 | cut -c 2-12`
 300 
  set date2=`\ls ${charp}[0-9_]*[0-9] | head -$n | tail -1 | cut -c 2-`
299 
  set timeinc=`${LAGRANTO}/goodies/gettidiff $date2 $date1`
 301 
  set timeinc=`${LAGRANTO}/goodies/gettidiff $date2 $date1`
300 
endif
 302 
endif
301 
if ( $timeinc == 0 ) then
 303 
if ( $timeinc == 0 ) then
302 
    echo " ERROR: cannot set the time increment between input files ***"
 304 
    echo " ERROR: cannot set the time increment between input files ***"
303 
    exit 1
 305 
    exit 1
304 
endif
 306 
endif
305 
 
 307 
 
- 
 
 308 
# Decide whether the time increment is smaller than 1 h; if so set the hires
- 
 
 309 
# flag and convert the timeincement into fractional hour
306 
# Search the first file to use
 310 
set flag = `echo ${timeinc} | grep -n '\.'`
- 
 
 311 
set hires = 0
- 
 
 312 
if ( "${flag}" != "" ) then
- 
 
 313 
  echo "Input files availability < 1 hour ($timeinc): switch to hires"
- 
 
 314 
  echo "mode. All input files must be of form: P{yyyymmdd_hhmm}"
- 
 
 315 
  set hires = 1
- 
 
 316 
  set timeinc_min = `echo "(${timeinc} * 100 + 0.5)/1" | bc`
- 
 
 317 
  set timeinc     = `echo "scale=2; ${timeinc} * 100 / 60" | bc`
- 
 
 318 
  goto hires
- 
 
 319 
endif
- 
 
 320 
 
- 
 
 321 
## Search the first file to use: We step through all P files and see whether they are
- 
 
 322 
# good P files. Let's first do the test for the first data file found. If it's ok, we
- 
 
 323 
# take it; if not, we step through all P files and find the good one
307 
set flag=0
 324 
set flag=0
308 
set td=
 325 
set td=
- 
 
 326 
 
- 
 
 327 
set date = `echo $file | cut -c 2-12`
- 
 
 328 
set td1  = `${LAGRANTO}/goodies/gettidiff ${startdate_ymdh} ${date}`
- 
 
 329 
set td2  = `${LAGRANTO}/goodies/gettidiff ${enddate_ymdh}   ${date}`
- 
 
 330 
 
- 
 
 331 
if (( $td1 < $timeinc || $td2 < $timeinc ) && ( $td1 >= 0 || $td2 >= 0 )) then
- 
 
 332 
   set datfiles=$date
- 
 
 333 
   if ( $td1 < $timeinc    ) set td=$td1
- 
 
 334 
   if ( $td2 < $timeinc    ) set td=$td2
- 
 
 335 
   if ( ( $startdate_min > 0 ) || ( $enddate_min > 0 ) ) @ td ++
- 
 
 336 
   goto label2
- 
 
 337 
endif
- 
 
 338 
 
309 
foreach i ( ${charp}????????_?? )
 339 
foreach i ( ${charp}????????_?? )
310 
 
 340 
 
311 
  set date = `echo $i | cut -c 2-12`
 341 
  set date = `echo $i | cut -c 2-12`
312 
  set td1  = `${LAGRANTO}/goodies/gettidiff ${startdate_ymdh} ${date}`
 342 
  set td1  = `${LAGRANTO}/goodies/gettidiff ${startdate_ymdh} ${date}`
313 
  set td2  = `${LAGRANTO}/goodies/gettidiff ${enddate_ymdh}   ${date}`
 343 
  set td2  = `${LAGRANTO}/goodies/gettidiff ${enddate_ymdh}   ${date}`
Line 320... Line 350...
320 
      goto label2
 350 
      goto label2
321 
  endif
 351 
  endif
322 
 
 352 
 
323 
end
 353 
end
324 
 
 354 
 
325 
# if no P/T-files are available for the specified time period, then $td is
 - 
 
326 
# still undefined
 - 
 
327 
if ( $td == "" ) then
 355 
if ( $td == "" ) then
328 
  echo " ERROR: no data files available for the specified time period"
 356 
  echo " ERROR: no data files available for the specified time period"
329 
  exit 1
 357 
  exit 1
330 
endif
 358 
endif
331 
 
 359 
 
- 
 
 360 
# If the input files come with time increment < 1 hour, handle it in this
- 
 
 361 
# special section. The aim is not to 'interfere' with the well established
- 
 
 362 
# code handling the standard case with time increment >= 1 hour. In this mode
- 
 
 363 
# we require that the start and end date are available as P files: no start is
- 
 
 364 
# allowed between the model output files
- 
 
 365 
hires:
- 
 
 366 
 
- 
 
 367 
set flag=0
- 
 
 368 
set td=0
- 
 
 369 
 
- 
 
 370 
if ( $direction == "f" ) then
- 
 
 371 
  set file=${charp}${startdate}
- 
 
 372 
else
- 
 
 373 
  set file=${charp}${enddate}
- 
 
 374 
endif
- 
 
 375 
 
- 
 
 376 
if ( -f $file ) then
- 
 
 377 
  set date     = `echo ${file} | cut -c 2-14`
- 
 
 378 
  set datfiles = $date
- 
 
 379 
  set td       = 0
- 
 
 380 
  goto label2
- 
 
 381 
else
- 
 
 382 
  echo
- 
 
 383 
  echo " ERROR: no data files available for the specified time period"
- 
 
 384 
  exit 1
- 
 
 385 
endif
- 
 
 386 
 
332 
# Everything is fine so far: proceed
 387 
# Everything is fine so far: proceed
333 
label2:
 388 
label2:
334 
 
 389 
 
335 
# Calculate the number of further files
 390 
# Calculate the number of further files; handle hires (timeinc<1 h) and
- 
 
 391 
# standard mode (timeinc >= 1 h) separately
- 
 
 392 
if ( ${hires} == 0 ) then
336 
@ num = ( $timediff_hh + $td ) / $timeinc + 1
 393 
 @ num = ( $timediff_hh + $td ) / $timeinc + 1
337 
@ dum1 = ( $num - 1 ) * $timeinc
 394 
 @ dum1 = ( $num - 1 ) * $timeinc
338 
@ dum2 = $timediff_hh + $td
 395 
 @ dum2 = $timediff_hh + $td
339 
if ( $dum1 != $dum2 ) @ num ++
 396 
 if ( $dum1 != $dum2 ) @ num ++
- 
 
 397 
else
- 
 
 398 
 set num = `echo "( ${timediff_hh}*60 + ${timediff_mm} ) / ( ${timeinc_min} ) + 1" | bc`
- 
 
 399 
endif
340 
 
 400 
 
341 
# Get a list of all needed files
 401 
# Get a list of all needed files
342 
set numfiles=$num
 402 
set numfiles=$num
343 
set sfiles=1
 403 
set sfiles=1
- 
 
 404 
echo $datfiles
344 
while ( $num > 1 )
 405 
while ( $num > 1 )
- 
 
 406 
  if ( $hires == 0 ) then
345 
  set date=`${LAGRANTO}goodies/newtime $date $timeinc`
 407 
     set date=`${LAGRANTO}/goodies/newtime $date $timeinc`
- 
 
 408 
  else
- 
 
 409 
     set date=`${LAGRANTO}/goodies/newtime $date $timeinc_min`
- 
 
 410 
  endif
346 
  if ( ! -f ${charp}${date} ) then
 411 
  if ( ! -f ${charp}${date} ) then
347 
    echo " ERROR: file with primary data is missing for $date"
 412 
    echo " ERROR: file with primary data is missing for $date"
348 
    exit 1
 413 
    exit 1
349 
  else if ( ! -f ${chars}${date} ) then
 414 
  else if ( ! -f ${chars}${date} ) then
350 
    set sfiles=0
 415 
    set sfiles=0
Line 366... Line 431...
366 
 
 431 
 
367 
 
 432 
 
368 
# Write some status information
 433 
# Write some status information
369 
echo "Primary file prefix               : ${charp}"
 434 
echo "Primary file prefix               : ${charp}"
370 
echo "Secondary file prefix             : ${chars}"
 435 
echo "Secondary file prefix             : ${chars}"
- 
 
 436 
if ( ${hires} == 0 ) then
371 
echo "Time increment for input files    : ${timeinc}"
 437 
  echo "Time increment for input files    : ${timeinc} h"
- 
 
 438 
else
- 
 
 439 
  echo "Time increment for input files    : ${timeinc_min} min (0$timeinc h) "
- 
 
 440 
endif
372 
echo "# input files                     : ${numfiles}"
 441 
echo "# input files                     : ${numfiles}"
373 
echo "First input file                  : $datfiles[1] "
442 
echo "First input file                  : $datfiles[1] "
374 
echo "Last input file                   : $datfiles[$numfiles] "
443 
echo "Last input file                   : $datfiles[$numfiles] "
375 
echo "${charp} files availability              : 1"
444 
echo "${charp} files availability              : 1"
376 
echo "${chars} files availability              : ${sfiles}"
445 
echo "${chars} files availability              : ${sfiles}"
Line 380... Line 449...
380 
else
 449 
else
381 
echo "Start time relative to last file  : $datfiles[$numfiles] - ${tstart} "
 450 
echo "Start time relative to last file  : $datfiles[$numfiles] - ${tstart} "
382 
echo "End time relative to first file   : $datfiles[1] + ${tend} "
 451 
echo "End time relative to first file   : $datfiles[1] + ${tend} "
383 
endif
 452 
endif
384 
 
 453 
 
- 
 
 454 
echo
- 
 
 455 
echo $datfiles
- 
 
 456 
 
385 
# ---------------------------------------------------------------------
 457 
# ---------------------------------------------------------------------
386 
# Check availability of input data
458 
# Check availability of input data
387 
# ---------------------------------------------------------------------
 459 
# ---------------------------------------------------------------------
388 
 
 460 
 
389 
echo
 461 
echo
Line 395... Line 467...
395 
 
 467 
 
396 
#   Preset values for <compfl> and <tprefix>
 468 
#   Preset values for <compfl> and <tprefix>
397 
    set tcompfl=1
 469 
    set tcompfl=1
398 
    set tprefix='P'
 470 
    set tprefix='P'
399
471 
 
- 
 
 472 
#   Split the variable name into its vectorial components, if necessary
- 
 
 473 
    set tvar0 = `echo ${tvar} | awk -F. '{print $1}' `
- 
 
 474 
    set tvar1 = `echo ${tvar} | awk -F. '{print $2}' `
- 
 
 475
400 
#   Check availability on P file
 476 
#   Check availability on P file - only for first vector component
401 
    foreach var ( `${LAGRANTO}/goodies/getvars ${charp}$datfiles[1]` )
 477 
    foreach var ( `${LAGRANTO}/goodies/getvars ${charp}$datfiles[1]` )
402 
       if ( "${var}" == "${tvar}" ) then
 478 
       if ( "${var}" == "${tvar0}" ) then
403 
          set tcompfl=0
 479 
          set tcompfl=0
404 
	  set tprefix="P"
 480 
              set tprefix="P"
405 
       endif
 481 
       endif
406 
    end
 482 
    end
407 
 
 483 
 
408 
#   Check availability on S file
484 
#   Check availability on S file
409 
    if ( ${sfiles} == 1 ) then
 485 
    if ( ${sfiles} == 1 ) then
410 
       foreach var ( `${LAGRANTO}/goodies/getvars ${chars}$datfiles[1]` )
 486 
       foreach var ( `${LAGRANTO}/goodies/getvars ${chars}$datfiles[1]` )
411 
         if ( "${var}" == "${tvar}" ) then
 487 
         if ( "${var}" == "${tvar0}" ) then
412 
            set tcompfl=0
 488 
            set tcompfl=0
413 
	    set tprefix="S"
 489 
                set tprefix="S"
414 
         endif
 490 
         endif
415 
       end
 491 
       end
416 
    endif
 492 
    endif
Line 420... Line 496...
420 
    echo "Temporary tracervar file <${tvfile}> created"
 496 
    echo "Temporary tracervar file <${tvfile}> created"
421 
    echo
 497 
    echo
422 
 
 498 
 
423 
endif
 499 
endif
424 
 
 500 
 
425 
 
 - 
 
426 
# Check if tracevars-file exists
 501 
# Check if tracevars-file exists
427 
if ( ! -f $tvfile ) then
 502 
if ( ! -f $tvfile ) then
428 
  echo  " ERROR:  file $tvfile was not found ***"
 503 
  echo  " ERROR:  file $tvfile was not found ***"
429 
  exit 1
 504 
  exit 1
430 
endif
 505 
endif
Line 439... Line 514...
439 
  echo " ERROR: tracevars-files must not contain empty lines ***"
 514 
  echo " ERROR: tracevars-files must not contain empty lines ***"
440 
  exit 1
 515 
  exit 1
441 
endif
 516 
endif
442 
set calf=`cat $tvfile | awk '{print $3}'`
 517 
set calf=`cat $tvfile | awk '{print $3}'`
443 
set tfil=`cat $tvfile | awk '{print $4}'`
 518 
set tfil=`cat $tvfile | awk '{print $4}'`
- 
 
 519 
 
- 
 
 520 
# Write some status information
- 
 
 521 
cat ${tvfile}
- 
 
 522 
echo
- 
 
 523 
echo "# Number of tracing variables : ${nlines}"
- 
 
 524 
echo "Fields are read from following files : ${tfil}"
- 
 
 525 
 
- 
 
 526 
# Loop over all variables - check availability
444 
foreach v ( $vars )
 527 
foreach v ( $vars )
445 
  if ( $calf[1] == 0 ) then
 528 
  if ( $calf[1] == 0 ) then
446 
    set v0 = `echo $v | awk 'BEGIN {FS = ":"}; {print $1}'`
 529 
    set v0 = `echo $v | awk 'BEGIN {FS = ":"}; {print $1}'`
447 
    set flag=`${LAGRANTO}goodies/getvars $tfil[1]$datfiles[1] | grep " $v0 " | wc -l`
 530 
    set flag=`${LAGRANTO}/goodies/getvars $tfil[1]$datfiles[1] | grep " $v0 " | wc -l`
448 
    set iscomment=`echo $v0 | cut -c 1`
531 
    set iscomment=`echo $v0 | cut -c 1`
449 
    if ( "${iscomment}" != "#" ) then
532 
    if ( "${iscomment}" != "#" ) then
450 
	 if ( $flag == 0 ) then
 533 
         if ( $flag == 0 ) then
451 
           echo " ERROR: variable $v listed in $tvfile is not on the $tfil[1]-files ***"
 534 
           echo " ERROR: variable $v listed in $tvfile is not on the $tfil[1]-files ***"
452 
           exit 1
 535 
           exit 1
Line 456... Line 539...
456 
  shift calf
 539 
  shift calf
457 
  shift tfil
 540 
  shift tfil
458 
end
 541 
end
459 
set ntrace=${nlines}
542 
set ntrace=${nlines}
460 
 
 543 
 
461 
# Write some status information
 - 
 
462 
cat ${tvfile}
 - 
 
463 
echo
 - 
 
464 
echo "# Number of tracing variables : ${ntrace}"
 - 
 
465 
 
 - 
 
466 
# ---------------------------------------------------------------------
 544 
# ---------------------------------------------------------------------
467 
# Prepare input file for trace and run it
 545 
# Prepare input file for trace and run it
468 
# ---------------------------------------------------------------------
 546 
# ---------------------------------------------------------------------
469 
 
 547 
 
470 
# Set times relative to the reference date
 548 
# Set times relative to the reference date
Line 481... Line 559...
481 
    end
 559 
    end
482 
  endif
 560 
  endif
483 
endif
 561 
endif
484 
 
 562 
 
485 
# ---------------------------------------------------------------------
 563 
# ---------------------------------------------------------------------
486 
# Prepare input file for caltra and run it
 564 
# Prepare input file for trace and run it
487 
# ---------------------------------------------------------------------
 565 
# ---------------------------------------------------------------------
488 
 
 566 
 
489 
# Write parameter file
 567 
# Write parameter file
490 
\rm -f ${parfile}
 568 
\rm -f ${parfile}
491 
touch ${parfile}
 569 
touch ${parfile}