========== 20191101: 1. Prepared meeting slides. 2. Met with Randy, Melanie and Eunjee. 3. Processed and analyzed the output: > cd ~/Catchment/CLM4.5 > grid_restore_hdeg clm4.5_DE720_nofire > grid_restore_hdeg clm4.5_DE720_nofire_PFT6Vcmax65 > grads -lc "run check_output.gs" 4. October report. 5. Set up a point experiment: double precipitation > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSlana_GridComp > nedit clsm_ensdrv_force_routines.F90 & Added to L375-378: ! double precipitation, fzeng, 1 Nov 2019 met_force_obs_tile_new%Rainf = met_force_obs_tile_new%Rainf * 2 met_force_obs_tile_new%Rainf_C = met_force_obs_tile_new%Rainf_C * 2 met_force_obs_tile_new%Snowf = met_force_obs_tile_new%Snowf * 2 > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules In case I didn't recompile after I restored compute_rc.F90 and CN_DriverMod.F90 last time: > cd ../GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp/ > gmake install > cd ../../../GEOSlana_GridComp > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_DoublePRCP > /bin/cp -pr Linux exec/clm4.5_Kn11_DoublePRCP/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_DoublePRCP/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DoublePRCP_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71267990 2019-11-01 15:59 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-01 16:01 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_200730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSlana_GridComp > cp clsm_ensdrv_force_routines.F90.restore clsm_ensdrv_force_routines.F90 > gmake install ========== 20191104: 1. Processed GMAO Nov forecast. 2. Processed and analyzed the output of clm4.5_DoublePRCP_point: > cd ~/Catchment/CLM4.5 > grads -lc "run check_output.gs" ========== 20191111: 1. Continued processing and analyzing the output of clm4.5_DoublePRCP_point: > cd ~/Catchment/CLM4.5 > grads -lc "run check_output.gs" 2. Read CLM4.5 Tech Note about photosynthesis. ========== 20191112: 1. Fire experiment: change PFT 6 Vcmax25top to 65.18 change the impact of surface soil moisture on fire occurrence: change exponent of f_theta from 2 to 1.5 use the original thresholds for the impact of btran on fire spread > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > nedit CN_DriverMod.F90 & Changed FLNR from 0.0716 to 0.1138 for type 6. > nedit CNFireMod.F90 & ! fire_m = exp(-SHR_CONST_PI *(m/0.69_r8)**2)*(1.0_r8 - max(0._r8, & fire_m = exp(-SHR_CONST_PI *(m/0.69_r8)**1.5)*(1.0_r8 - max(0._r8, & min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8))))* & min(1._r8,exp(SHR_CONST_PI*(forc_t(g)-SHR_CONST_TKFRZ)/10._r8)) Changed the btran lower threshold from 0.4 back to 0.3, and the btran upper threshold from 0.8 back to 0.7. spread_m = (1.0_r8 - max(0._r8,min(1._r8,(btran_col(c)/wtlf(c)-0.3_r8)/ & (0.7_r8-0.3_r8))))*(1.0-max(0._r8, & min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8)))) ! spread_m = (1.0_r8 - max(0._r8,min(1._r8,(btran_col(c)/wtlf(c)-0.4_r8)/ & ! (0.8_r8-0.4_r8))))*(1.0-max(0._r8, & ! min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8)))) > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_PFT6Vcmax65_origBtran_sfmc1 > /bin/cp -pr Linux exec/clm4.5_Kn11_PFT6Vcmax65_origBtran_sfmc1/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_PFT6Vcmax65_origBtran_sfmc1/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_PFT6Vcmax65_origBtran_sfmc1_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71263362 2019-11-12 12:17 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-12 12:27 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_000730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.restore CN_DriverMod.F90 > cp CNFireMod.F90.restore CNFireMod.F90 > gmake install 2. Fire experiment: change PFT 6 Vcmax25top to 65.18 change the impact of surface soil moisture on fire occurrence: change exponent of f_theta from 2 to 1.5, and change pi to 5 use the original thresholds for the impact of btran on fire spread > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > nedit CN_DriverMod.F90 & Changed FLNR from 0.0716 to 0.1138 for type 6. > nedit CNFireMod.F90 & ! fire_m = exp(-SHR_CONST_PI *(m/0.69_r8)**2)*(1.0_r8 - max(0._r8, & fire_m = exp(-5.0 *(m/0.69_r8)**1.5)*(1.0_r8 - max(0._r8, & min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8))))* & min(1._r8,exp(SHR_CONST_PI*(forc_t(g)-SHR_CONST_TKFRZ)/10._r8)) Changed the btran lower threshold from 0.4 back to 0.3, and the btran upper threshold from 0.8 back to 0.7. spread_m = (1.0_r8 - max(0._r8,min(1._r8,(btran_col(c)/wtlf(c)-0.3_r8)/ & (0.7_r8-0.3_r8))))*(1.0-max(0._r8, & min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8)))) ! spread_m = (1.0_r8 - max(0._r8,min(1._r8,(btran_col(c)/wtlf(c)-0.4_r8)/ & ! (0.8_r8-0.4_r8))))*(1.0-max(0._r8, & ! min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8)))) > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_PFT6Vcmax65_origBtran_sfmc2 > /bin/cp -pr Linux exec/clm4.5_Kn11_PFT6Vcmax65_origBtran_sfmc2/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_PFT6Vcmax65_origBtran_sfmc2/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_PFT6Vcmax65_origBtran_sfmc2_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71263362 2019-11-12 14:59 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-12 15:03 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_010730z. > sbatch lenkf.0.j 3. Fire experiment: change PFT 6 Vcmax25top to 65.18 change the impact of surface soil moisture on fire occurrence: change exponent of f_theta from 2 to 1.5, and change pi to 5 change the thresholds for the impact of btran on fire spread: lower threshold from 0.3 to 0.2, upper threshold from 0.7 to 0.6 > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > nedit CN_DriverMod.F90 & Changed FLNR from 0.0716 to 0.1138 for type 6. > nedit CNFireMod.F90 & ! fire_m = exp(-SHR_CONST_PI *(m/0.69_r8)**2)*(1.0_r8 - max(0._r8, & fire_m = exp(-5.0 *(m/0.69_r8)**1.5)*(1.0_r8 - max(0._r8, & min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8))))* & min(1._r8,exp(SHR_CONST_PI*(forc_t(g)-SHR_CONST_TKFRZ)/10._r8)) Changed the btran lower threshold from 0.4 back to 0.3, and the btran upper threshold from 0.8 back to 0.7. ! spread_m = (1.0_r8 - max(0._r8,min(1._r8,(btran_col(c)/wtlf(c)-0.3_r8)/ & ! (0.7_r8-0.3_r8))))*(1.0-max(0._r8, & ! min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8)))) ! spread_m = (1.0_r8 - max(0._r8,min(1._r8,(btran_col(c)/wtlf(c)-0.4_r8)/ & ! (0.8_r8-0.4_r8))))*(1.0-max(0._r8, & ! min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8)))) spread_m = (1.0_r8 - max(0._r8,min(1._r8,(btran_col(c)/wtlf(c)-0.2_r8)/ & (0.6_r8-0.2_r8))))*(1.0-max(0._r8, & min(1._r8,(forc_rh(g)-30._r8)/(80._r8-30._r8)))) > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_PFT6Vcmax65_Btran1_sfmc2 > /bin/cp -pr Linux exec/clm4.5_Kn11_PFT6Vcmax65_Btran1_sfmc2/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_PFT6Vcmax65_Btran1_sfmc2/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_PFT6Vcmax65_Btran1_sfmc2_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71263362 2019-11-12 15:47 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-12 15:50 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_030730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.restore CN_DriverMod.F90 > cp CNFireMod.F90.restore CNFireMod.F90 > gmake install ========== 20191113: 1. Processed and analyzed the output of clm4.5_PFT6Vcmax65_origBtran_sfmc1_point clm4.5_PFT6Vcmax65_origBtran_sfmc2_point clm4.5_PFT6Vcmax65_Btran1_sfmc2_point 2. Searched for and extracted CASA GFED4s NEP data for Huisheng. ========== 20191114: 1. Prepared meeting slides. 2. Met with Randy, Melanie and Eunjee. ========== 20191115: 1. Talked with Sarith about the C180 run for S2S forecast. 2. Further analyzed the output of clm4.5_PFT6Vcmax65_Btran1_sfmc2_point and the control run. Tile file: 100 125745 21.2500 4.7500 403 190 1.000000000000 125745 0.1799E+03 266364100000 15434 ========== 20191118: 1. Read literature trying to find a way to increase GPP while increasing burn and fire C loss. 2. GPP and fire experiment: Increase NDEP to a huge number so that there is no nitrogen limitation. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > nedit CN_DriverMod.F90 & L768-769: ! forc_ndep(:) = ndep(:) forc_ndep(:) = 1000./86400. * 10 ! gN/m2/s, so this is 10 kgN/m2/day, unlimited N supply, fzeng, 18 Nov 2019 > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_unlimitedN > /bin/cp -pr Linux exec/clm4.5_Kn11_unlimitedN/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_unlimitedN/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_unlimitedN_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71263362 2019-11-18 19:59 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-18 20:03 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_020730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.restore CN_DriverMod.F90 > gmake install ========== 20191119: 1. Analyzed the output of clm4.5_unlimitedN_point. 2. GPP and fire experiment: Modified model code so that stomatal conductance is also down-regulated by N limitation like what is done for GPP. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > nedit CN_DriverMod.F90 & > cp -p CN_DriverMod.F90 CN_DriverMod.F90.downreg > cd ../../../GEOSlana_GridComp/ > nedit process_cn.F90 & > cp -p process_cn.F90 process_cn.F90.downreg > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_downreg > /bin/cp -pr Linux exec/clm4.5_Kn11_downreg/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_downreg/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_downreg_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71267458 2019-11-19 20:08 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-19 20:11 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j Crashed: date_time_new 19800101_013000z forrtl: error (65): floating invalid Image PC Routine Line Source libirc.so 00002AAAAB590961 Unknown Unknown Unknown libirc.so 00002AAAAB58F0B7 Unknown Unknown Unknown libmpi_usempif08. 00002AAAAC2E3362 Unknown Unknown Unknown libmpi_usempif08. 00002AAAAC2E31B6 Unknown Unknown Unknown libmpi_usempif08. 00002AAAAC2D1A0C Unknown Unknown Unknown libmpi_usempif08. 00002AAAAC2B2E82 Unknown Unknown Unknown libc.so.6 00002AAAABA16910 Unknown Unknown Unknown LDASsaCN_mpi.x 00000000007AF9AC catchment_cn_mode 2124 catchmentCN.F90 LDASsaCN_mpi.x 00000000007AE738 catchment_cn_mode 1777 catchmentCN.F90 LDASsaCN_mpi.x 000000000079E740 catchment_cn_mode 727 catchmentCN.F90 LDASsaCN_mpi.x 0000000000655B44 process_cn_mp_pro 1901 process_cn.F90 LDASsaCN_mpi.x 00000000005C16C6 ens_driver_routin 375 ens_driver_routines.F90 LDASsaCN_mpi.x 00000000004AF170 MAIN__ 2101 cnlsm_ensdrv_main.F90 LDASsaCN_mpi.x 000000000049E6BE Unknown Unknown Unknown libc.so.6 00002AAAABA02C36 Unknown Unknown Unknown ========== 20191120: 1. Compiled the code with BOPT=g GPP and fire experiment: Modified model code so that stomatal conductance is also down-regulated by N limitation like what is done for GPP. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake clean > gmake install BOPT=g > cd ../../../GEOSlana_GridComp > gmake clean > gmake install BOPT=g > cd ../../Applications/LDAS_App/ > gmake clean > gmake install BOPT=g > cd ../../.. > /bin/cp -pr Linux exec/clm4.5_Kn11_downreg/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_downreg_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 99307805 2019-11-20 09:53 build/Linux/bin/LDASsaCN_mpi.x* > cd run > cp lenkf.0.j lenkf.0.i > nedit lenkf.0.i & > /bin/rm ../output/point/rc_out/* > ./lenkf.0.i Fixed the bug. It was able to run on Arm DDT. Stopped it at date_time_new 19800102_110730z. Re-compiled without BOPT=g. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > gmake clean > gmake install > cd ../../../GEOSlana_GridComp > gmake clean > gmake install BOPT=g > cd ../../Applications/LDAS_App/ > gmake clean > gmake install BOPT=g > cd ../../.. > /bin/cp -pr Linux exec/clm4.5_Kn11_downreg/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_downreg_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71267598 2019-11-20 11:13 build/Linux/bin/LDASsaCN_mpi.x* > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_010730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.restore CN_DriverMod.F90 > gmake install > cd ../../../GEOSlana_GridComp > cp process_cn.F90.restore process_cn.F90 > gmake install 2. GPP and fire experiment: Modified model code so that stomatal conductance is also down-regulated by N limitation like what is done for GPP. Increase PFT6 FLNR from 0.0716 to 0.096 so that Vcmax increase from 41 to 55. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.downreg CN_DriverMod.F90 > nedit CN_DriverMod.F90 & Change flnr of PFT6 from 0.0716 to 0.096 > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > cp process_cn.F90.downreg process_cn.F90 > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_downreg_PFT6Vcmax55 > /bin/cp -pr Linux exec/clm4.5_Kn11_downreg_PFT6Vcmax55/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_downreg_PFT6Vcmax55/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_downreg_PFT6Vcmax55_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71267598 2019-11-20 14:37 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-20 14:42 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_020730z. > sbatch lenkf.0.j 3. GPP and fire experiment: Modified model code so that stomatal conductance is also down-regulated by N limitation like what is done for GPP. Increase PFT6 FLNR from 0.0716 to 0.1138 so that Vcmax increase from 41 to 65.18. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.downreg CN_DriverMod.F90 > nedit CN_DriverMod.F90 & Change flnr of PFT6 from 0.0716 to 0.1138 > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > cp process_cn.F90.downreg process_cn.F90 > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_downreg_PFT6Vcmax65 > /bin/cp -pr Linux exec/clm4.5_Kn11_downreg_PFT6Vcmax65/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_downreg_PFT6Vcmax65/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_downreg_PFT6Vcmax65_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71267598 2019-11-20 14:57 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-20 15:01 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_040730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.restore CN_DriverMod.F90 > gmake install > cd ../../../GEOSlana_GridComp > cp process_cn.F90.restore process_cn.F90 > gmake install 4. Met with Huisheng for ~2.5 hours to take a look at the 3 NEE/NEP datasets. Finally we decided to you the one from Randy Kawa. 5. Plot our CatchCN4.5 monthly LAI climatology and compare it to the MODIS LAI monthly climatology plots Sarith made in the bcs directory. ========== 20191121: 1. Analyzed the output of clm4.5_downreg_point clm4.5_downreg_PFT6Vcmax55_point clm4.5_downreg_PFT6Vcmax65_point 2. GPP and fire experiment: Modified model code so that stomatal conductance is also down-regulated by N limitation like what is done for GPP. Increase PFT6 FLNR from 0.0716 to 0.0786 so that Vcmax increase from 41 to 45. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.downreg CN_DriverMod.F90 > nedit CN_DriverMod.F90 & Change flnr of PFT6 from 0.0716 to 0.0786 > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > cp process_cn.F90.downreg process_cn.F90 > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_downreg_PFT6Vcmax45 > /bin/cp -pr Linux exec/clm4.5_Kn11_downreg_PFT6Vcmax45/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_downreg_PFT6Vcmax45/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_downreg_PFT6Vcmax45_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71267598 2019-11-21 10:48 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-21 10:51 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_010730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.restore CN_DriverMod.F90 > gmake install > cd ../../../GEOSlana_GridComp > cp process_cn.F90.restore process_cn.F90 > gmake install 3. Prepared meeting slides. 4. Fire meeting. 5. NDEP experiment: Use NDEP = 36.8463 ng/m2/s for the SH Africa point (it is originally 15.6155). > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > nedit CN_DriverMod.F90 & ! forc_ndep(:) = ndep(:) forc_ndep(:) = ndep(:) * 36.8463 / 15.6155 > setenv ESMADIR /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05 > source $ESMADIR/src/g5_modules > gmake install > cd ../../../GEOSlana_GridComp > gmake install > cd ../../Applications/LDAS_App/ > gmake install > cd ../../.. > mkdir -p exec/clm4.5_Kn11_ndep1 > /bin/cp -pr Linux exec/clm4.5_Kn11_ndep1/. > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run > cp -p ../clm4.5_point/run/driver_inputs_0.5D_point.nml . > nedit D0.5_clm4.5_point.bat D0.5_clm4.5_point.exe & > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/exec/clm4.5_Kn11_ndep1/Linux/bin > source /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/g5_modules > ./ldsetup setup /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5 /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.exe /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/run/D0.5_clm4.5_point.bat --runmodel --monthsperjob 120 --landmodel catchCN > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_ndep1_point > ls -l build/Linux/bin/LDASsaCN_mpi.x -rwxr-xr-x 1 fzeng g0620 71325654 2019-11-21 16:04 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ lrwxrwxrwx 1 fzeng g0620 80 2019-11-21 16:07 output -> /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/tests/clm4.5/clm4.5_DE720_66/output/ > cd run > ./lenkf.0.j It's running. Stopped it at date_time_new 19800102_010730z. > sbatch lenkf.0.j Restore and compile: > cd /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/clm4.5/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3_V24_C05/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp > cp CN_DriverMod.F90.restore CN_DriverMod.F90 > gmake install ========== 20191122: 1. Analyzed the output of clm4.5_Kn11_ndep1 and emailed Randy, Melanie and Eunjee the plots. By increasing ndep, GPP, LAI, burned area, fire C loss and carbon pools all increase. ========== 20191129-20191202: 1. For IDS hindcast project: 80 simulations total. (1) Use 20160101 initial conditions and bias corrected hindcast meteorology from each of the 4 ensembles in each of the 10 years (40 simulations total) (2) Use the 4 ensembles of the 2016 bias corrected hindcast meteorology data and the initial conditions from each of the 10 years (40 simulations total). (1) Use 20160101 initial conditions and bias corrected hindcast meteorology from each of the 4 ensembles in each of the 10 years (40 simulations total) Take e0004s_transientCO2_20160101init_BChindcast_YYYYJan01_ens1 as an example. YYYY is each of the 10 years from 2007 to 2016. Check if the restart file is available: (one-time thing) > ls -l /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_05/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/Y2016/M01/e0004s_transientCO2.ens0000.catchcn_ldas_rst.20160101_0000z (1.1) Prepare initial conditions. (one-time thing for each of the 9 years -- not needed for 2016) > cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_05/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/Y2007/M01 > mv e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z.orig > ln -s ../../Y2016/M01/e0004s_transientCO2.ens0000.catchcn_ldas_rst.20160101_0000z e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z (1.2) Prepare files and run directory for experiment setup: > cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run > nedit M09_CN_e0004s_transientCO2_BChindcast.exe M09_e0004s_transientCO2_BChindcast.bat & M09_CN_e0004s_transientCO2_BChindcast.exe: exp_id = e0004s_transientCO2_20160101init_BChindcast_2007Jan01_ens1 exp_domain = SMAP_EASEv2_M09_GLOBAL N_ens = 1 start_time = 2007-01-01-00-00-00 end_time = 2007-10-01-00-00-00 force_dtstep = 21600 restart_path = /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_05/output/ restart_domain = SMAP_EASEv2_M09_GLOBAL restart_id = e0004s_transientCO2 met_tag = HindcastBC_netcdf met_path = /discover/nobackup/projects/gmao/geos_carb/fzeng/FORECASTS/GEOS5/CLIM/GEOS5v2/BCSD_Final/6-hourly/ M09_e0004s_transientCO2_BChindcast.bat: job-name = BChindcast_2007_ens1 (1.3) Run ldsetup: > cd /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/exec/e0004s_transientCO2_BChindcast/Linux/bin > source /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/src/g5_modules > ./ldsetup setup /discover/nobackup/projects/gmao/geos_carb/fzeng/Catchment/SMAP_EASEv2_M09 /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_CN_e0004s_transientCO2_BChindcast.exe /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_e0004s_transientCO2_BChindcast.bat --runmodel --monthsperjob 10 --landmodel catchCN (1.4) Check the executable, restart file and create year_co2.txt: > cd /discover/nobackup/projects/gmao/geos_carb/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_20160101init_BChindcast_2007Jan01_ens1 > ls -l build/Linux/bin/LDASsaCN_mpi.x (to make sure the executable is the right one) -rwxr-xr-x 1 fzeng g0620 69619628 2019-06-26 11:25 build/Linux/bin/LDASsaCN_mpi.x* > ls -l input/restart/ (to make sure the restart file is the right one) lrwxrwxrwx 1 fzeng s1460 80 2019-11-29 10:50 output -> /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_05/output/ (1.5) Create year_co2.txt: > cd run > echo 2007 > year_co2.txt > cat year_co2.txt (to make sure the year of CO2 is correct) 2007 [Correct!] (1.6) Modify the restart_id in the first job script: > nedit lenkf.0.j & [No change made!] -restart_path ../input/restart/output \ -restart_domain SMAP_EASEv2_M09_GLOBAL \ -restart_id e0004s_transientCO2 \ (1.7) Create LDAS.rc: > nedit LDAS.rc & INIT_MON: jan01 ENS_NUM: 1 NOTE: The space is important, so just copy and paste from /discover/nobackup/projects/gmao/geos_carb/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_BChindcast_2016Jan01_ens1/run/LDAS.rc (1.8) I moved the noaaCO2 data to a different directory a few days ago. Created a symbolic link so the existing executable still works. (one-time thing) /discover/nobackup/fzeng > ln -s data/noaaCO2 . (1.9) Did an interactive run to make sure it works: > interactive.py -A sp3 -n 140 -a g0620 -X --debug > cd /discover/nobackup/projects/gmao/geos_carb/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_20160101init_BChindcast_2007Jan01_ens1/run > ./lenkf.0.j It's running, and using the right restart file and meteorology data. Reading restart file ../input/restart/output/SMAP_EASEv2_M09_GLOBAL//rs/ens0000//Y2007/M01/e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z get_forcing(): bias-corrected GEOS-5 hindcast forcing data set opening../input/met_forcing/6-hourly//2007/jan01/ens1/GEOS5.200701.nc4 Stopped it at date_time_new 20070101_133000z. > ls -l ../input/restart/output/SMAP_EASEv2_M09_GLOBAL//rs/ens0000//Y2007/M01/ (one-time thing for each of the 10 years) total 7354336 lrwxrwxrwx 1 fzeng g0620 75 2019-11-29 07:40 e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z -> ../../Y2016/M01/e0004s_transientCO2.ens0000.catchcn_ldas_rst.20160101_0000z -rw-r--r-- 1 fzeng g0620 7547615847 2017-07-28 14:44 e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z.orig (1.10) Submit the job: > qsub lenkf.0.j (1.11) Process the four ensemble run output and check the output on GrADS. > cd ~/Catchment/SMAP_M09 > tile2grid_ease_BChindcast_init_met_monthly 2016 2007 1 > tile2grid_ease_BChindcast_init_met_monthly 2016 2007 2 > tile2grid_ease_BChindcast_init_met_monthly 2016 2007 3 > tile2grid_ease_BChindcast_init_met_monthly 2016 2007 4 Checked July GPP and NEE of the 4 ensembles using check_e0004s_transientCO2_BChindcast_monthly.gs on GrADS. Looks correct. (1.12) Restore the restart file when all the 2007 simulations are completed!! (not needed for 2016) > cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_05/output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/Y2007/M01 > mv e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z.20160101 > mv e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z.orig e0004s_transientCO2.ens0000.catchcn_ldas_rst.20070101_0000z