Changeset 4601 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

Timestamp:

Jul 1, 2023, 12:07:30 AM (11 months ago)

Author:

dcugnet

Message:

StratAer? commit, N. Lebas

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -126 +126 @@
 
 #ifdef CPP_StratAer
-    USE strataer_mod, ONLY: strataer_init
+    USE strataer_local_var_mod
+    USE strataer_nuc_mod, ONLY: strataer_nuc_init
+    USE strataer_emiss_mod, ONLY: strataer_emiss_init
 #endif
 
@@ -1414 +1416 @@
 #ifdef CPP_StratAer
        CALL strataer_init
+       CALL strataer_nuc_init
+       CALL strataer_emiss_init
 #endif
 
@@ -1599 +1603 @@
           piz_aero(:,:,:,:) = 0.
           cg_aero(:,:,:,:) = 0.
+          d_q_ch4(:,:) = 0.
 
           config_inca='none' ! default
@@ -1610 +1615 @@
        piz_aero(:,:,:,:) = 1.
        cg_aero(:,:,:,:)  = 0.
+       d_q_ch4(:,:) = 0.
 
        IF (aerosol_couple .AND. (config_inca /= "aero" &
@@ -6154 +6160 @@
     !
     !
+#ifdef CPP_StratAer
+    IF (ok_qemiss) THEN
+       flh2o=1
+       IF(flag_verbose_strataer) THEN
+          print *,'IN physiq_mod: ok_qemiss =yes (',ok_qemiss,'), flh2o=',flh2o
+          print *,'IN physiq_mod: flag_emit=',flag_emit,', nErupt=',nErupt
+          print *,'IN physiq_mod: nAerErupt=',nAerErupt
+       ENDIF
+       
+       SELECT CASE(flag_emit)
+       CASE(1) ! emission volc H2O dans LMDZ
+          DO ieru=1, nErupt
+             IF (year_cur==year_emit_vol(ieru).AND.&
+                  mth_cur==mth_emit_vol(ieru).AND.&
+                  day_cur>=day_emit_vol(ieru).AND.&
+                  day_cur<(day_emit_vol(ieru)+injdur)) THEN
+                
+                IF(flag_verbose_strataer) print *,'IN physiq_mod: date=',year_cur,mth_cur,day_cur
+                ! initialisation tendance q emission
+                d_q_emiss(:,:)=0.
+                ! daily injection mass emission - NL
+                m_H2O_emiss_vol_daily = m_H2O_emiss_vol(ieru)/(REAL(injdur)&
+                     *REAL(ponde_lonlat_vol(ieru)))
+                !
+                CALL STRATEMIT(pdtphys,pdtphys,latitude_deg,longitude_deg,t_seri,&
+                     pplay,paprs,tr_seri,&
+                     m_H2O_emiss_vol_daily,&
+                     xlat_min_vol(ieru),xlat_max_vol(ieru),&
+                     xlon_min_vol(ieru),xlon_max_vol(ieru),&
+                     altemiss_vol(ieru),&
+                     sigma_alt_vol(ieru),1,&
+                     1,nAerErupt+1,0)
+                
+                IF(flag_verbose_strataer) print *,'IN physiq_mod: min max d_q_emiss=',&
+                     minval(d_q_emiss),maxval(d_q_emiss)
+                
+                CALL add_phys_tend(du0, dv0, dt0, d_q_emiss, dql0, dqi0, dqbs0, paprs, &
+                     'q_emiss',abortphy,flag_inhib_tend,itap,0)
+                IF (abortphy==1) Print*,'ERROR ABORT TEND EMISS'
+             ENDIF
+          ENDDO
+          flh2o=0
+       END SELECT ! emission scenario (flag_emit)
+    ENDIF
+#endif
 
 !===============================================================
@@ -6744 +6795 @@
       qql2(:)=qql2(:)+(q_seri(:,k)+ql_seri(:,k)+qs_seri(:,k))*zmasse(:,k)
     ENDDO
+    
+#ifdef CPP_StratAer
+    IF (ok_qemiss) THEN
+       DO k = 1, klev
+          qql1(:) = qql1(:)+d_q_emiss(:,k)*zmasse(:,k)
+       ENDDO
+    ENDIF
+#endif
+    IF (ok_qch4) THEN
+       DO k = 1, klev
+          qql1(:) = qql1(:)+d_q_ch4_dtime(:,k)*zmasse(:,k)
+       ENDDO
+    ENDIF
+    
     DO i = 1, klon
       !--compute ratio of what q+ql should be with conservation to what it is