/[MITgcm]/MITgcm_contrib/jscott/igsm/src/radia_chem.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/igsm/src/radia_chem.F

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

revision 1.2 by jscott, Tue Aug 22 20:25:52 2006 UTC revision 1.3 by jscott, Mon Apr 23 21:20:18 2007 UTC
# Line 82  C Line 82  C
82  #if ( defined OCEAN_3D )  #if ( defined OCEAN_3D )
83  #include "AGRID.h"  #include "AGRID.h"
84  #endif  #endif
85           dimension SWNET(jm0,2),SWIN(jm0,2)
86    
87  #if ( defined CLM )  #if ( defined CLM )
88  #include "CLM.COM"  #include "CLM.h"
89  #endif  #endif
90  c  c
91        common/conprn/HPRNT        common/conprn/HPRNT
# Line 100  c Line 101  c
101       &,BSO4LAND(JM0),BSO4OCEAN(JM0),BSO4TOTAL(JM0)       &,BSO4LAND(JM0),BSO4OCEAN(JM0),BSO4TOTAL(JM0)
102        dimension DSWSRF(jm0),DLWSRF(jm0),DSWVIS(jm0),DSWNIR(jm0)        dimension DSWSRF(jm0),DLWSRF(jm0),DSWVIS(jm0),DSWNIR(jm0)
103        integer PCLOUD        integer PCLOUD
104        common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTSURF(JM0)  !     common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTSURF(JM0)
105       *,cfcld(JM0,3)  !    *,cfcld(JM0,3)
106    #include "TSRF.COM"
107        CHARACTER*4 JMNTHF,JMLAST        CHARACTER*4 JMNTHF,JMLAST
108        DATA JMLAST /'LAST'/        DATA JMLAST /'LAST'/
109        DATA TF/273.16/,TCIR/258.16/,STBO/.567257E-7/,IFIRST/1/,JDLAST/-9/5045.          DATA TF/273.16/,TCIR/258.16/,STBO/.567257E-7/,IFIRST/1/,JDLAST/-9/5045.  
# Line 335  c         cfcld(j,k)=1.+coefcl(k)*DTSURF Line 337  c         cfcld(j,k)=1.+coefcl(k)*DTSURF
337  c 925  continue  c 925  continue
338          DTSURFAV=0.          DTSURFAV=0.
339          do j=1,jm          do j=1,jm
340           DTSURFAV=DTSURFAV+DTSURF(J)*DXYP(j)           DTSURFAV=DTSURFAV+DT2MGL(J)*DXYP(j)
341          end do  !j          end do  !j
342           DTSURFAV=DTSURFAV/AREAG           DTSURFAV=DTSURFAV/AREAG
343          do j=1,jm          do j=1,jm
# Line 420  CE     END OF READING OF CLOUD Line 422  CE     END OF READING OF CLOUD
422          print *,'cfcld'          print *,'cfcld'
423          print 9456,cfcld          print 9456,cfcld
424          print *,' DTSURF'          print *,' DTSURF'
425          print 9456,DTSURF          print 9456,DT2MGL
426          print *,' DTSURFAV=',DTSURFAV          print *,' DTSURFAV=',DTSURFAV
427   9456   format(12f6.2)   9456   format(12f6.2)
428         endif         endif
# Line 487  c     RVOL=0.012 Line 489  c     RVOL=0.012
489    
490        JLAT=J                                                            5145.1          JLAT=J                                                            5145.1  
491        IF(MODRJ.EQ.0) CALL RCOMPJ                                        5146.          IF(MODRJ.EQ.0) CALL RCOMPJ                                        5146.  
492          SWIN(j,1)=0.0
493          SWNET(j,1)=0.0
494          SWIN(j,2)=0.0
495          SWNET(j,2)=0.0
496  C****                                                                   5147.    C****                                                                   5147.  
497  C**** MAIN I LOOP                                                       5148.    C**** MAIN I LOOP                                                       5148.  
498  C****                                                                   5149.    C****                                                                   5149.  
# Line 863  C        for TEM CLM Line 869  C        for TEM CLM
869            AJ(J,3)=AJ(J,3)+(SRNFLB(1+LM)*COSZ)*POCEAN            AJ(J,3)=AJ(J,3)+(SRNFLB(1+LM)*COSZ)*POCEAN
870            AJ(J,71)=AJ(J,71)-(TRNFLB(1+LM)-TRNFLB(1))*POCEAN            AJ(J,71)=AJ(J,71)-(TRNFLB(1+LM)-TRNFLB(1))*POCEAN
871  #if ( defined OCEAN_3D )  #if ( defined OCEAN_3D )
872            solarinc_ocean(J)=solarinc_ocean(J)+SRDFLB(1)*COSZ            SWIN(j,1)=SRDFLB(1)
873            solarnet_ocean(J)=solarnet_ocean(J)+SRNFLB(1)*COSZ            SWNET(j,1)=SRNFLB(1)
           navrado(j)=navrado(j)+1  
874  #endif  #endif
875  C  C
876            DO  K=2,9            DO  K=2,9
# Line 885  C Line 890  C
890            CJ(J,3)=CJ(J,3)+(SRNFLB(1+LM)*COSZ)*POICE            CJ(J,3)=CJ(J,3)+(SRNFLB(1+LM)*COSZ)*POICE
891            CJ(J,71)=CJ(J,71)-(TRNFLB(1+LM)-TRNFLB(1))*POICE            CJ(J,71)=CJ(J,71)-(TRNFLB(1+LM)-TRNFLB(1))*POICE
892  #if ( defined OCEAN_3D )  #if ( defined OCEAN_3D )
893            solarinc_ice(J)=solarinc_ice(J)+SRDFLB(1)*COSZ            SWIN(j,2)=SRDFLB(1)
894            solarnet_ice(J)=solarnet_ice(J)+SRNFLB(1)*COSZ            SWNET(j,2)=SRNFLB(1)
           navrad(j)=navrad(j)+1  
895  #endif  #endif
896  C  C
897            DO  K=2,9            DO  K=2,9
# Line 1146  C**** Line 1150  C****
1150       *  +TRHRS(I,J,LR))*COE(LR+LM)                                      5513.         *  +TRHRS(I,J,LR))*COE(LR+LM)                                      5513.  
1151    840 DO 860 J=1,JM                                                     5514.      840 DO 860 J=1,JM                                                     5514.  
1152  #if ( defined CLM )  #if ( defined CLM )
1153            dsw4clm(j)=DSWSRF(j)*COSZ1(1,j)            i=1
1154            dlw4clm(j)=DLWSRF(j)            dsw4clm(i,j)=DSWSRF(j)*COSZ1(1,j)
1155            swinr4clm(j)=DSWNIR(j)*COSZ1(1,j)            dlw4clm(i,j)=DLWSRF(j)
1156            swvis4clm(j)=DSWVIS(j)*COSZ1(1,j)            swinr4clm(i,j)=DSWNIR(j)*COSZ1(1,j)
1157              swvis4clm(i,j)=DSWVIS(j)*COSZ1(1,j)
1158  c    For TEM  c    For TEM
1159            swtd4tem(j)=swtd4tem(j)+S0*COSZ1(1,j)            swtd4tem(j)=swtd4tem(j)+S0*COSZ1(1,j)
1160            swsd4tem(j)=swsd4tem(j)+DSWSRF(j)*COSZ1(1,j)            swsd4tem(j)=swsd4tem(j)+DSWSRF(j)*COSZ1(1,j)
1161            nradd4tem(j)=nradd4tem(j)+1            nradd4tem(j)=nradd4tem(j)+1
1162  #endif  #endif
1163    #if ( defined OCEAN_3D )
1164              solarinc_ocean(J)=solarinc_ocean(J)+SWIN(j,1)*COSZ1(1,j)
1165              solarnet_ocean(J)=solarnet_ocean(J)+SWNET(j,1)*COSZ1(1,j)
1166              solarinc_ice(J)=solarinc_ice(J)+SWIN(j,2)*COSZ1(1,j)
1167              solarnet_ice(J)=solarnet_ice(J)+SWNET(j,2)*COSZ1(1,j)
1168              navrado(j)=navrado(j)+1
1169              navrad(j)=navrad(j)+1
1170    #endif
1171        IMAX=IM                                                           5515.          IMAX=IM                                                           5515.  
1172        IF(J.EQ.1.OR.J.EQ.JM) IMAX=1                                      5516.          IF(J.EQ.1.OR.J.EQ.JM) IMAX=1                                      5516.  
1173        DO 860 L=1,LM                                                     5517.          DO 860 L=1,LM                                                     5517.  
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22