[MITgcm-support] A question about gchem pkg

[Jean-Michel Campin]

Hi Haoran,

The model stops when PTRACERS_numInUse is reduced to 4 because the setting is not consistent:
pkg/ptracers holds the passive tracer array for pkg/gchem, and if there are more tracers
used in pkg/gchem than defined in pkg/ptracers it will just not fit.
Note that the opposite is fine (having more passive tracers than gchem tracers is OK).

So you might want to clarify why you want to have 4 ptracers (instead of 5 as used in gchem).
If you are currently using pkg/dic with 5 tracers (without iron, i.e., #undef ALLOW_FE) and 
want to drop the oxigen (changing #define ALLOW_O2 to #undef ALLOW_O2 in DIC_OPTIONS.h)
then you can make this change in dic/gchem first (might need to adjust some parameter files
"data.*") and leave ptracers unchanged. And when this is working then you can decide
to reduce PTRACERS_numInUse to just 4 (since then you will only be using 4 tracers in gchem).

Cheers,
Jean-Michel

On Fri, Jul 16, 2021 at 02:25:50PM +0800, Haoran Xu wrote:
> Hope you all going well!
> I am using the gchem pkg. One error occurs when i change the number of the ptracer to 4.It shows as following:
> 
> 
> (PID.TID 0062.0001) *** ERROR *** GCHEM_CHECK: Number of GCHEM tracers: gchem_Tracer_num =    5
> (PID.TID 0062.0001) *** ERROR *** GCHEM_CHECK:  exceeds number of pTr: PTRACERS_numInUse =    4
> (PID.TID 0062.0001) *** ERROR *** GCEHM_CHECK: detected  1 fatal error(s)
> (PID.TID 0062.0001) *** ERROR *** S/R ALL_PROC_DIE: ending the run
> 
> 
>  I donont know how to change the parameter  gchem_Tracer_num.
>  Is there any idea?

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support