20 | | It appears that the post-treatment procedes in two steps: step 1 is to combine the raw data into files containing all the years you want to average over (given by SeasonalFrequency in the config.card). Step 2 is to compare these files (located in JOBNAME/ATM/Analyse/SE/ in the IGCM directory tree...what comes before this depends on what you put in the config.card and what system you are running on) to a set of reference data. |
| 21 | http://wiki.ipsl.jussieu.fr/IGCMG/Outils/ferret/Atlas |
| 22 | |
| 23 | |
| 24 | It appears that the post-treatment proceeds in two steps: step 1 is to combine the raw data into files containing all the years you want to average over (given by SeasonalFrequency in the config.card). Step 2 is to compare these files (located in JOBNAME/ATM/Analyse/SE/ in the IGCM directory tree...what comes before this depends on what you put in the config.card and what system you are running on) to a set of reference data. Step 2 then compares this data to the reference simulation or observations, depending on what you are trying to do. |
| 25 | |
| 26 | As an example, let's say I have two simulations that I want to compare. |
| 27 | /ccc/store/cont003/dsm/p529grat/IGCM_OUT/LMDZOR/PROD/amip/MCOUP/ATM/Analyse/SE/MCOUP_SE_1980_1981_1M_histmth.nc |
| 28 | /ccc/store/cont003/dsm/p529grat/IGCM_OUT/LMDZOR/PROD/amip/TNV32/ATM/Analyse/SE/TNV32_SE_1980_1981_1M_histmth.nc |
| 29 | |
| 30 | Notice that both of these runs have already been processed by the seasonal frequency script, and therefore the monthly simulation data has been combined into a single files covering two years (the length of SeasonalFrequency in config.card). |
| 31 | |
| 32 | If you have loaded all the ipsl paths on, for example, Curie, you should be able to find all the ferret scripts (*.jnl) and ATLAS configuration files (*.cfg) in the $FER_ATLAS directory. |
| 33 | |
| 34 | {{{ |
| 35 | ls $FER_ATLAS/ |
| 36 | }}} |
| 37 | |
| 38 | We are interested in comparing two simulations, which according to the webpage above, means we need a diag2 ATLAS script. If you open the LMDZ*diag2.cfg files with a text editor, you can look at all the variables they compare and decide if you want to run all of them, or just one. We create an empty directory somewhere to hold the output, and go to that directory. |
| 39 | |
| 40 | {{{ |
| 41 | mkdir ATLAS |
| 42 | cd ATLAS |
| 43 | }}} |
| 44 | |
| 45 | Let's copy one of the ATLAS scripts here and modify it a bit. |
| 46 | |
| 47 | {{{ |
| 48 | cp $FER_ATLAS/atlas_p_LMDZ3.3_diag2_a.cfg test.cfg |
| 49 | ref=/ccc/store/cont003/dsm/p529grat/IGCM_OUT/LMDZOR/PROD/amip/TNV32/ATM/Analyse/SE/TNV32_SE_1980_1981_1M_histmth.nc |
| 50 | sed -e "s%your_reference_file=.*%your_reference_file=${ref}%" $FER_ATLAS/atlas_p_LMDZ3.3_diag2_a.cfg > atlas_p_LMDZ3.3_diag2_a.cfg |
| 51 | }}} |
| 52 | |
| 53 | Now my 'reference' simulation is TNV32. If you decide that you want to run multiple .cfg scripts, you will have to run the "sed" line for each script you want to use. |
| 54 | |
| 55 | {{{ |
| 56 | atlas -o atlas_model_vs_model atlas_p_LMDZ3.3_diag2_*.cfg /ccc/store/cont003/dsm/p529grat/IGCM_OUT/LMDZOR/PROD/amip/MCOUP/ATM/Analyse/SE/MCOUP_SE_1980_1981_1M_histmth.nc |
| 57 | }}} |
| 58 | |
| 59 | This command should be run in the same directory where I just created the modified .cfg files. With the "atlas_p_LMDZ3.3_diag2_*.cfg", it will run all the .cfg files that it finds matching that expression. The last argument is the name of the file I want to compare to my "reference" simulation, which I declared above. |
| 60 | |
| 61 | The output of this command is a series of folders all grouped into "atlas_model_vs_model/". In there you have pdf and html summaries of all the variables, plus individual directories for each variable if you just want to look at one of them. |