Changeset 3865 for LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

Timestamp:

Mar 23, 2021, 4:14:07 PM (3 years ago)

Author:

lmdz-users

Message:

Modifications from Thibaut to create an ESM with interactive CO2 + INCA aerosols

File:

• LMDZ6/trunk/libf/phylmd/infotrac_phy.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

@@ -20 +20 @@
   INTEGER, SAVE :: nbtr
 !$OMP THREADPRIVATE(nbtr)
+
+! ThL : number of tracers specific to INCA
+  INTEGER, SAVE :: nqINCA
+!$OMP THREADPRIVATE(nqINCA)
 
 #ifdef CPP_StratAer
@@ -33 +37 @@
   INTEGER, SAVE :: nqperes
 !$OMP THREADPRIVATE(nqperes)
+
+! ThL : nb de traceurs dans le traceur.def
+  INTEGER, SAVE :: nqexcl
+!$OMP THREADPRIVATE(nqexcl)
 
 ! Name variables
@@ -96 +104 @@
 CONTAINS
 
-  SUBROUTINE init_infotrac_phy(nqtot_,nqo_,nbtr_,tname_,ttext_,type_trac_,&
+  SUBROUTINE init_infotrac_phy(nqtot_,nqo_,nbtr_,nqINCA_,tname_,ttext_,type_trac_,&
                                niadv_,conv_flg_,pbl_flg_,solsym_,&
                                nqfils_,nqdesc_,nqdesc_tot_,iqfils_,iqpere_,&
@@ -118 +126 @@
     INTEGER,INTENT(IN) :: nqo_
     INTEGER,INTENT(IN) :: nbtr_
+    INTEGER,INTENT(IN) :: nqINCA_
 #ifdef CPP_StratAer
     INTEGER,INTENT(IN) :: nbtr_bin_
@@ -130 +139 @@
     CHARACTER(len=4),INTENT(IN) :: type_trac_
     INTEGER,INTENT(IN) :: niadv_ (nqtot_) ! equivalent dyn / physique
-    INTEGER,INTENT(IN) :: conv_flg_(nbtr_)
-    INTEGER,INTENT(IN) :: pbl_flg_(nbtr_)
-    CHARACTER(len=8),INTENT(IN) :: solsym_(nbtr_)
+    INTEGER,INTENT(IN) :: conv_flg_(nqINCA_)
+    INTEGER,INTENT(IN) :: pbl_flg_(nqINCA_)
+    CHARACTER(len=8),INTENT(IN) :: solsym_(nqINCA_)
     ! Isotopes:
     INTEGER,INTENT(IN) :: nqfils_(nqtot_)
@@ -163 +172 @@
     nqo=nqo_
     nbtr=nbtr_
+    nqINCA=nqINCA_
 #ifdef CPP_StratAer
     nbtr_bin=nbtr_bin_
@@ -179 +189 @@
     niadv(:)=niadv_(:)
     ALLOCATE(conv_flg(nbtr))
-    conv_flg(:)=conv_flg_(:)
+    IF (type_trac == 'inco') THEN
+      conv_flg(1)=1
+      conv_flg(2:nbtr)=conv_flg_(:)
+    ELSE
+      conv_flg(:)=conv_flg_(:)
+    ENDIF
     ALLOCATE(pbl_flg(nbtr))
-    pbl_flg(:)=pbl_flg_(:)
+    IF (type_trac == 'inco') THEN
+      pbl_flg(1)=1
+      pbl_flg(2:nbtr)=pbl_flg_(:)
+    ELSE
+      pbl_flg(:)=pbl_flg_(:)
+    ENDIF
     ALLOCATE(solsym(nbtr))
-    solsym(:)=solsym_(:)
-  
+    IF (type_trac == 'inco') THEN
+      solsym(1)='CO2'
+      solsym(2:nbtr)=solsym_(:)
+    ELSE
+      solsym(:)=solsym_(:)
+    ENDIF
+     
     IF(prt_level.ge.1) THEN
-      write(lunout,*) TRIM(modname)//": nqtot,nqo,nbtr",nqtot,nqo,nbtr
+      write(lunout,*) TRIM(modname)//": nqtot,nqo,nbtr,nqINCA",nqtot,nqo,nbtr,nqINCA
     ENDIF