/[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.10 by jmc, Fri Apr 27 15:50:40 2007 UTC revision 1.14 by dimitri, Fri May 23 09:56:25 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"
 #ifdef ALLOW_AUTODIFF_TAMC  
 # include "THSICE_2DYN.h"  
 #endif  
32  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
33  # include "AIM_FFIELDS.h"  # include "AIM_FFIELDS.h"
34  #endif  #endif
# Line 79  C-        fluxes : Line 77  C-        fluxes :
77              icFlxSW (i,j,bi,bj) = 0. _d 0              icFlxSW (i,j,bi,bj) = 0. _d 0
78              icFlxAtm(i,j,bi,bj) = 0. _d 0              icFlxAtm(i,j,bi,bj) = 0. _d 0
79              icFrwAtm(i,j,bi,bj) = 0. _d 0              icFrwAtm(i,j,bi,bj) = 0. _d 0
80  #ifdef ALLOW_AUTODIFF_TAMC  C-        needed when using advection/diffusion:
81              oceFWfx(i,j,bi,bj)  = 0. _d 0              oceFWfx(i,j,bi,bj)  = 0. _d 0
82              oceSflx(i,j,bi,bj)  = 0. _d 0              oceSflx(i,j,bi,bj)  = 0. _d 0
83              oceQnet(i,j,bi,bj)  = 0. _d 0              oceQnet(i,j,bi,bj)  = 0. _d 0
 #endif  
84             ENDDO             ENDDO
85            ENDDO            ENDDO
86  C-        oceanic mixed layer state :  C-        oceanic mixed layer state :
# Line 122  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                 IF (iceMask(i,j,bi,bj) .NE. 0. _d 0)
149         &            snowHeight(i,j,bi,bj) = 0.2 _d 0
150                ENDDO
151               ENDDO
152              ENDDO
153             ENDDO
154          ENDIF          ENDIF
155          IF ( thSIceSnowA_InitFile .NE. ' ' ) THEN          IF ( thSIceSnowA_InitFile .NE. ' ' ) THEN
156           CALL READ_REC_XY_RL(thSIceSnowA_InitFile,snowAge,1,0,myThid)           CALL READ_REC_XY_RL(thSIceSnowA_InitFile,snowAge,1,0,myThid)
157            ELSE
158    C     default to 3 days for snow age
159             DO bj = myByLo(myThid), myByHi(myThid)
160              DO bi = myBxLo(myThid), myBxHi(myThid)
161               DO j=1-OLy,sNy+OLy
162                DO i=1-OLx,sNx+OLx
163                 IF (snowHeight(i,j,bi,bj) .NE. 0. _d 0)
164         &            snowAge(i,j,bi,bj) = 259200. _d 0
165                ENDDO
166               ENDDO
167              ENDDO
168             ENDDO
169          ENDIF          ENDIF
170          IF ( thSIceEnthp_InitFile .NE. ' ' ) THEN          IF ( thSIceEnthp_InitFile .NE. ' ' ) THEN
171           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice1,1,0,myThid)           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice1,1,0,myThid)
172           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice2,2,0,myThid)           CALL READ_REC_XY_RL(thSIceEnthp_InitFile,Qice2,2,0,myThid)
173            ELSE
174    C     default to 3.4e5 J/kg for enthalpy where there is sea ice
175             DO bj = myByLo(myThid), myByHi(myThid)
176              DO bi = myBxLo(myThid), myBxHi(myThid)
177               DO j=1-OLy,sNy+OLy
178                DO i=1-OLx,sNx+OLx
179                 IF (iceMask(i,j,bi,bj) .NE. 0. _d 0) THEN
180                  Qice1(i,j,bi,bj) = 3.4 _d 5
181                  Qice2(i,j,bi,bj) = 3.4 _d 5
182                 ENDIF
183                ENDDO
184               ENDDO
185              ENDDO
186             ENDDO
187          ENDIF          ENDIF
188          IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN          IF ( thSIceTsurf_InitFile .NE. ' ' ) THEN
189           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.10  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22