Changeset 222

Timestamp:

Jul 14, 2011, 12:56:15 AM (13 years ago)

Author:

aslmd

Message:

MESSAGE DE JB MADELEINE

C'était juste pour vous dire que j'ai fait un peu de rangement dans le
DATAFILE concernant les propriétés de la poussière (directement dans le
répertoire /u/forget/WWW/datagcm/datafile/), car lorsque Ockert-Bell
dans le domaine solaire était utilisé avec Forget 1998 dans l'IR, le
rapport Solaire/IR "solsir" n'était pas rigoureusement égal à 2.

Du coup il faut changer les noms des fichiers dans suaer.F90

REMARQUE: Le rapport "solsir" est sorti dans les logs au firstcall, et
vous verrez qu'il n'est pas égal à 2 lorsque ces propriétés sont
utilisées. C'EST NORMAL! En effet, j'ai réglé la longueur d'onde de
référence de la poussière dans l'IR à 9.3um pour pouvoir comparer
l'opacité donnée en output avec TES. Or le fameux solsir=2 utilisé
partout date de l'époque où la longueur de référence était celle du
canal "T9" d'IRTM qui est 9um et non 9.3um!

En résumé, au firstcall, le GCM affiche Qext(0.67um)/Qext(9.3um), et non
l'ancien solsir qui lui vaut Qext(0.67um)/Qext(9um), d'où un résultat
qui n'est pas égal à deux. Mais rassurez-vous, Qext(0.67um)/Qext(9um)
dans les fichiers d'input vaut bien 2.

Également, j'ai un peu modifié aeropacity.F car il affichait le
"ccn_factor" à tout va dans les logs. Il est également attaché à ce mail.

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/aeropacity.F (modified) (3 diffs)
• libf/phymars/suaer.F90 (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -773 +773 @@
 == 04/07/2011 == AC
 - Minor setting modification to thermals
-- Added new flux optimization in thermcell_main_mars.F90, to run properly
-in 3D
+- Added new flux optimization in thermcell_main_mars.F90, to run properly in 3D
+
+== 14/07/2011 == JBM
+- Tidying up dust properties in DATAFILE for better consistency (cf. suaer.F90)
+- Cosmetic changes in aeropacity.F (changed comment and put a print inside a water flag)

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

@@ -172 +172 @@
         msolsir(1:nlayer,1:naerkind)=0
         mqextsqabs(1:nlayer,1:naerkind)=0
-        WRITE(*,*) "Typical profiles of solsir and Qext/Qabs(IR):"
+        WRITE(*,*) "Typical profiles of Qext(vis)/Qext(IR)"
+        WRITE(*,*) "  and Qext(IR)/Qabs(IR):"
         DO iaer = 1, naerkind ! Loop on aerosol kind
           WRITE(*,*) "Aerosol # ",iaer
@@ -193 +194 @@
 !       otherwise default value read from starfi.nc file will be used)
         call getin("tauvis",tauvis)
+
+!       Some information about the water cycle
+        IF (water) THEN
+          write(*,*) "water_param CCN reduc. fac. ", ccn_factor
+        ENDIF
 
         firstcall=.false.
@@ -490 +496 @@
 !         ccn(ig,l) = ccn(ig,l) / 8.
 c -----------------------------------------------------------------
-       write(*,*) "water_param CCN reduc. fac. ", ccn_factor
        DO l=1,nlayer
          DO ig=1,ngrid

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

@@ -100 +100 @@
 !       Visible domain:
         file_id(iaer,1) = 'optprop_dustvis_TM.dat'     !M.Wolff TM
+!       file_id(iaer,1) = 'optprop_dustvis_clancy.dat' !Clancy-Lee
 !       file_id(iaer,1) = 'optprop_dustvis_ockert.dat' !Ockert-Bell
-!       file_id(iaer,1) = 'optprop_dustvis.dat'        !Clancy-Lee
 !       Infrared domain:
         file_id(iaer,2) = 'optprop_dustir_TM.dat'      !M.Wolff TM
-!       file_id(iaer,2) = 'optprop_dustir_x0.5.dat'    !Toon-Forget
+!       Toon-Forget + solsir=2 using Clancy-Lee
+!       file_id(iaer,2) = 'optprop_dustir_clancy.dat'
+!       Toon-Forget + solsir=2 using Ockert-Bell
+!       file_id(iaer,2) = 'optprop_dustir_ockert.dat'
 !       Reference wavelength in the visible:
         longrefvis(iaer)=0.67E-6