Changeset 784 for trunk/LMDZ.MARS/libf


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/aeroptproperties.F

    r737 r784  
    179179
    180180      LOGICAL,SAVE :: out_qwg = .false.
    181       INTEGER, PARAMETER :: out_iaer = 1
     181      INTEGER, PARAMETER :: out_iaer = 2
     182      INTEGER :: out_ndim
    182183      REAL :: out_qext(ngridmx,nlayermx)
    183184      REAL :: out_omeg(ngridmx,nlayermx)
     
    11831184c=====Radiative properties - TESTS=================================
    11841185      IF (out_qwg) THEN
    1185       IF (ngrid.NE.1) THEN
     1186c     -------------------------------------------------------------
     1187        IF (ngrid.NE.1) THEN
     1188          out_ndim = 3
     1189        ELSE
     1190          out_ndim = 1
     1191        ENDIF
     1192c     -------------------------------------------------------------
    11861193        DO out_nchannel = 1, 2
    11871194c     -------------------------------------------------------------
     
    11981205          write(out_str(1:1),'(i1.1)') out_nchannel
    11991206          call WRITEDIAGFI(ngrid,'qextvis'//out_str,"Ext.efficiency","",
    1200      &                     3,out_qext)
     1207     &                     out_ndim,out_qext)
    12011208          call WRITEDIAGFI(ngrid,'omegvis'//out_str,"Sing.Scat.Alb.","",
    1202      &                     3,out_omeg)
     1209     &                     out_ndim,out_omeg)
    12031210          call WRITEDIAGFI(ngrid,'gvis'//out_str,"Asym.Factor","",
    1204      &                     3,out_g)
     1211     &                     out_ndim,out_g)
    12051212c     -------------------------------------------------------------
    12061213        ENDDO ! out_nchannel
     
    12191226          write(out_str(1:1),'(i1.1)') out_nchannel
    12201227       call WRITEDIAGFI(ngrid,'qextir'//out_str,"Ext.efficiency","",
    1221      &                     3,out_qext)
     1228     &                     out_ndim,out_qext)
    12221229       call WRITEDIAGFI(ngrid,'omegir'//out_str,"Sing.Scat.Alb.","",
    1223      &                     3,out_omeg)
     1230     &                     out_ndim,out_omeg)
    12241231       call WRITEDIAGFI(ngrid,'gir'//out_str,"Asym.Factor","",
    1225      &                     3,out_g)
     1232     &                     out_ndim,out_g)
    12261233c     -------------------------------------------------------------
    12271234        ENDDO ! out_nchannel
    12281235        call WRITEDIAGFI(ngrid,"omegvisref","Sing.Scat.Alb.","",
    1229      &                   3,omegaREFvis3d(1,1,out_iaer))
     1236     &                   out_ndim,omegaREFvis3d(1,1,out_iaer))
    12301237        call WRITEDIAGFI(ngrid,"omegirref","Sing.Scat.Alb.","",
    1231      &                   3,omegaREFir3d(1,1,out_iaer))
    1232       ENDIF ! ngrid.EQ.1
     1238     &                   out_ndim,omegaREFir3d(1,1,out_iaer))
    12331239      ENDIF ! out_qwg
    12341240c==================================================================
Note: See TracChangeset for help on using the changeset viewer.