/[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.12 by jahn, Fri Aug 24 19:45:36 2012 UTC revision 1.13 by stephd, Tue Oct 23 16:39:32 2012 UTC
# Line 92  c ANNA define variables for wavebands Line 92  c ANNA define variables for wavebands
92         _RL PARw_k(tlam,Nr)         _RL PARw_k(tlam,Nr)
93         _RL PARwup(tlam)         _RL PARwup(tlam)
94         _RL acdom_k(Nr,tlam)         _RL acdom_k(Nr,tlam)
95           _RL Ek_nll(npmax,tlam)
96           _RL EkoverE_nll(npmax,tlam)
97  #ifdef DAR_RADTRANS  #ifdef DAR_RADTRANS
98         integer iday,iyr,imon,isec,lp,wd,mydate(4)         integer iday,iyr,imon,isec,lp,wd,mydate(4)
99         _RL Edwsf(tlam),Eswsf(tlam)         _RL Edwsf(tlam),Eswsf(tlam)
# Line 131  C      always need for diagnostics Line 133  C      always need for diagnostics
133  #ifdef GEIDER  #ifdef GEIDER
134        _RL phychl(npmax)        _RL phychl(npmax)
135        _RL phychl_k(npmax,Nr)        _RL phychl_k(npmax,Nr)
136          _RL Ekl(npmax)
137          _RL EkoverEl(npmax)
138          _RL chl2cl(npmax)
139  #ifdef DYNAMIC_CHL  #ifdef DYNAMIC_CHL
140        _RL dphychl(npmax)        _RL dphychl(npmax)
141        _RL chlup(npmax)        _RL chlup(npmax)
142          _RL accliml(npmax)
143  #endif  #endif
144  #endif  #endif
145  #ifdef ALLOW_CDOM  #ifdef ALLOW_CDOM
# Line 964  c set other arguments to zero Line 970  c set other arguments to zero
970                  NfixPl(np)=0. _d 0                  NfixPl(np)=0. _d 0
971  #endif  #endif
972  #endif  #endif
973    #ifdef DAR_DIAG_PARW
974                    chl2cl(np)=0. _d 0
975    #endif
976    #ifdef DAR_DIAG_EK
977                    Ekl(np)=0. _d 0
978                    EkoverEl(np)=0. _d 0
979                    do ilam=1,tlam
980                      Ek_nll(np,ilam)=0. _d 0
981                      EkoverE_nll(np,ilam)=0. _d 0
982                    enddo
983    #endif
984               enddo               enddo
985    
986    
# Line 1033  c ANNA pass extra variables if WAVEBANDS Line 1050  c ANNA pass extra variables if WAVEBANDS
1050  #endif  #endif
1051  #ifdef GEIDER  #ifdef GEIDER
1052       O                       phychl,       O                       phychl,
1053    #ifdef DAR_DIAG_EK
1054         I                       Ekl, EkoverEl,
1055    #endif
1056    #ifdef DAR_DIAG_PARW
1057         I                       chl2cl,
1058    #endif
1059  #ifdef DYNAMIC_CHL  #ifdef DYNAMIC_CHL
1060       I                       dphychl,       I                       dphychl,
1061       I                       chlup,       I                       chlup,
1062    #ifdef DAR_DIAG_EK
1063         O                       accliml,
1064    #endif
1065  #endif  #endif
1066  #ifdef ALLOW_CDOM  #ifdef ALLOW_CDOM
1067       O                       dcdoml,       O                       dcdoml,
# Line 1043  c ANNA pass extra variables if WAVEBANDS Line 1069  c ANNA pass extra variables if WAVEBANDS
1069  #endif  #endif
1070  #ifdef WAVEBANDS  #ifdef WAVEBANDS
1071       I                       PARw_k(1,k),       I                       PARw_k(1,k),
1072    #ifdef DAR_DIAG_EK
1073         I                       Ek_nll, EkoverE_nll,
1074    #endif
1075  #endif  #endif
1076  #endif  #endif
1077  #ifdef ALLOW_PAR_DAY  #ifdef ALLOW_PAR_DAY
# Line 1381  c    &                       deltaTclock Line 1410  c    &                       deltaTclock
1410       &                           phychl(np)*dtplankton       &                           phychl(np)*dtplankton
1411               enddo               enddo
1412  #endif  #endif
1413    #ifdef DAR_DIAG_PARW
1414                do ilam=1,tlam
1415                   PARwave(i,j,k,bi,bj,ilam)=PARwave(i,j,k,bi,bj,ilam)+
1416         &                           PARw_k(ilam,k)*dtplankton
1417                enddo
1418                do np=1,npmax
1419                  chl2cave(i,j,k,bi,bj,np)=chl2cave(i,j,k,bi,bj,np)+
1420         &                          chl2cl(np)*dtplankton
1421                enddo
1422    #endif
1423  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
1424  c            print*,'acdom',k,acdom_k(k,darwin_diag_acdom_ilam)  c            print*,'acdom',k,acdom_k(k,darwin_diag_acdom_ilam)
1425               aCDOMave(i,j,k,bi,bj)=aCDOMave(i,j,k,bi,bj)+               aCDOMave(i,j,k,bi,bj)=aCDOMave(i,j,k,bi,bj)+
# Line 1446  Coj            no Eu at surface (yet) Line 1485  Coj            no Eu at surface (yet)
1485       &                                 rmud*dtplankton       &                                 rmud*dtplankton
1486               endif               endif
1487  #endif  #endif
1488    #ifdef DAR_DIAG_EK
1489                do np=1,npmax
1490                 Ekave(i,j,k,bi,bj,np)=Ekave(i,j,k,bi,bj,np)+
1491         &                        Ekl(np)*dtplankton
1492                 EkoverEave(i,j,k,bi,bj,np)=EkoverEave(i,j,k,bi,bj,np)+
1493         &                        EkoverEl(np)*dtplankton
1494                 acclimave(i,j,k,bi,bj,np)=acclimave(i,j,k,bi,bj,np)+
1495         &                        accliml(np)*dtplankton
1496                 do ilam=1,tlam
1497                    Ek_nlave(i,j,k,bi,bj,np,ilam)=
1498         &                        Ek_nlave(i,j,k,bi,bj,np,ilam)+
1499         &                        Ek_nll(np,ilam)*dtplankton
1500                    EkoverE_nlave(i,j,k,bi,bj,np,ilam)=
1501         &                        EkoverE_nlave(i,j,k,bi,bj,np,ilam)+
1502         &                        EkoverE_nll(np,ilam)*dtplankton
1503                 enddo
1504                enddo
1505    #endif
1506  #ifdef DAR_DIAG_RSTAR  #ifdef DAR_DIAG_RSTAR
1507               do np=1,npmax               do np=1,npmax
1508                 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.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22