Ignore:
Timestamp:
Jul 17, 2018, 4:26:24 PM (6 years ago)
Author:
mvals
Message:

Mars GCM:
Cosmetic changes:

  • variables suffixes for sub-grid scale cloud scheme changed from --1 to --clf

MV

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

    r1969 r1972  
    401401      logical clearsky
    402402      ! 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)
    413412
    414413c=======================================================================
     
    744743               CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,
    745744     &              albedo,emis,mu0,zplev,zplay,pt,tsurf,fract,
    746      &              dist_sol,igout,zdtlw1,zdtsw1,fluxsurf_lw1,
    747      &              fluxsurf_sw1,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,
    749748     &              rice,nuice,co2ice, clearsky, totcloudfrac)
    750749               clearsky = .false.  ! just in case.
     
    754753                  tf=totcloudfrac(ig)
    755754                  ntf=1.-tf
    756                   fluxsurf_lw(ig) = ntf*fluxsurf_lw1(ig)
     755                  fluxsurf_lw(ig) = ntf*fluxsurf_lwclf(ig)
    757756     &                                          + tf*fluxsurf_lw(ig)
    758                   fluxsurf_sw(ig,1) = ntf*fluxsurf_sw1(ig,1)
     757                  fluxsurf_sw(ig,1) = ntf*fluxsurf_swclf(ig,1)
    759758     &                                          + tf*fluxsurf_sw(ig,1)
    760                   fluxsurf_sw(ig,2) = ntf*fluxsurf_sw1(ig,2)
     759                  fluxsurf_sw(ig,2) = ntf*fluxsurf_swclf(ig,2)
    761760     &                                          + tf*fluxsurf_sw(ig,2)
    762                   fluxtop_lw(ig)  = ntf*fluxtop_lw1(ig)
     761                  fluxtop_lw(ig)  = ntf*fluxtop_lwclf(ig)
    763762     &                                          + tf*fluxtop_lw(ig)
    764                   fluxtop_sw(ig,1)  = ntf*fluxtop_sw1(ig,1) 
     763                  fluxtop_sw(ig,1)  = ntf*fluxtop_swclf(ig,1) 
    765764     &                                          + tf*fluxtop_sw(ig,1)
    766                   fluxtop_sw(ig,2)  = ntf*fluxtop_sw1(ig,2) 
     765                  fluxtop_sw(ig,2)  = ntf*fluxtop_swclf(ig,2) 
    767766     &                                          + tf*fluxtop_sw(ig,2)
    768                   taucloudtes(ig) = ntf*taucloudtes1(ig)
     767                  taucloudtes(ig) = ntf*taucloudtesclf(ig)
    769768     &                                          + tf*taucloudtes(ig)
    770                   zdtlw(ig,1:nlayer) = ntf*zdtlw1(ig,1:nlayer)
     769                  zdtlw(ig,1:nlayer) = ntf*zdtlwclf(ig,1:nlayer)
    771770     &                                         + tf*zdtlw(ig,1:nlayer)
    772                   zdtsw(ig,1:nlayer) = ntf*zdtsw1(ig,1:nlayer)
     771                  zdtsw(ig,1:nlayer) = ntf*zdtswclf(ig,1:nlayer)
    773772     &                                         + tf*zdtsw(ig,1:nlayer)
    774773               ENDDO
Note: See TracChangeset for help on using the changeset viewer.