/[MITgcm]/MITgcm/pkg/exf/exf_zenithangle_table.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_zenithangle_table.F

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

revision 1.2 by gforget, Tue Apr 13 16:24:50 2010 UTC revision 1.3 by gforget, Wed Apr 14 23:02:18 2010 UTC
# Line 47  C     == local variables == Line 47  C     == local variables ==
47    
48  C     == end of interface ==  C     == end of interface ==
49    
50          _BEGIN_MASTER( myThid )
51    
52  c solar constant  c solar constant
53  c --------------  c --------------
54        SOLC   = 1368. _d 0        SOLC   = 1368. _d 0
55    
        DO bj = myByLo(myThid),myByHi(myThid)  
         DO bi = myBxLo(myThid),myBxHi(myThid)  
56           DO iLat=1,181           DO iLat=1,181
57            DO iTyear=1,366            DO iTyear=1,366
58    
# Line 114  c ------------------- Line 113  c -------------------
113  c compute weighted average of albedo  c compute weighted average of albedo
114  c ----------------------------------  c ----------------------------------
115           if ( 0.5 _d 0 * tmpINT2 .GT. tmpINT1) then           if ( 0.5 _d 0 * tmpINT2 .GT. tmpINT1) then
116              zen_albedo_table(iTyear,iLat,bi,bj)=tmpINT1/tmpINT2              zen_albedo_table(iTyear,iLat)=tmpINT1/tmpINT2
117           else           else
118              zen_albedo_table(iTyear,iLat,bi,bj)=0.5 _d 0              zen_albedo_table(iTyear,iLat)=0.5 _d 0
119           endif           endif
120    
121            ENDDO            ENDDO
122           ENDDO           ENDDO
         ENDDO  
        ENDDO  
123    
124          _END_MASTER( myThid )
125          _BARRIER
126    
127    
128  c determine interpolation coefficient for each grid point  c determine interpolation coefficient for each grid point
# Line 133  c determine interpolation coefficient fo Line 132  c determine interpolation coefficient fo
132            DO i = 1,sNx            DO i = 1,sNx
133             LLLAT=yC(i,j,bi,bj)+91. _d 0             LLLAT=yC(i,j,bi,bj)+91. _d 0
134  c ensure that it is in valid range    c ensure that it is in valid range  
135             LLLAT=max(LLLAT, 1._d 0)             LLLAT=max(LLLAT, 1. _d 0)
136             LLLAT=min(LLLAT, 181._d 0)             LLLAT=min(LLLAT, 181. _d 0)
137  c store  c store
138             zen_albedo_pointer(i,j,bi,bj)=LLLAT             zen_albedo_pointer(i,j,bi,bj)=LLLAT
139            ENDDO            ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22