derpycode / cgenie.muffin

cGENIE.muffin Earth system model

Installation instructions and a series of tutorials for running the model can be found in the muffindoc repository. (Also see: muffin3.install.Ubuntu.pdf and muffin3.install.MacOS.pdf for a simpler installation processing when using the master.python3 branch.)

cgenie.muffin's Issues

Update to Python3

Seem that Python code in cgenie.muffin is python27 code.
Python3 is newer and more popular.

Python interface to create cgenie input map files?

There is need of create Python/r or some free language program to create input files for
cgenie, because all users does not have Matlab.

licence

At the moment cGENIE is not free. Is there a licence intended but missing?

Docker container for cgenie?

Hi

Is there a Linux virtual machine or Docker containerised version of cgenie available, or a set of build instructions for it?

I'm trying to evaluate how Docker containers might be used to support various applications in a distance education setting and a colleague tipped me off to this application as one that would be handy to try out for a course they're involved with...

How do I simulate adding alkalinity into a specific grid index compared a baseline scenario?

Hello! cgenie.muffin novice here. I've managed to cobble together a simulation for RCP 2.6. I'd like to see what would happen to the global DIC if I injected around 10,000 tons of DIC and Ca into a specific grid index for a year. Do you know how I could go about doing that?

What I've done so far:

• Simulate RCP 2.6
• Create a copy of the configuration files and forcing files for RCP 2.6 (hence known as RCP_MOD)
• Add into RCP_MOD biogem_force_flux_ocn_Ca_sig.dat and biogem_force_flux_ocn_DIC_sig.dat, both of which have the following:

-START-OF-DATA-
0.0      0.0
2030.0   0.0
2030.0   10000
2031.0   10000
2031.0   0
999999.9 0.0
-END-OF-DATA-

• Add configure_forcings_ocn.dat into RCP_MOD with the following:

  00 01 02 03   04 05 06  07  08  09
  \/ \/ \/ \/   \/ \/ \/  \/  \/  \/

-START-OF-DATA-
 03  f  f  0.1  T  T   0  22  11  16  '[dissolved inorganic carbon (DIC) (mol kg-1)]'
 35  f  f  0.1  T  T   0  22  11  16  '[Ca (mol kg-1)]'
-END-OF-DATA-

  /\ /\ /\ /\   /\ /\ /\  /\  /\  /\
  00 01 02 03   04 05 06  07  08  09

DATA FORMAT AND ORDER
---------------------

COLUMN #00: TRACER NUMBER
COLUMN #01: include restoring forcing of tracer? [DATA FORMAT: 'T'/'F' ('t'/'f')]
COLUMN #02: restrict restoring forcing to surface? [DATA FORMAT: 'T'/'F' ('t'/'f')]
COLUMN #03: time constant of restoring forcing (years) [DATA FORMAT: real]
COLUMN #04: include flux forcing of tracer? [DATA FORMAT: 'T'/'F' ('t'/'f')]
COLUMN #05: scale flux forcing of tracer? [DATA FORMAT: 'T'/'F' ('t'/'f')]
COLUMN #06: make forcing uniform over this dimension (3 == 3D, 2 == 2D, 0 == point, -1 == SURFACE, -2 == BENTHIC, ELSE spatially explicit forcing) [DATA FORMAT: interger]
COLUMN #07: i grid location of point forcing [DATA FORMAT: interger]
COLUMN #08: j grid location of point forcing [DATA FORMAT: interger]
COLUMN #09: k grid location of point forcing [DATA FORMAT: interger]
COLUMN END: TRACER DESCRIPTION

• Modify the configuration file to have bg_par_ocn_force_scale_val_3=22722.1086116792 and bg_par_ocn_force_scale_val_35=22722.1086116792

However, after doing all of this and running RCP_MOD, I strangely get less global DIC when I look at biogem_series_ocn_DIC.res. Do you know why this would be the case? Is this because adding a forcing flux will override the existing flux calculated from the simulation instead of adding onto the existing flux? Any advice on how I should go about trying to accomplish my goal?

I'll attach all of the files I reference. If you made it this far, thanks for even taking the time to read this over!

rcp2p6_config_file.SPIN.txt
rcp2p6_forcing_files.zip
rcp2p6_mod_config_file.SPIN.txt
rcp2p6_mod_forcing_files.zip
rcp2p6_mod_simulation_results.SPIN.zip
rcp2p6_simulation_results.SPIN.zip