1 | #!/bin/ksh |
---|
2 | |
---|
3 | set -vx |
---|
4 | |
---|
5 | ########################################################################################### |
---|
6 | # ajoute_bord.py |
---|
7 | ########################################################################################### |
---|
8 | # usage : |
---|
9 | # prompt> script_regrid restart_file_to_keep_the_grid.nc restart_file_to_get_the_values.nc |
---|
10 | |
---|
11 | # example in ksh shell : |
---|
12 | # > script_regrid CM1414A_20001231_sechiba_rest.nc CM5PIRC11_29901231_sechiba_rest.nc > shell_test_script_regrid_SECHIBA 2>&1 |
---|
13 | ########################################################################################### |
---|
14 | |
---|
15 | |
---|
16 | PATH_FILE_IN_TO_GRID=$1 |
---|
17 | FILE_IN_TO_GRID=$( basename $1 ) |
---|
18 | DIR_IN_TO_GRID=$( dirname $1 ) |
---|
19 | |
---|
20 | PATH_FILE_IN_VALUES=$2 |
---|
21 | FILE_IN_VALUES=$( basename $2 ) |
---|
22 | DIR_IN_VALUES=$( dirname $2 ) |
---|
23 | |
---|
24 | if [ X${DIR_IN_TO_GRID} != X. ] ; then |
---|
25 | cp ${1} . |
---|
26 | fi |
---|
27 | if [ X${DIR_IN_VALUES} != X. ] ; then |
---|
28 | cp ${2} . |
---|
29 | fi |
---|
30 | |
---|
31 | # Add true CF convention compatible with CDAT cdms library |
---|
32 | # -------------------------------------------------------- |
---|
33 | CF_restart.py ${FILE_IN_TO_GRID} |
---|
34 | CF_FILE_IN_TO_GRID=CF_${FILE_IN_TO_GRID} |
---|
35 | |
---|
36 | CF_restart.py ${FILE_IN_VALUES} |
---|
37 | CF_FILE_IN_VALUES=CF_${FILE_IN_VALUES} |
---|
38 | |
---|
39 | |
---|
40 | |
---|
41 | # Interpol the CF restart |
---|
42 | # ----------------------- |
---|
43 | # For SECHIBA restart files, we delete routing variables : |
---|
44 | # 1) with restget calls order in the code : |
---|
45 | # routingcounter,routingarea,routetogrid,routetobasin,basinid,topoindex,fastres,slowres,streamres,lakeres, |
---|
46 | # previous_outflow,lakeinflow,returnflow,riverflow,coastalflow,hydrographs, |
---|
47 | # runoff_route,drainage_route,evapot_route,precip_route,humrel_route,totnobio_route,vegtot_route |
---|
48 | # 2) with alphabetic order : |
---|
49 | # basinid, coastalflow, |
---|
50 | # drainage_route, evapot_route, fastres, humrel_route, hydrographs, lakeinflow, lakeres, |
---|
51 | # precip_route, previous_outflow, returnflow, riverflow, routetobasin, routetogrid, routingarea, |
---|
52 | # runoff_route, slowres, streamres, topoindex, totnobio_route, vegtot_route |
---|
53 | |
---|
54 | restart_interpol.py ${CF_FILE_IN_TO_GRID} ${CF_FILE_IN_VALUES} |
---|
55 | regrid_FILE_OUT=regrid_${CF_FILE_IN_VALUES} |
---|
56 | |
---|
57 | # Change CF convention to restcom axes |
---|
58 | # ------------------------------------ |
---|
59 | ncrename -d lat,y -d lon,x -d lev,z ${regrid_FILE_OUT} |
---|
60 | ncks -O -a -x -v lat,lon,lev,bounds_lat,bounds_lon ${regrid_FILE_OUT} ${regrid_FILE_OUT} |
---|
61 | |
---|
62 | # Name of the output file |
---|
63 | # ----------------------- |
---|
64 | FILE_OUT=grille_${regrid_FILE_OUT} |
---|
65 | |
---|
66 | # Correct mask of the output file |
---|
67 | # ------------------------------- |
---|
68 | ajoute_bords.py ${FILE_IN_TO_GRID} ${regrid_FILE_OUT} ${FILE_OUT} |
---|
69 | |
---|