/[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.11 by stephd, Mon Dec 3 20:49:55 2007 UTC revision 1.18 by gforget, Thu Apr 24 21:35:53 2008 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"  
5    
6  CBOP  CBOP
7  C !ROUTINE: GCHEM_FORCING_INT  C !ROUTINE: DIC_INI_FORCING
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE DIC_INI_FORCING( myThid )        SUBROUTINE DIC_INI_FORCING( myThid )
# Line 21  C !USES: =============================== Line 20  C !USES: ===============================
20  #include "EEPARAMS.h"  #include "EEPARAMS.h"
21  #include "PARAMS.h"  #include "PARAMS.h"
22  #include "GRID.h"  #include "GRID.h"
23  c#include "PTRACERS_SIZE.h"  #include "DIC_VARS.h"
 c#include "PTRACERS_PARAMS.h"  
 c#include "PTRACERS_FIELDS.h"  
 #include "GCHEM.h"  
 #include "DIC_ABIOTIC.h"  
 #include "DIC_BIOTIC.h"  
24  #include "DIC_LOAD.h"  #include "DIC_LOAD.h"
25    
26  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
# Line 41  c !LOCAL VARIABLES: ==================== Line 35  c !LOCAL VARIABLES: ====================
35         CHARACTER*(MAX_LEN_MBUF) msgBuf         CHARACTER*(MAX_LEN_MBUF) msgBuf
36  c  c
37  C First call requires that we initialize everything to zero for safety  C First call requires that we initialize everything to zero for safety
38        CALL LEF_ZERO( wspeed0,myThid )        CALL LEF_ZERO( dicwind0,myThid )
39        CALL LEF_ZERO( wspeed1,myThid )        CALL LEF_ZERO( dicwind1,myThid )
40        CALL LEF_ZERO( atmosp0,myThid )        CALL LEF_ZERO( atmosp0,myThid )
41        CALL LEF_ZERO( atmosp1,myThid )        CALL LEF_ZERO( atmosp1,myThid )
42        CALL LEF_ZERO( silica0,myThid )        CALL LEF_ZERO( silica0,myThid )
# Line 58  C First call requires that we initialize Line 52  C First call requires that we initialize
52        CALL LEF_ZERO( feinput1,myThid )        CALL LEF_ZERO( feinput1,myThid )
53  # endif  # endif
54    
55         _EXCH_XY_R4(wspeed0, myThid )         _EXCH_XY_R4(dicwind0, myThid )
56         _EXCH_XY_R4(wspeed1, myThid )         _EXCH_XY_R4(dicwind1, myThid )
57         _EXCH_XY_R4(atmosp0, myThid )         _EXCH_XY_R4(atmosp0, myThid )
58         _EXCH_XY_R4(atmosp1, myThid )         _EXCH_XY_R4(atmosp1, myThid )
59         _EXCH_XY_R4(silica0, myThid )         _EXCH_XY_R4(silica0, myThid )
# Line 85  C First call requires that we initialize Line 79  C First call requires that we initialize
79  #endif  #endif
80  #endif  #endif
81    
           QSW_underice=.FALSE.  
 #ifdef USE_QSW_UNDERICE  
           QSW_underice=.TRUE.  
 #endif  
 #if (defined (ALLOW_SEAICE) && defined (USE_QSW))  
 c if using Qsw and seaice, then ice fraction is already  
 c taken into account  
           QSW_underice=.TRUE.  
 #endif  
 #if (defined (ALLOW_THSICE) && defined (USE_QSW))  
           QSW_underice=.TRUE.  
 #endif  
   
   
82  c set reasonable values to those that need at least something  c set reasonable values to those that need at least something
83         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
84          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22