/[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.11 by jmc, Thu Aug 30 15:02:51 2007 UTC revision 1.12 by dimitri, Thu May 22 01:59:51 2008 UTC
# Line 26  C     === Global variables === Line 26  C     === Global variables ===
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
27  #include "PARAMS.h"  #include "PARAMS.h"
28  #include "FFIELDS.h"  #include "FFIELDS.h"
29    #include "DYNVARS.h"
30  #include "THSICE_PARAMS.h"  #include "THSICE_PARAMS.h"
31  #include "THSICE_VARS.h"  #include "THSICE_VARS.h"
32  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
# Line 118  C--     Read ice pickup fields Line 119  C--     Read ice pickup fields
119        ELSE        ELSE
120  C--     Read initial conditions:  C--     Read initial conditions:
121          _BEGIN_MASTER( myThid )          _BEGIN_MASTER( myThid )
         IF ( thSIceFract_InitFile .NE. ' ' ) THEN  
          CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid)  
         ENDIF  
122          IF ( thSIceThick_InitFile .NE. ' ' ) THEN          IF ( thSIceThick_InitFile .NE. ' ' ) THEN
123           CALL READ_REC_XY_RL(thSIceThick_InitFile,iceHeight,1,0,myThid)           CALL READ_REC_XY_RL(thSIceThick_InitFile,iceHeight,1,0,myThid)
124          ENDIF          ENDIF
125            IF ( thSIceFract_InitFile .NE. ' ' ) THEN
126             CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid)
127            ELSE
128             DO bj = myByLo(myThid), myByHi(myThid)
129              DO bi = myBxLo(myThid), myBxHi(myThid)
130               DO j=1-OLy,sNy+OLy
131                DO i=1-OLx,sNx+OLx
132                 IF (iceHeight(i,j,bi,bj) .NE. 0. _d 0)
133         &            iceMask(i,j,bi,bj) = 1. _d 0
134                ENDDO
135               ENDDO
136              ENDDO
137             ENDDO
138            ENDIF
139          IF ( thSIceSnowH_InitFile .NE. ' ' ) THEN          IF ( thSIceSnowH_InitFile .NE. ' ' ) THEN
140           CALL READ_REC_XY_RL(thSIceSnowH_InitFile,snowHeight,1,0,myThid)           CALL READ_REC_XY_RL(thSIceSnowH_InitFile,snowHeight,1,0,myThid)
141            ELSE
142             DO bj = myByLo(myThid), myByHi(myThid)
143              DO bi = myBxLo(myThid), myBxHi(myThid)
144               DO j=1-OLy,sNy+OLy
145                DO i=1-OLx,sNx+OLx
146                 snowHeight(i,j,bi,bj) = 0.2 _d 0 * iceMask(i,j,bi,bj)
147                ENDDO
148               ENDDO
149              ENDDO
150             ENDDO
151          ENDIF          ENDIF
152          IF ( thSIceSnowA_InitFile .NE. ' ' ) THEN          IF ( thSIceSnowA_InitFile .NE. ' ' ) THEN
153           CALL READ_REC_XY_RL(thSIceSnowA_InitFile,snowAge,1,0,myThid)           CALL READ_REC_XY_RL(thSIceSnowA_InitFile,snowAge,1,0,myThid)
154            ELSE
155             DO bj = myByLo(myThid), myByHi(myThid)
156              DO bi = myBxLo(myThid), myBxHi(myThid)
157               DO j=1-OLy,sNy+OLy
158                DO i=1-OLx,sNx+OLx
159                 IF (snowHeight(i,j,bi,bj) .NE. 0. _d 0)
160         &            snowAge(i,j,bi,bj) = 259200. _d 0
161                ENDDO
162               ENDDO
163              ENDDO
164             ENDDO
165          ENDIF          ENDIF
166          IF ( thSIceEnthp_InitFile .NE. ' ' ) THEN          IF ( thSIceEnthp_InitFile .NE. ' ' ) THEN
167           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice1,1,0,myThid)           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice1,1,0,myThid)
168           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice2,2,0,myThid)           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice2,2,0,myThid)
169            ELSE
170             DO bj = myByLo(myThid), myByHi(myThid)
171              DO bi = myBxLo(myThid), myBxHi(myThid)
172               DO j=1-OLy,sNy+OLy
173                DO i=1-OLx,sNx+OLx
174                 IF (iceMask(i,j,bi,bj) .NE. 0. _d 0) THEN
175                  Qice1(i,j,bi,bj) = 3.4 _d 5
176                  Qice2(i,j,bi,bj) = 3.4 _d 5
177                 ENDIF
178                ENDDO
179               ENDDO
180              ENDDO
181             ENDDO
182          ENDIF          ENDIF
183          IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN          IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN
184           CALL READ_REC_XY_RS(thSIceTsurf_InitFile,Tsrf,1,0,myThid)           CALL READ_REC_XY_RS(thSIceTsurf_InitFile,Tsrf,1,0,myThid)
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22