Changeset 2169 for LMDZ5/trunk/libf/phylmd

Timestamp:

Dec 15, 2014, 4:47:34 PM (10 years ago)

Author:

acozic

Message:

1- correct a missing in commit 2146 in infotrac
2- modify the declaration of solsym to use the same variable in differents routines

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phytrac_mod.F90 (modified) (4 diffs)
• rrtm/readaerosol_optic_rrtm.F90 (modified) (2 diffs)
• tracinca_mod.F90 (modified) (4 diffs)
• traclmdz_mod.F90 (modified) (2 diffs)
• tracreprobus_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phytrac_mod.F90

@@ -308 +308 @@
 !$OMP THREADPRIVATE(lessivage)
 
-    CHARACTER(len=8),DIMENSION(nbtr) :: solsym
     !RomP >>>
     INTEGER,SAVE  :: iflag_lscav_omp,iflag_lscav
@@ -557 +556 @@
             cdragh,  coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon,iflag_vdf_trac>=0,sh, &
             rh, pphi, ustar, wstar, ale_bl, ale_wake,  u10m, v10m, &
-            tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
+            tr_seri, source, d_tr_cl,d_tr_dec, zmasse)               !RomP
 
     CASE('inca')
@@ -572 +571 @@
             tau_aero, piz_aero, cg_aero,        ccm,       &
             rfname,                                        &
-            tr_seri,  source,   solsym)      
+            tr_seri,  source)      
 
     CASE('repr')
@@ -580 +579 @@
             presnivs, xlat, xlon, pphis, pphi, &
             t_seri, pplay, paprs, sh , &
-            tr_seri, solsym)
+            tr_seri)
 
     END SELECT

LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

@@ -84 +84 @@
   REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer
   REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer_pi !RAF  
-  !  REAL, DIMENSION(klon,naero_tot)      :: fractnat_allaer !RAF delete??
-  character(len=8), dimension(nbtr) :: tracname
+
   integer :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM
   integer :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M
@@ -104 +103 @@
      !--convert to ug m-3 unit for consistency with offline fields
      !
-#ifdef INCA
-     call tracinca_name(tracname)
-#endif
-
      do i=1,nbtr
-        select case(trim(tracname(i)))
+        select case(trim(solsym(i)))
            case ("ASBCM")
               id_ASBCM = i

LMDZ5/trunk/libf/phylmd/tracinca_mod.F90

@@ -35 +35 @@
        tau_aero, piz_aero, cg_aero,        ccm,       &
        rfname,                                        &
-       tr_seri,  source,   solsym)      
+       tr_seri,  source)      
 
 !========================================================
@@ -113 +113 @@
   ! Output arguments
     REAL,DIMENSION(klon,nbtr), INTENT(OUT)        :: source  ! a voir lorsque le flux de surface est prescrit 
-    CHARACTER(len=8),DIMENSION(nbtr), INTENT(OUT) :: solsym
 
 !=======================================================================================
@@ -134 +133 @@
        pdel(:,k) = paprs(:,k) - paprs (:,k+1)
     END DO
-    
-    zpmfu(:,:)=pmfu(:,:)
+  
+#ifdef INCA
+    IF (config_inca == 'aero') THEN 
+       zpmfu(:,:)=pmfu(:,:)       
+    ELSE IF (config_inca == 'aeNP') THEN
+       zpmfu(:,:)=upwd(:,:)
+    ENDIF
 
-    IF (config_inca == 'aero') THEN
-#ifdef INCA
-       CALL aerosolmain(                    &
-            aerosol_couple,tr_seri,pdtphys, &
-            pplay,pdel,prfl,pmflxr,psfl,    &
-            pmflxs,zpmfu,itop_con,ibas_con,  &
-            pphi,airephy,nstep,rneb,t_seri, &      
-            rh,tau_aero,piz_aero,cg_aero,   &
-            rfname,ccm,lafin)
+    CALL aerosolmain(                    &
+         aerosol_couple,tr_seri,pdtphys, &
+         pplay,pdel,prfl,pmflxr,psfl,    &
+         pmflxs,zpmfu,itop_con,ibas_con,  &
+         pphi,airephy,nstep,rneb,t_seri, &      
+         rh,tau_aero,piz_aero,cg_aero,   &
+         rfname,ccm,lafin, config_inca)
 #endif
-    END IF
-
-    IF (config_inca == 'aeNP') THEN
-#ifdef INCA
-       zpmfu(:,:)=upwd(:,:)
-       CALL aerosolmainNP(                  &
-            aerosol_couple,tr_seri,pdtphys, &
-            pplay,pdel,prfl,pmflxr,psfl,    &
-            pmflxs,zpmfu,itop_con,ibas_con,  &
-            pphi,airephy,nstep,rneb,t_seri, &      
-            rh,lafin)
-#endif
-    END IF
 
 
@@ -196 +185 @@
          iip1,       & !nx
          jjp1,       & !ny
-         source,     &
-         solsym)
+         source )
 #endif
     

LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90

@@ -338 +338 @@
        cdragh,  coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon, couchelimite, sh, &
        rh, pphi, ustar, wstar, ale_bl, ale_wake,  zu10m, zv10m, &
-!!          tr_seri, source, solsym, d_tr_cl, zmasse)                      !RomP
-          tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
+       tr_seri, source, d_tr_cl,d_tr_dec, zmasse)               !RomP
     
     USE dimphy
@@ -397 +396 @@
 
 ! Output argument
-    CHARACTER(len=8),DIMENSION(nbtr), INTENT(OUT) :: solsym
     REAL,DIMENSION(klon,nbtr), INTENT(OUT)        :: source  ! a voir lorsque le flux de surface est prescrit 
     REAL,DIMENSION(klon,klev,nbtr), INTENT(OUT)   :: d_tr_cl ! Td couche limite/traceur

LMDZ5/trunk/libf/phylmd/tracreprobus_mod.F90

@@ -9 +9 @@
        presnivs, xlat, xlon, pphis, pphi, &
        t_seri, pplay, paprs, sh , &
-       tr_seri, solsym)
+       tr_seri)
 
     USE dimphy
@@ -42 +42 @@
 !----------------
     REAL,DIMENSION(klon,klev,nbtr),INTENT(INOUT)  :: tr_seri ! Concentration Traceur [U/KgA]  
-    CHARACTER(len=8),DIMENSION(nbtr), INTENT(OUT) :: solsym