Changeset 5453 for LMDZ6/branches/contrails/libf/phylmd/physiq_mod.F90

Timestamp:

Dec 23, 2024, 8:19:39 PM (37 hours ago)

Author:

aborella

Message:

Added water emissions and IO routines for contrails.
Work to be done: contrails initial cross section and radiative transfer

File:

• LMDZ6/branches/contrails/libf/phylmd/physiq_mod.F90 (modified) (5 diffs)

LMDZ6/branches/contrails/libf/phylmd/physiq_mod.F90

@@ -75 +75 @@
     USE write_field_phy
     use wxios_mod, ONLY: g_ctx, wxios_set_context
+    USE lmdz_aviation, ONLY : read_aviation_emissions, aviation_water_emissions
     USE lmdz_lscp, ONLY : lscp
     USE lmdz_call_cloud_optics_prop, ONLY : call_cloud_optics_prop
@@ -327 +328 @@
        dqvc_adj, dqvc_sub, dqvc_con, dqvc_mix, qsatliq, qsatice, &
        !-- LSCP - aviation and contrails variables
-       rcont_seri, d_rcont_dyn, flight_dist, flight_h2o, &
+       d_q_avi, rcont_seri, d_rcont_dyn, flight_dist, flight_h2o, &
        Tcritcont, qcritcont, potcontfraP, potcontfraNP, &
        dcf_avi, dqi_avi, dqvc_avi, &
@@ -3863 +3864 @@
          ratqs,ratqsc,ratqs_inter_,sigma_qtherm)
 
+    !--Read the aviation emissions
+    IF ( ok_plane_h2o .OR. ok_plane_contrail ) THEN
+      CALL read_aviation_emissions(klon, klev, latitude_deg, longitude_deg, pplay, &
+                                   flight_dist, flight_h2o)
+    ENDIF
+
+    !--Add the water emissions from aviation
+    IF ( ok_plane_h2o ) THEN
+       CALL aviation_water_emissions(klon, klev, phys_tstep, paprs, pplay, &
+            t_seri, q_seri, cell_area, flight_h2o, d_q_avi)
+       CALL add_phys_tend(du0, dv0, dt0, d_q_avi, dql0, dqi0, dqbs0, paprs, &
+            'avi', abortphy, flag_inhib_tend, itap, 0)
+       d_q_avi = d_q_avi / phys_tstep
+    ENDIF
+
     !
     ! Appeler le processus de condensation a grande echelle
@@ -3875 +3891 @@
 
     IF (ok_new_lscp) THEN
-
  
     DO k = 1, klev
@@ -3883 +3898 @@
       ENDDO
     ENDDO
-
-
-    !--mise à jour de flight_m et flight_h2o dans leur module
-    !IF (ok_plane_h2o .OR. ok_plane_contrail) THEN
-    !  CALL airplane(debut,pphis,pplay,paprs,t_seri)
-    !ENDIF
 
     CALL lscp(klon,klev,phys_tstep,missing_val,paprs,pplay,omega, &