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

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

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

revision 1.11 by jmc, Wed Dec 28 19:18:54 2011 UTC revision 1.12 by jmc, Fri Jan 13 21:11:30 2012 UTC
# Line 34  C     i,j,k,bi,bj - Loop counters Line 34  C     i,j,k,bi,bj - Loop counters
34  #ifndef SEAICE_CGRID  #ifndef SEAICE_CGRID
35        _RS  mask_uice        _RS  mask_uice
36  #endif  #endif
37  cif(  #ifdef SHORTWAVE_HEATING
38  cif   Helper variable for determining the fraction of sw radiation  cif   Helper variable for determining the fraction of sw radiation
39  cif   penetrating the model shallowest layer  cif   penetrating the model shallowest layer
       INTEGER dummyIter  
40        _RL dummyTime        _RL dummyTime
41        _RL swfracba(2)        _RL swfracba(2)
42        _RL FACTORM        _RL tmpFac
43        INTEGER IMAX  #endif /* SHORTWAVE_HEATING */
 cif)  
44    
45        IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN        IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
46         kSurface        = Nr         kSurface        = Nr
# Line 57  C     Initialize MNC variable informatio Line 55  C     Initialize MNC variable informatio
55          CALL SEAICE_MNC_INIT( myThid )          CALL SEAICE_MNC_INIT( myThid )
56        ENDIF        ENDIF
57    
58  cif(        _BEGIN_MASTER(myThid)
59  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
60        IMAX      = 2         tmpFac    = -1.0
61        FACTORM   = -1.0         dummyTime = 1.0
62        dummyTime = 1.0         swfracba(1) = ABS(rF(1))
63        dummyIter = 0         swfracba(2) = ABS(rF(2))
64        swfracba(1) = abs(rF(1))         CALL SWFRAC(
65        swfracba(2) = abs(rF(2))       I       2, tmpFac,
       CALL SWFRAC(  
      I       IMAX,FACTORM,  
66       U       swfracba,       U       swfracba,
67       I       dummyTime,dummyIter,myThid)       I       dummyTime, 0, myThid )
68        SWFRACB = swfracba(2)         SWFracB = swfracba(2)
69  #endif  #else /* SHORTWAVE_HEATING */
70  cif)         SWFracB = 0. _d 0
71    #endif /* SHORTWAVE_HEATING */
72          _END_MASTER(myThid)
73    
74  C--   Initialize grid info  C--   Initialize grid info
75        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22