Evaluation of LMDZOR6 with aerosols defined with 1 and 3 parameters

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 2471 and ORCHIDEE to trunk 3013 in mod.def
3. Install the model: cd modipsl/util ; ./model LMDZOR_v6

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. Activate RRTM mode in modipsl/modeles/LMDZ/makelmdz_fcm: rrtm=true

5. Compile with L=79: cd modipsl/config/LMDZOR ; gmake LMD144142-L79

To use RRTM

1. Edit COMP/lmdz.card and a choose New Physic (NP) version:

   [UserChoices]
   # Physics package to use :
   # LMDZ_Physics=AP for standard/old physics(defalut), can be used with LMDZ4_AR5 or LMDZ5/trunk sources
   # LMDZ_Physics=NPv5.17h for new physics, to be used with LMDZ5/trunk revision 1554 or later
   
   LMDZ_Physics=NPv5.17h
   

2. Activate RRTM flag with correct number of short wave bands into the corresponding physiq.def file, physiq.def_L39_NPv3.2 in my case:

   # Rayonnement
   #*********************
   # activation du nouveau code de rayonnement RRTM
   # 0 : Ancien code et 1 : RRTM (D=0)
   iflag_rrtm=1
   
   ## NSW need to be 6 when iflag_rrtm=1
   NSW=6
   

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. Add the variable output in PARAM/file_def_histmth_lmdz.xml, PARAM/file_def_histday_lmdz.xml... (i.e for each ouput frequency):
   

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

CREATE LIMIT : create boundaries based on HadISST 1870-2007 (ELIHAD07-144x142x79)

• Physics: NPv5.17h with RRTM (iflag_rrtm=1 and NSW=6).
• Config: EXPERIMENT/amip/config.card
• COMP/lmdz.card: LMDZ_Physics=NPv5.17h 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

TESTS

testLMDZ6.2.NP5.17h: new physic with "official" forcings (no RRTM)

• Description: Test this version model "basically" to check that everything is ok and nb CPU.
• Model compile with RRTM and strato forcing.
• Options: no RRTM and strato forcing
• Restart: use limits from ELIHAD07-144x142x79
• Boundary: use limit.ncfrom ELIHAD07-144x142x79
• COMP/lmdz.card: LMDZ_Physics=NP5.17h and flag_aerosol=6
• config.card:

  CalendarType=360d
  
  StartDate=1991-01-01
  EndDate=1994-12-30
  
  PeriodLength=1M
  
  RebuildFrequency=NONE
  PackFrequency=4Y
  

Results: running

testLMDZ6rrtm.2.NP5.17h: new physic+RRTM with "official" forcings

• Description: Test RRTM with default config. Evaluate time to process 4Y.
• Model compile with RRTM and strato forcing.
• Options: activate RRTM and no strato forcing
• Restart: use limits from ELIHAD07-144x142x79
• Boundary: use limit.ncfrom ELIHAD07-144x142x79
• COMP/lmdz.card: LMDZ_Physics=NP5.17h and flag_aerosol=6
• config.card:

  CalendarType=360d
  
  StartDate=1991-01-01
  EndDate=1994-12-30
  
  PeriodLength=1M
  
  RebuildFrequency=NONE
  PackFrequency=4Y
  

Results: running

testLMDZ6rrtm.2.SATO: new physic, RRTM and Sato AOD

• Description: Test RRTM + SATO AOD.
• Model compile with RRTM and strato forcing.
• Options: activate RRTM and strato forcing (AOD)
• Restart: use limits from ELIHAD07-144x142x79
• Boundary: use limit.ncfrom ELIHAD07-144x142x79
• COMP/lmdz.card: LMDZ_Physics=NP5.17h and flag_aerosol=6
• config.card:

  CalendarType=360d
  
  StartDate=1991-01-01
  EndDate=1994-12-30
  
  PeriodLength=1M
  
  RebuildFrequency=NONE
  PackFrequency=4Y
  

Results: -

TESTS

In dir:

Name	#MPI	OMP	flag_aer	PerLen?	PackFq?	Physic	RRTM	Out	Forcings	Results	Comment
TestLMDZ6.2.NP5.17h	47	8CPU	False	1M	4Y	AP	No	DEVT	Offical
TestLMDZ6rrtm.2.NP5.17h	47	8CPU	False	1M	4Y	NP5.17h	Yes	DEVT	Offical