#===============================================================# # # prepetm.inp # # LOSCAR Model: Long-term Ocean-atmosphere-Sediment # CArbon cycle Reservoir Model # # *** LOSCAR comes with ABSOLUTELY NO WARRANTY *** # *** Use at your own risk. DO NOT DISTRIBUTE *** # # Model parameter (control) input file # example: pre-PETM default # # When using LOSCAR, cite as: # # Zeebe, R. E., Zachos, J. C., and Dickens, G. R. Carbon dioxide # forcing alone insufficient to explain Paleocene-Eocene Thermal # Maximum warming. Nature Geoscience, doi:10.1038/NGEO578 (2009) # # Richard E. Zeebe # School of Ocean and Earth # Science and Technology # Department of Oceanography # University of Hawaii at Manoa # 1000 Pope Road, MSB 504 # Honolulu, HI 96822, USA # email: loscar.model@gmail.com # # # NOTES # # Input line format: # # PARAMETER VALUE1 VALUE2 ... # (delimited by spaces) # # *** Do not change input format of a single line # (yields errors) # # - Parameter sequence (order of lines) can be random. # - If parameters are omitted, internal default values # are used. # - Lines starting with '#' are ignored. # - Generally the following convention is used: # TRUE = ON = 1 # FALSE = OFF = 0 # # *** Change one parameter in small increments at a time # (not all at once in big steps). # # - New steady-state: If a large parameter change is needed # for a new steady-state (but the model aborts), you may # try applying a series of small changes and running the # model into a new steady-state between each increment # (restart from there). # # # parameter def.value unit note # --------------------------------------------------------- # RESTART dat/y0prepetm-2.x.dat load restart file # SVSTART dat/y0prepetm-2.x.dat save file for restart # EMSFILE dat/Emss/1000_0500.dat load fossil fuel emission file # TSTART 0. years >= 0.0 # TFINAL 200.e3 years > TSTART # CINP 1000. Pg C Carbon input (into atmosphere) # D13CIN -55. per mil d13C Carbon input # TCIN0 0. years time start C input # TCSPAN 6.e3 years time span C input (duration) # EPSLV 1.e-4 - solver accuracy (1.e-4 ~sloppy) # FCSML 0.05 - numerics: linear f_calcite drop # during dissolution if fc < FCSML. # KMAX 1000 - approximate # output values. # (range 2 to MAXSTP+1, see below) # TSNS 0 - switch Temp sensitivity to pCO2 ON/OFF (1/0) # SCLIM 3.0 deg C T-sensitivity (2.0-? C), requires TSNS 1 # TEMP 25. 25. 25. 16. 16. 16. 12. 12. 12. 12. 18. 14. 12. # deg C temp. boxes (3Low,3Int,3Deep,1High,3Tethys) # SAL 34.72 ... - salinity boxes # PCO2SI 1000.0 ppmv Final steady-state pCO2, balancing # long-term Si weath fluxes (> 200 ppmv) # FDICI 00.00 % change init. DIC all boxes (range +- x%) # FALKI 00.00 % change init. ALK all boxes (range +- x%) # FPO4I 00.00 % change init. PO4 all boxes (range +- x%) # THC 25. Sv conveyor transport # FBIOL 0.80 - Efficiency Low-Lat Corg Biopump # (useful range 0.70 to 0.95) # CBIOH 1.8 mol C/m2/y Export High-Lat Corg Biopump # RRAIN 6.7 - rain ratio (Corg/CaCO3) # FSHLF 4.5 - raise shelf relative deep rain, see below # (arbitrary factor, useful range: 1.0-6.0) # FINC 15.8e12 mol C/y Initial CaCO3 riverine flux # CALC 20.0e-3 mol/kg Seawater [Ca2+] # MAGN 30.0e-3 mol/kg Seawater [Mg2+] # --------------------------------------------------------- # # More notes # # KMAX: Solver stores the solution during integration at # approximate time intervals (TFINAL-TSTART)/KMAX. # KMAX may be set larger than nstep, the actual #steps # taken, yet #output <= nstep. MAXSTP is the max number # of steps allowed. MAXSTP=100000 (V 2.0.3). # # TSNS: default = 0 (recommended). Set to 1 to obtain results of # Zeebe et al., Science (2008). For other runs, be careful # when setting TSNS=1 (can lead to surprises e.g. when pCO2 # drops below initial pCO2 value). # # FINC: Initial (!) CaCO3 riverine flux. The riverine flux # during the run can change depending on pCO2. # # FSHLF: Raising FSHLF increases shelf CaCO3 rain to shallow # water sediments relative to pelagic CaCO3 rain to deep # sea sediments (basin to shelf partitioning at constant # total rain). For example, FSHLF = 4.5 in the Paleo setup # reduces pelagic rain by 22% (relative to FSHLF = 1.0) # on average in Atl, Ind, Pac. # # #===============================================================# # input starts below this line # declare Restart files to be used and saved #RESTART dat/y0prepetm-2.0.4.dat #RESTART dat/restart_jan_equi_1500ppm.dat RESTART dat/restart_jan_1Myr.dat #SVSTART dat/y0prepetm-2.0.4.dat #SVSTART dat/restart_jan_equi_1500ppm.dat #SVSTART dat/restart_jan_1Myr.dat #EMSFILE dat/Emss/1000_0500.dat #EMSFILE dat/Emss/carbon_emission_1Myr_22100GtC.dat EMSFILE dat/Emss/carbon_emission_1Myr_22100GtC_3rd4thPulseLong.dat # End year of Run TFINAL-TSTART is length of run # To equilibrate at new pCO2 needs approx. 5 M.yr. #TSTART 0. TSTART 1000.e3 #TFINAL 200.e3 TFINAL 4000.e3 #TFINAL 5000.e3 # Carbon input can be declared here or in EMSFILE above. Here, there will be constant C emission over TCSPAN years # Has to be commented out if no carbon input is desired #CINP 1000. D13CIN -20.0 #TCIN0 0. #TCSPAN 6.e3 #For my emission scenarios EPSLV has to be lowered to approx 1.e-5 which increases computing time significantly #(10-15min) #EPSLV 1.e-4 #EPSLV 5.e-5 EPSLV 1.e-4 FCSML 0.05 KMAX 1000 #KMAX 10000 # Sensitivity of Ocean Temp to pCO2; 0 or 1 TSNS 1 # Climate Sensitivity, 3.0 , Climber exhibits sensitivity of 2.86 in Late Triassic runs SCLIM 2.86 #TEMP 25. 25. 25. 16. 16. 16. 12. 12. 12. 12. 18. 14. 12. # Try to use ocean temperatures simulated by Climber # Average low-latitudes (55S:55N) over 0:100m, 100:1000m, 1000:4000m -> ~23, 12, 2 degC #TEMP 23. 23. 23. 12. 12. 12. 2. 2. 2. 5. 27. 14. 3. TEMP 23. 23. 23. 12. 12. 12. 2. 2. 2. 5. 23. 12. 2. # Average Salinity in Climber is also 34.73 SAL 34.72 34.72 34.72 34.72 34.72 34.72 34.72 34.72 34.72 34.72 34.72 34.72 34.72 # jan #PCO2SI 1500.0 PCO2SI 1500.0 FDICI 00.00 FALKI 00.00 FPO4I 00.00 # Default for PETM was 25, for pre-Industrial 20; maybe weaker case more suitable for end-Triassic #THC 25. THC 20. FBIOL 0.80 CBIOH 1.8 RRAIN 6.7 # Increasing FSHLF from 4.5 to 6 leads to a slightly slower drawdown of CO2 #FSHLF 4.5 FSHLF 6.0 FINC 15.83409493e12 CALC 20.0e-3 MAGN 30.0e-3