/[MITgcm]/MITgcm/pkg/dic/dic_ini_forcing.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_ini_forcing.F

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

revision 1.7 by stephd, Mon May 7 15:58:19 2007 UTC revision 1.11 by stephd, Mon Dec 3 20:49:55 2007 UTC
# Line 21  C !USES: =============================== Line 21  C !USES: ===============================
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  #include "GRID.h"  #include "GRID.h"
24  #include "PTRACERS_SIZE.h"  c#include "PTRACERS_SIZE.h"
25  #include "PTRACERS.h"  c#include "PTRACERS_PARAMS.h"
26    c#include "PTRACERS_FIELDS.h"
27  #include "GCHEM.h"  #include "GCHEM.h"
28  #include "DIC_ABIOTIC.h"  #include "DIC_ABIOTIC.h"
29  #include "DIC_BIOTIC.h"  #include "DIC_BIOTIC.h"
# Line 84  C First call requires that we initialize Line 85  C First call requires that we initialize
85  #endif  #endif
86  #endif  #endif
87    
88              QSW_underice=.FALSE.
89    #ifdef USE_QSW_UNDERICE
90              QSW_underice=.TRUE.
91    #endif
92    #if (defined (ALLOW_SEAICE) && defined (USE_QSW))
93    c if using Qsw and seaice, then ice fraction is already
94    c taken into account
95              QSW_underice=.TRUE.
96    #endif
97    #if (defined (ALLOW_THSICE) && defined (USE_QSW))
98              QSW_underice=.TRUE.
99    #endif
100    
101    
102  c set reasonable values to those that need at least something  c set reasonable values to those that need at least something
103         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
104          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
105           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
106            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
107               WIND(i,j,bi,bj) = 5.d0*maskC(i,j,1,bi,bj)               WIND(i,j,bi,bj)   = 5. _d 0*maskC(i,j,1,bi,bj)
108               ATMOSP(i,j,bi,bj)   =1.d0*maskC(i,j,1,bi,bj)               AtmosP(i,j,bi,bj) = 1. _d 0*maskC(i,j,1,bi,bj)
109               SILICA(i,j,bi,bj)   =7.6838e-3*maskC(i,j,1,bi,bj)               SILICA(i,j,bi,bj) = 7.6838 _d -3*maskC(i,j,1,bi,bj)
110               FIce(i,j,bi,bj) =0.d0               FIce(i,j,bi,bj)   = 0. _d 0
111  #ifdef READ_PAR  #ifdef READ_PAR
112               PAR(i,j,bi,bj) =100.d0*maskC(i,j,1,bi,bj)               PAR(i,j,bi,bj)    = 100. _d 0*maskC(i,j,1,bi,bj)
113  #endif  #endif
114  #ifdef ALLOW_FE  #ifdef ALLOW_FE
115               InputFe(i,j,bi,bj)  = 1e-11*maskC(i,j,1,bi,bj)               InputFe(i,j,bi,bj)= 1. _d -11*maskC(i,j,1,bi,bj)
116  #endif  #endif
117            ENDDO            ENDDO
118           ENDDO           ENDDO
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22