/[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.3 by jmc, Tue Jul 13 18:03:31 2004 UTC revision 1.5 by stephd, Tue May 1 22:12:12 2007 UTC
# Line 12  C !INTERFACE: ========================== Line 12  C !INTERFACE: ==========================
12    
13  C !DESCRIPTION:  C !DESCRIPTION:
14  c initializes forcing fields to zero  c initializes forcing fields to zero
15    c (or to reasonable values)
16    
17  C !USES: ===============================================================  C !USES: ===============================================================
18        IMPLICIT NONE        IMPLICIT NONE
# Line 24  C !USES: =============================== Line 25  C !USES: ===============================
25  #include "PTRACERS.h"  #include "PTRACERS.h"
26  #include "GCHEM.h"  #include "GCHEM.h"
27  #include "DIC_ABIOTIC.h"  #include "DIC_ABIOTIC.h"
 #ifdef DIC_BIOTIC  
28  #include "DIC_BIOTIC.h"  #include "DIC_BIOTIC.h"
29  #include "DIC_LOAD.h"  #include "DIC_LOAD.h"
 #endif  
30    
31  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
32  C  myThid               :: thread number  C  myThid               :: thread number
# Line 35  C  myThid               :: thread number Line 34  C  myThid               :: thread number
34  CEOP  CEOP
35    
36  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
37    
38    c !LOCAL VARIABLES: ===================================================
39           INTEGER bi,bj,i,j
40  c  c
41  C First call requires that we initialize everything to zero for safety  C First call requires that we initialize everything to zero for safety
42        CALL LEF_ZERO( wspeed0,myThid )        CALL LEF_ZERO( wspeed0,myThid )
# Line 45  C First call requires that we initialize Line 47  C First call requires that we initialize
47        CALL LEF_ZERO( silica1,myThid )        CALL LEF_ZERO( silica1,myThid )
48        CALL LEF_ZERO( ice0,myThid )        CALL LEF_ZERO( ice0,myThid )
49        CALL LEF_ZERO( ice1,myThid )        CALL LEF_ZERO( ice1,myThid )
50    #ifdef READ_PAR
51          CALL LEF_ZERO( par0,myThid )
52          CALL LEF_ZERO( par1,myThid )
53    #endif
54  # ifdef ALLOW_FE  # ifdef ALLOW_FE
55        CALL LEF_ZERO( feinput0,myThid )        CALL LEF_ZERO( feinput0,myThid )
56        CALL LEF_ZERO( feinput1,myThid )        CALL LEF_ZERO( feinput1,myThid )
# Line 58  C First call requires that we initialize Line 64  C First call requires that we initialize
64         _EXCH_XY_R4(silica1, myThid )         _EXCH_XY_R4(silica1, myThid )
65         _EXCH_XY_R4(ice0, myThid )         _EXCH_XY_R4(ice0, myThid )
66         _EXCH_XY_R4(ice1, myThid )         _EXCH_XY_R4(ice1, myThid )
67    #ifdef READ_PAR
68           _EXCH_XY_R4(par0, myThid )
69           _EXCH_XY_R4(par1, myThid )
70    #endif
71  # ifdef ALLOW_FE  # ifdef ALLOW_FE
72         _EXCH_XY_R4(feinput0, myThid )         _EXCH_XY_R4(feinput0, myThid )
73         _EXCH_XY_R4(feinput1, myThid )         _EXCH_XY_R4(feinput1, myThid )
74  # endif  # endif
75    
76    c set reasonable values to those that need at least something
77           DO bj = myByLo(myThid), myByHi(myThid)
78            DO bi = myBxLo(myThid), myBxHi(myThid)
79             DO j=1-Oly,sNy+Oly
80              DO i=1-Olx,sNx+Olx
81                 WIND(i,j,bi,bj) = 5.d0*maskC(i,j,1,bi,bj)
82                 ATMOSP(i,j,bi,bj)   =1.d0*maskC(i,j,1,bi,bj)
83                 SILICA(i,j,bi,bj)   =7.6838e-3*maskC(i,j,1,bi,bj)
84                 FIce(i,j,bi,bj) =0.d0
85    #ifdef READ_PAR
86                 PAR(i,j,bi,bj) =100.d0*maskC(i,j,1,bi,bj)
87    #endif
88    #ifdef ALLOW_FE
89                 InputFe(i,j,bi,bj)  = 1e-11*maskC(i,j,1,bi,bj)
90    #endif
91              ENDDO
92             ENDDO
93            ENDDO
94           ENDDO
95    
96    
97  #endif  #endif
98        RETURN        RETURN
99        END        END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22