/[MITgcm]/MITgcm/pkg/seaice/seaice_budget_ocean.F

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

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

-revision 1.1 by mlosch,
Thu Dec 14 08:36:20 2006 UTC
+revision 1.8 by mlosch,
Wed Jun 24 08:25:05 2009 UTC
 Line 8 
 CStartOfInterface
       I     UG,
       U     TSURF,
       O     netHeatFlux, SWHeatFlux,
-      I     bi, bj )
+      I     bi, bj, myTime, myIter, myThid )
  C     /================================================================\
  C     | SUBROUTINE seaice_budget_ocean                                 |
  C     | o Calculate surface heat fluxes over open ocean                |
 Line 24 
 C     === Global variables ===
  #include "EEPARAMS.h"
  #include "FFIELDS.h"
  #include "SEAICE_PARAMS.h"
- #include "SEAICE_FFIELDS.h"
  #ifdef SEAICE_VARIABLE_FREEZING_POINT
- #include "DYNVARS.h"
+ # include "DYNVARS.h"
- #endif /* SEAICE_VARIABLE_FREEZING_POINT */
+ #endif
+ #ifdef ALLOW_EXF
+ # include "EXF_OPTIONS.h"
+ # include "EXF_FIELDS.h"
+ #endif
+ #ifdef SEAICE_CLIM_AIR
+       COMMON/SEAICE_DYNVARS_1/AREA
+       _RL AREA       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,3,nSx,nSy)
+ #endif
  C     === Routine arguments ===
  C     INPUT:
  C     UG      :: thermal wind of atmosphere
  C     TSURF   :: surface temperature of ocean in Kelvin
  C     bi,bj   :: loop indices
+ C     myTime  :: Simulation time
+ C     myIter  :: Simulation timestep number
+ C     myThid  :: Thread no. that called this routine.
  C     OUTPUT:
  C     netHeatFlux :: net surface heat flux over open water or under ice
  C     SWHeatFlux  :: short wave heat flux over open water or under ice
-       _RL UG         (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL UG         (1:sNx,1:sNy)
        _RL TSURF      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
-       _RL netHeatFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL netHeatFlux(1:sNx,1:sNy)
-       _RL SWHeatFlux (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL SWHeatFlux (1:sNx,1:sNy)
-       INTEGER bi, bj
+       _RL myTime
+       INTEGER bi, bj, myIter, myThid
  CEndOfInterface
  C     === Local variables ===
-Line 49 
 C     i,j - Loop counters
+Line 60 
 C     i,j - Loop counters
        INTEGER i, j
  #ifndef SEAICE_EXTERNAL_FLUXES
        INTEGER ITER
-       _RL  QS1, C1, C2, C3, C4, C5, TB, D1, D1W, D1I, D3
+       _RL  QS1, TB, D1, D1W, D3, TMELT
-       _RL  TMELT, TMELTP, XKI, XKS, HCUT, ASNOW, XIO
  C     effective conductivity of combined ice and snow
        _RL  effConduct
  C     specific humidity at ice surface
-Line 59 
 C     powers of temperature
+Line 69 
 C     powers of temperature
        _RL  t1, t2, t3, t4
  C     local copies of global variables
-       _RL tsurfLoc   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL tsurfLoc   (1:sNx,1:sNy)
-       _RL atempLoc   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL atempLoc   (1:sNx,1:sNy)
-       _RL lwdownLoc  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL lwdownLoc  (1:sNx,1:sNy)
-       _RL ALB        (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL ALB        (1:sNx,1:sNy)
  C     coefficients of Hibler (1980), appendix B
-       _RL A1         (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL A1         (1:sNx,1:sNy)
-       _RL A2         (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL A2         (1:sNx,1:sNy)
-       _RL A3         (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
  C     auxiliary variable
-       _RL B          (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RL B          (1:sNx,1:sNy)
+ #ifdef SEAICE_CLIM_AIR
+       _RL aqhLoc     (1:sNx,1:sNy)
+       _RL fac
+       logical first, changed
+       integer count0, count1
+ C--   Compute indices and weights for seasonal interpolation
+       call cal_GetMonthsRec(
+      O     fac, first, changed,
+      O     count0, count1,
+      I     mytime, myiter, mythid
+      &     )
+ #endif /* SEAICE_CLIM_AIR */
  C NOW DEFINE ASSORTED CONSTANTS
  C SATURATION VAPOR PRESSURE CONSTANT
        QS1=0.622 _d +00/1013.0 _d +00
- C MAYKUTS CONSTANTS FOR SAT. VAP. PRESSURE TEMP. POLYNOMIAL
-       C1=2.7798202 _d -06
-       C2=-2.6913393 _d -03
-       C3=0.97920849 _d +00
-       C4=-158.63779 _d +00
-       C5=9653.1925 _d +00
  C FREEZING TEMPERATURE OF SEAWATER
        TB=271.2 _d +00
  C SENSIBLE HEAT CONSTANT
        D1=SEAICE_sensHeat
  C WATER LATENT HEAT CONSTANT
        D1W=SEAICE_latentWater
- C ICE LATENT HEAT CONSTANT
-       D1I=SEAICE_latentIce
  C STEFAN BOLTZMAN CONSTANT TIMES 0.97 EMISSIVITY
        D3=SEAICE_emissivity
  C MELTING TEMPERATURE OF ICE
        TMELT=273.16 _d +00
-       TMELTP=273.159 _d +00
- C ICE CONDUCTIVITY
-       XKI=SEAICE_iceConduct
- C SNOW CONDUCTIVITY
-       XKS=SEAICE_snowConduct
- C CUTOFF SNOW THICKNESS
-       HCUT=SEAICE_snowThick
- C PENETRATION SHORTWAVE RADIATION FACTOR
-       XIO=SEAICE_shortwave
        DO J=1,sNy
         DO I=1,sNx
-Line 107 
 C PENETRATION SHORTWAVE RADIATION FACTOR
+Line 113 
 C PENETRATION SHORTWAVE RADIATION FACTOR
          SWHeatFlux (I,J) = 0. _d 0
  C
          tsurfLoc (I,J) = MIN(273.16 _d 0+MAX_TICE,TSURF(I,J,bi,bj))
+ # ifdef ALLOW_ATM_TEMP
  C     Is this necessary?
          atempLoc (I,J) = MAX(273.16 _d 0+MIN_ATEMP,ATEMP(I,J,bi,bj))
+ # endif
+ #ifdef SEAICE_CLIM_AIR
+         atempLoc (I,J) = AREA(I,J,bi,bj) *
+      &       ( fac * SEAICE_clim_atemp(count0) + (1-fac) *
+      &       SEAICE_clim_atemp(count1) ) +
+      &       (1-AREA(I,J,bi,bj)) * atempLoc(I,J)
+         aqhLoc (I,J)   = AREA(I,J,bi,bj) *
+      &       ( fac * SEAICE_clim_aqh(count0) + (1-fac) *
+      &       SEAICE_clim_aqh(count1) ) +
+      &       (1-AREA(I,J,bi,bj)) * aqh(I,J,bi,bj)
+ #endif /* SEAICE_CLIM_AIR */
+ # ifdef ALLOW_DOWNWARD_RADIATION
          lwdownLoc(I,J) = MAX(MIN_LWDOWN,LWDOWN(I,J,bi,bj))
+ # endif
         ENDDO
        ENDDO
  #endif /* SEAICE_EXTERNAL_FLUXES */
-Line 123 
 C WATER ALBEDO IS ASSUMED TO BE THE CONS
+Line 143 
 C WATER ALBEDO IS ASSUMED TO BE THE CONS
          SWHeatFlux (I,J) =  Qsw(I,J,bi,bj)
  #else /* SEAICE_EXTERNAL_FLUXES undefined */
          ALB(I,J)=SEAICE_waterAlbedo
+ # ifdef ALLOW_DOWNWARD_RADIATION
+ # ifdef SEAICE_CLIM_AIR
+         A1(I,J)=(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj)
+      &       +lwdownLoc(I,J)*0.97 _d 0
+      &       +D1*UG(I,J)*atempLoc(I,J)+D1W*UG(I,J)*aqhLoc (I,J)
+ #else
          A1(I,J)=(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj)
       &       +lwdownLoc(I,J)*0.97 _d 0
       &       +D1*UG(I,J)*atempLoc(I,J)+D1W*UG(I,J)*AQH(I,J,bi,bj)
+ #endif /* SEAICE_CLIM_AIR */
          B(I,J)=QS1*6.11 _d +00*EXP(17.2694 _d +00
       &       *(tsurfLoc(I,J)-TMELT)
       &       /(tsurfLoc(I,J)-TMELT+237.3 _d +00))
-Line 133 
 C WATER ALBEDO IS ASSUMED TO BE THE CONS
+Line 160 
 C WATER ALBEDO IS ASSUMED TO BE THE CONS
       &       -D3*(tsurfLoc(I,J)**4)
          netHeatFlux(I,J)=-A1(I,J)-A2(I,J)
          SWHeatFlux (I,J)=-(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj)
+ # endif /* ALLOW_DOWNWARD_RADIATION */
  #endif /* SEAICE_EXTERNAL_FLUXES */
         ENDDO
        ENDDO

 Legend:



Removed from v.1.1
 


changed lines


 
Added in v.1.8
 Legend:



Removed from v.1.1
 


changed lines


 
Added in v.1.8
-Removed from v.1.1
+Added in v.1.8

	ViewVC Help
Powered by ViewVC 1.1.22