Evaluation of LMDZOR6

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:

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

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):
   

        <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="-" />
   

4. Activate the variable output in PARAM/file_def_hist*_lmdz.xml:
   

   <field field_ref="OD_10um_STRAT" level="2" />
   

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

  [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:

  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.ncfrom ELIHAD07-9695x39
• COMP/lmdz.card: LMDZ_Physics=AP and flag_aerosol=6
• config.card:

  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.ncfrom ELIHAD07-9695x39
• COMP/lmdz.card: LMDZ_Physics=AP and flag_aerosol=6

  [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

  flag_aerosol_strat=n
  

• config.card:

  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.ncfrom ELIHAD07-9695x39
• COMP/lmdz.card: LMDZ_Physics=AP and flag_aerosol=6

  [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

  flag_aerosol_strat=y
  

• config.card:

  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