Changeset 1353 for trunk/LMDZ.MARS/libf/phymars/conf_phys.F

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.

File:

• trunk/LMDZ.MARS/libf/phymars/conf_phys.F (modified) (1 diff)

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 ?"