Test IPSLCM6.0.10 RRTM with CMIP6 volcanic forcings

Model installation

1. Start to import modIPSL: svn co http://forge.ipsl.jussieu.fr/igcmg/svn/modipsl/trunk modipsl
2. Install the model: cd modipsl/util ; ./model IPSLCM6.0.10-LR
3. Add in code some diagnostics: procedure is describded ​here.
4. Compile: cd modipsl/config/LMDZOR ; gmake IPSLCM6-LR

NOTE: file_def_XXX.xml are taken from the model and not in PARAM subdir in default config. Change directory in lmdz.card.

CM6010.LOWTROPO.TOPSW

This model version has got a 100mb troupopause shifted.
To do so, modify the LMDZ/libf/phylmd/stratosphere_mask.F90 file as below around line 63 (code between !NL):

  call twmo(klev, t, p, psrf, zsrf, plimu, pliml, gamma, ptrp, ttrp, ztrp)
!NL: Shift tropopause lower (10.000Pa)
  tp(i)=ptrp+10000.0 ! test to decrease tropopause of 100mb
!NL: END
  ttp(i)=ttrp
  ztp(i)=ztrp
!NL: Shift tropopause lower (10.000Pa)
!  IF (ptrp.lt.0.0) THEN
  IF (tp(i).lt.0.0) THEN
!NL: END
    tperr = tperr+1
  ENDIF

CM6010.NOSTMASK.TOPSW

This model version has no troupopause mask.
To do so, modify the LMDZ/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90 file as below around line 284 (code between !NL):

!--total vertical aod at the 5 SW wavelengths
!--for now use band 3 AOD into all 5 wavelengths
!--it is only a reasonable approximation for 550 nm (wave=2)
    band=3
    DO i=1, klon
    DO k=1, klev
!NL: Remove stratomask
!      IF (stratomask(i,k).GT.0.999999) THEN ! To comment
!NL: END
        DO wave=1, nwave_sw
          tausum_aero(i,wave,id_STRAT_phy)=tausum_aero(i,wave,id_STRAT_phy)+tau_aer_strat(i,k,band)
        ENDDO
!NL: Remove stratomask
!      ENDIF !comment
!NL: END
    ENDDO
    ENDDO

!MK
    DO i=1, klon
       DO k=1, klev
!NL: Remove stratomask
!          IF (stratomask(i,k).LT.0.999999) THEN !comment
!             tau_aer_strat(i,k,band)=1.e-15 !comment
!          ENDIF !comment
!NL: END
          tau_aer_stratoo(i,k)=tau_aer_strat(i,k,band)
       ENDDO
    ENDDO
    print *,'readaerosolstrato2_rrtm: tau_aer_stratoo(',klon,',',klev,')=', tau_aer_stratoo(klon,klev)

TESTS

All outputs were made with IPSLCM6.0.10.1, physic=NPv6.0.10. All simulations are with actuel type of atmosphere and restart is from 1989-12-31 of the reference pdControl simulation CM6010.SnwCnd?-LR-pdctrl-04. Anomalies are performed by substracting control run values (ie CM6010.TOPSW.CTL) in 1991-1992 period.

Name of simulations and Nomenclature:

• CM6010.TOPSW3: CMIP6 forcing and double call to physic to get only volcanic contribution (topswad).
  
• CM6010.TOPSW.CTL: corresponds to CM6010.TOPSW3 with CMIP6 volcanic forcing set 1960 value.
  
• CM6010.TOPSW.pdCTL: corresponds to CM6010.TOPSW3 with CMIP6 volcanic forcing historical average (tau*wstrat.2D.ave.nc).
  
• CM6010.LOWTROPO.TOPSW: corresponds to CM6010.TOPSW3 with a traupopause shift of 100mb.
  
• CM6010.LOWTROPO.TOPSW.OZ: corresponds to CM6010.LOWTROPO.TOPSW with annual Ozone values instead of climato.
  
• CM6010.LOWTROPO.pdCTRL: corresponds to CM6010.TOPSW.pdCTRL with a traupopause shift of 100mb.
  
• CM6010.NOSTMASK.TOPSW: corresponds to CM6010.TOPSW3 with no stratomask.
  
• CM6010.NOSTMASK.pdCTRL: corresponds to CM6010.TOPSW.pdCTRL with no stratomask.
  
• CM6010.PinaOB.OZ*: corresponds to CM6010.TOPSW3 with O. Boucher new strato forcings (5 members).
  

Name	#CPU	flag_aer	PerLen?	PackFq?	Physic	Restart	RRTM	Out	Volc forcing	Results	Comment
CM6010.TOPSW.CTL	1049	2	1M	1Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	DEVT	1960 CMIP6	ok	-
CM6010.TOPSW.pdCTL*	1049	2	1M	5Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	DEVT	CMIP6 histAVE		-
CM6010.TOPSW3	1049	2	1M	1Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	TEST	CMIP6	ok	-
CM6010.TOPSW4	1049	2	1Y	1Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	TEST	CMIP6	ok	-
CM6010.LOWTROPO.pdCTRL	1049	2	1M	5Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	DEVT	CMIP6 histAVE		tropo+10000Pa
CM6010.LOWTROPO.TOPSW	1049	2	1M	5Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	DEVT	CMIP6		tropo+10000Pa
CM6010.LOWTROPO.TOPSW.OZ	1049	2	1M	5Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	DEVT	CMIP6		(tropo+10000Pa)&annualOzone
CM6010.NOSTMASK.pdCTRL	697	2	1M	5Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	DEVT	CMIP6 histAVE		no stratomask
CM6010.NOSTMASK.TOPSW	1049	2	1M	5Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	DEVT	CMIP6		no stratomask
CM6010.PinaOB.OZ*	1049	2	1M	1Y	NPv6.0.10	CM6010.SnwCnd?-LR-pdctrl-04	Yes	TEST	CMIP6		OB new tau*lstrat forcings

Members restart infos

Member nb	restart simulation	restart year
	CM6010.SnwCnd?-LR-pdctrl-04	1989-12-31
2	CM6010.SnwCnd?-LR-pdctrl-04	1994-12-31
3	CM6010.SnwCnd?-LR-pdctrl-04	1999-12-31
4	CM6010.SnwCnd?-LR-pdctrl-04	2004-12-31
5	CM6010.SnwCnd?-LR-pdctrl-04	2009-12-31

CM6010.TOPSW3

CM6010.LOWTROPO.TOPSW

Ozone climato

Ozone annual values

CM6010.NOSTMASK.TOPSW