Back to home page

MITgcm
 
 

    

File indexing completed on 2018-03-02 18:37:22 UTC

view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
65007c221b Jean*0001 #include "AIM_OPTIONS.h"
                0002 
2a80e4d00e Jean*0003 CBOP
                0004 C     !ROUTINE: AIM_SURF_BC
                0005 C     !INTERFACE:
d0a9461855 Jean*0006       SUBROUTINE AIM_SURF_BC(
                0007      U                        tYear,
                0008      O                        aim_sWght0, aim_sWght1,
                0009      I                        bi, bj, myTime, myIter, myThid )
2a80e4d00e Jean*0010 
                0011 C     !DESCRIPTION: \bv
65007c221b Jean*0012 C     *================================================================*
                0013 C     | S/R AIM_SURF_BC
d0a9461855 Jean*0014 C     | Set surface Boundary Conditions
65007c221b Jean*0015 C     |  for the atmospheric physics package
                0016 C     *================================================================*
                0017 c     | was part of S/R FORDATE in Franco Molteni SPEEDY code (ver23).
                0018 C     | For now, surface BC are loaded from files (S/R AIM_FIELDS_LOAD)
                0019 C     |  and imposed (= surf. forcing).
d0a9461855 Jean*0020 C     | In the future, will add
65007c221b Jean*0021 C     |  a land model and a coupling interface with an ocean GCM
                0022 C     *================================================================*
2a80e4d00e Jean*0023 C     \ev
                0024 
                0025 C     !USES:
65007c221b Jean*0026       IMPLICIT NONE
                0027 
                0028 C     -------------- Global variables --------------
                0029 C-- size for MITgcm & Physics package :
                0030 #include "AIM_SIZE.h"
                0031 
                0032 C-- MITgcm
                0033 #include "EEPARAMS.h"
                0034 #include "PARAMS.h"
72b4e12f25 Dani*0035 #include "GRID.h"
6fb2765f24 Jean*0036 c #include "DYNVARS.h"
65007c221b Jean*0037 c #include "SURFACE.h"
                0038 
                0039 C-- Physics package
a27dd2281d Jean*0040 #include "AIM_PARAMS.h"
65007c221b Jean*0041 #include "AIM_FFIELDS.h"
                0042 c #include "AIM_GRID.h"
                0043 #include "com_forcon.h"
                0044 #include "com_forcing.h"
                0045 c #include "com_physvar.h"
4c105e462e Davi*0046 #include "AIM_CO2.h"
65007c221b Jean*0047 
e19eef4074 Jean*0048 C-- Coupled to the Ocean :
                0049 #ifdef COMPONENT_MODULE
                0050 #include "CPL_PARAMS.h"
                0051 #include "ATMCPL.h"
                0052 #endif
                0053 
2a80e4d00e Jean*0054 C     !INPUT/OUTPUT PARAMETERS:
65007c221b Jean*0055 C     == Routine arguments ==
d0a9461855 Jean*0056 C     tYear      :: Fraction into year
                0057 C     aim_sWght0 :: weight for time interpolation of surface BC
                0058 C     aim_sWght1 :: 0/1 = time period before/after the current time
                0059 C     bi,bj      :: Tile indices
                0060 C     myTime     :: Current time of simulation ( s )
                0061 C     myIter     :: Current iteration number in simulation
                0062 C     myThid     :: my Thread number Id.
                0063       _RL     tYear
                0064       _RL     aim_sWght0, aim_sWght1
                0065       INTEGER bi, bj
                0066       _RL     myTime
                0067       INTEGER myIter, myThid
2a80e4d00e Jean*0068 CEOP
65007c221b Jean*0069 
                0070 #ifdef ALLOW_AIM
2a80e4d00e Jean*0071 C     !FUNCTIONS:
                0072 C     !LOCAL VARIABLES:
                0073 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0074 C--   Local Variables originally (Speedy) in common bloc (com_forcing.h):
                0075 C--   COMMON /FORDAY/ Daily forcing fields (updated in FORDATE)
                0076 C     oice1      :: sea ice fraction
                0077 C     snow1      :: snow depth (mm water)
                0078       _RL     oice1(NGP)
                0079       _RL     snow1(NGP)
                0080 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
65007c221b Jean*0081 C     == Local variables ==
d0a9461855 Jean*0082 C     i,j,k,I2,k   :: Loop counters
                0083       INTEGER i,j,I2,k, nm0
                0084       _RL t0prd, tNcyc, tmprd, dTprd
65007c221b Jean*0085       _RL SDEP1, IDEP2, SDEP2, SWWIL2, RSW, soilw_0, soilw_1
72b4e12f25 Dani*0086       _RL RSD, alb_land, oceTfreez, ALBSEA1, ALPHA, CZEN, CZEN2
972eaa29a0 Dani*0087       _RL RZEN, ZS, ZC, SJ, CJ, TMPA, TMPB, TMPL, hlim
b3097ed02d Jean*0088 c     _RL DALB, alb_sea
2a80e4d00e Jean*0089 #ifdef ALLOW_AIM_CO2
870deba1f6 Davi*0090 #ifdef ALLOW_DIAGNOSTICS
fb7fcd6036 Davi*0091       _RL pCO2scl
870deba1f6 Davi*0092 #endif
2a80e4d00e Jean*0093 #endif /* ALLOW_AIM_CO2 */
65007c221b Jean*0094 
4b996cbf32 Jean*0095 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0096 C-    Set Land-sea mask (in [0,1]) from aim_landFr to fMask1:
                0097       DO j=1,sNy
                0098         DO i=1,sNx
                0099           I2 = i+(j-1)*sNx
                0100           fMask1(I2,1,myThid) = aim_landFr(i,j,bi,bj)
                0101         ENDDO
                0102       ENDDO
                0103 
65007c221b Jean*0104       IF (aim_useFMsurfBC) THEN
                0105 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0106 
d0a9461855 Jean*0107 C     aim_surfForc_TimePeriod :: Length of forcing time period (e.g. 1 month)
                0108 C     aim_surfForc_NppCycle   :: Number of time period per Cycle (e.g. 12)
                0109 C     aim_surfForc_TransRatio ::
                0110 C-     define how fast the (linear) transition is from one month to the next
                0111 C       = 1                 -> linear between 2 midle month
                0112 C       > TimePeriod/deltaT -> jump from one month to the next one
                0113 
                0114 C--   Calculate weight for linear interpolation between 2 month centers
                0115         t0prd = myTime / aim_surfForc_TimePeriod
                0116         tNcyc = aim_surfForc_NppCycle
                0117         tmprd = t0prd - 0.5 _d 0 + tNcyc
                0118         tmprd = MOD(tmprd,tNcyc)
                0119 C-     indices of previous month (nm0) and next month (nm1):
                0120         nm0 = 1 + INT(tmprd)
                0121 c       nm1 = 1 + MOD(nm0,aim_surfForc_NppCycle)
                0122 C-     interpolation weight:
                0123         dTprd = tmprd - (nm0 - 1)
                0124         aim_sWght1 = 0.5 _d 0+(dTprd-0.5 _d 0)*aim_surfForc_TransRatio
                0125         aim_sWght1 = MAX( 0. _d 0, MIN(1. _d 0, aim_sWght1) )
                0126         aim_sWght0 = 1. _d 0 - aim_sWght1
                0127 
65007c221b Jean*0128 C--   Compute surface forcing at present time (linear Interp in time)
                0129 C     using F.Molteni surface BC form ; fields needed are:
4b996cbf32 Jean*0130 C     1. Sea  Surface temperatures  (in situ Temp. [K])
                0131 C     2. Land Surface temperatures  (in situ Temp. [K])
                0132 C     3. Soil moisture         (between 0-1)
                0133 C     4. Snow depth, Sea Ice : used to compute albedo (=> local arrays)
                0134 C     5. Albedo                (between 0-1)
65007c221b Jean*0135 
d0a9461855 Jean*0136 C-    Surface Temperature:
65007c221b Jean*0137         DO j=1,sNy
                0138          DO i=1,sNx
                0139           I2 = i+(j-1)*sNx
d0a9461855 Jean*0140           sst1(I2,myThid) = aim_sWght0*aim_sst0(i,j,bi,bj)
65007c221b Jean*0141      &                    + aim_sWght1*aim_sst1(i,j,bi,bj)
                0142           stl1(I2,myThid) = aim_sWght0*aim_lst0(i,j,bi,bj)
                0143      &                    + aim_sWght1*aim_lst1(i,j,bi,bj)
                0144          ENDDO
                0145         ENDDO
                0146 
                0147 C-    Soil Water availability : (from F.M. INFORC S/R)
                0148         SDEP1 = 70. _d 0
                0149         IDEP2 =  3. _d 0
                0150         SDEP2 = IDEP2*SDEP1
                0151 
                0152         SWWIL2= SDEP2*SWWIL
                0153         RSW   = 1. _d 0/(SDEP1*SWCAP+SDEP2*(SWCAP-SWWIL))
d0a9461855 Jean*0154 
65007c221b Jean*0155         DO j=1,sNy
                0156          DO i=1,sNx
                0157           I2 = i+(j-1)*sNx
d0a9461855 Jean*0158           soilw_0 = ( aim_sw10(i,j,bi,bj)
65007c221b Jean*0159      &     +aim_veget(i,j,bi,bj)*
                0160      &      MAX(IDEP2*aim_sw20(i,j,bi,bj)-SWWIL2, 0. _d 0)
d0a9461855 Jean*0161      &              )*RSW
                0162           soilw_1 = ( aim_sw11(i,j,bi,bj)
65007c221b Jean*0163      &     +aim_veget(i,j,bi,bj)*
                0164      &      MAX(IDEP2*aim_sw21(i,j,bi,bj)-SWWIL2, 0. _d 0)
d0a9461855 Jean*0165      &              )*RSW
                0166           soilw1(I2,myThid) = aim_sWght0*soilw_0
65007c221b Jean*0167      &                      + aim_sWght1*soilw_1
                0168           soilw1(I2,myThid) = MIN(1. _d 0, soilw1(I2,myThid) )
                0169          ENDDO
                0170         ENDDO
                0171 
                0172 C-    Set snow depth & sea-ice fraction :
                0173         DO j=1,sNy
                0174          DO i=1,sNx
                0175           I2 = i+(j-1)*sNx
                0176           snow1(I2) = aim_sWght0*aim_snw0(i,j,bi,bj)
d0a9461855 Jean*0177      &              + aim_sWght1*aim_snw1(i,j,bi,bj)
65007c221b Jean*0178           oice1(I2) = aim_sWght0*aim_oic0(i,j,bi,bj)
d0a9461855 Jean*0179      &              + aim_sWght1*aim_oic1(i,j,bi,bj)
65007c221b Jean*0180          ENDDO
                0181         ENDDO
                0182 
b3097ed02d Jean*0183         IF (aim_splitSIOsFx) THEN
                0184 C-    Split Ocean and Sea-Ice surf. temp. ; remove ice-fraction < 1 %
cdcb187d4c Jean*0185 c        oceTfreez = tFreeze - 1.9 _d 0
                0186          oceTfreez = celsius2K - 1.9 _d 0
b3097ed02d Jean*0187          DO J=1,NGP
d0a9461855 Jean*0188           sti1(J,myThid) = sst1(J,myThid)
b3097ed02d Jean*0189           IF ( oice1(J) .GT. 1. _d -2 ) THEN
5ef3ea4448 Ed H*0190             sst1(J,myThid) = MAX(sst1(J,myThid),oceTfreez)
d0a9461855 Jean*0191             sti1(J,myThid) = sst1(J,myThid)
b3097ed02d Jean*0192      &                     +(sti1(J,myThid)-sst1(J,myThid))/oice1(J)
                0193           ELSE
                0194             oice1(J) = 0. _d 0
                0195           ENDIF
                0196          ENDDO
                0197         ELSE
                0198          DO J=1,NGP
d0a9461855 Jean*0199           sti1(J,myThid) = sst1(J,myThid)
b3097ed02d Jean*0200          ENDDO
                0201         ENDIF
                0202 
65007c221b Jean*0203 C-    Surface Albedo : (from F.M. FORDATE S/R)
b3097ed02d Jean*0204 c_FM    DALB=ALBICE-ALBSEA
65007c221b Jean*0205         RSD=1. _d 0/SDALB
972eaa29a0 Dani*0206         ALPHA= 2. _d 0*PI*(TYEAR+10. _d 0/365. _d 0)
6fb2765f24 Jean*0207 #ifdef ALLOW_INSOLATION
7f98c35e47 Davi*0208         ZS = - SIN(OBLIQ * deg2rad) * COS(ALPHA)
                0209         ZC =   ASIN( ZS )
                0210         ZC =   COS(ZC)
6fb2765f24 Jean*0211 #else /* ALLOW_INSOLATION */
                0212         RZEN = COS(ALPHA) * ( -23.45 _d 0 * deg2rad)
                0213         ZC = COS(RZEN)
                0214         ZS = SIN(RZEN)
                0215 #endif /* ALLOW_INSOLATION */
65007c221b Jean*0216         DO j=1,sNy
                0217          DO i=1,sNx
                0218 c_FM      SNOWC=MIN(1.,RSD*SNOW1(I,J))
                0219 c_FM      ALBL=ALB0(I,J)+MAX(ALBSN-ALB0(I,J),0.0)*SNOWC
                0220 c_FM      ALBS=ALBSEA+DALB*OICE1(I,J)
                0221 c_FM      ALB1(I,J)=FMASK1(I,J)*ALBL+FMASK0(I,J)*ALBS
                0222           I2 = i+(j-1)*sNx
                0223           alb_land = aim_albedo(i,j,bi,bj)
                0224      &       + MAX( 0. _d 0, ALBSN-aim_albedo(i,j,bi,bj) )
                0225      &        *MIN( 1. _d 0, RSD*snow1(I2))
b3097ed02d Jean*0226 c         alb_sea  = ALBSEA + DALB*oice1(I2)
d0a9461855 Jean*0227 c         alb1(I2,0,myThid) = alb_sea
b3097ed02d Jean*0228 c    &        + (alb_land - alb_sea)*fMask1(I2,1,myThid)
72b4e12f25 Dani*0229           ALBSEA1 = ALBSEA
                0230           IF ( aim_selectOceAlbedo .EQ. 1) THEN
                0231            SJ = SIN(yC(i,j,bi,bj) * deg2rad)
972eaa29a0 Dani*0232            CJ = COS(yC(i,j,bi,bj) * deg2rad)
72b4e12f25 Dani*0233            TMPA = SJ*ZS
972eaa29a0 Dani*0234            TMPB = CJ*ZC
72b4e12f25 Dani*0235            TMPL = -TMPA/TMPB
972eaa29a0 Dani*0236            IF (TMPL .GE. 1.0 _d 0) THEN
72b4e12f25 Dani*0237             CZEN = 0.0 _d 0
972eaa29a0 Dani*0238            ELSEIF (TMPL .LE. -1.0 _d 0) THEN
72b4e12f25 Dani*0239             CZEN = (2.0 _d 0)*TMPA*PI
972eaa29a0 Dani*0240             CZEN2= PI*((2.0 _d 0)*TMPA*TMPA + TMPB*TMPB)
72b4e12f25 Dani*0241             CZEN = CZEN2/CZEN
                0242            ELSE
972eaa29a0 Dani*0243             hlim = ACOS(TMPL)
                0244             CZEN = 2.0 _d 0*(TMPA*hlim + TMPB*SIN(hlim))
                0245             CZEN2= 2.0 _d 0*TMPA*TMPA*hlim
                0246      &          + 4.0 _d 0*TMPA*TMPB*SIN(hlim)
                0247      &          + TMPB*TMPB*( hlim + 0.5 _d 0*SIN(2.0 _d 0*hlim) )
72b4e12f25 Dani*0248             CZEN = CZEN2/CZEN
                0249            ENDIF
                0250            ALBSEA1 = ( ( 2.6 _d 0 / (CZEN**(1.7 _d 0) + 0.065 _d 0) )
                0251      &          + ( 15. _d 0 * (CZEN-0.1 _d 0) * (CZEN-0.5 _d 0)
                0252      &          * (CZEN-1.0 _d 0) ) ) / 100.0 _d 0
                0253           ENDIF
b3097ed02d Jean*0254           alb1(I2,1,myThid) = alb_land
72b4e12f25 Dani*0255 C_DE      alb1(I2,2,myThid) = ALBSEA
2a80e4d00e Jean*0256           alb1(I2,2,myThid) = 0.5 _d 0 * ALBSEA
72b4e12f25 Dani*0257      &        + 0.5 _d 0 * ALBSEA1
b3097ed02d Jean*0258           alb1(I2,3,myThid) = ALBICE
65007c221b Jean*0259          ENDDO
                0260         ENDDO
                0261 
                0262 C-- else aim_useFMsurfBC
                0263       ELSE
                0264 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0265 
d0a9461855 Jean*0266 C-    safer to initialise output argument aim_sWght0,1
                0267 C     even if they are not used when aim_useFMsurfBC=F
                0268         aim_sWght1 = 0. _d 0
                0269         aim_sWght0 = 1. _d 0
                0270 
65007c221b Jean*0271 C-    Set surface forcing fields needed by atmos. physics package
                0272 C     1. Albedo                (between 0-1)
                0273 C     2. Sea  Surface temperatures  (in situ Temp. [K])
                0274 C     3. Land Surface temperatures  (in situ Temp. [K])
                0275 C     4. Soil moisture         (between 0-1)
d0a9461855 Jean*0276 C        Snow depth, Sea Ice (<- no need for now)
65007c221b Jean*0277 
                0278 C      Set surface albedo data (in [0,1]) from aim_albedo to alb1 :
                0279        IF (aim_useMMsurfFc) THEN
                0280         DO j=1,sNy
                0281          DO i=1,sNx
                0282           I2 = i+(j-1)*sNx
b3097ed02d Jean*0283           alb1(I2,1,myThid) = aim_albedo(i,j,bi,bj)
                0284           alb1(I2,2,myThid) = aim_albedo(i,j,bi,bj)
                0285           alb1(I2,3,myThid) = aim_albedo(i,j,bi,bj)
65007c221b Jean*0286          ENDDO
                0287         ENDDO
                0288        ELSE
                0289         DO j=1,sNy
                0290          DO i=1,sNx
                0291           I2 = i+(j-1)*sNx
b3097ed02d Jean*0292           alb1(I2,1,myThid) = 0.
                0293           alb1(I2,2,myThid) = 0.
                0294           alb1(I2,3,myThid) = 0.
65007c221b Jean*0295          ENDDO
                0296         ENDDO
                0297        ENDIF
                0298 C      Set surface temperature data from aim_S/LSurfTemp to sst1 & stl1 :
                0299        IF (aim_useMMsurfFc) THEN
                0300         DO j=1,sNy
                0301          DO i=1,sNx
                0302           I2 = i+(j-1)*sNx
d0a9461855 Jean*0303           sst1(I2,myThid) = aim_sst0(i,j,bi,bj)
                0304           stl1(I2,myThid) = aim_sst0(i,j,bi,bj)
                0305           sti1(I2,myThid) = aim_sst0(i,j,bi,bj)
65007c221b Jean*0306          ENDDO
                0307         ENDDO
                0308        ELSE
                0309         DO j=1,sNy
                0310          DO i=1,sNx
                0311           I2 = i+(j-1)*sNx
                0312           sst1(I2,myThid) = 300.
                0313           stl1(I2,myThid) = 300.
b3097ed02d Jean*0314           sti1(I2,myThid) = 300.
65007c221b Jean*0315          ENDDO
                0316         ENDDO
                0317        ENDIF
                0318 
d0a9461855 Jean*0319 C-     Set soil water availability (in [0,1]) from aim_sw10 to soilw1 :
65007c221b Jean*0320        IF (aim_useMMsurfFc) THEN
                0321         DO j=1,sNy
                0322          DO i=1,sNx
                0323           I2 = i+(j-1)*sNx
d0a9461855 Jean*0324           soilw1(I2,myThid) = aim_sw10(i,j,bi,bj)
65007c221b Jean*0325          ENDDO
                0326         ENDDO
                0327        ELSE
                0328         DO j=1,sNy
                0329          DO i=1,sNx
                0330           I2 = i+(j-1)*sNx
                0331           soilw1(I2,myThid) = 0.
                0332          ENDDO
                0333         ENDDO
                0334        ENDIF
                0335 
d0a9461855 Jean*0336 C-     Set Snow depth and Sea Ice
65007c221b Jean*0337 C      (not needed here since albedo is loaded from file)
b3097ed02d Jean*0338         DO j=1,sNy
                0339          DO i=1,sNx
                0340           I2 = i+(j-1)*sNx
                0341           oice1(I2) = 0.
                0342           snow1(I2) = 0.
                0343          ENDDO
                0344         ENDDO
65007c221b Jean*0345 
                0346 C-- endif/else aim_useFMsurfBC
                0347       ENDIF
                0348 
e19eef4074 Jean*0349 #ifdef COMPONENT_MODULE
                0350       IF ( useCoupler ) THEN
d0a9461855 Jean*0351 C--   take surface data from the ocean component
97ecf03a3f Jean*0352 C     to replace MxL fields (if use sea-ice) or directly AIM SST
                0353         CALL ATM_APPLY_IMPORT(
                0354      I           aim_landFr,
d0a9461855 Jean*0355      U           sst1(1,myThid), oice1,
                0356      I           myTime, myIter, bi, bj, myThid )
e19eef4074 Jean*0357       ENDIF
870deba1f6 Davi*0358 #endif /* COMPONENT_MODULE */
4c105e462e Davi*0359 
                0360 #ifdef ALLOW_AIM_CO2
2b1795828b Jean*0361       DO j=1-OLy,sNy+OLy
                0362         DO i=1-OLx,sNx+OLx
                0363            aim_CO2(i,j,bi,bj)= atm_pCO2
4c105e462e Davi*0364         ENDDO
2a80e4d00e Jean*0365       ENDDO
870deba1f6 Davi*0366 #ifdef ALLOW_DIAGNOSTICS
fb7fcd6036 Davi*0367       IF ( useDiagnostics ) THEN
                0368          pCO2scl = 1. _d 6
2b1795828b Jean*0369          CALL DIAGNOSTICS_SCALE_FILL( aim_CO2(1-OLx,1-OLy,bi,bj),
                0370      &               pCO2scl, 1, 'aim_pCO2', 1, 1, 2, bi, bj, myThid )
4c105e462e Davi*0371       ENDIF
870deba1f6 Davi*0372 #endif /* ALLOW_DIAGNOSTICS */
2a80e4d00e Jean*0373 #endif /* ALLOW_AIM_CO2 */
4c105e462e Davi*0374 
a27dd2281d Jean*0375 #ifdef ALLOW_LAND
                0376       IF (useLand) THEN
                0377 C-    Use land model output instead of prescribed Temp & moisture
d0a9461855 Jean*0378         CALL AIM_LAND2AIM(
b3097ed02d Jean*0379      I           aim_landFr, aim_veget, aim_albedo, snow1,
d0a9461855 Jean*0380      U           stl1(1,myThid), soilw1(1,myThid), alb1(1,1,myThid),
                0381      I           myTime, myIter, bi, bj, myThid )
a27dd2281d Jean*0382       ENDIF
                0383 #endif /* ALLOW_LAND */
                0384 
cdcb187d4c Jean*0385 #ifdef ALLOW_THSICE
                0386       IF (useThSIce) THEN
                0387 C-    Use thermo. sea-ice model output instead of prescribed Temp & albedo
d0a9461855 Jean*0388         CALL AIM_SICE2AIM(
97ecf03a3f Jean*0389      I           aim_landFr,
d0a9461855 Jean*0390      U           sst1(1,myThid), oice1,
                0391      O           sti1(1,myThid), alb1(1,3,myThid),
                0392      I           myTime, myIter, bi, bj, myThid )
cdcb187d4c Jean*0393       ENDIF
                0394 #endif /* ALLOW_THSICE */
                0395 
b3097ed02d Jean*0396 C-- set the sea-ice & open ocean fraction :
                0397         DO J=1,NGP
                0398           fMask1(J,3,myThid) =(1. _d 0 - fMask1(J,1,myThid))
                0399      &                        *oice1(J)
d0a9461855 Jean*0400           fMask1(J,2,myThid) = 1. _d 0 - fMask1(J,1,myThid)
b3097ed02d Jean*0401      &                                 - fMask1(J,3,myThid)
                0402         ENDDO
                0403 
                0404 C-- set the mean albedo :
                0405         DO J=1,NGP
                0406           alb1(J,0,myThid) = fMask1(J,1,myThid)*alb1(J,1,myThid)
                0407      &                     + fMask1(J,2,myThid)*alb1(J,2,myThid)
                0408      &                     + fMask1(J,3,myThid)*alb1(J,3,myThid)
                0409         ENDDO
                0410 
e749d70ece Jean*0411 C-- initialize surf. temp. change to zero:
                0412         DO k=1,3
                0413          DO J=1,NGP
                0414           dTsurf(J,k,myThid) = 0.
                0415          ENDDO
                0416         ENDDO
                0417 
b3097ed02d Jean*0418         IF (.NOT.aim_splitSIOsFx) THEN
                0419          DO J=1,NGP
                0420           fMask1(J,3,myThid) = 0. _d 0
d0a9461855 Jean*0421           fMask1(J,2,myThid) = 1. _d 0 - fMask1(J,1,myThid)
b3097ed02d Jean*0422          ENDDO
                0423         ENDIF
                0424 
65007c221b Jean*0425 #endif /* ALLOW_AIM */
                0426 
                0427       RETURN
                0428       END
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated from https://github.com/MITgcm/MITgcm by the 2.2.1-MITgcm-0.1 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!