Changeset 2388


Ignore:
Timestamp:
Jul 2, 2020, 3:42:33 PM (4 years ago)
Author:
cmathe
Message:

mprovedco2clouds_mod: case satu < 1 removed, debug in 4.2 (write)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/improvedco2clouds_mod.F90

    r2387 r2388  
    482482          end if ! of if co2useh2o
    483483        end if ! of if No > 1e-30
    484 
    485         ! WIP: cas où satu < 1 et donc si on a de la glace de co2 et que T > Tcond, alors on relache notre glace en vap
    486       else
    487           if (tcond(ig,l) <= zt(ig,l) .and. zq(ig,l,igcm_co2_ice) >= threshold) then
    488             zq(ig,l,igcm_co2) = zq(ig,l,igcm_co2) + zq(ig,l,igcm_co2_ice)
    489             zq(ig,l,igcm_co2_ice) = 0.
    490           end if
    491484      end if ! of is satu > 1
    492485!----------------------------------------------------------------------------------------------------------------------!
     
    499492      ! we trigger crystal growth
    500493      if (zq(ig,l,igcm_ccnco2_number) * tauscaling(ig) + threshold >= 1.) then
     494        write(*,*)'yes I am in', zq(ig, l, igcm_ccnco2_number), zq(ig, l, igcm_co2_ice), satu
    501495
    502496        call updaterice_microco2(zq(ig,l,igcm_co2_ice), zq(ig,l,igcm_ccnco2_mass), zq(ig,l,igcm_ccnco2_number), &
     
    510504        facteurmax = abs(tcond(ig,l)-zt(ig,l)) * (cpp/lw)
    511505
    512         ! call scheme of microphys. mass growth for CO2
     506        ! call scheme of microphys. mass growth for CO2 (evaporation/condensation)
    513507        call massflowrateCO2(pplay(ig,l), zt(ig,l), satu, riceco2(ig,l), mmean(ig,l), Ic_rice)
    514508
Note: See TracChangeset for help on using the changeset viewer.