Changeset 3418 for LMDZ6/trunk

Timestamp:

Dec 5, 2018, 3:21:39 PM (6 years ago)

Author:

acozic

Message:

Add modification to take into account value read in INCA restart file

Location:

LMDZ6/trunk/libf/phylmd

Files:

• physiq_mod.F90 (modified) (3 diffs)
• phytrac_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1144 +1144 @@
     REAL, dimension(klon, klev) :: cldfrarad   ! fraction nuageuse
 
+#ifdef INCA
+    ! set de variables utilisees pour l'initialisation des valeurs provenant de INCA
+    REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_tauinca
+    REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_pizinca
+    REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_cginca
+    REAL, DIMENSION(klon,klev,nbands) :: init_ccminca
+#endif
+    REAL, DIMENSION(klon,nbtr) :: init_source
+
     !lwoff=y : offset LW CRE for radiation code and other schemes
     REAL, SAVE :: betalwoff 
@@ -1763 +1772 @@
                io_lon, &
                io_lat, &
-               chemistry_couple)
+               chemistry_couple, &
+               init_source, &
+               init_tauinca, &
+               init_pizinca, &
+               init_cginca, &
+               init_ccminca)
+
+
+          ! initialisation des variables depuis le restart de inca
+          ccm(:,:,:) = init_ccminca
+          tau_aero(:,:,:,:) = init_tauinca
+          piz_aero(:,:,:,:) = init_pizinca
+          cg_aero(:,:,:,:) = init_cginca
+!          
+
 
           CALL VTe(VTinca)
@@ -4545 +4568 @@
          rfname, &
          d_tr_dyn, &                                 !<<RomP
-         tr_seri)
+         tr_seri, init_source)
 #endif
     ENDIF    ! (iflag_phytrac=1)

LMDZ6/trunk/libf/phylmd/phytrac_mod.F90

@@ -77 +77 @@
        rfname,                                        &
        d_tr_dyn,                                      &   ! RomP
-       tr_seri)         
+       tr_seri, init_source)         
     !
     !======================================================================
@@ -256 +256 @@
     REAL,DIMENSION(klon,klev,nbtr),INTENT(INOUT) :: tr_seri ! Concentration Traceur [U/KgA]
     REAL,DIMENSION(klon,klev)                    :: sourceBE 
+    REAL,DIMENSION(klon,nbtr), INTENT(IN) :: init_source
 
     !=======================================================================================
@@ -472 +473 @@
           CALL traclmdz_init(pctsrf,xlat,xlon,ftsol,tr_seri,t_seri,pplay,sh,pdtphys,aerosol,lessivage)
        CASE('inca')
-          source(:,:)=0.
+          source(:,:)=init_source(:,:)
           CALL tracinca_init(aerosol,lessivage)
        CASE('repr')