Changes between Version 5 and Version 6 of Documentation/UserGuide/VerifyClimate


Ignore:
Timestamp:
07/21/14 15:17:26 (7 years ago)
Author:
mmcgrath
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Documentation/UserGuide/VerifyClimate

    v5 v6  
    8888 
    8989{{{ 
    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]' 
     90echo $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=-20E:40E,Y=30N:75N]' 
    9191}}} 
    9292 
    9393If you then run ATLAS as above with this .cfg, you should get the maps and zonal differences only over Europe for the two simulations. 
    9494 
     95NOTE: Do not use 20W; use -20E instead to mean 20 degrees west longitude.  While it is true that 20W is the same as -20E, Ferret seems to have trouble with some of the commands if you use 20W, and the maps come out strange.  I assume the same is true for 20S and -20N. 
     96 
    9597== Reference Data == 
    9698 
    9799The script which I've found works best for reference data are atlas_p_LMDZ4_diag1_a.cfg and atlas_p_LMDZ4_diag1_e.cfg, since they include comparing the temperature at 2 m (t2m) against CRU.  You can use these scripts in exactly the same way as the scripts above, except that you do not need to change your_reference_file; the reference file is always the observational data.  It seems that these scripts already include 2variables1diffandzonal_map_zoomed_new.jnl, so if you want to zoom you can skip that step.  You will still have to add the two arguments about shade and the zoom region, though. 
     100 
     101