/[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.3 by mlosch,
Mon Apr 16 22:43:01 2007 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, myThid )
+      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, myThid
+       _RL myTime
+       INTEGER bi, bj, myIter, myThid
  CEndOfInterface
  C     === Local variables ===
-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
-Line 90 
 C MELTING TEMPERATURE OF ICE
+Line 113 
 C MELTING TEMPERATURE OF ICE
          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 106 
 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 116 
 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.3
 


changed lines


 
Added in v.1.8
 Legend:



Removed from v.1.3
 


changed lines


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

	ViewVC Help
Powered by ViewVC 1.1.22