========== 20170103: 1. Updated higgins rain data through Dec 31, 2016. 2. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170104: 1. Prepared for meeting with Randy this afternoon. 2. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. 3. Processed GMAO GEOS5 seasonal forecast for January. ========== 20170105: 1. The CF0180 run finished the 13th cycle. Showed Randy the carbon state at the end of the 12th cycle yesterday. Randy said FPAR looks good enough and FPAR is most important for Yehui's forecast, so we can stop the CF0180 run for now when the 13th cycle is completed. Checked and processed the output: discover15:/discover/nobackup/fzeng/Catchment/M2n5P/CF0180x6C/CF0180x6C/output/CF0180x6C_GLOBAL/rs/ens0000 > ls -l Y2016/M01/ total 2999168 -rw-r--r-- 1 fzeng g0620 3075334887 2017-01-04 17:42 CF0180x6C.ens0000.catchcn_ldas_rst.20160101_0000z discover15:/discover/nobackup/fzeng/Catchment/M2n5P/CF0180x6C/CF0180x6C/output/CF0180x6C_GLOBAL/rs/ens0000 > ls -l ../../cat/ens_avg/Y2015/M12 total 214528 -rw-r--r-- 1 fzeng g0620 219664000 2017-01-04 17:42 CF0180x6C.ens_avg.ldas_tile_monthly_out.201512.bin Process the output (covert from tile-space to grid-space): ~/Catchment/CF0180x6C > tile2grid_CF0180x6C_monthly CF0180x6C Check the output on GrADS. Looks good. cd /discover/nobackup/fzeng/Catchment/M2n5P/CF0180x6C/ mv CF0180x6C CF0180x6C_13 cd CF0180x6C_13 rm -rf lenkf_orig 2. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170106: 1. For Eunjee: computed this following map which is the fraction of land in each grid cell for the 144x91 tile system used in our CO2 experiments /discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb/land_fraction program created for this purpose: /home/fzeng/Catchment/M2n5P_m0001/calculate_land_frac.m The read_tilecoord.m, read_tile_data.m, read_tilegrids.m, tile2grid.m used in calculate_land_frac.m were copied from /discover/nobackup/fzeng/offline_code/LDASsa_m3-15_2-CN/src/Components/GEOSlana_GridComp with two minor modifications. 2. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. 3. The e0004s_preindustrialCO2 finished the 6th cycle. /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l Y2016/M01 total 7354336 -rw-r--r-- 1 fzeng g0620 7547615847 2017-01-06 01:00 e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z dali12:/discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l ../../cat/ens_avg/Y2015/M12 total 526496 -rw-r--r-- 1 fzeng g0620 539112640 2017-01-06 01:00 e0004s_preindustrialCO2.ens_avg.ldas_tile_monthly_out.201512.bin Process the output (covert from tile-space to grid-space): ~/Catchment/SMAP_M09 > tile2grid_ease_spinup_monthly e0004s_preindustrialCO2 Check the output on GrADS. Looks good. cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/ mv e0004s_preindustrialCO2 e0004s_preindustrialCO2_06 cd e0004s_preindustrialCO2_06 rm -rf lenkf_orig cd output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/ /bin/rm -rf Y200[1-9] Y201[0-5] mkdir -p Y2001/M01 cd Y2001/M01 ln -s ../../Y2016/M01/e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20010101_0000z cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run Make sure that in M09_CN_e0004s_preindustrialCO2.exe: "restart_path = /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2_06/output/" "restart_id = e0004s_preindustrialCO2" Run ldsetup to set up the next cycle: cd /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/exec/e0004s_preindustrialCO2/Linux/bin source /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/src/g5_modules ./ldsetup setup /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09 /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_CN_e0004s_preindustrialCO2.exe /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_e0004s_preindustrialCO2.bat --runmodel --monthsperjob 10 --landmodel catchCN --submit Double-checked the executable and the restart file. They are both correct. ========== 20170109: 1. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170110: 1. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170111: 1. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. 2. Investigated why the carbon states increase in the high latitudes in the pre-industrial CO2 spinup cycles. ========== 20170112: 1. Investigated why the carbon states increase in the high latitudes in the pre-industrial CO2 spinup cycles. Updated Randy on what I found about this. ========== 20170113: 1. Met with Randy briefly to talk about why the carbon states increase in the high latitudes in the pre-industrial CO2 spinup cycles. 2. Worked on compute_rc.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170115: 1. The e0004s_preindustrialCO2 finished the 7th cycle. dali11:/discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l Y2016/M01 total 7354336 -rw-r--r-- 1 fzeng g0620 7547615847 2017-01-13 15:13 e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z dali11:/discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l ../../cat/ens_avg/Y2015/M12 total 526496 -rw-r--r-- 1 fzeng g0620 539112640 2017-01-13 15:13 e0004s_preindustrialCO2.ens_avg.ldas_tile_monthly_out.201512.bin Process the output (covert from tile-space to grid-space): ~/Catchment/SMAP_M09 > tile2grid_ease_spinup_monthly e0004s_preindustrialCO2 Check the output on GrADS. Looks good. cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/ mv e0004s_preindustrialCO2 e0004s_preindustrialCO2_07 cd e0004s_preindustrialCO2_07/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/ /bin/rm -rf Y200[1-9] Y201[0-5] mkdir -p Y2001/M01 cd Y2001/M01 ln -s ../../Y2016/M01/e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20010101_0000z cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run Make sure that in M09_CN_e0004s_preindustrialCO2.exe: "restart_path = /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2_07/output/" "restart_id = e0004s_preindustrialCO2" Run ldsetup to set up the next cycle: cd /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/exec/e0004s_preindustrialCO2/Linux/bin source /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/src/g5_modules ./ldsetup setup /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09 /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_CN_e0004s_preindustrialCO2.exe /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_e0004s_preindustrialCO2.bat --runmodel --monthsperjob 10 --landmodel catchCN --submit Double-checked the executable and the restart file. They are both correct. ========== 20170117: 1. Worked on catch_types.F90 and clsm_ensdrv_init_routines.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. 2. Updated higgins rain data through Jan 16, 2017. 3. Read the code to see how to connect rc to fpar. This is for future use of MODIS FPAR in Catchment-CN to calculate rc. Checked if parsun and parsha are used in anywhere else. Subroutine stomata: ppf = 4.6 * apar(n,nv) ! gkw: used in 8.3 j = ppf * qe25(ivt) ! gkw: used in 8.3 wj = max(ci-cp,0.)*j/(ci+2.*cp)*c3psn(ivt) + j*(1.-c3psn(ivt)) ! gkw: 8.3 psn(n,nv) = min(wj,wc,we) xn=1.-je/j ! gkw: 0.8 factor removed 20141108 (email from Jung-Eun) xn=max(xn,0.) ! gkw: added 1/31/14 if (wj <= 0.) xn=0. call fluorescence(xn,fs) sif(n,nv) = fs*ppf APAR affects psn and sif other than rs. It's only used in computer_rc.F90 and the subroutines therein, not used anywhere else. Modified compute_rc.F90 (see /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp/compute_rc.F90.MODISfpar) to use MODIS FPAR. ========== 20170118: 1. Worked on process_cn.F90 for upgrading to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170119: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170120: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170123: 1. The e0004s_preindustrialCO2 finished the 8th cycle. dali17:/discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l Y2016/M01 total 7354336 -rw-r--r-- 1 fzeng g0620 7547615847 2017-01-22 15:54 e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z dali17:/discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l ../../cat/ens_avg/Y2015/M12 total 526496 -rw-r--r-- 1 fzeng g0620 539112640 2017-01-22 15:54 e0004s_preindustrialCO2.ens_avg.ldas_tile_monthly_out.201512.bin Process the output (covert from tile-space to grid-space): ~/Catchment/SMAP_M09 > tile2grid_ease_spinup_monthly e0004s_preindustrialCO2 Check the output on GrADS. Looks good. cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/ mv e0004s_preindustrialCO2 e0004s_preindustrialCO2_08 cd e0004s_preindustrialCO2_08/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/ /bin/rm -rf Y200[1-9] Y201[0-5] mkdir -p Y2001/M01 cd Y2001/M01 ln -s ../../Y2016/M01/e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20010101_0000z cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run Make sure that in M09_CN_e0004s_preindustrialCO2.exe: "restart_path = /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2_08/output/" "restart_id = e0004s_preindustrialCO2" Run ldsetup to set up the next cycle: cd /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/exec/e0004s_preindustrialCO2/Linux/bin source /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/src/g5_modules ./ldsetup setup /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09 /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_CN_e0004s_preindustrialCO2.exe /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_e0004s_preindustrialCO2.bat --runmodel --monthsperjob 10 --landmodel catchCN --submit Double-checked the executable and the restart file. They are both correct. 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170124: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170125: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170126: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170130: 1. Updated higgins rain data through Jan 29, 2017. 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170131: 1. The e0004s_preindustrialCO2 finished the 9th cycle. dali10:/discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l Y2016/M01 total 7354336 -rw-r--r-- 1 fzeng g0620 7547615847 2017-01-30 14:31 e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z dali10:/discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000 > ls -l ../../cat/ens_avg/Y2015/M12 total 526496 -rw-r--r-- 1 fzeng g0620 539112640 2017-01-30 14:32 e0004s_preindustrialCO2.ens_avg.ldas_tile_monthly_out.201512.bin Process the output (covert from tile-space to grid-space): ~/Catchment/SMAP_M09 > tile2grid_ease_spinup_monthly e0004s_preindustrialCO2 Check the output on GrADS. Looks good. cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/ mv e0004s_preindustrialCO2 e0004s_preindustrialCO2_09 cd e0004s_preindustrialCO2_09/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/ /bin/rm -rf Y200[1-9] Y201[0-5] mkdir -p Y2001/M01 cd Y2001/M01 ln -s ../../Y2016/M01/e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20160101_0000z e0004s_preindustrialCO2.ens0000.catchcn_ldas_rst.20010101_0000z cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run Make sure that in M09_CN_e0004s_preindustrialCO2.exe: "restart_path = /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_preindustrialCO2_09/output/" "restart_id = e0004s_preindustrialCO2" Run ldsetup to set up the next cycle: cd /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/exec/e0004s_preindustrialCO2/Linux/bin source /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/src/g5_modules ./ldsetup setup /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09 /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_CN_e0004s_preindustrialCO2.exe /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_e0004s_preindustrialCO2.bat --runmodel --monthsperjob 10 --landmodel catchCN --submit Double-checked the executable and the restart file. They are both correct. 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily.