Dealing with /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl//TEST-CM72-SIMPLE-ROUTING.10 The analysis relies on files from the following model output directories /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/ATM/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/ATM/Output/MO/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/OCE/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/OCE/Output/MO/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/ICE/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/ICE/Output/MO/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/SRF/Output/DA/ /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.10/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 1869 , or 10 years WARNING: at least 30 years are needed for a good accuracy 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.330024653e+18 OCE vol end = 1.329990018e+18 These estimates of d OCE (converted to Sv) are approximate d OCE vol (Sv) = -0.112792189878 d SIC ice (Sv) = 0.0211135327394 d SIC snow (Sv) = 0.00180724224058 ------------------------------------------------------------------------------------ --NEMO fluxes (for the records) EMP_OCE= 2.1894105 EMP_OCE + CALVING= 2.231965568 EMP_SIC= -0.16968247 EMP_OCE + EMP_SIC + CALVING = 2.062283098 ICEBERG = 0.029660049 ICESHELF= 0.029660049 CALVING = 0.042555068 FRIVER = 1.8640181 ------------------------------------------------------------------------------------ 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= 0.00651841710198 ------------------------------------------------------------------------------------ --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.000131488394216 d ATM from prlw (Sv) = 2.58783717022e-06 d ATM from prsw (Sv) = 8.68861954592e-07 ------------------------------------------------------------------------------------ --LMDZ fluxes (for the records) PRECIP - EVAP should be < 0.0001 Sv PRECIP - EVAP (Sv) = 17.481745 - 17.482133 = -0.00038818458 Next three lines are optional Numbers come in this order OCE + SIC + TER + LIC = TOTAL RAIN by tile and total (Sv) = 12.439897 + 0.011122027 + 3.603294 + 0.0061488786 = 16.0604619056 SNOW by tile and total (Sv) = 0.43288859 + 0.22097882 + 0.60874447 + 0.15866986 = 1.42128174 EVAP by tile and total (Sv) = 15.139333 + 0.027685598 + 2.3081109 + 0.0070036605 = 17.4821331585 Net water budget (evap-precip) by land surface type in LMDz BILO TER (Sv)= -1.9039275 BILO OCE (Sv)= 2.2665482 BILO SIC (Sv)= -0.20441525 BILO LIC (Sv)= -0.15781508 The BILO tot diag below should be similar (+/-5e-6) to PRECIP-EVAP above but with opposite sign BILO tot (Sv)= 0.00039037 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)= 2.06213295 ------------------------------------------------------------------------------------ Water created in LMDZ should be < 0.001 Sv WATER_CREATED_IN_ATM (Sv)= -0.000516216275091 ------------------------------------------------------------------------------------ --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.10190774 d LIC (Sv) = 0.05590734 d LIC diagnosed from msnow (Sv) = 0.0231983713534 ------------------------------------------------------------------------------------ --SECHIBA change in stores dSTOCK Soil Hum (Sv) = -0.001108905 dSTOCK Intercept (Sv) = -2.1608219e-07 dSTOCK Snow (Sv) = 0.0037836357 dSTOCK Stream (Sv) = 0.0 dSTOCK FastR (Sv) = 0.0 dSTOCK SlowR (Sv) = 0.0 dSTOCK Lake (Sv) = 0.0 dSTOCK Pond (Sv) = 0.0 dSTOCK Flood (Sv) = 0.0 dSTOCK Total (Sv) = 0.00267451461781 ------------------------------------------------------------------------------------ --SECHIBA fluxes EVAP SRF (Sv)= 2.9291907 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.0976643 + 0 + 0.10026724 + 1.1613162 + 0.56994296 EVAP TOTAL SRF= 2.9291907 RAINFALL SRF= 3.603294 SNOWFALL SRF= 0.60874447 The SECH NET SRF (precip - evap over land) should be the same as BILO TER above but opposite sign SECH NET SRF= 1.9039275 FLOW to the ocean should the same (<0.0001 Sv) as FRIVER above FLOW to ocean (Sv)= 1.8640114 ------------------------------------------------------------------------------------ Created water by ORCHIDEE should be < 0.0001 Sv WATER CREATED BY ORCHIDEE= 1.8640114 - 1.9039275 + 0.00267451461781 = -0.0372415853822 ------------------------------------------------------------------------------------ SUMMARY d ATM (Sv) = -0.000128031695091 d LIC (Sv) = 0.05590734 d SRF (Sv) = 0.00267451461781 d OCE (Sv) = -0.112792189878 d SIC ice (Sv) = 0.0211135327394 d SIC snow (Sv) = 0.00180724224058 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.0314175919753 WATER CREATED TOTAL (Sv) = -0.0312393845553 ------------------------------------------------------------------------------------