= IPSLCM6-MR =

Sur ada, juillet 2017. Mise à jour 6 septembre 2017. Nouvelle simulation : 30 juillet 2018.

== Juillet 2018 : IPSLCM6.1.5 en MR ==

== Simulation ==

Faire 60 ans ou plus si ada est vide.

== Mise en oeuvre ==
 * IPSLCM6.1.5-LR
 * Résolution : LMDZ 256x256x279 eORCA1
  * Modification dans AA_make :
{{{
+++ AA_make     (working copy)

+IPSLCM6-MR :
+       $(M_K) ORCA1LIM3xLMD256256-L79

+ORCA1LIM3xLMD256256-L79 : libioipsl oasis3-mct xios liborchidee orca1lim3 lmdz256x256x79 verif
+       echo "ORCA1LIM3xLMD256256-L79" >.resol
+       echo "RESOL_ATM_3D=256x256x79" >>.resol
+       echo "$(LIB_MPI)" >.libmpi

+lmdz256x256x79:
+       $(M_K) lmdz RESOL_LMDZ=256x256x79
}}}
 * pas d'aérosols strato : config.def_actuel:flag_aerosol_strat=0
 * start/startphy d'une simu précédente peu importe laquelle
 * pdControl 
 * différences liés au MR :
{{{
diff  gcm.def_256x256_NPv6.0.11  gcm.def_144x142_NPv6.0.11
< day_step=1344 > day_step=672
< iphysiq=14 > iphysiq=7
< ngroup=5 > ngroup=4
< tetagdiv=1800.
< tetagrot=10800.
< tetatemp=5400.
}}}
   * LMDZ : physiq.def.6.1 avec comme seuls changements :
{{{
-cld_lc_lsc=0.000065
-cld_lc_con=0.000065
+cld_lc_lsc=0.000055
+cld_lc_con=0.000055
}}}
et aussi :
   * NEMO : comme 6.1 
 * ozone analytique (read_climoz=0)
 * aerosols nat et 1980 de Frédéric  flag_aerosol=6, ok_ade=y, ok_aie=y, ok_cdnc=y
 * départ depuis pdControl LR CM61  pour l’océan
 * Equilibrage à la louche : 
   * 220 coeurs pour NEMO, 128x8 (1 024) pour LMDZ et 1 pour XIOS. Cela prend en tout 1 248 coeurs.
   * 1 200 s par mois. 14 400 s par an. 40h pour 10 ans. 6 ans par jour. 50 000 h pour 10 ans. 
   * Voir aussi : http://forge.ipsl.jussieu.fr/igcmg/wiki/PerformancesIPSLCM6

== Simulations ==

 * essai orographie :  CM6011.3-MR-pd-split-A-A03 
   * restart 2239-12-31 de CM6011.3-MR-pd-split-A-A02 soit 50 ans avant à la fin de cette simu
   * 50 ans OK, 11 mars 2018, rsync et dods_cp fait.
   * modif : 
{{{
< ok_orolf=n
> ok_orolf=y
< sso_gkdrag=1.2
> sso_gkdrag=0.6
}}}

 * 2ème simulation longue : CM6011.3-MR-pd-split-A-A02 
  * monitoring : http://prodn.idris.fr/thredds/fileServer/ipsl_public/rgzi003/IPSLCM6-MR/DEVT/pdControl/CM6011.3-MR-pd-split-A-A02/MONITORING/index.html
  * cld_lc*=0.000205 jusque 2049 (10 ans)
  * cld_lc*=0.00018  jusque 2069 (20 ans)
  * cld_lc*= 0.00016
  * en 2121 : retour à une simulation de type split ie sans boost albedo glace.
  {{{
  rn_alb_sdry  =  0.85
  rn_alb_smlt  =  0.75
  rn_alb_idry  =  0.60
  rn_alb_imlt  =  0.50
}}}
  * Attention : rsync à faire pour rapatriement des résultats sur curie pour les atlas CLIMAF. Fait

 * Simulation longue :  CM6011.3-MR-pd-cvoro-A01. Voir [http://prodn.idris.fr/thredds/fileServer/ipsl_public/rgzi003/IPSLCM6-MR/DEVT/pdControl/CM6011.3-MR-pd-cvoro-A01/MONITORING/index.html monitoring CM6011.3-MR-pd-cvoro-A01] ; Atlas d'évaluation CESMEP: https://vesg.ipsl.upmc.fr/thredds/fileServer/IPSLFS/jservon/C-ESM-EP/CM6-MR_FirstResults_20170906_jservon/C-ESM-EP_CM6-MR_FirstResults_20170906.html
  * physiq.def de NPv6.0.11cvoro
  * A réchauffer un peu, après 2024 :
{{{
< cld_lc_lsc=0.000145
< cld_lc_con=0.000145
> cld_lc_lsc=0.000115
> cld_lc_con=0.000115
}}}
  * Intermonitoring comparatif LR MR : [http://webservices2017.ipsl.fr/interMonitoring_fromHermes/tmp/interMonitoring_plot01_B9aJEK_prod/ Intermonitoring CM6011.3-MR-pd-cvoro-A01 vs CM6011.2-pd-cvoro-01]
  * Plantage en 211708
  * Tourner par mois pour avoir facilement le point de plantage, mais pas de reproductibilité 1M+1M /1Y
 * Simulation 2 ans pour essai/validation : CM6011.3-MR-pd-A03 . Voir [http://prodn.idris.fr/thredds/fileServer/ipsl_public/rgzi003/IPSLCM6-MR/DEVT/pdControl/CM6011.3-MR-pd-A03/MONITORING/index.html monitoring CM6011.3-MR-pd-A03]  ===> Faire simulation longue avec tuning NPv6.0.11cvoro
 * Forcé LMDZOR : [http://prodn.idris.fr/thredds/fileServer/ipsl_public/rgzi003/LMDZOR/DEVT/clim/LMDZOR-6011-MR-01/MONITORING/index.html monitoring LMDZOR-6011-MR-01]

== Mise en oeuvre ==
 * IPSLCM6.0.11-LR
 * Résolution : LMDZ 256x256x279 eORCA1
  * Modification dans AA_make :
{{{
+++ AA_make     (working copy)

+IPSLCM6-MR :
+       $(M_K) ORCA1LIM3xLMD256256-L79

+ORCA1LIM3xLMD256256-L79 : libioipsl oasis3-mct xios liborchidee orca1lim3 lmdz256x256x79 verif
+       echo "ORCA1LIM3xLMD256256-L79" >.resol
+       echo "RESOL_ATM_3D=256x256x79" >>.resol
+       echo "$(LIB_MPI)" >.libmpi

+lmdz256x256x79:
+       $(M_K) lmdz RESOL_LMDZ=256x256x79
}}}
 * Poids OASIS : fait
  * Vérification poids OASIS : souci pour TAUU sur le point 0° 180°E
  * Voir : wiki:IPSLCM6/Coupling/BugBilinear
 * pas d'aérosols strato : config.def_actuel:flag_aerosol_strat=0
 * start/startphy créés par Frédéric. 
 * essai avec physique 6.0.11 figé le 12 juillet 2017
 * différences liés au MR :
{{{
diff  gcm.def_256x256_NPv6.0.11  gcm.def_144x142_NPv6.0.11
< day_step=1344 > day_step=672
< iphysiq=14 > iphysiq=7
< ngroup=5 > ngroup=4
}}}
  * explosion au Sud, changement des teta divisés par 2 :
{{{
tetagdiv=1800.
tetagrot=10800.
tetatemp=5400.
}}}
 * et les paramètres de 6011.3 :
   * LMDZ : physiq.def.6.0.11 avec comme seuls changements par rapport à 6.0.11.2 :
{{{
-cld_lc_lsc=0.000075
-cld_lc_con=0.000075
+cld_lc_lsc=0.00007
+cld_lc_con=0.00007
}}}
   * NEMO : comme 6.0.11.2 avec comme seuls changements :
{{{
-  nn_mxl0     = 1         ! type of scaling under sea-ice
+ nn_mxl0     = 2         ! type of scaling under sea-ice
+ ln_lc      = .true.     !  Langmuir cell parameterisation (Axell 2002)
+ rn_lc      =   0.20     !  coef. associated to Langmuir cells
}}}
 * ozone analytique (read_climoz=0)
 * aerosols nat et 1980 de Frédéric  flag_aerosol=6, ok_ade=y, ok_aie=y, ok_cdnc=y
 * départ depuis repos/Levitus pour l’océan
 * Equilibrage à la louche : 
   * 220 coeurs pour NEMO, 128x8 (1 024) pour LMDZ et 1 pour XIOS. Cela prend en tout 1 248 coeurs.
   * 1 200 s par mois. 14 400 s par an. 40h pour 10 ans. 6 ans par jour. 50 000 h pour 10 ans. 
   * Voir aussi : http://forge.ipsl.jussieu.fr/igcmg/wiki/PerformancesIPSLCM6
 * Si on veut éliminer les processeurs terre, suivre les explications de C. :
{{{
- cd NEMOGCM/TOOLS
- ./maketools -n MPP_PREP -m X64_CURIE(ADA) -j8
- cd MPP_PREP
- ln -sf /ccc/work/cont003/igcmg/igcmg/IGCM/OCE/NEMO/ORCA1_LIM3_PISCES/v3.6_stable/eORCA1.2_bathy_meter.nc  bathymetry.nc
- vi namelist

 &NAMSPACE
   jpk=75
   jpiglo =  362
   jpjglo =  332
   jpidta =  362
   jpjdta =  332
   nizoom = 1
   njzoom = 1

 &NAMPROC
   jprocx= 220

 &NAMFILE
   cbathy='bathymetry.nc'

- ./mpp_optimiz_zoom_nc.exe

- creation d'un fichier processor.layout - fichier texte- que tu édites. T'as tous les découpages possibles

Pour le même découpage que tu as, tu peux prendre : jpni=11 / jpnj=20 / jpnij=174
Mais il te propose comme choix optimum            : jpni=22  jpnj=10  / jpnij=171
}}}

1er atlas :
http://prodn.idris.fr/thredds/fileServer/ipsl_public/rgzi003/IPSLCM6-MR/DEVT/pdControl/CM6011.2-MR-pd-A03/MONITORING/index.html

Celui-ci existe, mais avec le bug d’interpolation de TAUU sur le point à l’equateur 180E, sans aerosols du tout, avec les performances de 3 bandes par tâche MPI et la physique collée à 6.0.11.2 :
http://prodn.idris.fr/thredds/fileServer/ipsl_public/rgzi003/IPSLCM6-MR/DEVT/pdControl/CM6011.2-MR-pd-A01/MONITORING/index.html