= Evaluation of LMDZOR6 with RRTM (test 2435) =


== Model installation ==

1. Start to import modIPSL: `svn co http://forge.ipsl.jussieu.fr/igcmg/svn/modipsl/trunk modipsl`
2. Change LMDZ revision to `testing 2435` and ORCHIDEE to `r3013` in `mod.def`
3. Install the model: `cd modipsl/util ; ./model LMDZOR6`

3b. To enable debug version: edit `modipsl/config/LMDZOR/Makefile` and add ''debug'' option:
{{{#!bash
lmdz: ../../modeles/LMDZ
   (cd ../../modeles/LMDZ ; ./makelmdz_fcm -debug -d $(RESOL_LMDZ) [...]
   (cd ../../modeles/LMDZ ; ./makelmdz_fcm -debug -d $(RESOL_LMDZ) [...]
}}}

4. Compile: `cd modipsl/config/LMDZOR ; gmake`


[[br]]

== To use aerosol strato ==

1. Edit corresponding `PARAM/*.def` file and add `flag_aerosol_strat=y`
2. Edit `COMP/lmdz.card` and configure the access to the AOD file containing a ''TAUSTRAT'' variable.
For example: 
{{{
[BoundaryFiles]
List=   (${STORE_DIR}/INPUT/SATO/taustrat.${year}.nc, taustrat.nc)
}}}
3. Check that output variable `OD_10um_STRAT` is defined and `OD550_STRAT` is activated in LMDZ (file `modipsl/modeles/LMDZ/DefLists/field_def_lmdz.xml`):[[br]]
{{{#!xml
     <field id="OD550_STRAT" long_name="Aerosol Optical depth at 550 nm " unit="1" enabled="TRUE" />
     <field id="OD_10um_STRAT" long_name="Stratospheric Aerosol Optical depth at 10 um" unit="-" />
}}}
[[br]]
4. Activate the variable output in `PARAM/file_def_hist*_lmdz.xml`:[[br]]
{{{#!xml
<field field_ref="OD_10um_STRAT" level="2" />
}}}
[[br]]

== Tests ==

=== run to create boundaries based on HadISST 1870-2007: ELIHAD07-9695x39 ===

- No RRTM and debug.
- Config: `EXPERIMENT/amip/config.card`
- ''COMP/lmdz.card'': `LMDZ_Physics=AP` and `flag_aerosol=6`
{{{#!bash
[BoundaryFiles]
List=   (/ccc/store/cont003/dsm/p25khod/dmf_import/INPUT/amipbc_sst_360x180_${year}.nc, amipbc_sst_1x1.nc), \
        (/ccc/store/cont003/dsm/p25khod/dmf_import/INPUT/amipbc_sic_360x180_${year}.nc, amipbc_sic_1x1.nc), \
}}}

- ''config.card'': 
{{{#!bash
CalendarType=360d

StartDate=1870-01-01
EndDate=2007-12-30

PeriodLength=1Y
}}}

- >100 years

Results: '''OK'''


=== run AP with "official" forcings: testLMDZ6 ===

- Model compile without RRTM and debug.
- Restart: `start.nc` and `restartphy.nc` from ''ELIHAD07-9695x39''
- Boundary: use `limit.nc`from ''ELIHAD07-9695x39''
- ''COMP/lmdz.card'': `LMDZ_Physics=AP` and `flag_aerosol=6`
- ''config.card'':
{{{#!bash
CalendarType=360d

StartDate=1991-01-01
EndDate=1994-12-30

PeriodLength=1Y

RebuildFrequency=NONE
PackFrequency=4Y
}}}

- 4 years in DEVT

Results: '''OK'''

=== run AP with "TSI" forcing only + Sato AOD copy: testLMDZ6satoNoVolc ===

- Model compile without RRTM and debug.
- Restart: `start.nc` and `restartphy.nc` from ''ELIHAD07-9695x39''
- Boundary: use `limit.nc`from ''ELIHAD07-9695x39''
- ''COMP/lmdz.card'': `LMDZ_Physics=AP` and `flag_aerosol=6`
{{{#!bash
[BoundaryFiles]
List=   (/ccc/store/cont003/dsm/p25khod/INPUT/SATO2012/taustrat.${year}.nc, taustrat.nc), \ [...]

ListNonDel= (/ccc/store/cont003/dsm/p25khod/FORC_LM_CMIP5/Formatted/TSI_Wm2_pmpip3_Kriv_WLS_0800_2000.txt, SOLARANDVOLCANOES.txt),\ (...)
}}}
- ''PARAM/physiq.def_L39-AP'': no activation of the readaerostrato
{{{#!bash
flag_aerosol_strat=n
}}}

- ''config.card'':
{{{#!bash
CalendarType=360d

StartDate=1991-01-01
EndDate=1994-12-30

PeriodLength=1M

RebuildFrequency=NONE
PackFrequency=4Y
}}}

- 4 years in DEVT

Results: ''Running''

=== run AP with TSI + Sato AOD: testLMDZ6sato ===

- Model compile without RRTM and debug.
- Restart: `start.nc` and `restartphy.nc` from ''ELIHAD07-9695x39''
- Boundary: use `limit.nc`from ''ELIHAD07-9695x39''
- ''COMP/lmdz.card'': `LMDZ_Physics=AP` and `flag_aerosol=6`
{{{#!bash
[BoundaryFiles]
List=   (/ccc/store/cont003/dsm/p25khod/INPUT/SATO2012/taustrat.${year}.nc, taustrat.nc), \ [...]

ListNonDel= (/ccc/store/cont003/dsm/p25khod/FORC_LM_CMIP5/Formatted/TSI_Wm2_pmpip3_Kriv_WLS_0800_2000.txt, SOLARANDVOLCANOES.txt),\ (...)
}}}
- ''PARAM/physiq.def_L39-AP'': no activation of the readaerostrato
{{{#!bash
flag_aerosol_strat=y
}}}

- ''config.card'':
{{{#!bash
CalendarType=360d

StartDate=1991-01-01
EndDate=1994-12-30

PeriodLength=1M

RebuildFrequency=NONE
PackFrequency=4Y
}}}

- 4 years in DEVT

Results: '''ERROR''' (in first month)


=== run AP with "official" forcings with debug mode: debugLMDZ6 ===

- Same setup as ''run AP with "official" forcings'' with compilation -debug in compilation option

Results: '''ERROR''' in first month (no debug message seems to appear in out_execution?!)
         '''ERROR''' with last LMDZ trunk (r2436)


=== TESTS ===

In dir: `/ccc/work/cont003/gen2211/lebasn/LMDZOR6_testing/modipsl.testing2435/config/LMDZOR_v6/`

||Name          ||#CPU||OMP||Debug||flag_aer||PerLen||PackFq||RebFq||Out ||Forcings||Results ||Comment                 ||
||TestLMDZ6     ||32  ||Yes||off  ||False   ||1Y    ||4Y    ||NONE ||DEVT||Offical ||Ok 4Y   ||                        ||
||TestLMDZ6ter  ||32  ||Yes||off  ||False   ||1Y    ||4Y    ||NONE ||DEVT||Offical ||   ||OD550_STRAT                        ||
||Te6.SATOnoVolc||32  ||Yes||off  ||False   ||1M    ||4Y    ||NONE ||DEVT||TSI+sato||Ok 4Y   ||Only TSI + aod file copy||
||TestLMDZ6.SATO||32  ||Yes||off  ||True    ||1M    ||4Y    ||NONE ||DEVT||TSI+sato||Crash 1M||                        ||
||TestLMDZ6.SATO||2   ||No ||off  ||True    ||1M    ||4Y    ||NONE ||DEVT||TSI+sato||5M+Crash||No aer_flag output_write||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1M    ||NONE ||TEST||TSI+sato||Crash 1M||No aer_flag output_write||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok2M+stop||Only conf_phys ||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok2M+stop||Only conf+open_phys ||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok2M+stop||conf+open_phys+readaer||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok2M+stop||conf+open_phys+readaer+lwsw-aeroar4||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok13M+stop||conf+open_phys+readaer+lwsw-aeroar4||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok3M+stop||conf+open_phys+readaer+lwsw||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Crash 1M || ||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok2M+stop||out_write-histwrite_phy()||
||TestLMDZ6.SATO||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok 4Y    ||OD_10um_STRAT in PARAM/xml||
||TesLMDZ6.SATO2||32  ||No ||off  ||True    ||1M    ||1Y    ||NONE ||TEST||TSI+sato||Ok 1Y    ||OD550_STRAT true in PARAM/xml||
||TesLMDZ6.SATO2||32  ||No ||off  ||True    ||1Y    ||4Y    ||NONE ||DEVT||TSI+sato||Ok 4Y    ||OD550_STRAT+LMDZ_NbPeriod_adjust=3||

==== Outputs with TesLMDZ6.SATO2 ====

[[Image(https://pagesperso.locean-ipsl.upmc.fr/nillod/img/LMDZ6/LMDZ6_Tops0_Anom.gif)]]

[[Image(https://pagesperso.locean-ipsl.upmc.fr/nillod/img/LMDZ6/LMDZ6_Tops0_Topl0_glbAveAnom.gif)]]