Many minor updates to my version of the model in Code_withCellModel
There are 2 changes that are also relevant to other versions of the OCIM-BGC model.
in neglogpost.m: fix bug in save DOPl output
In the model output structure, both the DOP and DOPl fields contained the same DOP output due to an index error when extracting eqPcycle model solution in neglogpost
changed:
DOPl(iwet) = P(1+2*nwet:3*nwet) ;
to:
DOPl(iwet) = P(1+3*nwet:4*nwet) ;
in make_datafile_24layer/bin_WOA_nut.m: convert WOA data from umol/kg to mmol/m^3 using the decadal average seawater density from WOA, because OCIM-BGC assumes nutrient data units are in mmol/m^3.
I saved the converted and regridded data to new mat file in DATA/BGC_24layer: fileName = O2_Nut_WOA_mmolperm3_91x180x24.mat
(FYI, the commit messages starting with "(LOOK)" include changes that are not specific to the cell model/ my C2P function project)
Many minor updates to my version of the model in
Code_withCellModelThere are 2 changes that are also relevant to other versions of the OCIM-BGC model.
in
neglogpost.m: fix bug in save DOPl output In the model output structure, both the DOP and DOPl fields contained the same DOP output due to an index error when extracting eqPcycle model solution in neglogpost changed:DOPl(iwet) = P(1+2*nwet:3*nwet) ;to:DOPl(iwet) = P(1+3*nwet:4*nwet) ;in
make_datafile_24layer/bin_WOA_nut.m: convert WOA data from umol/kg to mmol/m^3 using the decadal average seawater density from WOA, because OCIM-BGC assumes nutrient data units are in mmol/m^3. I saved the converted and regridded data to new mat file inDATA/BGC_24layer: fileName =O2_Nut_WOA_mmolperm3_91x180x24.mat(FYI, the commit messages starting with "(LOOK)" include changes that are not specific to the cell model/ my C2P function project)