/[MITgcm]/MITgcm/pkg/atm2d/calc_1dto2d.F
ViewVC logotype

Diff of /MITgcm/pkg/atm2d/calc_1dto2d.F

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

revision 1.1 by jscott, Wed Sep 6 15:32:39 2006 UTC revision 1.2 by jscott, Fri Apr 6 19:56:55 2007 UTC
# Line 42  C     Accumulate runoff into bands (runo Line 42  C     Accumulate runoff into bands (runo
42          IF (ib.LT.numBands) ibj2= rband(ib)          IF (ib.LT.numBands) ibj2= rband(ib)
43          run_b(ib)=0. _d 0          run_b(ib)=0. _d 0
44          DO j=ibj1,ibj2          DO j=ibj1,ibj2
45            run_b(ib)=run_b(ib)+atm_runoff(atm_oc_ind(j))*atm_oc_wgt(j) +            run_b(ib)=run_b(ib) +
46       &              atm_runoff(atm_oc_ind(j)+1)*(1. _d 0-atm_oc_wgt(j))       &              atm_runoff(atm_oc_ind(j))*atm_oc_frac1(j) +
47         &              atm_runoff(atm_oc_ind(j)+1)*atm_oc_frac2(j)
48          ENDDO          ENDDO
49        ENDDO        ENDDO
50        
51        DO j=1,sNy        DO j=1,sNy
52          DO i=1,sNx          DO i=1,sNx
53    
# Line 73  C  Tabulate following diagnostic fluxes Line 74  C  Tabulate following diagnostic fluxes
74       &           precipo_2D(i,j)*dtatmo*(1. _d 0-iceMask(i,j,1,1))       &           precipo_2D(i,j)*dtatmo*(1. _d 0-iceMask(i,j,1,1))
75              runoff_atm(i,j)= runoff_atm(i,j) +              runoff_atm(i,j)= runoff_atm(i,j) +
76       &           runoff_2D(i,j)*dtatmo       &           runoff_2D(i,j)*dtatmo
77            ENDIF           ENDIF
78    
79          ENDDO          ENDDO
80        ENDDO        ENDDO
81    
 C      PRINT *,'*** bottom calc_1to2d; evapo_2D',evapo_2D(JBUGI,JBUGJ)  
 C      PRINT *,'*** bottom calc_1to2d; precipo_2D',precipo_2D(JBUGI,JBUGJ)  
 C      PRINT *,'*** bottom calc_1to2d; runoff_2D',runoff_2D(JBUGI,JBUGJ)  
 C      PRINT *,'*** bottom calc_1to2d; qneto_2D',qneto_2D(JBUGI,JBUGJ)  
   
82        RETURN        RETURN
83        END        END
84    
# Line 152  C     wght  - weight of this atmos cell Line 148  C     wght  - weight of this atmos cell
148    
149        IF (useAltDeriv) THEN        IF (useAltDeriv) THEN
150          qneto_2D(i,j)= qneto_2D(i,j) + atm_dFdt_ocnq(ind)*          qneto_2D(i,j)= qneto_2D(i,j) + atm_dFdt_ocnq(ind)*
151       &             (sstFromOcn(i,j)-ctocn(ind)*wgt)       &             (sstFromOcn(i,j)-ctocn(ind))*wgt
152          evapo_2D(i,j)= evapo_2D(i,j) + atm_dLdt_ocnq(ind)*          evapo_2D(i,j)= evapo_2D(i,j) + atm_dLdt_ocnq(ind)*
153       &             (sstFromOcn(i,j)-ctocn(ind)*wgt)       &             (sstFromOcn(i,j)-ctocn(ind))*wgt
154          IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN          IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN
155            qneti_2D(i,j)=qneti_2D(i,j)+atm_dFdt_iceq(ind)*            qneti_2D(i,j)=qneti_2D(i,j)+atm_dFdt_iceq(ind)*
156       &             (Tsrf(i,j,1,1)-ctice(ind)*wgt)       &             (Tsrf(i,j,1,1)-ctice(ind))*wgt
157            evapi_2D(i,j)=evapi_2D(i,j)+atm_dLdt_iceq(ind)*            evapi_2D(i,j)=evapi_2D(i,j)+atm_dLdt_iceq(ind)*
158       &             (Tsrf(i,j,1,1)-ctice(ind)*wgt)       &             (Tsrf(i,j,1,1)-ctice(ind))*wgt
159          ENDIF          ENDIF
160        ELSE        ELSE
161          qneto_2D(i,j)= qneto_2D(i,j) + atm_dFdt_ocn(ind)*          qneto_2D(i,j)= qneto_2D(i,j) + atm_dFdt_ocn(ind)*
162       &             (sstFromOcn(i,j)-ctocn(ind)*wgt)       &             (sstFromOcn(i,j)-ctocn(ind))*wgt
163          evapo_2D(i,j)= evapo_2D(i,j) + atm_dLdt_ocn(ind)*          evapo_2D(i,j)= evapo_2D(i,j) + atm_dLdt_ocn(ind)*
164       &             (sstFromOcn(i,j)-ctocn(ind)*wgt)       &             (sstFromOcn(i,j)-ctocn(ind))*wgt
165          IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN          IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN
166            qneti_2D(i,j)= qneti_2D(i,j) + atm_dFdt_ice(ind)*            qneti_2D(i,j)= qneti_2D(i,j) + atm_dFdt_ice(ind)*
167       &             (Tsrf(i,j,1,1)-ctice(ind)*wgt)       &             (Tsrf(i,j,1,1)-ctice(ind))*wgt
168            evapi_2D(i,j)= evapi_2D(i,j)+atm_dLdt_ice(ind)*            evapi_2D(i,j)= evapi_2D(i,j)+atm_dLdt_ice(ind)*
169       &             (Tsrf(i,j,1,1)-ctice(ind)*wgt)       &             (Tsrf(i,j,1,1)-ctice(ind))*wgt
170          ENDIF          ENDIF
171        ENDIF        ENDIF
172    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22