/[MITgcm]/MITgcm/pkg/seaice/seaice_solve4temp.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_solve4temp.F

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

revision 1.34 by gforget, Mon Mar 5 15:21:45 2012 UTC revision 1.35 by jmc, Tue Mar 6 01:28:11 2012 UTC
# Line 217  c     SurfMeltTemp = 273.159 _d +00 Line 217  c     SurfMeltTemp = 273.159 _d +00
217  C     Initialize variables  C     Initialize variables
218        DO J=1,sNy        DO J=1,sNy
219         DO I=1,sNx         DO I=1,sNx
220    C     initialise output arrays:
221            TSURFout (I,J) = TSURFin(I,J)
222            F_ia     (I,J) = 0. _d 0
223            IcePenetSW(I,J)= 0. _d 0
224            FWsublim (I,J) = 0. _d 0
225  C     HICE_ACTUAL is modified in this routine, but at the same time  C     HICE_ACTUAL is modified in this routine, but at the same time
226  C     used to decided where there is ice, therefore we save this information  C     used to decided where there is ice, therefore we save this information
227  C     here in a separate array  C     here in a separate array
228          iceOrNot  (I,J) = HICE_ACTUAL(I,J) .GT. 0. _d 0          iceOrNot  (I,J) = HICE_ACTUAL(I,J) .GT. 0. _d 0
         IcePenetSW(I,J) = 0. _d 0  
229          absorbedSW(I,J) = 0. _d 0          absorbedSW(I,J) = 0. _d 0
230          qhice    (I,J) = 0. _d 0          qhice    (I,J) = 0. _d 0
231          dqh_dTs  (I,J) = 0. _d 0          dqh_dTs  (I,J) = 0. _d 0
         F_ia     (I,J) = 0. _d 0  
232          F_lh     (I,J) = 0. _d 0          F_lh     (I,J) = 0. _d 0
233          F_lwu    (I,J) = 0. _d 0          F_lwu    (I,J) = 0. _d 0
234          F_sens   (I,J) = 0. _d 0          F_sens   (I,J) = 0. _d 0
235  C     Make a local copy of LW, surface & atmospheric temperatures  C     Make a local copy of LW, surface & atmospheric temperatures
236          tsurfLoc (I,J) = TSURFin(I,J)          tsurfLoc (I,J) = TSURFin(I,J)
         TSURFout (I,J) = TSURFin(I,J)  
237  c       tsurfLoc (I,J) = MIN( celsius2K+MAX_TICE, TSURFin(I,J) )  c       tsurfLoc (I,J) = MIN( celsius2K+MAX_TICE, TSURFin(I,J) )
238          lwdownLoc(I,J) = MAX( MIN_LWDOWN, LWDOWN(I,J,bi,bj) )          lwdownLoc(I,J) = MAX( MIN_LWDOWN, LWDOWN(I,J,bi,bj) )
239          atempLoc (I,J) = MAX( celsius2K+MIN_ATEMP, ATEMP(I,J,bi,bj) )          atempLoc (I,J) = MAX( celsius2K+MIN_ATEMP, ATEMP(I,J,bi,bj) )
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22