/[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.13 by dimitri, Thu May 22 15:54:37 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    C     default to 100% ice concentration when thickness is non-zero
129             DO bj = myByLo(myThid), myByHi(myThid)
130              DO bi = myBxLo(myThid), myBxHi(myThid)
131               DO j=1-OLy,sNy+OLy
132                DO i=1-OLx,sNx+OLx
133                 IF (iceHeight(i,j,bi,bj) .NE. 0. _d 0)
134         &            iceMask(i,j,bi,bj) = 1. _d 0
135                ENDDO
136               ENDDO
137              ENDDO
138             ENDDO
139            ENDIF
140          IF ( thSIceSnowH_InitFile .NE. ' ' ) THEN          IF ( thSIceSnowH_InitFile .NE. ' ' ) THEN
141           CALL READ_REC_XY_RL(thSIceSnowH_InitFile,snowHeight,1,0,myThid)           CALL READ_REC_XY_RL(thSIceSnowH_InitFile,snowHeight,1,0,myThid)
142            ELSE
143    C     default to <= 20 cm snow thickness where there is ice
144             DO bj = myByLo(myThid), myByHi(myThid)
145              DO bi = myBxLo(myThid), myBxHi(myThid)
146               DO j=1-OLy,sNy+OLy
147                DO i=1-OLx,sNx+OLx
148                 snowHeight(i,j,bi,bj) = 0.2 _d 0
149                ENDDO
150               ENDDO
151              ENDDO
152             ENDDO
153          ENDIF          ENDIF
154          IF ( thSIceSnowA_InitFile .NE. ' ' ) THEN          IF ( thSIceSnowA_InitFile .NE. ' ' ) THEN
155           CALL READ_REC_XY_RL(thSIceSnowA_InitFile,snowAge,1,0,myThid)           CALL READ_REC_XY_RL(thSIceSnowA_InitFile,snowAge,1,0,myThid)
156            ELSE
157    C     default to 3 days for snow age
158             DO bj = myByLo(myThid), myByHi(myThid)
159              DO bi = myBxLo(myThid), myBxHi(myThid)
160               DO j=1-OLy,sNy+OLy
161                DO i=1-OLx,sNx+OLx
162                 IF (snowHeight(i,j,bi,bj) .NE. 0. _d 0)
163         &            snowAge(i,j,bi,bj) = 259200. _d 0
164                ENDDO
165               ENDDO
166              ENDDO
167             ENDDO
168          ENDIF          ENDIF
169          IF ( thSIceEnthp_InitFile .NE. ' ' ) THEN          IF ( thSIceEnthp_InitFile .NE. ' ' ) THEN
170           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice1,1,0,myThid)           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice1,1,0,myThid)
171           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice2,2,0,myThid)           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice2,2,0,myThid)
172            ELSE
173    C     default to 3.4e5 J/kg for enthalpy where there is sea ice
174             DO bj = myByLo(myThid), myByHi(myThid)
175              DO bi = myBxLo(myThid), myBxHi(myThid)
176               DO j=1-OLy,sNy+OLy
177                DO i=1-OLx,sNx+OLx
178                 IF (iceMask(i,j,bi,bj) .NE. 0. _d 0) THEN
179                  Qice1(i,j,bi,bj) = 3.4 _d 5
180                  Qice2(i,j,bi,bj) = 3.4 _d 5
181                 ENDIF
182                ENDDO
183               ENDDO
184              ENDDO
185             ENDDO
186          ENDIF          ENDIF
187          IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN          IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN
188           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.13
  ViewVC Help
Powered by ViewVC 1.1.22