Changeset 1972 for trunk/LMDZ.MARS/libf/phymars
- Timestamp:
- Jul 17, 2018, 4:26:24 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
r1969 r1972 401 401 logical clearsky 402 402 ! flux for the part without clouds 403 real zdtsw1(ngrid,nlayer) 404 real zdtlw1(ngrid,nlayer) 405 real fluxsurf_lw1(ngrid) 406 real fluxsurf_sw1(ngrid,2) 407 real fluxtop_lw1(ngrid) 408 real fluxtop_sw1(ngrid,2) 409 REAL taucloudtes1(ngrid) 410 ! tf: fraction of clouds, ntf: fraction without clouds 411 real tf, ntf 412 REAL rave2(ngrid), totrave2(ngrid) ! Mean water ice mean radius (m) 403 real zdtswclf(ngrid,nlayer) 404 real zdtlwclf(ngrid,nlayer) 405 real fluxsurf_lwclf(ngrid) 406 real fluxsurf_swclf(ngrid,2) 407 real fluxtop_lwclf(ngrid) 408 real fluxtop_swclf(ngrid,2) 409 real taucloudtesclf(ngrid) 410 real tf, ntf ! tf: fraction of clouds, ntf: fraction without clouds 411 real rave2(ngrid), totrave2(ngrid) ! Mean water ice mean radius (m) 413 412 414 413 c======================================================================= … … 744 743 CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq, 745 744 & albedo,emis,mu0,zplev,zplay,pt,tsurf,fract, 746 & dist_sol,igout,zdtlw 1,zdtsw1,fluxsurf_lw1,747 & fluxsurf_sw 1,fluxtop_lw1,fluxtop_sw1,tauref,tau,748 & aerosol,dsodust,tauscaling,taucloudtes1,rdust,745 & dist_sol,igout,zdtlwclf,zdtswclf,fluxsurf_lwclf, 746 & fluxsurf_swclf,fluxtop_lwclf,fluxtop_swclf,tauref, 747 & tau,aerosol,dsodust,tauscaling,taucloudtesclf,rdust, 749 748 & rice,nuice,co2ice, clearsky, totcloudfrac) 750 749 clearsky = .false. ! just in case. … … 754 753 tf=totcloudfrac(ig) 755 754 ntf=1.-tf 756 fluxsurf_lw(ig) = ntf*fluxsurf_lw 1(ig)755 fluxsurf_lw(ig) = ntf*fluxsurf_lwclf(ig) 757 756 & + tf*fluxsurf_lw(ig) 758 fluxsurf_sw(ig,1) = ntf*fluxsurf_sw 1(ig,1)757 fluxsurf_sw(ig,1) = ntf*fluxsurf_swclf(ig,1) 759 758 & + tf*fluxsurf_sw(ig,1) 760 fluxsurf_sw(ig,2) = ntf*fluxsurf_sw 1(ig,2)759 fluxsurf_sw(ig,2) = ntf*fluxsurf_swclf(ig,2) 761 760 & + tf*fluxsurf_sw(ig,2) 762 fluxtop_lw(ig) = ntf*fluxtop_lw 1(ig)761 fluxtop_lw(ig) = ntf*fluxtop_lwclf(ig) 763 762 & + tf*fluxtop_lw(ig) 764 fluxtop_sw(ig,1) = ntf*fluxtop_sw 1(ig,1)763 fluxtop_sw(ig,1) = ntf*fluxtop_swclf(ig,1) 765 764 & + tf*fluxtop_sw(ig,1) 766 fluxtop_sw(ig,2) = ntf*fluxtop_sw 1(ig,2)765 fluxtop_sw(ig,2) = ntf*fluxtop_swclf(ig,2) 767 766 & + tf*fluxtop_sw(ig,2) 768 taucloudtes(ig) = ntf*taucloudtes 1(ig)767 taucloudtes(ig) = ntf*taucloudtesclf(ig) 769 768 & + tf*taucloudtes(ig) 770 zdtlw(ig,1:nlayer) = ntf*zdtlw 1(ig,1:nlayer)769 zdtlw(ig,1:nlayer) = ntf*zdtlwclf(ig,1:nlayer) 771 770 & + tf*zdtlw(ig,1:nlayer) 772 zdtsw(ig,1:nlayer) = ntf*zdtsw 1(ig,1:nlayer)771 zdtsw(ig,1:nlayer) = ntf*zdtswclf(ig,1:nlayer) 773 772 & + tf*zdtsw(ig,1:nlayer) 774 773 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.