Changeset 77 for trunk/mars/libf

Timestamp:

Feb 24, 2011, 4:14:33 PM (14 years ago)

Author:

aslmd

Message:

LMD_MM_MARS: nuages et poussiere radiativement actifs

--> tests concluants
VERSION DE REFERENCE pour NOUVELLE PHYSIQUE
-- il reste le nest a implementer (cf. lignes trop longues)

M 76 mars/libf/phymars/dimradmars.h
M 76 mars/libf/phymars/callradite.F
M 76 mesoscale/LMDZ.MARS.new/myGCM/callphys.def
M 76 mesoscale/TESTS/newphys_tracers/namelist.input
reglages pour nuages et poussiere radiativement actifs

M 76 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
M 76 mars/libf/phymars/meso_physiq.F
correction pour bien sortir les depots de glace d'eau au sol [QSURFICE]
correction pour que pd_scalar soit egal a .true. par defaut

M 76 mars/libf/phymars/vdifc.F
M 76 mesoscale/NOTES.txt
M 76 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
commentaires et notes

Location:

trunk/mars/libf/phymars

Files:

• callradite.F (modified) (1 diff)
• dimradmars.h (modified) (1 diff)
• meso_physiq.F (modified) (2 diffs)
• vdifc.F (modified) (1 diff)

trunk/mars/libf/phymars/callradite.F

@@ -263 +263 @@
 c        PLEASE MAKE SURE that you set up the right number of
 c          scatterers in dimradmars.h (naerkind);
-         name_iaer(1) = "dust_conrath"   !! poussiere classique
-c          name_iaer(1) = "dust_doubleq"
+c         name_iaer(1) = "dust_conrath"   !! poussiere classique
+          name_iaer(1) = "dust_doubleq"
 cc        name_iaer(2) = "dust_submicron" !! JB: experimental
-c          name_iaer(2) = "h2o_ice"
+          name_iaer(2) = "h2o_ice"
 c        ----------------------------------------------------------
 

trunk/mars/libf/phymars/dimradmars.h

@@ -27 +27 @@
 ! (ex: naerkind=2 if you use one dust mode and active ice ...)
       integer naerkind
-      parameter (naerkind=1)
-!      parameter (naerkind=2)
+!      parameter (naerkind=1)
+      parameter (naerkind=2)
 
 ! Reference wavelengths used to compute reference optical depth (m)

trunk/mars/libf/phymars/meso_physiq.F

@@ -192 +192 @@
       REAL ust(ngridmx)    !! pour LES avec isfflx!=0
       LOGICAL flag_LES     !! pour LES avec isfflx!=0
-      REAL qsurflast(ngridmx) !! pour diagnostics
+      REAL qsurfice(ngridmx) !! pour diagnostics
       real alpha,lay1 ! coefficients for building layers
       integer iloop
@@ -1577 +1577 @@
          wco2ice(ig) = co2ice(ig)  !! co2 ice 
 
-         !!! TEMP TEMP TEMP TEMP TEMP TEMP TEMP
          !!! specific to WRF WRF WRF
          !!! just to output water ice on surface
-         !!! [it might not be water ice on surface but OK]
          !!! uncomment the Registry entry
-         qsurflast(ig) = qsurf(ig,nqmx)
+         IF (igcm_h2o_ice .ne. 0) qsurfice(ig) = qsurf(ig,igcm_h2o_ice)
 
       enddo

trunk/mars/libf/phymars/vdifc.F

@@ -474 +474 @@
            else if (doubleq) then
              do ig=1,ngrid
+           !!! soulevement constant
                  pdqsdif(ig,igcm_dust_mass) =
      &             -alpha_lift(igcm_dust_mass)