/[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.12 by dimitri, Thu May 22 01:59:51 2008 UTC revision 1.13 by dimitri, Thu May 22 15:54:37 2008 UTC
# Line 125  C--     Read initial conditions: Line 125  C--     Read initial conditions:
125          IF ( thSIceFract_InitFile .NE. ' ' ) THEN          IF ( thSIceFract_InitFile .NE. ' ' ) THEN
126           CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid)           CALL READ_REC_XY_RL(thSIceFract_InitFile,iceMask,1,0,myThid)
127          ELSE          ELSE
128    C     default to 100% ice concentration when thickness is non-zero
129           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
130            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
131             DO j=1-OLy,sNy+OLy             DO j=1-OLy,sNy+OLy
# Line 139  C--     Read initial conditions: Line 140  C--     Read initial conditions:
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          ELSE
143    C     default to <= 20 cm snow thickness where there is ice
144           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
145            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
146             DO j=1-OLy,sNy+OLy             DO j=1-OLy,sNy+OLy
147              DO i=1-OLx,sNx+OLx              DO i=1-OLx,sNx+OLx
148               snowHeight(i,j,bi,bj) = 0.2 _d 0 * iceMask(i,j,bi,bj)               snowHeight(i,j,bi,bj) = 0.2 _d 0
149              ENDDO              ENDDO
150             ENDDO             ENDDO
151            ENDDO            ENDDO
# Line 152  C--     Read initial conditions: Line 154  C--     Read initial conditions:
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          ELSE
157    C     default to 3 days for snow age
158           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
159            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
160             DO j=1-OLy,sNy+OLy             DO j=1-OLy,sNy+OLy
# Line 167  C--     Read initial conditions: Line 170  C--     Read initial conditions:
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          ELSE
173    C     default to 3.4e5 J/kg for enthalpy where there is sea ice
174           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
175            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
176             DO j=1-OLy,sNy+OLy             DO j=1-OLy,sNy+OLy
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22