/[MITgcm]/MITgcm/pkg/thsice/thsice_ini_vars.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_ini_vars.F

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

revision 1.7 by jmc, Sun Jun 25 22:35:17 2006 UTC revision 1.8 by jmc, Wed Apr 4 01:55:45 2007 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "THSICE_OPTIONS.h"  #include "THSICE_OPTIONS.h"
5    #ifdef ALLOW_AIM
6    # include "AIM_OPTIONS.h"
7    #endif
8    
9  CBOP  CBOP
10  C     !ROUTINE: THSICE_INI_VARS  C     !ROUTINE: THSICE_INI_VARS
# Line 22  C     === Global variables === Line 25  C     === Global variables ===
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
27  #include "PARAMS.h"  #include "PARAMS.h"
 #include "GRID.h"  
28  #include "FFIELDS.h"  #include "FFIELDS.h"
29  #include "THSICE_PARAMS.h"  #include "THSICE_PARAMS.h"
30  #include "THSICE_VARS.h"  #include "THSICE_VARS.h"
31    #ifdef ALLOW_AIM
32    # include "AIM_FFIELDS.h"
33    #endif
34    
35  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
36  C     == Routine arguments ==  C     == Routine arguments ==
37  C     myThid -  Number of this instance of INI_RBCS  C     myThid :: My Thread Id. number
38        INTEGER myThid        INTEGER myThid
39  CEOP  CEOP
40    
# Line 38  C     == Local variables == Line 43  C     == Local variables ==
43  C     bi,bj  :: Loop counters  C     bi,bj  :: Loop counters
44  C     i,j    :: Loop counters  C     i,j    :: Loop counters
45        INTEGER bi, bj        INTEGER bi, bj
46        INTEGER I, J        INTEGER i, j
47        INTEGER prec        INTEGER prec
48  c     CHARACTER*(MAX_LEN_FNAM) fn  c     CHARACTER*(MAX_LEN_FNAM) fn
49        _RL v2Loc        _RL v2Loc
50    
51    
52  c     set up ice arrays to zero if starting ice  c     set up ice arrays to zero if starting ice
53          DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
54           DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
55  C-        state variables :  C-        state variables :
56            DO j=1-OLy,sNy+OLy            DO j=1-OLy,sNy+OLy
57             DO i=1-OLx,sNx+OLx             DO i=1-OLx,sNx+OLx
# Line 68  C-        fluxes : Line 73  C-        fluxes :
73              flxCndBt(i,j,bi,bj) = 0. _d 0              flxCndBt(i,j,bi,bj) = 0. _d 0
74              snowPrc(i,j,bi,bj)  = 0. _d 0              snowPrc(i,j,bi,bj)  = 0. _d 0
75              siceAlb(i,j,bi,bj)  = 0. _d 0              siceAlb(i,j,bi,bj)  = 0. _d 0
 #ifdef COUPLE_MODEL  
             dFdT(i,j,bi,bj)     = 0. _d 0  
 #endif  
76              icFlxSW (i,j,bi,bj) = 0. _d 0              icFlxSW (i,j,bi,bj) = 0. _d 0
77              icFlxAtm(i,j,bi,bj) = 0. _d 0              icFlxAtm(i,j,bi,bj) = 0. _d 0
78              icFrwAtm(i,j,bi,bj) = 0. _d 0              icFrwAtm(i,j,bi,bj) = 0. _d 0
# Line 86  C-        oceanic mixed layer state : Line 88  C-        oceanic mixed layer state :
88               v2ocMxL(i,j,bi,bj) = v2Loc               v2ocMxL(i,j,bi,bj) = v2Loc
89             ENDDO             ENDDO
90            ENDDO            ENDDO
91           ENDDO  #ifdef ALLOW_AIM
92              IF ( useAIM ) THEN
93    C-        Mask mixed layer depth : depth is used in thsice slab_ocean
94    C         and this mask is used in thsice_advdiff and if coupled
95               DO j=1-OLy,sNy+OLy
96                DO i=1-OLx,sNx+OLx
97                 IF ( aim_landFr(i,j,bi,bj).EQ.1. _d 0 )
98         &       hOceMxL(i,j,bi,bj) = 0.
99                ENDDO
100               ENDDO
101              ENDIF
102    #endif /* ALLOW_AIM */
103          ENDDO          ENDDO
104          ENDDO
105    
106  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
107    
# Line 131  C--     Read initial conditions: Line 145  C--     Read initial conditions:
145          _EXCH_XY_R8(Qice1,  myThid)          _EXCH_XY_R8(Qice1,  myThid)
146          _EXCH_XY_R8(Qice2,  myThid)          _EXCH_XY_R8(Qice2,  myThid)
147          _EXCH_XY_R8(snowAge,myThid)          _EXCH_XY_R8(snowAge,myThid)
 #ifdef COUPLE_MODEL  
         _EXCH_XY_R8(dFdT, myThid)  
 #endif  
148    
149  #ifdef ATMOSPHERIC_LOADING  #ifdef ATMOSPHERIC_LOADING
150  C--     Initialise Sea-Ice Loading for SeaIce-Dynamics :  C--     Initialise Sea-Ice Loading for SeaIce-Dynamics :
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22