/[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.7 by stephd, Fri Jun 29 20:41:59 2012 UTC revision 1.10 by jahn, Thu Aug 23 21:48:24 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    
 c choice which field to take pCO2 from for pCO2limit  
 c this assumes we use Ttendency from offline  
 #include "FFIELDS.h"  
   
47    
48  C     === Global variables ===  C     === Global variables ===
49  c tracers  c tracers
# Line 98  c ANNA define variables for wavebands Line 95  c ANNA define variables for wavebands
95  #ifdef DAR_RADTRANS  #ifdef DAR_RADTRANS
96         integer iday,iyr,imon,isec,lp,wd,mydate(4)         integer iday,iyr,imon,isec,lp,wd,mydate(4)
97         _RL Edwsf(tlam),Eswsf(tlam)         _RL Edwsf(tlam),Eswsf(tlam)
98         _RL Edz(tlam,Nr),Esz(tlam,Nr),Euz(tlam,Nr),Eutop(tlam,Nr)         _RL Edz(tlam,Nr),Esz(tlam,Nr),Euz(tlam,Nr)
99           _RL Estop(tlam,Nr),Eutop(tlam,Nr)
100         _RL tirrq(nr)         _RL tirrq(nr)
101         _RL tirrwq(tlam,nr)         _RL tirrwq(tlam,nr)
102           _RL amp1(tlam,nr), amp2(tlam,nr)
103         _RL solz         _RL solz
104         _RL rmud         _RL rmud
105         _RL actot,bctot,bbctot         _RL actot,bctot,bbctot
# Line 627  c           add water and CDOM Line 626  c           add water and CDOM
626              bt_k(k,ilam) = bw(ilam) + bctot + bpart_k(k,ilam)              bt_k(k,ilam) = bw(ilam) + bctot + bpart_k(k,ilam)
627              bb_k(k,ilam) = darwin_bbw*bw(ilam)+bbctot+bbpart_k(k,ilam)              bb_k(k,ilam) = darwin_bbw*bw(ilam)+bbctot+bbpart_k(k,ilam)
628              bb_k(k,ilam) = MAX(darwin_bbmin, bb_k(k,ilam))              bb_k(k,ilam) = MAX(darwin_bbmin, bb_k(k,ilam))
629    c           initialize output variables
630                Edz(ilam,k) = 0.0
631                Esz(ilam,k) = 0.0
632                Euz(ilam,k) = 0.0
633                Estop(ilam,k) = 0.0
634                Eutop(ilam,k) = 0.0
635                amp1(ilam,k) = 0.0
636                amp2(ilam,k) = 0.0
637            ENDDO            ENDDO
638           ENDDO           ENDDO
639    
# Line 636  c           add water and CDOM Line 643  c           add water and CDOM
643       I                darwin_radtrans_kmax,darwin_radtrans_niter,       I                darwin_radtrans_kmax,darwin_radtrans_niter,
644       O                Edz,Esz,Euz,Eutop,       O                Edz,Esz,Euz,Eutop,
645       O                tirrq,tirrwq,       O                tirrq,tirrwq,
646         O                amp1,amp2,
647       I                myThid)       I                myThid)
648  #else  #else
649  c dzlocal ?????  c dzlocal ?????
# Line 1338  Coj            no Eu at surface (yet) Line 1346  Coj            no Eu at surface (yet)
1346                 Euave(i,j,k,bi,bj,ilam)=Euave(i,j,k,bi,bj,ilam)+                 Euave(i,j,k,bi,bj,ilam)=Euave(i,j,k,bi,bj,ilam)+
1347       &                                 Euz(ilam,k-1)*dtplankton       &                                 Euz(ilam,k-1)*dtplankton
1348                endif                endif
1349                  Estave(i,j,k,bi,bj,ilam)=Estave(i,j,k,bi,bj,ilam)+
1350         &                                 Estop(ilam,k)*dtplankton
1351                Eutave(i,j,k,bi,bj,ilam)=Eutave(i,j,k,bi,bj,ilam)+                Eutave(i,j,k,bi,bj,ilam)=Eutave(i,j,k,bi,bj,ilam)+
1352       &                                 Eutop(ilam,k)*dtplankton       &                                 Eutop(ilam,k)*dtplankton
1353               enddo               enddo
1354  #endif  #endif
1355    #ifdef DAR_DIAG_IRR_AMPS
1356                 do ilam = 1,tlam
1357                   amp1ave(i,j,k,bi,bj,ilam)=amp1ave(i,j,k,bi,bj,ilam)+
1358         &                                 amp1(ilam,k)*dtplankton
1359                   amp2ave(i,j,k,bi,bj,ilam)=amp2ave(i,j,k,bi,bj,ilam)+
1360         &                                 amp2(ilam,k)*dtplankton
1361                 enddo
1362    #endif
1363  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
1364               do ilam = 1,tlam               do ilam = 1,tlam
1365                 aave(i,j,k,bi,bj,ilam)=aave(i,j,k,bi,bj,ilam)+                 aave(i,j,k,bi,bj,ilam)=aave(i,j,k,bi,bj,ilam)+
# Line 1366  Coj            no Eu at surface (yet) Line 1384  Coj            no Eu at surface (yet)
1384       &                                 bbpart_k(k,ilam)*dtplankton       &                                 bbpart_k(k,ilam)*dtplankton
1385               enddo               enddo
1386  #endif  #endif
1387    #ifdef DAR_RADTRANS
1388                 if (k.eq.1) then
1389                   rmudave(i,j,bi,bj)=rmudave(i,j,bi,bj)+
1390         &                                 rmud*dtplankton
1391                 endif
1392    #endif
1393  #ifdef DAR_DIAG_RSTAR  #ifdef DAR_DIAG_RSTAR
1394               do np=1,npmax               do np=1,npmax
1395                 Rstarave(i,j,k,bi,bj,np)=Rstarave(i,j,k,bi,bj,np)+                 Rstarave(i,j,k,bi,bj,np)=Rstarave(i,j,k,bi,bj,np)+
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22