/[MITgcm]/MITgcm/pkg/fizhi/fizhi_swrad.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_swrad.F

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

revision 1.13 by molod, Mon Jul 26 19:51:08 2004 UTC revision 1.14 by molod, Wed Jul 28 01:25:07 2004 UTC
# Line 120  C ------------------------------ Line 120  C ------------------------------
120        enddo        enddo
121        enddo        enddo
122    
123        if (first .and. myid.eq.0 ) then        if (first .and. myid.eq.1 ) then
124        print *        print *
125        print *,'Low-Level Clouds are Grouped between levels: ',        print *,'Low-Level Clouds are Grouped between levels: ',
126       .         lm,' and ',low_level       .         lm,' and ',low_level
# Line 150  C ************************************** Line 150  C **************************************
150        ENDDO        ENDDO
151        ENDDO        ENDDO
152    
   
153  C **********************************************************************  C **********************************************************************
154  c ****        Compute Two-Dimension Total Cloud Fraction (0-1)      ****  c ****        Compute Two-Dimension Total Cloud Fraction (0-1)      ****
155  C **********************************************************************  C **********************************************************************
# Line 1404  c-----integration over the k-distributio Line 1403  c-----integration over the k-distributio
1403              do i= 1, m              do i= 1, m
1404    
1405               tauwv=xk(ik)*wh(i,j,k)               tauwv=xk(ik)*wh(i,j,k)
1406    
1407  c-----compute total optical thickness, single scattering albedo,  c-----compute total optical thickness, single scattering albedo,
1408  c     and asymmetry factor.  c     and asymmetry factor.
1409    
# Line 1530  c     in certain parallel processors. Line 1529  c     in certain parallel processors.
1529    
1530  c-----flux calculations  c-----flux calculations
1531    
1532          call cldflx (m,n,np,ict,icb,cc,rr,tt,td,rs,ts,         do k= 1, np+1
1533       *               fclr,fall,fsdir,fsdif)          do j= 1, n
1534             do i= 1, m
1535              fclr(i,j,k) = 0.
1536              fall(i,j,k) = 0.
1537             enddo
1538            enddo
1539           enddo
1540           do j= 1, n
1541            do i= 1, m
1542             fsdir(i,j) = 0.
1543             fsdif(i,j) = 0.
1544            enddo
1545           enddo
1546    
1547    c       call cldflx (m,n,np,ict,icb,cc,rr,tt,td,rs,ts,
1548    c    *               fclr,fall,fsdir,fsdif)
1549    
1550         do k= 1, np+1         do k= 1, np+1
1551          do j= 1, n          do j= 1, n
# Line 1887  c-----compute reflectance and transmitta Line 1901  c-----compute reflectance and transmitta
1901    
1902  c-----flux calculations  c-----flux calculations
1903    
1904          call cldflx (m,n,np,ict,icb,cc,rr,tt,td,rs,ts,         do k= 1, np+1
1905       *               fclr,fall,fsdir,fsdif)          do j= 1, n
1906             do i= 1, m
1907              fclr(i,j,k) = 0.
1908              fall(i,j,k) = 0.
1909             enddo
1910            enddo
1911           enddo
1912           do j= 1, n
1913            do i= 1, m
1914             fsdir(i,j) = 0.
1915             fsdif(i,j) = 0.
1916            enddo
1917           enddo
1918    c       call cldflx (m,n,np,ict,icb,cc,rr,tt,td,rs,ts,
1919    c    *               fclr,fall,fsdir,fsdif)
1920    
1921         do k= 1, np+1         do k= 1, np+1
1922          do j= 1, n          do j= 1, n
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22