Changeset 784 for trunk/LMDZ.MARS/libf/phymars/aeroptproperties.F

Timestamp:

Sep 17, 2012, 4:13:22 PM (12 years ago)

Author:

jbmadeleine

Message:

Minor change: added optional outputs for the 1D model in
aeroptproperties.F.

File:

• trunk/LMDZ.MARS/libf/phymars/aeroptproperties.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/aeroptproperties.F

@@ -179 +179 @@
 
       LOGICAL,SAVE :: out_qwg = .false.
-      INTEGER, PARAMETER :: out_iaer = 1
+      INTEGER, PARAMETER :: out_iaer = 2
+      INTEGER :: out_ndim
       REAL :: out_qext(ngridmx,nlayermx)
       REAL :: out_omeg(ngridmx,nlayermx)
@@ -1183 +1184 @@
 c=====Radiative properties - TESTS=================================
       IF (out_qwg) THEN
-      IF (ngrid.NE.1) THEN
+c     -------------------------------------------------------------
+        IF (ngrid.NE.1) THEN
+          out_ndim = 3
+        ELSE
+          out_ndim = 1
+        ENDIF
+c     -------------------------------------------------------------
         DO out_nchannel = 1, 2
 c     -------------------------------------------------------------
@@ -1198 +1205 @@
           write(out_str(1:1),'(i1.1)') out_nchannel
           call WRITEDIAGFI(ngrid,'qextvis'//out_str,"Ext.efficiency","",
-     &                     3,out_qext)
+     &                     out_ndim,out_qext)
           call WRITEDIAGFI(ngrid,'omegvis'//out_str,"Sing.Scat.Alb.","",
-     &                     3,out_omeg)
+     &                     out_ndim,out_omeg)
           call WRITEDIAGFI(ngrid,'gvis'//out_str,"Asym.Factor","",
-     &                     3,out_g)
+     &                     out_ndim,out_g)
 c     -------------------------------------------------------------
         ENDDO ! out_nchannel 
@@ -1219 +1226 @@
           write(out_str(1:1),'(i1.1)') out_nchannel
        call WRITEDIAGFI(ngrid,'qextir'//out_str,"Ext.efficiency","",
-     &                     3,out_qext)
+     &                     out_ndim,out_qext)
        call WRITEDIAGFI(ngrid,'omegir'//out_str,"Sing.Scat.Alb.","",
-     &                     3,out_omeg)
+     &                     out_ndim,out_omeg)
        call WRITEDIAGFI(ngrid,'gir'//out_str,"Asym.Factor","",
-     &                     3,out_g)
+     &                     out_ndim,out_g)
 c     -------------------------------------------------------------
         ENDDO ! out_nchannel 
         call WRITEDIAGFI(ngrid,"omegvisref","Sing.Scat.Alb.","",
-     &                   3,omegaREFvis3d(1,1,out_iaer))
+     &                   out_ndim,omegaREFvis3d(1,1,out_iaer))
         call WRITEDIAGFI(ngrid,"omegirref","Sing.Scat.Alb.","",
-     &                   3,omegaREFir3d(1,1,out_iaer))
-      ENDIF ! ngrid.EQ.1
+     &                   out_ndim,omegaREFir3d(1,1,out_iaer))
       ENDIF ! out_qwg
 c==================================================================