| 62 | |
| 63 | == Zooming == |
| 64 | |
| 65 | Sometimes you are only interested in a certain region of the world. This is certainly the case when you are using a zoomed grid to run the simulation. Is there a way to get ATLAS to only print out the results over a single region? |
| 66 | |
| 67 | The answer is yes. There is a script, 2variables1diffandzonal_map_zoomed_new.jnl, which calculates both the maps and the zonal differences for two variables over just a particular region. In the .cfg files, the script 2variables1diffandzonal_map.jnl is currently used. The first step is to change 2variables1diffandzonal_map.jnl to 2variables1diffandzonal_map_zoomed_new.jnl in every line. So, for example, if before there was |
| 68 | |
| 69 | {{{ |
| 70 | echo $input | ferretloop - $paper -r $your_reference_file -i "Precip (mm/day)" -s $frequency -o precip 2variables1diffandzonal_map.jnl '"Precip (mm/day)" inputfile (precip[l=sloop]*86400) referencefile (precip[l=sloop]*86400) /pal=precip/lev=(-INF)(0.5)(1,6,1)(8,14,2)(20)(INF) /pal=pal_01/lev=(-INF)(-5)(-2,2,1)(-0.5)(-0.2)(0.2)(0.5)(5)(INF) /hlim=0:10 /hlim=-3:3' |
| 71 | }}} |
| 72 | |
| 73 | afterwards there should be |
| 74 | |
| 75 | {{{ |
| 76 | echo $input | ferretloop - $paper -r $your_reference_file -i "Precip (mm/day)" -s $frequency -o precip 2variables1diffandzonal_map_zoomed_new.jnl '"Precip (mm/day)" inputfile (precip[l=sloop]*86400) referencefile (precip[l=sloop]*86400) /pal=precip/lev=(-INF)(0.5)(1,6,1)(8,14,2)(20)(INF) /pal=pal_01/lev=(-INF)(-5)(-2,2,1)(-0.5)(-0.2)(0.2)(0.5)(5)(INF) /hlim=0:10 /hlim=-3:3' |
| 77 | }}} |
| 78 | |
| 79 | The command |
| 80 | |
| 81 | {{{ |
| 82 | sed -i.bak -e "s%2variables1diffandzonal_map.jnl%2variables1diffandzonal_map_zoomed_new.jnl%" atlas_p_LMDZ3.3_diag2_a.cfg |
| 83 | }}} |
| 84 | |
| 85 | will do that in one step, copying the file first to something.bak and the overwriting the existing file. |
| 86 | |
| 87 | The new script also takes two additional arguments. The first is the ferret command used to create the map (either "shade" or "fill"), and the second is the region of the zoom. So if we want to look at Europe, we can use a line like the following: |
| 88 | |
| 89 | {{{ |
| 90 | echo $input | ferretloop - $paper -r $your_reference_file -i "Precip (mm/day)" -s $frequency -o precip 2variables1diffandzonal_map_zoomed_new.jnl '"Precip (mm/day)" inputfile (precip[l=sloop]*86400) referencefile (precip[l=sloop]*86400) /pal=precip/lev=(-INF)(0.5)(1,6,1)(8,14,2)(20)(INF) /pal=pal_01/lev=(-INF)(-5)(-2,2,1)(-0.5)(-0.2)(0.2)(0.5)(5)(INF) /hlim=0:10 /hlim=-3:3 shade [X=10W:40E,Y=30N:75N]' |
| 91 | }}} |
| 92 | |
| 93 | If you then run ATLAS as above with this .cfg, you should get the maps and zonal differences only over Europe for the two simulations. |