REAL fluxsurf_lw(ngrid)      !incident LW (IR) surface flux (W.m-2)
      REAL fluxsurf_dn_sw(ngrid,2) ! Incident SW (solar) surface flux (W.m-2)
      REAL fluxsurf_up_sw(ngrid,2) ! Reflected SW (solar) surface flux (W.m-2)
      REAL fluxtop_lw(ngrid)       !Outgoing LW (IR) flux to space (W.m-2)
      REAL fluxtop_dn_sw(ngrid,2) ! Incoming SW (solar) flux from space (W.m-2)
      REAL fluxtop_up_sw(ngrid,2) ! Outgoing SW (solar) flux to space (W.m-2)

辐射计算中水汽的影响

The north (NPH) and south (SPH) polar hood clouds of Mars are associated with very cold atmospheric temperatures at mid to high latitudes during fall–winter–spring seasons.(Haberle et al., 2017)

为什么NPH水冰云和水汽输送的季节相反？

    seq  = exp(2.*sig(zt(ig,l))*mh2o / (rho_ice*rgp*zt(ig,l)*
     &             max(rice(ig,l),1.e-7)))
! sig ~ 0.1  mh2o=0.018  rho_ice=920  rgp=8.31
! rice < 5e-4  1e-5~1e-4
    
        dMice =
     & (zq(ig,l,igcm_h2o_vap)-seq*zqsat(ig,l))
     &   /(res*zqsat(ig,l)/(cste*No*rice(ig,l)) + 1.)
    
! cste = 4*pi*rho_ice*microtimestep ~1.42e6   
! microtimestep = ptimestep/15 ~ 123
! zqsat ~ 1.5e-4
! res ~ 2e11
! No   = zq(ig,l,igcm_ccn_number)* tauscaling(ig) + 1.e-30
! tauscaling ~ 0.003
     
       zq(ig,l,igcm_h2o_ice) = zq(ig,l,igcm_h2o_ice)+dMice
       zq(ig,l,igcm_h2o_vap) = zq(ig,l,igcm_h2o_vap)-dMice