Changeset 1353

Timestamp:

Oct 6, 2014, 6:31:06 PM (10 years ago)

Author:

tnavarro

Message:

New option dustiropacity in callphys.def to change the reference IR opacity of dust + New output dsodust (density-scaled opacity). Without the use of this option, nothing changes for the uninformed user.

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/aeropacity.F (modified) (4 diffs)
• libf/phymars/callkeys.h (modified) (3 diffs)
• libf/phymars/callradite.F (modified) (3 diffs)
• libf/phymars/conf_phys.F (modified) (1 diff)
• libf/phymars/physiq.F (modified) (3 diffs)
• libf/phymars/suaer.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -2140 +2140 @@
 - Bug fix: Enforce recomputation of reaction rates at every call of
   routine prodsandlosses (in paramfoto_compact.F)
+
+== 06/10/2014 == TN
+- New option dustiropacity in callphys.def to change the reference IR opacity
+  of dust. Default is 'tes' (like before), the new one is 'mcs'. A new output,
+  dsodust, is the density scaled opacity of dust, to be compared with MCS
+  measurements, along with dustiropacity='mcs'. Future evolutions of this option
+  could include other instruments, or be used for water ice.

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

@@ -1 +1 @@
       SUBROUTINE aeropacity(ngrid,nlayer,nq,zday,pplay,pplev,ls,
-     &    pq,tauscaling,tauref,tau,taucloudtes,aerosol,reffrad,nueffrad,
-     &    QREFvis3d,QREFir3d,omegaREFvis3d,omegaREFir3d)
+     &    pq,tauscaling,tauref,tau,taucloudtes,aerosol,dsodust,reffrad,
+     &    nueffrad,QREFvis3d,QREFir3d,omegaREFvis3d,omegaREFir3d)
                                                    
 ! to use  'getin'
@@ -340 +340 @@
      &          pq(ig,cstdustlevel,igcm_dust_mass) *
      &          ( pplev(ig,l) - pplev(ig,l+1) ) / g
+                ! DENSITY SCALED OPACITY IN INFRARED:
+                dsodust(ig,l) =
+     &          (  0.75 * QREFir3d(ig,cstdustlevel,iaer) /
+     &          ( rho_dust * reffrad(ig,cstdustlevel,iaer) )  ) *
+     &          pq(ig,cstdustlevel,igcm_dust_mass)
               ENDDO
             ELSE
@@ -348 +353 @@
      &          pq(ig,l,igcm_dust_mass) *
      &          ( pplev(ig,l) - pplev(ig,l+1) ) / g
+                ! DENSITY SCALED OPACITY IN INFRARED:
+                dsodust(ig,l) =
+     &          (  0.75 * QREFir3d(ig,l,iaer) /
+     &          ( rho_dust * reffrad(ig,l,iaer) )  ) *
+     &          pq(ig,l,igcm_dust_mass)
               ENDDO
             ENDIF
@@ -458 +468 @@
      &                aerosol(ig,l,iaerdust(iaer))* tauscaling(ig))
           ENDDO
+        ENDDO
+      ENDDO
+
+      DO l=1,nlayer
+        DO ig=1,ngrid
+          dsodust(ig,l) = max(1E-20,dsodust(ig,l)* tauscaling(ig))
         ENDDO
       ENDDO

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -19 +19 @@
      
       COMMON/callkeys_r/topdustref,solarcondate,semi,alphan,euveff,     &
-     &   tke_heat_flux
+     &   tke_heat_flux,dustrefir
      
       LOGICAL callrad,calldifv,calladj,callcond,callsoil,               &
@@ -28 +28 @@
      &   ,calltherm,callrichsl,callslope,tituscap,callyamada4
 
+      COMMON/aeroutput/dustiropacity
 
       logical callemis
@@ -61 +62 @@
       integer nircorr
 
+      character(len=100) dustiropacity
+      real               dustrefir 
+  
       integer swrtype ! type of short wave (solar wavelength) radiative
       ! transfer to use 1: Fouquart 2: Toon.

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

@@ -2 +2 @@
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
      $     dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
-     &     tauref,tau,aerosol,tauscaling,taucloudtes,rdust,rice,
+     &     tauref,tau,aerosol,dsodust,tauscaling,taucloudtes,rdust,rice,
      &     nuice,co2ice)
 
@@ -180 +180 @@
                                !   (direct comparison with TES)
       REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind)
+      REAL,INTENT(OUT) :: dsodust(ngrid,nlayer)
       REAL,INTENT(OUT) :: rdust(ngrid,nlayer)  ! Dust geometric mean radius (m)
       REAL,INTENT(OUT) :: rice(ngrid,nlayer)   ! Ice geometric mean radius (m)
@@ -373 +374 @@
 c     Computing aerosol optical depth in each layer:
       CALL aeropacity(ngrid,nlayer,nq,zday,pplay,pplev,ls,
-     &      pq,tauscaling,tauref,tau,taucloudtes,aerosol,reffrad,
-     &      nueffrad,QREFvis3d,QREFir3d,omegaREFvis3d,omegaREFir3d)
+     &     pq,tauscaling,tauref,tau,taucloudtes,aerosol,dsodust,reffrad,
+     &     nueffrad,QREFvis3d,QREFir3d,omegaREFvis3d,omegaREFir3d)
 
 c     Starting loop on sub-domain

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

@@ -371 +371 @@
 #endif
 
+! Dust IR opacity
+         write(*,*)" Wavelength for infrared opacity of dust ?"
+         write(*,*)" Choices are:"
+         write(*,*)" tes  --- > 9.3 microns  [default]"
+         write(*,*)" mcs  --- > 21.6 microns"
+         !
+         ! WARNING WARNING WARNING WARNING WARNING WARNING
+         !
+         ! BEFORE ADDING A NEW VALUE, BE SURE THAT THE
+         ! CORRESPONDING WAVELENGTH IS IN THE LOOKUP TABLE,
+         ! OR AT LEAST NO TO FAR, TO AVOID FALLACIOUS INTERPOLATIONS.
+         !
+         dustiropacity="tes" !default value - is expected to shift to mcs one day
+         call getin("dustiropacity",dustiropacity)
+         write(*,*)" dustiropacity = ",trim(dustiropacity)
+         select case (trim(dustiropacity))
+           case ("tes")
+             dustrefir = 9.3E-6
+           case ("mcs")
+             dustrefir = 21.6E-6
+           case default
+              write(*,*) trim(dustiropacity),
+     &                  " is not a valid option for dustiropacity"
+             stop
+         end select
+
 ! callddevil
          write(*,*)" dust lifted by dust devils ?"

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

@@ -231 +231 @@
       REAL tau(ngrid,naerkind)     ! Column dust optical depth at each point
                                    ! AS: TBD: this one should be in a module !
+      REAL dsodust(ngrid,nlayer)   ! density-scaled opacity (in infrared) 
       REAL zls                       !  solar longitude (rad)
       REAL zday                      ! date (time since Ls=0, in martian days)
@@ -669 +670 @@
      $     emis,mu0,zplev,zplay,pt,tsurf,fract,dist_sol,igout,
      $     zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
-     $     tauref,tau,aerosol,tauscaling,taucloudtes,rdust,rice,
+     $     tauref,tau,aerosol,dsodust,tauscaling,taucloudtes,rdust,rice,
      $     nuice,co2ice)
 
@@ -2183 +2184 @@
              call WRITEDIAGFI(ngrid,'dustN','Dust number',
      &                        'part/kg',3,ndust)
+             call WRITEDIAGFI(ngrid,'dsodust',
+     &                        'dust density scaled opacity',
+     &                        'm2.kg-1',3,dsodust)
 c             call WRITEDIAGFI(ngrid,"tauscaling",
 c     &                    "dust conversion factor"," ",2,tauscaling)

trunk/LMDZ.MARS/libf/phymars/suaer.F90

@@ -114 +114 @@
 !                     For dust: change readtesassim accordingly;
 !       Reference wavelength in the infrared:
-        longrefir(iaer)=9.3E-6
+        longrefir(iaer)=dustrefir
 !==================================================================
         CASE("dust_doubleq") aerkind! Two-moment scheme for dust
@@ -127 +127 @@
 !       If not equal to 0.67e-6 -> change readtesassim accordingly;
 !       Reference wavelength in the infrared:
-        longrefir(iaer)=9.3E-6
+        longrefir(iaer)=dustrefir
 !==================================================================
         CASE("dust_submicron") aerkind   ! Small dust population
@@ -139 +139 @@
 !       If not equal to 0.67e-6 -> change readtesassim accordingly;
 !       Reference wavelength in the infrared:
-        longrefir(iaer)=9.3E-6
+        longrefir(iaer)=dustrefir
 !==================================================================
         CASE("h2o_ice") aerkind             ! Water ice crystals