/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.6 by stephd, Thu May 31 21:08:25 2012 UTC revision 1.8 by jahn, Sat Jun 30 19:23:04 2012 UTC
# Line 21  C======================================= Line 21  C=======================================
21  #include "PARAMS.h"  #include "PARAMS.h"
22  #include "GRID.h"  #include "GRID.h"
23  #include "DYNVARS.h"  #include "DYNVARS.h"
24  #ifdef USE_QSW  c for Qsw and/or surfaceForcingT
25    c choice which field to take pCO2 from for pCO2limit
26    c this assumes we use Ttendency from offline
27  #include "FFIELDS.h"  #include "FFIELDS.h"
 #endif  
28  #ifdef ALLOW_LONGSTEP  #ifdef ALLOW_LONGSTEP
29  #include "LONGSTEP.h"  #include "LONGSTEP.h"
30  #endif  #endif
# Line 43  c ANNA include wavebands_params.h Line 44  c ANNA include wavebands_params.h
44  #include "WAVEBANDS_PARAMS.h"  #include "WAVEBANDS_PARAMS.h"
45  #endif  #endif
46    
47    
48  C     === Global variables ===  C     === Global variables ===
49  c tracers  c tracers
50         _RL Ptr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,nDarwin)         _RL Ptr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,nDarwin)
# Line 194  c some local variables Line 196  c some local variables
196         _RL PSiupl         _RL PSiupl
197         _RL Tlocal         _RL Tlocal
198         _RL Slocal         _RL Slocal
199           _RL pCO2local
200         _RL Qswlocal         _RL Qswlocal
201         _RL NH4l         _RL NH4l
202         _RL NO2l         _RL NO2l
# Line 810  c for explicit sinking of particulate ma Line 813  c for explicit sinking of particulate ma
813               Slocal = salt(i,j,k,bi,bj)               Slocal = salt(i,j,k,bi,bj)
814  #endif  #endif
815    
816    c choice where to get pCO2 from
817    c taking from igsm dic run - fed through Tflux array
818    c               pCO2local=surfaceForcingT(i,j,bi,bj)
819    c or from darwin carbon module
820    #ifdef ALLOW_CARBON
821                   pCO2local=pCO2(i,j,bi,bj)
822    #else
823                   pCO2local=280. _d -6
824    #endif
825    
826               freefu = max(freefe(i,j,k),0. _d 0)               freefu = max(freefe(i,j,k),0. _d 0)
827               if (k.eq.1) then               if (k.eq.1) then
828                 inputFel = inputFe(i,j,bi,bj)                 inputFel = inputFe(i,j,bi,bj)
# Line 926  c ANNA pass extra variables if WAVEBANDS Line 939  c ANNA pass extra variables if WAVEBANDS
939       I                       pofeupl, psiupl,       I                       pofeupl, psiupl,
940       I                       PARl,       I                       PARl,
941       I                       Tlocal, Slocal,       I                       Tlocal, Slocal,
942         I                       pCO2local,
943       I                       freefu, inputFel,       I                       freefu, inputFel,
944       I                       bottom, dzlocal,       I                       bottom, dzlocal,
945       O                       Rstarl, RNstarl,       O                       Rstarl, RNstarl,
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22