adding tracers

[gbolzon]

Adding 2 tracers to:

• BFMcoupler_LOAD.h
• BFMcoupler_VARS.h
• BFMcoupler_calc_tendency.F
• BFMcoupler_fields_load.F
• BFMcoupler_ini_forcing.F
• BFMcoupler_readparams.F

[gbolzon]

diff_apply deve fare su GAD_OPTIONS.h

-#undef GAD_ALLOW_TS_SOM_ADV
+#define GAD_ALLOW_TS_SOM_ADV

E su OBCS_PARAMS.h
PARAMETER ( OBCS_tideCompSize = 8 )

[gbolzon]

Su EXF_OPTIONS.h

#define  ALLOW_RUNOFTEMP
# define EXF_CALC_ATMRHO
# define ALLOW_ZENITHANGLE
#define USE_EXF_INTERPOLATION

[gbolzon]

@stefanoquerin I thing the code is ready

• diff_apply covers all files we need (excluding swfrac.F)
• CONC01, CONC02 ... are the new names

[gbolzon]

The mechanism:

LINES, position_line = strings_and_position(inputfile, searchstring)
diff_apply.insert_lines(...new_lines,, position_line)

has this behaviour:
If searchstring corresponds to line 77 of inputfile, ==> position_line=76
No, inside insert_lines(), the firsts 76 lines are replicated.
Then we add the new lines.
So, they are applied just before line 77, while the docstring says 'after'

[gbolzon]

We are correcting here some bugs in migration 66j --> 68z

1. gchem_apply_tendency needs 2 new variables defined in gchem_tr_register.F, which needs BFMcoupler_tr_register.F

We need MITgcm/pkg/gchem/gchem_tr_register.F with this patch

+#ifdef ALLOW_BFMCOUPLER
+      IF ( useBFMcoupler ) THEN
+        CALL BFMcoupler_TR_REGISTER(
+     U                gchem_Tracer_num, gchem_sepFTr_num,
+     I                myThid )
+      ENDIF
+#endif

applied before line 90 WRITE(msgBuf,'(2A,2I5)') 'GCHEM_TR_REGISTER: ',

2. In order to read the added tracers, the main cycle of

   MITgcmBFM_coupler/BFMcoupler_fields_load.F

   Lines 256 to 257 in ba41186

   	DO bj = myByLo(myThid), myByHi(myThid)
   	DO bi = myBxLo(myThid), myBxHi(myThid)

   
   must be closed after the added tracers, not here!
   

   MITgcmBFM_coupler/BFMcoupler_fields_load.F

   Lines 383 to 388 in ba41186

   	ENDDO
   	ENDDO
   	C Adding tracers
   	IF ( BFMcoupler_CONC01BotForcFile.NE. ' ' ) THEN
   	DO j=1-OLy,sNy+OLy

3. In BFMcoupler_calc_tendency.F decay of passive tracers is working only with BFMv2