Changes between Version 8 and Version 9 of IPSLCM6/calving

Timestamp:

01/09/21 16:52:39 (3 years ago)

Author:

oboucher

Comment:

--

IPSLCM6/calving

@@ -9 +9 @@
 Détail : 1 mois de simu 
 
-Script cdo : on pondère par tmaskutil et l'aire des mailles pour l'océan. On pondère par l'aire des mailles et pourc_oce+pourc_sic pour l'atmosphère (moyenne mensuelle mais la somme des deux est un invariant). 
+Script cdo : on pondère par maskutil_T (eORCA025.1) et l'aire des mailles pour l'océan. On pondère par l'aire des mailles et pourc_oce+pourc_sic pour l'atmosphère (moyenne mensuelle mais la somme des deux est un invariant). 
 
 {{{ 
@@ -25 +25 @@
 filesnow2=${dir}/CPL/Output/MO/${exp}_${yr1}0101_${yr2}0131_OTotSnow.nc
 fileoce=${dir}/OCE/Output/MO/${exp}_${yr1}0101_${yr2}0131_1M_grid_T.nc
-maskoce=/ccc/work/cont003/igcmg/igcmg/IGCM/OCE/NEMO/eORCA025.4/GRIDS/eORCA025_mesh_mask.nc
+maskoce=/ccc/work/cont003/igcmg/igcmg/IGCM/OCE/NEMO/eORCA025.1/eORCA025.1_coordinates_mask.nc
 
 sumevap1=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -seltimestep,1/10 -selname,COTOTEVA ${fileevap1} -seltimestep,1 -expr,'toto=aire/100.*(pourc_oce+pourc_sic)' ${fileatm}`
-sumevap2=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -mul -seltimestep,1/10 -selname,OTotEvap ${fileevap2} -selvar,tmaskutil ${maskoce} -gridarea ${fileoce}`
+sumevap2=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -mul -seltimestep,1/10 -selname,OTotEvap ${fileevap2} -selvar,maskutil_T ${maskoce} -gridarea ${fileoce}`
 sumrain1=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -seltimestep,1/10 -selname,COTOTRAI ${filerain1} -seltimestep,1 -expr,'toto=aire/100.*(pourc_oce+pourc_sic)' ${fileatm}`
-sumrain2=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -mul -seltimestep,1/10 -selname,OTotRain ${filerain2} -selvar,tmaskutil ${maskoce} -gridarea ${fileoce}`
+sumrain2=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -mul -seltimestep,1/10 -selname,OTotRain ${filerain2} -selvar,maskutil_T ${maskoce} -gridarea ${fileoce}`
 sumsnow1=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -seltimestep,1/10 -selname,COTOTSNO ${filesnow1} -seltimestep,1 -expr,'toto=aire/100.*(pourc_oce+pourc_sic)' ${fileatm}`
-sumsnow2=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -mul -seltimestep,1/10 -selname,OTotSnow ${filesnow2} -selvar,tmaskutil ${maskoce} -gridarea ${fileoce}`
+sumsnow2=`cdo outputf,%12.10g,10 -divc,1.e9 -fldsum -mul -mul -seltimestep,1/10 -selname,OTotSnow ${filesnow2} -selvar,maskutil_T ${maskoce} -gridarea ${fileoce}`
 
 echo 'evap in=' $sumevap1
@@ -46 +46 @@
 {{{ 
 evap in= 16.58377264 16.5842436 16.73348733 16.70182978 16.66627324 16.63719809 16.59929387 16.55878592 16.52697635 16.4860888
-evap out= 16.55351427 16.55350019 16.70269787 16.67088149 16.63496352 16.60587634 16.56772969 16.52697626 16.49513478 16.45415434
+evap out= 16.58377268 16.58424364 16.73348737 16.70182981 16.66627328 16.63719813 16.59929391 16.55878595 16.52697638 16.48608884
 rain in= 12.76418216 11.79569265 11.90206466 12.29730812 12.41094174 12.47357032 12.53198387 12.55976865 12.5649098 12.58444113
-rain out= 12.74270543 11.77684503 11.88402995 12.27952764 12.395365 12.45955016 12.51887749 12.5477753 12.55303025 12.57239723
+rain out= 12.76418217 11.79569266 11.90206466 12.29730811 12.41094173 12.47357031 12.53198387 12.55976866 12.56490981 12.58444114
 snow in= 0.6307305996 0.613739933 0.6188579544 0.6443643258 0.6604737915 0.6693078788 0.6729603416 0.6815849762 0.684608096 0.6879445765
-snow out= 0.6216787337 0.6042128991 0.608490469 0.6332710421 0.6475660592 0.6557678674 0.6592571085 0.6672780528 0.6703978744 0.6741682585
-}}}
-
-Les écarts sont très importants : 4e décimale et 0.2% pour rain et evap, 2e décimale et 2% pour snow.  
-Bien sûr, le biais est moindre pour P+S-E car il y a compensation d'erreur. 
+snow out= 0.6307305919 0.6137399242 0.6188579491 0.6443643249 0.6604737919 0.6693078792 0.672960342 0.6815849766 0.6846080966 0.6879445774
+
+}}}
+
+Les écarts sont très faibles. 
+Attention ça ne marche pas si on prend le tmaskutil de eORCA025.4. 
 
 == Exemple avec le VLR - script python ==