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


Ignore:
Timestamp:
2017-03-16T11:39:11+01:00 (4 years ago)
Author:
frrh
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/1851/General

    v62 v63  
    513513 
    514514 
    515  
    516  
    517  
     515Poking 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.  
     516 
     517The get values from: 
     518{{{ 
     519               zpho = max(0.,trn(ji,jj,jk,jpdin)) / 16.0 !! phosphate via DIN and Redfiel               
     520               zdic = max(0.,trn(ji,jj,jk,jpdic)) !! dissolved inorganic carbon 
     521               zalk = max(0.,trn(ji,jj,jk,jpalk)) !! alkalinity 
     522}}} 
     523 
     524 
     525So the trn array is implicated - but interesting only some fields in it.  
     526 
     527 
     528 
     529