Changeset 2033 for trunk

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

Files:

• README (modified) (1 diff)
• libf/phystd/aeroptproperties.F90 (modified) (7 diffs)
• libf/phystd/radcommon_h.F90 (modified) (1 diff)
• libf/phystd/suaer_corrk.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/README

@@ -1404 +1404 @@
   once the value of L_NGAUSS has been set.
 - turned optci, optcv and callcorrk into modules in the process.
+
+== 30/10/2018 == EM
+- 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.

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