= Conservation de l'évaporation : où est l'erreur ? = 

== Exemple avec le VLR == 

Expérience : /ccc/work/cont003/ra5424/p25khod/IPSLCM6.5_work_15DEC/modipsl/config/IPSLCM6/CM65-VLR-pd-tn-01-test

Sorties : /ccc/store/cont003/gencmip6/p25khod/IGCM_OUT/IPSLCM6/DEVT/pdControl/CM65-VLR-pd-tn-01-test

Détails : 12 mois de simu, 18500101_18501231, cpl_old_calving=no, bug fixé dans LMDZ pour le calcul du flux cumulé par bande de latitude 

Problème : le flux de évap qui repart du coupleur diffère un peu trop du flux de évap qui arrive, alors que les flux de pluie et neige sont très bien conservés 

Fichier d'aires NEMO : ceux du fichier de sortie grid_T de NEMO 

Fichier masque NEMO : /ccc/cont003/home/lmd/oboucher/BILAN_EAU/maskutil_T.nc

Diagnostic 1 - script python 

{{{
import xarray as xr
import numpy as np

exp="CM65-VLR-pd-tn-01-test"
dir="/ccc/store/cont003/gencmip6/p25khod/IGCM_OUT/IPSLCM6/DEVT/pdControl/"+exp+"/"

file_atm=dir+"ATM/Output/MO/"+exp+"_18500101_18501231_1M_histmth.nc"

file_calving_1=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_COCALVIN.nc"
file_calving_2=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_OCalving.nc"

file_evap_1=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_COTOTEVA.nc"
file_evap_2=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_OTotEvap.nc"

file_rain_1=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_COTOTRAI.nc"
file_rain_2=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_OTotRain.nc"

file_snow_1=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_COTOTSNO.nc"
file_snow_2=dir+"CPL/Output/MO/"+exp+"_18500101_18501231_OTotSnow.nc"

fileoce=dir+"OCE/Output/MO/"+exp+"_18500101_18501231_1M_grid_T.nc"
maskoce="/ccc/cont003/home/lmd/oboucher/BILAN_EAU/maskutil_T.nc"

#--ouverture des xarray
xr_atm=xr.open_dataset(file_atm)
xr_calving_1=xr.open_dataset(file_calving_1)
xr_calving_2=xr.open_dataset(file_calving_2)
xr_evap_1=xr.open_dataset(file_evap_1)
xr_evap_2=xr.open_dataset(file_evap_2)
xr_rain_1=xr.open_dataset(file_rain_1)
xr_rain_2=xr.open_dataset(file_rain_2)
xr_snow_1=xr.open_dataset(file_snow_1)
xr_snow_2=xr.open_dataset(file_snow_2)
xr_ocean=xr.open_dataset(fileoce)
xr_mask=xr.open_dataset(maskoce)

#--extraction des array numpy - conversion en longdouble pour les sommes
atm_aire0=np.longdouble(xr_atm['aire'].values)
atm_oce=np.longdouble(xr_atm['pourc_oce'].values)
atm_sic=np.longdouble(xr_atm['pourc_sic'].values)

wbilo_oce=np.longdouble(xr_atm['wbilo_oce'].values)
wbilo_sic=np.longdouble(xr_atm['wbilo_sic'].values)

atm_aire=atm_aire0*(atm_oce+atm_sic)/100.

calving_1=np.longdouble(xr_calving_1['COCALVIN'].values)
calving_2=np.longdouble(xr_calving_2['OCalving'].values)

evap_1=np.longdouble(xr_evap_1['COTOTEVA'].values)
evap_2=np.longdouble(xr_evap_2['OTotEvap'].values)

rain_1=np.longdouble(xr_rain_1['COTOTRAI'].values)
rain_2=np.longdouble(xr_rain_2['OTotRain'].values)

snow_1=np.longdouble(xr_snow_1['COTOTSNO'].values)
snow_2=np.longdouble(xr_snow_2['OTotSnow'].values)

oce_area=np.longdouble(xr_ocean['area'].values)

mask=np.longdouble(xr_mask['maskutil_T'][:,:].values)

#--totaux sur la grille - conversion en Sv
tot_calving_1=np.sum(calving_1,axis=(1,2))/1.e9
tot_calving_2=np.sum(calving_2*oce_area*mask,axis=(1,2))/1.e9
print('calving in  Sv =',tot_calving_1[0:10])
print('calving out Sv =',tot_calving_2[0:10])

}}}

Résultats pour le calving : 
{{{

('calving in  Sv =', array([ 0.0,  0.025110938,  0.07871894,  0.060956317,  0.066686016,
        0.032709929,  0.049833877,  0.033137503,  0.032376687,  0.059238586], dtype=float128))
('calving out Sv =', array([ 0.0,  0.025110938,  0.07871894,  0.060956317,  0.066686016,
        0.032709929,  0.049833878,  0.033137503,  0.032376687,  0.059238586], dtype=float128))

}}}

Tout va bien, c'est la même chose en input et en output 

{{{
tot_evap_1=np.sum(evap_1*atm_aire[0,:,:],axis=(1,2))/1.e9
tot_evap_2=np.sum(evap_2*oce_area*mask,axis=(1,2))/1.e9
tot_rain_1=np.sum(rain_1*atm_aire[0,:,:],axis=(1,2))/1.e9
tot_rain_2=np.sum(rain_2*oce_area*mask,axis=(1,2))/1.e9
tot_snow_1=np.sum(snow_1*atm_aire[0,:,:],axis=(1,2))/1.e9
tot_snow_2=np.sum(snow_2*oce_area*mask,axis=(1,2))/1.e9

print('oceanic area from grid_T file')
print('rain in  Sv =',tot_rain_1[0:10])
print('rain out Sv =',tot_rain_2[0:10])
print('snow in  Sv =',tot_snow_1[0:10])
print('snow out Sv =',tot_snow_2[0:10])
print('evap in  Sv =',tot_evap_1[0:10])
print('evap out Sv =',tot_evap_2[0:10])
print('evap diff mSv=',(tot_evap_1[1:10]-tot_evap_2[1:10])*1000.)

}}}

Résultats : 
{{{
('rain in  Sv =', array([ 0.0,  3.4585411,  8.8186534,  11.318086,  11.570535,  11.836867,
        11.357751,  11.295551,  11.242278,  11.49349], dtype=float128))
('rain out Sv =', array([ 0.0,  3.4585411,  8.8186534,  11.317966,  11.570457,  11.836863,
        11.35775,  11.295479,  11.242256,  11.493486], dtype=float128))
('snow in  Sv =', array([ 0.0,  0.43405517,  0.48368972,  0.41251261,  0.36736733,
        0.39909007,  0.46146092,  0.3871201,  0.33174664,  0.45059974], dtype=float128))
('snow out Sv =', array([ 0.0,  0.43405517,  0.48368972,  0.41251261,  0.36736733,
        0.39909007,  0.46146092,  0.38712011,  0.33174664,  0.45059974], dtype=float128))
('evap in  Sv =', array([ 0.0,  12.638699,  12.227281,  13.689618,  14.567473,  14.940661,
        14.904491,  14.177904,  13.691065,  13.653616], dtype=float128))
('evap out Sv =', array([ 0.0,  12.638585,  12.226845,  13.689247,  14.566987,  14.940099,
        14.903804,  14.177131,  13.690155,  13.652478], dtype=float128))
('evap diff mSv=', array([ 0.11373978,  0.43656301,  0.37129745,  0.48537692,  0.5611987,
        0.68672948,  0.77286054,  0.91008484,  1.1375948], dtype=float128))
}}}

Rain et snow sont très bien conservés, mais ce n'est pas le cas de Evap avec des différences de 1 mSv.