Dealing with /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl//VALID-CM622-LR.01 The analysis relies on files from the following model output directories /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/ATM/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/ATM/Output/MO/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/OCE/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/OCE/Output/MO/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/ICE/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/ICE/Output/MO/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/SRF/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01/SRF/Output/MO/ NOTE all numbers are in Sv (1 Sv = 10^6 m3/s) of freshwater NOTE 0.001 Sv = 8.8 mm sea level rise per century We are dealing with years from 1860 to 1899 , or 40 years This is good, the longer the more acurate because of some approximations made TER=land OCE=ocean SIC=sea-ice LIC=land ice ------------------------------------------------------------------------------------ --NEMO change in stores (for the records) OCE vol beg = 1.330052416e+18 OCE vol end = 1.330052141e+18 These estimates of d OCE (converted to Sv) are approximate d OCE vol (Sv) = -0.000223890950025 d SIC ice (Sv) = 0.00309059347129 d SIC snow (Sv) = 0.000385300903943 ------------------------------------------------------------------------------------ --NEMO fluxes (for the records) EMP_OCE= -2.0328228 EMP_OCE + CALVING= -2.000451146 EMP_SIC= 0.24161038 EMP_OCE + EMP_SIC + CALVING = -1.758840766 ICEBERG = 0.042934212 ICESHELF= 0.042934212 CALVING = 0.032371654 FRIVER = 1.707276 ------------------------------------------------------------------------------------ WATER CREATED BY OCE+SIC is computed as dOCE_vol + dSIC_volice + dSIC_volsno + EMP_OCE + EMP_SIC - FRIVER - ICEBERG - ICESHELF and should be as small as possible, ideally ~ 0.001 Sv over a 50 to 100 year period WATER CREATED BY OCE+SIC= -3.58110484057 ------------------------------------------------------------------------------------ --LMDZ changes in stores (for the records) should be less than 0.001 Sv for prw and even smaller for cloud water d ATM from prw (Sv) = -0.000233175458206 d ATM from prlw (Sv) = 3.1581551243e-07 d ATM from prsw (Sv) = -1.46509082477e-07 ------------------------------------------------------------------------------------ --LMDZ fluxes (for the records) PRECIP - EVAP should be < 0.0001 Sv PRECIP - EVAP (Sv) = 17.659639 - 17.658873 = 0.00076679739 Next three lines are optional Numbers come in this order OCE + SIC + TER + LIC = TOTAL RAIN by tile and total (Sv) = 12.561464 + 0.020011039 + 3.6451891 + 0.0043768446 = 16.2310409836 SNOW by tile and total (Sv) = 0.48797176 + 0.30725748 + 0.51328031 + 0.12008871 = 1.42859826 EVAP by tile and total (Sv) = 15.168288 + 0.032327558 + 2.4520315 + 0.0062251774 = 17.6588722354 Net water budget (evap-precip) by land surface type in LMDz BILO TER (Sv)= -1.706438 BILO OCE (Sv)= 2.118853 BILO SIC (Sv)= -0.29494096 BILO LIC (Sv)= -0.11824037 The BILO tot diag below should be similar (+/-5e-6) to PRECIP-EVAP above but with opposite sign BILO tot (Sv)= -0.00076633 The BILO OCE+SIC diag below should be the same (within <0.001 Sv) as EMP_OCE+EMP_ICE+CALVING above BILO OCE+SIC (Sv)= 1.82391204 ------------------------------------------------------------------------------------ Water created in LMDZ should be < 0.001 Sv WATER_CREATED_IN_ATM (Sv)= 0.000533791238224 ------------------------------------------------------------------------------------ --LIC change in store and fluxes Atmospheric calving can be different from that reaching the ocean because there is a buffer CALVING atm (Sv)= 0.1182194 d LIC (Sv) = 2.097e-05 d LIC diagnosed from msnow (Sv) = 9.97734890284e-05 ------------------------------------------------------------------------------------ --SECHIBA change in stores dSTOCK Soil Hum (Sv) = -0.0010406314 dSTOCK Intercept (Sv) = 4.4341497e-09 dSTOCK Snow (Sv) = 0.00017943763 dSTOCK Stream (Sv) = -5.5093285e-06 dSTOCK FastR (Sv) = -2.3736166e-06 dSTOCK SlowR (Sv) = -3.8532624e-05 dSTOCK Lake (Sv) = 0.00011069203 dSTOCK Pond (Sv) = 0 dSTOCK Flood (Sv) = 0 dSTOCK Total (Sv) = -0.00079691287495 ------------------------------------------------------------------------------------ --SECHIBA fluxes EVAP SRF (Sv)= 3.0928526 This diagnoses how different is EVAP from the sum of known terms in SECHIBA The terms are evap bare soil, evap flooding plains, sublimation, transpiration, interception The difference with EVAP SRF should be small <0.0001 Sv EVAP TERMS SRF= 1.1426879 + 0 + 0.10058734 + 1.238021 + 0.61155639 EVAP TOTAL SRF= 3.09285263 RAINFALL SRF= 3.6451891 SNOWFALL SRF= 0.51328031 The SECH NET SRF (precip - evap over land) should be the same as BILO TER above but opposite sign SECH NET SRF= 1.706438 FLOW to the ocean should the same (<0.0001 Sv) as FRIVER above FLOW to ocean (Sv)= 1.7072795 ------------------------------------------------------------------------------------ Created water by ORCHIDEE should be < 0.0001 Sv WATER CREATED BY ORCHIDEE= 1.7072795 - 1.706438 + -0.00079691287495 = 4.45871250502e-05 ------------------------------------------------------------------------------------ SUMMARY d ATM (Sv) = -0.000233006151776 d LIC (Sv) = 2.097e-05 d SRF (Sv) = -0.00079691287495 d OCE (Sv) = -0.000223890950025 d SIC ice (Sv) = 0.00309059347129 d SIC snow (Sv) = 0.000385300903943 Two different diagnostics of created water rate in the coupled model should be < 0.003 Sv in the current state of affairs WATER CREATED TOTAL (Sv) = 0.00224305439848 WATER CREATED TOTAL (Sv) = -3.58052646221 ------------------------------------------------------------------------------------