Changeset 2033 for trunk/LMDZ.GENERIC/libf

Timestamp:

Oct 30, 2018, 2:38:18 PM (6 years ago)

Author:

emillour

Message:

Generic GCM:

• Cleanup around aeroptproperties.90, remove arrays that were not initialized but still used, probably for nothing since removing them leads to no change in results. Possibly to revisit and further clean up later.

EM

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• aeroptproperties.F90 (modified) (7 diffs)
• radcommon_h.F90 (modified) (1 diff)
• suaer_corrk.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/aeroptproperties.F90

@@ -7 +7 @@
       use radinc_h,    only: L_NSPECTI,L_NSPECTV,nsizemax,naerkind
       use radcommon_h, only: QVISsQREF,omegavis,gvis,QIRsQREF,omegair,gir
-      use radcommon_h, only: qrefvis,qrefir,omegarefvis,omegarefir
+      use radcommon_h, only: qrefvis,qrefir,omegarefir !,omegarefvis
       use radcommon_h, only: radiustab,nsize
 
@@ -189 +189 @@
       REAL,INTENT(OUT) :: QREFir3d(ngrid,nlayer,naerkind)
 
-      REAL :: omegaREFvis3d(ngrid,nlayer,naerkind)
-      REAL :: omegaREFir3d(ngrid,nlayer,naerkind)
+!      REAL :: omegaREFvis3d(ngrid,nlayer,naerkind)
+!      REAL :: omegaREFir3d(ngrid,nlayer,naerkind)
 
       DO iaer = 1, naerkind ! Loop on aerosol kind
@@ -211 +211 @@
               QREFvis3d(ig,lg,iaer)=QREFvis(iaer,1)
               QREFir3d(ig,lg,iaer)=QREFir(iaer,1)
-              omegaREFvis3d(ig,lg,iaer)=omegaREFvis(iaer,1)
-              omegaREFir3d(ig,lg,iaer)=omegaREFir(iaer,1)
+!              omegaREFvis3d(ig,lg,iaer)=omegaREFvis(iaer,1)
+!              omegaREFir3d(ig,lg,iaer)=omegaREFir(iaer,1)
             ENDDO
           ENDDO
@@ -296 +296 @@
                     (1-kint)*QREFvis(iaer,radius_id) +          &
                     kint*QREFvis(iaer,radius_id+1)
-            omegrefVISa(gausind,iaer)=                          &
-                    (1-kint)*omegaREFvis(iaer,radius_id) +      &
-                    kint*omegaREFvis(iaer,radius_id+1)
+            omegrefVISa(gausind,iaer)= 0
+!            omegrefVISa(gausind,iaer)=                          &
+!                    (1-kint)*omegaREFvis(iaer,radius_id) +      &
+!                    kint*omegaREFvis(iaer,radius_id+1)
           ELSE ! INFRARED DOMAIN ----------------------------------
             DO m=1,L_NSPECTI
@@ -357 +358 @@
                     (1-kint)*QREFvis(iaer,radius_id) +          &
                     kint*QREFvis(iaer,radius_id+1)
-            omegrefVISb(gausind,iaer)=                          &
-                    (1-kint)*omegaREFvis(iaer,radius_id) +      &
-                    kint*omegaREFvis(iaer,radius_id+1)
+            omegrefVISb(gausind,iaer)= 0
+!            omegrefVISb(gausind,iaer)=                          &
+!                    (1-kint)*omegaREFvis(iaer,radius_id) +      &
+!                    kint*omegaREFvis(iaer,radius_id+1)
           ELSE ! INFRARED DOMAIN ----------------------------------
             DO m=1,L_NSPECTI
@@ -712 +714 @@
                         k1*qrefVISgrid(grid_i,1,iaer) +            &
                         k2*qrefVISgrid(grid_i+1,1,iaer)
-          omegaREFvis3d(ig,lg,iaer) =                              &
-                        k1*omegrefVISgrid(grid_i,1,iaer) +         &
-                        k2*omegrefVISgrid(grid_i+1,1,iaer)
+!          omegaREFvis3d(ig,lg,iaer) =                              &
+!                        k1*omegrefVISgrid(grid_i,1,iaer) +         &
+!                        k2*omegrefVISgrid(grid_i+1,1,iaer)
           ELSE                   ! INFRARED -----------------------
           DO m=1,L_NSPECTI
@@ -730 +732 @@
                         k1*qrefIRgrid(grid_i,1,iaer) +             &
                         k2*qrefIRgrid(grid_i+1,1,iaer)
-          omegaREFir3d(ig,lg,iaer) =                               &
-                        k1*omegrefIRgrid(grid_i,1,iaer) +          &
-                        k2*omegrefIRgrid(grid_i+1,1,iaer)
+!          omegaREFir3d(ig,lg,iaer) =                               &
+!                        k1*omegrefIRgrid(grid_i,1,iaer) +          &
+!                        k2*omegrefIRgrid(grid_i+1,1,iaer)
           ENDIF                  ! --------------------------------
         ENDDO !nlayer

trunk/LMDZ.GENERIC/libf/phystd/radcommon_h.F90

@@ -118 +118 @@
       REAL :: QREFvis(naerkind,nsizemax)
       REAL :: QREFir(naerkind,nsizemax)
-      REAL :: omegaREFvis(naerkind,nsizemax)
+!      REAL :: omegaREFvis(naerkind,nsizemax)
       REAL :: omegaREFir(naerkind,nsizemax)
 

trunk/LMDZ.GENERIC/libf/phystd/suaer_corrk.F90

@@ -9 +9 @@
       use radcommon_h, only: QVISsQREF,omegavis,gvis,QIRsQREF,omegair,gir
       use radcommon_h, only: radiustab,nsize,tstellar
-      use radcommon_h, only: qrefvis,qrefir,omegarefvis,omegarefir
+      use radcommon_h, only: qrefvis,qrefir,omegarefir !,omegarefvis
       use aerosol_mod
       use callkeys_mod, only: tplanet