Changes between Version 62 and Version 63 of ticket/1851/General

Timestamp:

2017-03-16T11:39:11+01:00 (7 years ago)

Author:

frrh

Comment:

--

ticket/1851/General

@@ -513 +513 @@
 
 
-
-
-
+Poking into the MEDUSA code with glob_sum on key variables leads us to CO2FLUX_out_cpl then zn_co2_flx then fgco2 which diverges on TS 34. This takes us into the massive trcbio_medusa where calculations are done in an enormous loop on a cel by cell basis - so we don't gobal arrays to examine and the numbers of fields involved renders it impossible to examine each scalar. So I have to introduce special new arrays to store scalars of interest and sum them outside the loop for comparison. This leads to the call to mocsy_interface at around line 1800. Examining all declared "input" fields shows diffs in zpho, zdic and zalk. 
+
+The get values from:
+{{{
+               zpho = max(0.,trn(ji,jj,jk,jpdin)) / 16.0 !! phosphate via DIN and Redfiel              
+               zdic = max(0.,trn(ji,jj,jk,jpdic)) !! dissolved inorganic carbon
+               zalk = max(0.,trn(ji,jj,jk,jpalk)) !! alkalinity
+}}}
+
+
+So the trn array is implicated - but interesting only some fields in it. 
+
+
+
+