6,11d5
< ! Modified by Tarassova 2015
< !  Climate aerosol (Kinne, 2013) coarse mode included
< !  Fine aerosol mode 2000 included
< !  Modifications (5) are marked by 
< !  !tar begin   and  !tar end
< !
85c79
<        ifalb  , ifsst ,ifco2flx, &
---
>        ifalb  , ifsst , &
88,91d81
< !tar begin
< !climate aerosol selection parameter
<        ifaeros, &
< !tar end
140d129
<       gco2flx        , &
144,158d132
< !tar begin  
< !climate aerosol optical parameters of coarse mode   
<       aod            , &
<       asy            , &
<       ssa            , &
<       z_aer          , &    
< !tar end
< !
< !tar begin  
< !climate aerosol optical parameters of coarse mode   
<       aodF            , &
<       asyF            , &
<       ssaF            , &
<       z_aerF          , &    
< !tar end 
200c174
<       sfc,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy,poda,tmin2m   ,tmax2m 
---
>       sfc,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy
218a193
>   REAL(KIND=r8)   :: rbyg
326c301
<                        kMax          ,path          ,fNameSibVeg_in, &
---
>                        kMax          ,delsig        ,path          ,fNameSibVeg_in, &
329c304
<                        idate         ,idatec        , &
---
>                        idate         ,idatec        ,si            ,sl            , &
335a311
>     REAL(KIND=r8)   , INTENT(IN   ) :: delsig(kMax)
346a323,324
>     REAL(KIND=r8)   , INTENT(IN   ) :: si(kMax+1)
>     REAL(KIND=r8)   , INTENT(IN   ) :: sl(kMax)
349a328
>     rbyg  =gasr/grav*delsig(1)*0.5_r8
414c393
<                               tod       ,idate     ,&
---
>                               tod       ,si(1)     ,sl(1)     ,idate     ,&
427a407,408
>                                si1       ,&
>                                sl1       ,&
443a425,426
>    REAL(KIND=r8), INTENT(IN   ) :: si1
>    REAL(KIND=r8), INTENT(IN   ) :: sl1
487c470
<    CALL vegin ()
---
>    CALL vegin (si1   , sl1)
490c473
<    OPEN (UNIT=nfsibt, FILE=TRIM(fNameSibmsk),FORM='UNFORMATTED', ACCESS='DIRECT', RECL=LRecIN,&
---
>    OPEN (UNIT=nfsibt, FILE=TRIM(fNameSibmsk),FORM='FORMATTED', ACCESS='SEQUENTIAL',&
499c482
<    OPEN (UNIT=nftgz0,FILE=TRIM(fNameTg3zrl), FORM='UNFORMATTED', ACCESS='DIRECT', RECL=LRecIN, &
---
>    OPEN (UNIT=nftgz0,FILE=TRIM(fNameTg3zrl), FORM='FORMATTED', ACCESS='SEQUENTIAL', &
508c491
<    OPEN (UNIT=nfzol,FILE=TRIM(fNameRouLen),FORM='UNFORMATTED', ACCESS='DIRECT', RECL=LRecIN, &
---
>    OPEN (UNIT=nfzol,FILE=TRIM(fNameRouLen),FORM='FORMATTED', ACCESS='SEQUENTIAL',&
516c499
<    READ(UNIT=nfsibt, REC=1) ibuf
---
>    READ(nfsibt, *) ibuf
555c538
<       CALL getsbc (iMax ,jMax  ,kMax, AlbVisDiff,gtsea,gco2flx,gndvi,soilm,sheleg,o3mix,tracermix,wsib3d,&
---
>       CALL getsbc (iMax ,jMax  ,kMax, AlbVisDiff,gtsea,gndvi,soilm,sheleg,o3mix,tracermix,wsib3d,&
557c540
<            ifalb,ifsst,ifco2flx,ifndvi,ifslm ,ifslmSib2,ifsnw,ifozone,iftracer, &
---
>            ifalb,ifsst,ifndvi,ifslm ,ifslmSib2,ifsnw,ifozone,iftracer, &
561,562c544,545
<       CALL ReadGetNFTGZ(nftgz0,irec,buf(1:iMax,1:jMax,1),buf(1:iMax,1:jMax,2),buf(1:iMax,1:jMax,3))
<       READ (UNIT=nfzol, REC=1) brf
---
>       CALL ReadGetNFTGZ(nftgz0,irec,buf(1:iMax,1:jMax,1),buf(1:iMax,1:jMax,2),buf(1:iMax,1:jMax,3),'txt')
>       READ (nfzol, *) brf
689d671
<             IF(iMaskSSiB(i,j) >= 1_i8)gco2flx(i,j)=0.0_r8
766,767c748,749
<       READ(UNIT=nfsibi) gl0   ,zorl  ,gtsea ,gco2flx,tseam,qsfc0,tsfc0,qsfcm,tsfcm,HML,HUML,HVML,TSK,z0sea,&
<       TC_SeaIce,TGS_SeaIce,TD_SeaIce,TA_SeaIce,SNOA_SeaIce,SNOB_SeaIce,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy,poda,tmin2m   ,tmax2m 
---
>       READ(UNIT=nfsibi) gl0   ,zorl  ,gtsea ,tseam,qsfc0,tsfc0,qsfcm,tsfcm,HML,HUML,HVML,TSK,z0sea,&
>       TC_SeaIce,TGS_SeaIce,TD_SeaIce,TA_SeaIce,SNOA_SeaIce,SNOB_SeaIce,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy
775d756
<          IF(ifco2flx == -2 .or. ifco2flx == -3 )ifco2flx=-1
777c758
<          CALL getsbc (iMax ,jMax  ,kMax, AlbVisDiff,gtsea,gco2flx,gndvi,soilm,sheleg,o3mix,tracermix,wsib3d,&
---
>          CALL getsbc (iMax ,jMax  ,kMax, AlbVisDiff,gtsea,gndvi,soilm,sheleg,o3mix,tracermix,wsib3d,&
779c760
<            ifalb,ifsst,ifco2flx,ifndvi,ifslm ,ifslmSib2,ifsnw,ifozone,iftracer, &
---
>            ifalb,ifsst,ifndvi,ifslm ,ifslmSib2,ifsnw,ifozone,iftracer, &
788d768
<                   IF(iMaskSSiB(i,j) >= 1_i8) gco2flx(i,j)=0.0_r8
849d828
<        prsi               ,prsl                 ,phii             ,phil             ,&
856,857c835,836
<        tseam              ,tsurf                ,qsurf            ,&
<        imask              ,itype                ,tg               ,&
---
>        tseam              ,tsurf                ,qsurf            ,sigki            ,&
>        delsig             ,imask                ,itype            ,tg               ,&
896,899d874
<     REAL(KIND=r8), INTENT(IN   ) :: prsi     (nCols,kMax+1)  !     prsi     - real, pressure at layer interfaces             ix,levs+1  Pa
<     REAL(KIND=r8), INTENT(IN   ) :: prsl     (nCols,kMax)    !     prsl     - real, mean layer presure                       ix,levs   Pa
<     REAL(KIND=r8), INTENT(IN   ) :: phii     (nCols,kMax+1)  !===>  PHIH(K+1)  INPUT GEOPOTENTIAL @ EDGES  IN MKS units (m)
<     REAL(KIND=r8), INTENT(IN   ) :: phil     (nCols,kMax)    !===>  PHIL(K)    INPUT GEOPOTENTIAL @ LAYERS IN MKS units (m)
932a908,910
> 
>     REAL(KIND=r8), INTENT(IN   ) :: sigki    (1)
>     REAL(KIND=r8), INTENT(IN   ) :: delsig   (1:kMax)
1262,1263c1240
<           !dzm   (ncount)=rbyg*tm(ncount)           
<           dzm   (ncount)=0.5_r8*MAX((phii(i,2) - phii(i,1)),0.5_r8)
---
>           dzm   (ncount)=rbyg*tm(ncount)
1271,1272c1248
<           !bps   (ncount)=sigki(1)
<           bps   (ncount)= (prsi(i,1)/(prsl(i,1)))**(gasr/cp)
---
>           bps   (ncount)=sigki(1)
1276c1252
<           psb   (ncount)=psur(ncount)*(prsi(i,1)/prsi(i,1)) -(prsi(i,2)/prsi(i,1) )
---
>           psb   (ncount)=psur(ncount)*delsig(1)
1652,1653c1628,1629
<            prsi(1:ncols,1:kMax+1),prsl(1:ncols,1:kMax),phii(1:nCols,1:kMax+1),phil(1:nCols,1:kMax),&
<            xsea  (1:nCols)           ,dtc3x                       ,SIN(colrad(1:nCols))        ,&
---
>            gps   (1:nCols)             ,xsea  (1:nCols)           ,dtc3x                       ,&
>            SIN(colrad(1:nCols))        ,sigki (1)                 ,delsig(1:kMax)              ,&
7136c7112
<    REAL(KIND=r8) :: IceOceanAlb (nCols,2,2)  
---
>    REAL(KIND=r8) :: IceOceanAlb (2,2)  
7142,7143c7118,7119
<     REAL(KIND=r8) ::         tau(ncols)             !aerosol/cloud optical depth
<     REAL(KIND=r8) ::         chl (ncols)            !chlorophyll concentration in mg/m3
---
>     REAL(KIND=r8) ::         tau             !aerosol/cloud optical depth
>     REAL(KIND=r8) ::         chl             !chlorophyll concentration in mg/m3
7232,7233c7208,7209
<     IceOceanAlb=0.0_r8 
<     chl = 0.10_r8              !chlorophyll concentration in mg/m3
---
> 
>     chl = 0.10              !chlorophyll concentration in mg/m3
7239,7248c7215
<     tau=0.0_r8
<     DO k=1,kMax
<        DO i=1,ncols
<           tau(i)=tau(i)+ taud(i,k) ! tau = SUM(taud(i,1:kMax))              !aerosol/cloud optical depth
<        END DO
<     END DO   
<     DO i=1,ncols
<         tau(i)=MIN(MAX(tau(i),0.0_r8),25.0_r8)      !aerosol/cloud optical depth
<     END DO
< 
---
>   
7262c7229,7230
<              avisb(i)=GetOceanAlb(i,tau(i),f,wind(i),chl(i),0.2_r8,0.7_r8) !   s (shortwave)  = 0.2-0.7 micro-meters
---
>              tau = SUM(taud(i,1:kMax))              !aerosol/cloud optical depth
>              avisb(i)=GetOceanAlb(tau,f,wind(i),chl,0.2_r8,0.7_r8) !   s (shortwave)  = 0.2-0.7 micro-meters
7264c7232
<              anirb(i)=GetOceanAlb(i,tau(i),f,wind(i),chl(i),0.7_r8,5.0_r8) !   l (longwave)   = 0.7-5.0 micro-meters
---
>              anirb(i)=GetOceanAlb(tau,f,wind(i),chl,0.7_r8,5.0_r8) !   l (longwave)   = 0.7-5.0 micro-meters
7269c7237
<                 IceOceanAlb(i,:,:)=GetIceOceanAlb(i,latco,month(i),xVisBeam(i),xVisDiff(i),&
---
>                 IceOceanAlb=GetIceOceanAlb(i,latco,month(i),xVisBeam(i),xVisDiff(i),&
7271,7274c7239,7242
<                 avisb(i)=IceOceanAlb(i,1,1)
<                 avisd(i)=IceOceanAlb(i,1,2)
<                 anirb(i)=IceOceanAlb(i,2,1)
<                 anird(i)=IceOceanAlb(i,2,2)
---
>                 avisb(i)=IceOceanAlb(1,1)
>                 avisd(i)=IceOceanAlb(1,2)
>                 anirb(i)=IceOceanAlb(2,1)
>                 anird(i)=IceOceanAlb(2,2)
7305c7273
<                 IceOceanAlb(i,:,:)=GetIceOceanAlb(i,latco,month(i),xVisBeam(i),xVisDiff(i),&
---
>                 IceOceanAlb=GetIceOceanAlb(i,latco,month(i),xVisBeam(i),xVisDiff(i),&
7307,7310c7275,7278
<                 avisb(i)=IceOceanAlb(i,1,1)
<                 avisd(i)=IceOceanAlb(i,1,2)
<                 anirb(i)=IceOceanAlb(i,2,1)
<                 anird(i)=IceOceanAlb(i,2,2)
---
>                 avisb(i)=IceOceanAlb(1,1)
>                 avisd(i)=IceOceanAlb(1,2)
>                 anirb(i)=IceOceanAlb(2,1)
>                 anird(i)=IceOceanAlb(2,2)
7752c7720
<   SUBROUTINE vegin ( )
---
>   SUBROUTINE vegin ( si1   , sl1)
7754,7755c7722,7723
<     REAL(KIND=r8), PARAMETER :: si1=1
<     REAL(KIND=r8), PARAMETER :: sl1=0.8888
---
>     REAL(KIND=r8),  INTENT(in   ) :: si1
>     REAL(KIND=r8),  INTENT(in   ) :: sl1
9527,9528c9495,9496
<        READ(UNIT=nfsibi) gl0   ,zorl  ,gtsea,gco2flx,tseam,qsfc0,tsfc0,qsfcm,tsfcm,HML,HUML,HVML,TSK,z0sea,&
<       TC_SeaIce,TGS_SeaIce,TD_SeaIce,TA_SeaIce,SNOA_SeaIce,SNOB_SeaIce,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy,poda,tmin2m   ,tmax2m 
---
>        READ(UNIT=nfsibi) gl0   ,zorl  ,gtsea,tseam,qsfc0,tsfc0,qsfcm,tsfcm,HML,HUML,HVML,TSK,z0sea,&
>       TC_SeaIce,TGS_SeaIce,TD_SeaIce,TA_SeaIce,SNOA_SeaIce,SNOB_SeaIce,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy
9533c9501
<        CALL getsbc (iMax ,jMax  ,kMax, AlbVisDiff,gtsea,gco2flx,gndvi,soilm,sheleg,o3mix,tracermix,wsib3d,&
---
>        CALL getsbc (iMax ,jMax  ,kMax, AlbVisDiff,gtsea,gndvi,soilm,sheleg,o3mix,tracermix,wsib3d,&
9535c9503
<            ifalb,ifsst,ifco2flx,ifndvi,ifslm ,ifslmSib2,ifsnw,ifozone,iftracer, &
---
>            ifalb,ifsst,ifndvi,ifslm ,ifslmSib2,ifsnw,ifozone,iftracer, &
9641,9642c9609,9610
<        WRITE(UNIT=nfsibo) gl0 ,zorl,gtsea,gco2flx,tseam,qsfc0,tsfc0,qsfcm,tsfcm,HML,HUML,HVML,TSK,z0sea,&
<       TC_SeaIce,TGS_SeaIce,TD_SeaIce,TA_SeaIce,SNOA_SeaIce,SNOB_SeaIce,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy,poda,tmin2m   ,tmax2m 
---
>        WRITE(UNIT=nfsibo) gl0 ,zorl,gtsea,tseam,qsfc0,tsfc0,qsfcm,tsfcm,HML,HUML,HVML,TSK,z0sea,&
>       TC_SeaIce,TGS_SeaIce,TD_SeaIce,TA_SeaIce,SNOA_SeaIce,SNOB_SeaIce,PBL_CoefKm, PBL_CoefKh,tauresx,tauresy
