========== 20170601: 1. Continued reading Eunjee's manuscript and provided comments. 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170602: 1. May report 2. Made time series plots of 3hrly, daily and monthly CT CO2 for Eunjee. (1) [NH] Boreal, North America : CA-Qfo Site name: QC-Eastern Old Black Spruce (EOBS) https://fluxnet.ornl.gov/site/269 Lat: 49.69247, Lon: -74.34204 (2) [NH] Boreal, Europe: RU-Fyo Site name: Fyodorovskoye http://sites.fluxdata.org/RU-Fyo/ Lat: 56.4615, Lon: 32.9221 (3) [SH] Tropical, South America: TNF Site name: Tapajos National Forest Lat: -2.51, Lon: -54.58 (4) [SH] Temperate, Australia: AU-Tum Site name: Tumbarumba http://sites.fluxdata.org/AU-Tum/ Lat: -35.6566, Lon: 148.1517 3. Processed GMAO GEOS5 seasonal forecast for June. ========== 20170605: 1. Updated higgins rain data through Jun 4. 2. Sarith created new bcs: /discover/nobackup/smahanam/bcs/Icarus-NL/Icarus-NL_Reynolds More GEOS-5 test runs: Catchment-CN (1) /discover07:/discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm/src/Applications/GEOSgcm_App > gcm_setup Enter the Experiment ID: I10p1C90Reynolds Enter a 1-line Experiment Description: I10p1C90Reynolds Do you wish to CLONE an old experiment? (Default: NO or FALSE) Enter the Atmospheric Horizontal Resolution code: ----------------------------------------------------------- Lat/Lon Cubed-Sphere ----------------------------------------------------------- b -- 2 deg c48 -- 2 deg c -- 1 deg c90 -- 1 deg d -- 1/2 deg c180 -- 1/2 deg (56-km) e -- 1/4 deg (35-km) c360 -- 1/4 deg (28-km) c720 -- 1/8 deg (14-km) c1440 - 1/16 deg ( 7-km) c90 Enter the Atmospheric Model Vertical Resolution: LM (Default: 72) Do you wish to run the COUPLED Ocean/Sea-Ice Model? (Default: NO or FALSE) Enter the Data_Ocean Horizontal Resolution code: o1 (1 -deg, 360x180 Reynolds) Default o2 (1/4-deg, 1440x720 MERRA-2) o3 (1/8-deg, 2880x1440 OSTIA) CS (Cubed-Sphere OSTIA) o1 Enter the choice of Land Surface Model: use 1 (Default: Catchment Model), 2 (CatchmentCN Model) 2 Do you wish to run the Runoff Routing Model? (Default: NO or FALSE) Do you wish to run GOCART with Actual or Climatological Aerosols? (Enter: A (Default) or C) Enter the GOCART Emission Files to use: MERRA2 (Default), PIESA, CMIP, NR, MERRA2-DD or OPS: Enter the tag or directory (/filename) of the HISTORY.AGCM.rc.tmpl to use (To use HISTORY.AGCM.rc.tmpl from current build, Type: Current ) ------------------------------------------------------------------------- Hit ENTER to use Default Tag/Location: (Current) Enter Desired Location for the HOME Directory (to contain scripts and RC files) Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds Enter Desired Location for the EXPERIMENT Directory (to contain model output and restart files) Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds Enter Location for Build directory containing: src/ Linux/ etc... Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.catchcn_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/catchcn_internal_rst /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.fvcore_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/fvcore_internal_rst /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.gocart_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/gocart_internal_rst /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.lake_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/lake_internal_rst /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.landice_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/landice_internal_rst /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.moist_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/moist_internal_rst /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.pchem_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/pchem_internal_rst /bin/cp -p /discover/nobackup/fzeng/geos5_sims/cf90_catchcn_INL/M2_rst/OUTPUT/C90c_d5124_m2_jan10.saltwater_internal_rst.20121231_21z.bin /discover/nobackup/fzeng/geos5_sims/I10p1C90Reynolds/saltwater_internal_rst Create cap_restart 20030101 210000 Edit CAP.rc END_DATE: 20030103 210000 Edit HISTORY.rc In COLLECTIONS, comment out everything other than "geosgcm_surf" [L57] geosgcm_surf.frequency: 240000 [Lines 206-242] Uncomment geosgcm_surf collection for carbon variables (e.g., CNLAI, ICESOI) Edit gcm_run.j [Line 152] setenv BCSDIR /discover/nobackup/ltakacs/bcs/Icarus/Icarus_Reynolds => setenv BCSDIR /discover/nobackup/smahanam/bcs/Icarus-NL/Icarus-NL_Reynolds Complained about saltwater_internal_rst: numbers of tiles don't match. (2) (2.1) Set up experiment directory: /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm/src/Applications/GEOSgcm_App > ./gcm_setup --nocvs Enter the Experiment ID: I10p1_CatchCN Enter a 1-line Experiment Description: I10p1_CatchCN Do you wish to CLONE an old experiment? (Default: NO or FALSE) Enter the Atmospheric Horizontal Resolution code: ----------------------------------------------------------- Lat/Lon Cubed-Sphere ----------------------------------------------------------- b -- 2 deg c48 -- 2 deg c -- 1 deg c90 -- 1 deg d -- 1/2 deg c180 -- 1/2 deg (56-km) e -- 1/4 deg (35-km) c360 -- 1/4 deg (28-km) c720 -- 1/8 deg (14-km) c1440 - 1/16 deg ( 7-km) c90 Enter the Atmospheric Model Vertical Resolution: LM (Default: 72) Do you wish to run the COUPLED Ocean/Sea-Ice Model? (Default: NO or FALSE) Enter the Data_Ocean Horizontal Resolution code: o1 (1 -deg, 360x180 Reynolds) Default o2 (1/4-deg, 1440x720 MERRA-2) o3 (1/8-deg, 2880x1440 OSTIA) CS (Cubed-Sphere OSTIA) o1 Enter the choice of Land Surface Model: use 1 (Default: Catchment Model), 2 (CatchmentCN Model) 2 Do you wish to run the Runoff Routing Model? (Default: NO or FALSE) Do you wish to run GOCART with Actual or Climatological Aerosols? (Enter: A (Default) or C) Enter the GOCART Emission Files to use: MERRA2 (Default), PIESA, CMIP, NR, MERRA2-DD or OPS: Enter the tag or directory (/filename) of the HISTORY.AGCM.rc.tmpl to use (To use HISTORY.AGCM.rc.tmpl from current build, Type: Current ) ------------------------------------------------------------------------- Hit ENTER to use Default Tag/Location: (Current) Enter Desired Location for the HOME Directory (to contain scripts and RC files) Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN Enter Desired Location for the EXPERIMENT Directory (to contain model output and restart files) Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN Enter Location for Build directory containing: src/ Linux/ etc... Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm (2.2) Create restart files: /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN > mkdir M2_rst /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN/M2_rst > ln -s /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm/Linux/bin /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN/M2_rst > bin/regrid.pl -catchcn g5_modules_perl_wrapper: Setting BASEDIR and modules for discover18 Do you want script to select MERRA-2 input restarts for you (y/n)? [n] y Enter OUTPUT directory for new restarts: OUTPUT mkdir OUTPUT Enter ymd (yyyymmdd): 20021231 Enter hour (hh): 21 ---------------- Atmosphere Grids ---------------- Lat/Lon Grids Cubed-Sphere Grids ------------- ------------------ a = 4 deg C12 C180 C1000 b = 2 deg C24 C360 C1440 c = 1 deg C48 C500 C2880 d = 1/2 deg C90 C720 e = 1/4 deg FVCORE: /archive/users/gmao_ops/MERRA2/gmao_ops/GEOSadas-5_12_4/d5124_m2_jan00/rs/Y2002/M12/d5124_m2_jan00.fvcore_internal_rst.20021231_21z.bin Getting atmosphere grid resolution from INPUT fvcore file ... DONE INPUT atmosphere grid: C180 Enter OUTPUT atmosphere grid: C90 Atmosphere levels ----------------- 72 132 137 144 INPUT atmosphere levels: 72 Enter OUTPUT atmosphere levels: [72] Ocean Grids ----------- c = 360x180 (Reynolds) e = 1440x720 (MERRA-2) f = 2880x1440 (OSTIA) CS = same as atmosphere (OSTIA cubed-sphere) INPUT ocean grid: e Enter OUTPUT ocean grid: [c] OUTPUT ocean grid: c Sample GCM tags --------------- F14 : Fortuna-1_4 ......... Fortuna-1_4_p1 F20 : Fortuna-2_0 ......... Fortuna-2_0 F21 : Fortuna-2_1 ......... Fortuna-2_5_pp2 G10 : Ganymed-1_0 ......... Ganymed-1_0_BETA4 G10p : Ganymed-1_0_p1 ...... Ganymed-1_0_p6 G20 : Ganymed-2_0 ......... Ganymed-2_1_p6 G30 : Ganymed-3_0 ......... Ganymed-3_0_p1 G40 : Ganymed-4_0 ......... Heracles-5_4_p3 HNL : Icarus-NL ........... Icarus-NL ICA : Icarus .............. Icarus Sample DAS tags --------------- 214 : GEOSdas-2_1_4 ....... GEOSdas-2_1_4-m4 540 : GEOSadas-5_4_0 ...... GEOSadas-5_5_3 561 : GEOSadas-5_6_1 ...... GEOSadas-5_7_3_p2 580 : GEOSadas-5_8_0 ...... GEOSadas-5_9_1 591p : GEOSadas-5_9_1_p1 ... GEOSadas-5_9_1_p9 5A0 : GEOSadas-5_10_0 ..... GEOSadas-5_10_0_p1 5B0 : GEOSadas-5_10_0_p2 .. GEOSadas-5_11_0 512 : GEOSadas-5_12_2 ..... GEOSadas-5_16_5 INPUT tag: Ganymed-4_0_MERRA-2 Type 'bcs' to see BCS tags or Enter GCM or DAS tag for outputs: [Ganymed-4_0] Icarus-NL =============== REGRID RESTARTS =============== ESMABIN: /gpfsm/dnb31/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm/Linux/bin ESMATAG: Icarus-1_0_p1 logfile: /gpfsm/dnb31/fzeng/geos5_sims/I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.20021231_21z.log upper-air: yes surface: yes - route: no - catch: yes - catchCN: yes bkg/satb: yes rst.lcv: yes namelabel: no # inputs #------- . CVS tag: Ganymed-4_0 . expid: d5124_m2_jan00 . date: 20021231 . hour: 21 . atmos grid: 180x1080 (C180) . atmos levs: 72 . ocean grid: 1440x720 (e) . BCS tag: Ganymed-4_0_MERRA-2 . surflay: 50 . rstdir: /archive/users/gmao_ops/MERRA2/gmao_ops/GEOSadas-5_12_4/d5124_m2_jan00/rs/Y2002/M12 # outputs #-------- . CVS tag: Icarus-NL . expid: C90c_d5124_m2_jan00 . date: 20021231 . hour: 21 . atmos grid: 90x540 (C90) . atmos levs: 72 . ocean grid: 360x180 (c) . BCS tag: Icarus-NL_Reynolds . surflay: 50 . rescale: yes . outdir: OUTPUT . workdir: OUTPUT/20021231_29007 =============== (2.3) Copy restart files: /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN/M2_rst: /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.catchcn_internal_rst.20021231_21z.bin ../catchcn_internal_rst /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.fvcore_internal_rst.20021231_21z.bin ../fvcore_internal_rst /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.gocart_internal_rst.20021231_21z.bin ../gocart_internal_rst /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.lake_internal_rst.20021231_21z.bin ../lake_internal_rst /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.landice_internal_rst.20021231_21z.bin ../landice_internal_rst /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.moist_internal_rst.20021231_21z.bin ../moist_internal_rst /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.pchem_internal_rst.20021231_21z.bin ../pchem_internal_rst /bin/cp -p OUTPUT/C90c_d5124_m2_jan00.saltwater_internal_rst.20021231_21z.bin ../saltwater_internal_rst (2.4) Modify settings: Create cap_restart 20021231 210000 Edit CAP.rc END_DATE: 20030103 210000 Edit HISTORY.rc In COLLECTIONS, comment out everything other than "geosgcm_surf" [L57] geosgcm_surf.frequency: 240000 [Lines 206-242] Uncomment geosgcm_surf collection for carbon variables (e.g., CNLAI, ICESOI) Edit gcm_run.j [Line 152] setenv BCSDIR /discover/nobackup/ltakacs/bcs/Icarus/Icarus_Reynolds => setenv BCSDIR /discover/nobackup/smahanam/bcs/Icarus-NL/Icarus-NL_Reynolds (2.5) Run interactively: ssh -XY discover-sp3 interactive.py -A sp3 -n 96 -a g0620 -X --debug setenv ESMADIR /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm/ source $ESMADIR/src/g5_modules cd /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN ./gcm_run.j Successful! And output carbon fluxes look correct!! ========== 20170606: 1. More GEOS-5 test run: Catchment-CN (1) Repeat (2.3) through (2.5) above. The only difference is now cap_restart is set to: 20030101 210000 Crashed: Error! Found tiles in saltwater. Expect to find 168075 tiles. Your restarts are probably for a different ocean. So the date in cap_restart has to match the date set when running regrid.pl? No. The run also crashed when I set cap_restart to 20021231 210000, identical to what I did yesterday, giving the same error message about ocean. (2) /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm/src/Applications/GEOSgcm_App > ./gcm_setup --nocvs Enter the Experiment ID: I10p1_CatchCN2 Enter a 1-line Experiment Description: I10p1_CatchCN2 Do you wish to CLONE an old experiment? (Default: NO or FALSE) Enter the Atmospheric Horizontal Resolution code: ----------------------------------------------------------- Lat/Lon Cubed-Sphere ----------------------------------------------------------- b -- 2 deg c48 -- 2 deg c -- 1 deg c90 -- 1 deg d -- 1/2 deg c180 -- 1/2 deg (56-km) e -- 1/4 deg (35-km) c360 -- 1/4 deg (28-km) c720 -- 1/8 deg (14-km) c1440 - 1/16 deg ( 7-km) c90 Enter the Atmospheric Model Vertical Resolution: LM (Default: 72) Do you wish to run the COUPLED Ocean/Sea-Ice Model? (Default: NO or FALSE) Enter the Data_Ocean Horizontal Resolution code: o1 (1 -deg, 360x180 Reynolds) Default o2 (1/4-deg, 1440x720 MERRA-2) o3 (1/8-deg, 2880x1440 OSTIA) CS (Cubed-Sphere OSTIA) o1 Enter the choice of Land Surface Model: use 1 (Default: Catchment Model), 2 (CatchmentCN Model) 2 Do you wish to run the Runoff Routing Model? (Default: NO or FALSE) Do you wish to run GOCART with Actual or Climatological Aerosols? (Enter: A (Default) or C) Enter the GOCART Emission Files to use: MERRA2 (Default), PIESA, CMIP, NR, MERRA2-DD or OPS: Enter the tag or directory (/filename) of the HISTORY.AGCM.rc.tmpl to use (To use HISTORY.AGCM.rc.tmpl from current build, Type: Current ) ------------------------------------------------------------------------- Hit ENTER to use Default Tag/Location: (Current) Enter Desired Location for the HOME Directory (to contain scripts and RC files) Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN2 Enter Desired Location for the EXPERIMENT Directory (to contain model output and restart files) Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN2 Enter Location for Build directory containing: src/ Linux/ etc... Hit ENTER to use Default Location: ---------------------------------- Default: /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN2 > ls AGCM.rc archive/ CAP.rc* convert/ forecasts/ fvcore_layout.rc* gcm_run.j* GEOSgcm.x* HISTORY.rc plot/ post/ RC/ regress/ /discover/nobackup/fzeng/geos5_sims/I10p1_CatchCN2: /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.catchcn_internal_rst.20021231_21z.bin catchcn_internal_rst /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.fvcore_internal_rst.20021231_21z.bin fvcore_internal_rst /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.gocart_internal_rst.20021231_21z.bin gocart_internal_rst /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.lake_internal_rst.20021231_21z.bin lake_internal_rst /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.landice_internal_rst.20021231_21z.bin landice_internal_rst /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.moist_internal_rst.20021231_21z.bin moist_internal_rst /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.pchem_internal_rst.20021231_21z.bin pchem_internal_rst /bin/cp -p ../I10p1_CatchCN/M2_rst/OUTPUT/C90c_d5124_m2_jan00.saltwater_internal_rst.20021231_21z.bin saltwater_internal_rst Create cap_restart 20030101 210000 Edit CAP.rc END_DATE: 20030103 210000 Edit HISTORY.rc In COLLECTIONS, comment out everything other than "geosgcm_surf" [L57] geosgcm_surf.frequency: 240000 [Lines 206-242] Uncomment geosgcm_surf collection for carbon variables (e.g., CNLAI, ICESOI) Edit gcm_run.j [Line 152] setenv BCSDIR /discover/nobackup/ltakacs/bcs/Icarus/Icarus_Reynolds => setenv BCSDIR /discover/nobackup/smahanam/bcs/Icarus-NL/Icarus-NL_Reynolds Crashed and gave the same error message about ocean. Setting cap_restart to 20021231 210000 gave the same error message. Tried a few times but all failed. But then suddenly it worked, for both 20021231 210000 and 20030101 210000 in cap_restart. Don't know why it didn't work earlier. 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170607: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170608: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170609: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170612: 1. Moved and renamed the matlab script that plots variables in catchcn_internal_rst. ~ > mkdir matlab ~ > mv geos5/check_rst_vars.m matlab/ ~/matlab > mv check_rst_vars.m plot_catchcn_internal_rst_vars.m 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes_daily. ========== 20170613-20170615: 1. Manuscript review ========== 20170616: 1. /discover/nobackup/fzeng/geos5_code/Icarus-1_0_p1/GEOSagcm/src/GEOSgcs_GridComp/GEOSgcm_GridComp/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp: cvs upd -r 1.2 CNDecompMod.F90 cvs upd -r 1.3.4.1 GEOS_CatchCNGridComp.F90 Huisheng, I just learned from Sarith that one change to CNDecompMod.F90 that we made a while ago is not included in Icarus (this change is included in Heracles). In your Icarus-1_0_p1, under GEOSphysics_GridComp/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatchCN_GridComp, please update CNDecompMod.F90 by: cvs upd -r 1.2 CNDecompMod.F90 (Note: you may do "~ltakacs/cvstools/cvsxdiff -r 1.2 CNDecompMod.F90" to take a look at the difference first) Also, Sarith updated DTCN in GEOS_CatchCNGridComp.F90 under the same directory (this is what he said in his email this morning), so please update this file in your source code as well: cvs upd -r 1.3.4.1 GEOS_CatchCNGridComp.F90 (Note: you may do "~ltakacs/cvstools/cvsxdiff -r 1.3.4.1 GEOS_CatchCNGridComp.F90" to take a look at the difference first) After you update GEOS_CatchCNGridComp.F90, you don't need this block below in your AGCM.rc any more. The default value of DTCN will be 5400 (set in this updated GEOS_CatchCNGridComp.F90) # Time step for catchment-CN (CN) # ----------------------------------------- DTCN: 5400 # 1.5-hourly Please feel free to let me know if you have any questions. Fanwei 2. /discover/nobackup/fzeng/clm4-to-clm4.5/LDAS/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatch_GridComp/mk_restarts: cvs upd -r 1.4.2.1 mk_LDASsaRestarts.F90 NOTE: use the catchcn_internal_rst in OutData2 (not the final one) for now 3. Manuscript review ========== 20170619: 1. Manuscript review 2. Set up an offline run using 3hrly CO2 and GR=0.25*GPP (originally GR=0.30*GPP) with Eunjee, using process_cat.F90.3hrlyCO2, compute_rc.F90 and CNGRespMod.F90.GRpt25 in /discover/nobackup/fzeng/offline_code/LDASsa_m3-15_2-CN/src/Components/GEOSlana_GridComp/M2n5P_m0001CTco2. 3hrly CO2, GR=0.25GPP: cd /discover/nobackup/fzeng/Catchment/M2n5P/m0001 mkdir -p M2.n5P.3hco2.glb.GRpt25/RUN/rs/ens0000/Y2001/M01 cd M2.n5P.3hco2.glb.GRpt25 cp -p /discover/nobackup/elee15/offline/sims/M2.n5P.3htCO2.c5/CN_restart . cp -p /discover/nobackup/elee15/offline/sims/M2.n5P.3htCO2.c5/RUN/rs/ens0000/Y2001/M01/mM2.n5P.3htCO2.ens0000.catch_ldas_rst.20010101_0000z.bin RUN/rs/ens0000/Y2001/M01/mM2.n5P.3hco2.glb.GRpt25.ens0000.catch_ldas_rst.20010101_0000z.bin cp /discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb/lenkf.j . /discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb.GRpt25 > ls -l CN_restart -rw-r--r-- 1 fzeng g0620 1530820456 2016-11-08 19:34 CN_restart /discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb.GRpt25 > ls -l RUN/rs/ens0000/Y2001/M01/ total 40512 -rw-r--r-- 1 fzeng g0620 41453776 2016-11-08 19:34 mM2.n5P.3hco2.glb.GRpt25.ens0000.catch_ldas_rst.20010101_0000z.bin Processed the available output: ~/Catchment/M2n5P_m0001 > tile2grid_FV_144x91_monthly M2.n5P.3hco2.glb.GRpt25 Code change in this experiment: src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatch_GridComp/CNGRespMod.F90 Line 153: grperc = 0.3 Changed to grperc = 0.25 ========== 20170620: 1. The SMAP M09 transient CO2 run finished the 3rd cycle. CO2 years: 1909-1944; meteorology years: 1981-2016. Process the output (covert from tile-space to grid-space): ~/Catchment/SMAP_M09 > tile2grid_ease_spinup_monthly e0004s_transientCO2 Check the July 2016 CO2, LAI, GPP, NPP and NEE output on GrADS. Looks good. Verify the output using ~/Catchment/SMAP_M09/verify_e0004s_transientCO2.gs cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/ mv e0004s_transientCO2 e0004s_transientCO2_03 cd e0004s_transientCO2_03 rm -rf lenkf_orig cd output/SMAP_EASEv2_M09_GLOBAL/rs/ens0000/ /bin/rm -rf Y198? Y199? Y200? Y201[0-6] (NOTE: DO NOT do this for the last, i.e. 5th, cycle!!! Keep all the restart files!!!) mkdir -p Y1981/M01 cd Y1981/M01 ln -s ../../Y2017/M01/e0004s_transientCO2.ens0000.catchcn_ldas_rst.20170101_0000z e0004s_transientCO2.ens0000.catchcn_ldas_rst.19810101_0000z cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run nedit M09_CN_e0004s_transientCO2.exe & restart_path = /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_03/output/ Run ldsetup to set up a new experiment: cd /discover/nobackup/fzeng/LDASsa_m3-16_0_p2_CatchCatchCN_for_MERRA3/exec/e0004s_transientCO2/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_transientCO2.exe /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/run/M09_e0004s_transientCO2.bat --runmodel --monthsperjob 10 --landmodel catchCN (NOTE: don't use the "--submit" option because we need to modify the job scripts, see below) Check the executable and restart file: cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2 ls -l build/Linux/bin/LDASsaCN_mpi.x (to make sure the executable is the right one) -rwxr-xr-x 1 fzeng g0620 69819515 2017-04-27 10:10 build/Linux/bin/LDASsaCN_mpi.x* ls -l input/restart/ (to make sure the restart file is the right one) lrwxrwxrwx 1 fzeng g0620 80 2017-06-20 10:39 output -> /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_03/output/ Copy year_co2.txt and the *all.txt files from previous cycle: cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2_03/run cp -p year_co2.txt *all.txt /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2/run/. cd /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2/run Check if these files are correct, particularly year_co2.txt! cat year_co2.txt (to make sure the year of CO2 is correct) 1945 (correct!) nedit *all.txt & Modify the job scripts before submission: Under /discover/nobackup/fzeng/Catchment/SMAP_EASEv2_M09/e0004s_transientCO2/run: mkdir ../lenkf_orig (Important: lenkf_orig has to be one level up. Otherwise the lenkf.*.j files in lenkf_orig will be counted by edit_lenkf.csh) cp -p lenkf.*.j ../lenkf_orig/. cp ../../e0004s_transientCO2_03/run/edit_lenkf.csh . ./edit_lenkf.csh Check the lenkf.*.j files and they look correct. qsub lenkf.0.j NOTE: for transientCO2 run, (1) Need to use edit_lenkf.csh to add these lines below to the end of each job file lenkf.N.j (as above): cat year_ts.txt >> year_ts_all.txt cat EEAco2g_ts.txt >> EEAco2g_ts_all.txt cat co2_scalar_ts.txt >> co2_scalar_ts_all.txt qsub lenkf.N+1.j (2) Need to copy edit_lenkf.csh, year_co2.txt, year_ts_all.txt, EEAco2g_ts_all.txt and co2_scalar_ts_all.txt from previous cycle. 2. 3hrly CO2, GR=0.25GPP run: Processed the available output: ~/Catchment/M2n5P_m0001 > tile2grid_FV_144x91_monthly M2.n5P.3hco2.glb.GRpt25 3. Submitted my comments for manuscript review. 4. Emailed Randy the results from the SMAP M09 transient CO2 run: GPP is not increasing everywhere. 5. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes/notes_daily. ========== 20170621: 1. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes/notes_daily. ========== 20170622: 1. The M2.n5P.3hco2.glb.GRpt25 run is done. Processed the output: ~/Catchment/M2n5P_m0001 > tile2grid_FV_144x91_monthly M2.n5P.3hco2.glb.GRpt25 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes/notes_daily. ========== 20170623: 1. Set up two experiments: /discover/nobackup/fzeng/offline_code/LDASsa_m3-15_2-CN/src/Components/GEOSlana_GridComp/M2n5P_m0001CTco2: process_cat.F90.3hrlyCO2 process_cat.F90.3hrlyCO2+200ppm compute_rc.F90 CNGRespMod.F90.GRpt27 CNAllocationMod.F90.GRpt27 Code change in these experiments: src/Components/GEOSsurface_GridComp/GEOSland_GridComp/GEOScatch_GridComp/CNGRespMod.F90 & CNAllocationMod.F90 Line 153: grperc = 0.3 Changed to grperc = 0.27 (1) 3hrly CO2, GR=0.27GPP: cd /discover/nobackup/fzeng/Catchment/M2n5P/m0001 mkdir -p M2.n5P.3hco2.glb.GRpt27/RUN/rs/ens0000/Y2001/M01 cd M2.n5P.3hco2.glb.GRpt27 cp -p /discover/nobackup/elee15/offline/sims/M2.n5P.3htCO2.c5/CN_restart . cp -p /discover/nobackup/elee15/offline/sims/M2.n5P.3htCO2.c5/RUN/rs/ens0000/Y2001/M01/mM2.n5P.3htCO2.ens0000.catch_ldas_rst.20010101_0000z.bin RUN/rs/ens0000/Y2001/M01/mM2.n5P.3hco2.glb.GRpt27.ens0000.catch_ldas_rst.20010101_0000z.bin cp /discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb.GRpt25/lenkf.j . discover08:/discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb.GRpt27 > ls -l CN_restart -rw-r--r-- 1 fzeng g0620 1530820456 2016-11-08 19:34 CN_restart discover08:/discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb.GRpt27 > ls -l RUN/rs/ens0000/Y2001/M01/ total 40512 -rw-r--r-- 1 fzeng g0620 41453776 2016-11-08 19:34 mM2.n5P.3hco2.glb.GRpt27.ens0000.catch_ldas_rst.20010101_0000z.bin (2) 3hrly+200ppm CO2, GR=0.27GPP: cd /discover/nobackup/fzeng/Catchment/M2n5P/m0001 mkdir -p M2.n5P.3hco2+200ppm.glb.GRpt27/RUN/rs/ens0000/Y2001/M01 cd M2.n5P.3hco2+200ppm.glb.GRpt27 cp -p /discover/nobackup/elee15/offline/sims/M2.n5P.3htCO2.c5/CN_restart . cp -p /discover/nobackup/elee15/offline/sims/M2.n5P.3htCO2.c5/RUN/rs/ens0000/Y2001/M01/mM2.n5P.3htCO2.ens0000.catch_ldas_rst.20010101_0000z.bin RUN/rs/ens0000/Y2001/M01/mM2.n5P.3hco2+200ppm.glb.GRpt27.ens0000.catch_ldas_rst.20010101_0000z.bin cp /discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2.glb.GRpt27/lenkf.j . discover08:/discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2+200ppm.glb.GRpt27 > ls -l CN_restart -rw-r--r-- 1 fzeng g0620 1530820456 2016-11-08 19:34 CN_restart discover08:/discover/nobackup/fzeng/Catchment/M2n5P/m0001/M2.n5P.3hco2+200ppm.glb.GRpt27 > ls -l RUN/rs/ens0000/Y2001/M01/ total 40512 -rw-r--r-- 1 fzeng g0620 41453776 2016-11-08 19:34 mM2.n5P.3hco2+200ppm.glb.GRpt27.ens0000.catch_ldas_rst.20010101_0000z.bin 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes/notes_daily. ========== 20170626: 1. Processed the output: ~/Catchment/M2n5P_m0001 > tile2grid_FV_144x91_monthly M2.n5P.3hco2.glb.GRpt27 ~/Catchment/M2n5P_m0001 > tile2grid_FV_144x91_monthly M2.n5P.3hco2+200ppm.glb.GRpt27 2. Worked on upgrading the CN model in Catchment-CN to CLM4.5. See notes /discover/nobackup/fzeng/clm4-to-clm4.5/notes/notes_daily.