Ignore:
Timestamp:
Apr 2, 2025, 4:04:40 PM (3 months ago)
Author:
aborella
Message:

Multiple changes:

  • tracers which were ratios are now absolute quantities. This is needed because when the ratio

is not defined, some aberrations may occur

  • added a new tracer for total specific humidity in contrails
  • rework of the mixing process for cirrus clouds (and contrails)
  • changed the numerical integration of ice crystals' sublimation
  • subroutines do not take real inputs anymore (at least klon tables)
  • added more radiative diagnostics for contrails
Location:
LMDZ6/branches/contrails/libf/phylmd/dyn1d
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/contrails/libf/phylmd/dyn1d/old_lmdz1d.f90

    r5368 r5601  
    2222       zgam, zmax0, zmea, zpic, zsig, &
    2323       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, &
    24        ql_ancien, qs_ancien, qbs_ancien, cf_ancien, rvc_ancien, &
     24       ql_ancien, qs_ancien, qbs_ancien, &
     25       cf_ancien, qvc_ancien, cfa_ancien, qta_ancien, &
    2526       prlw_ancien, prsw_ancien, prbsw_ancien, prw_ancien, &
    2627       u10m,v10m,ale_wake,ale_bl_stat
    … …  
    872873        IF ( ok_ice_supersat ) THEN
    873874          cf_ancien = 0.
    874           rvc_ancien = 0.
     875          qvc_ancien = 0.
     876        ENDIF
     877        IF ( ok_plane_contrail ) THEN
     878          cfa_ancien = 0.
     879          qta_ancien = 0.
    875880        ENDIF
    876881!jyg<

  • LMDZ6/branches/contrails/libf/phylmd/dyn1d/scm.f90

    r5450 r5601  
    1818       zgam, zmax0, zmea, zpic, zsig, &
    1919       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, &
    20        ql_ancien, qs_ancien, qbs_ancien, cf_ancien, rvc_ancien, &
     20       ql_ancien, qs_ancien, qbs_ancien, &
     21       cf_ancien, qvc_ancien, cfa_ancien, qta_ancien, &
    2122       prlw_ancien, prsw_ancien, prbsw_ancien, prw_ancien, &
    2223       u10m,v10m,ale_wake,ale_bl_stat, ratqs_inter_
    … …  
    614615        IF ( ok_ice_supersat ) THEN
    615616          cf_ancien = 0.
    616           rvc_ancien = 0.
     617          qvc_ancien = 0.
     618        ENDIF
     619        IF ( ok_plane_contrail ) THEN
     620          cfa_ancien = 0.
     621          qta_ancien = 0.
    617622        ENDIF
    618623        rain_fall=0.
Note: See TracChangeset for help on using the changeset viewer.